A Time-Calibrated Road Map of Brassicaceae Species Radiation and Evolutionary History

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.

Transition of survival strategies under global climate shifts in the grape family

Faced with environmental changes, plants may either move to track their ancestral niches or evolve to adapt to new niches. Vitaceae, the grape family, has evolved diverse adaptive traits facilitating a global expansion in wide-ranging habitats, making it ideal for investigating transition between move and evolve strategies and exploring the underlying mechanisms. Here we inferred the patterns of biogeographic diversification and trait evolution in Vitaceae based on a robust phylogeny with dense sampling including 495 species (~52% of Vitaceae species). Vitaceae probably originated from Asia-the diversity centre of extant genera and the major source of dispersals. Boundaries of the Eocene, Oligocene and Miocene were identified as turning points in shifting strategies. A significant decrease in move strategy was identified during the Oligocene, followed by increases in move and evolve. After the Miocene, evolve began to dominate, during which increased niche opportunities and key trait innovations played important roles.

Impact of whole-genome duplications on structural variant evolution in Cochlearia

Polyploidy, the result of whole-genome duplication (WGD), is a major driver of eukaryote evolution. Yet WGDs are hugely disruptive mutations, and we still lack a clear understanding of their fitness consequences. Here, we study whether WGDs result in greater diversity of genomic structural variants (SVs) and how they influence evolutionary dynamics in a plant genus, Cochlearia (Brassicaceae). By using long-read sequencing and a graph-based pangenome, we find both negative and positive interactions between WGDs and SVs. Masking of recessive mutations due to WGDs leads to a progressive accumulation of deleterious SVs across four ploidal levels (from diploids to octoploids), likely reducing the adaptive potential of polyploid populations. However, we also discover putative benefits arising from SV accumulation, as more ploidy-specific SVs harbor signals of local adaptation in polyploids than in diploids. Together, our results suggest that SVs play diverse and contrasting roles in the evolutionary trajectories of young polyploids.

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.

IPD3, a master regulator of arbuscular mycorrhizal symbiosis, affects genes for immunity and metabolism of non-host Arabidopsis when restored long after its evolutionary loss

Arbuscular mycorrhizal symbiosis (AM) is a beneficial trait originating with the first land plants, which has subsequently been lost by species scattered throughout the radiation of plant diversity to the present day, including the model Arabidopsis thaliana. To explore if elements of this apparently beneficial trait are still present and could be reactivated we generated Arabidopsis plants expressing a constitutively active form of Interacting Protein of DMI3, a key transcription factor that enables AM within the Common Symbiosis Pathway, which was lost from Arabidopsis along with the AM host trait. We characterize the transcriptomic effect of expressing IPD3 in Arabidopsis with and without exposure to the AM fungus (AMF) Rhizophagus irregularis, and compare these results to the AM model Lotus japonicus and its ipd3 knockout mutant cyclops-4. Despite its long history as a non-AM species, restoring IPD3 in the form of its constitutively active DNA-binding domain to Arabidopsis altered expression of specific gene networks. Surprisingly, the effect of expressing IPD3 in Arabidopsis and knocking it out in Lotus was strongest in plants not exposed to AMF, which is revealed to be due to changes in IPD3 genotype causing a transcriptional state, which partially mimics AMF exposure in non-inoculated plants. Our results indicate that molecular connections to symbiosis machinery remain in place in this nonAM species, with implications for both basic science and the prospect of engineering this trait for agriculture.

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.

Research progress on Brassicaceae plants: a bibliometrics analysis

The Brassicaceae is a worldwide family that produces ornamental flowers, edible vegetables, and oilseed plants, with high economic value in agriculture, horticulture, and landscaping. This study used the Web of Science core dataset and the CiteSpace bibliometric tool to quantitatively visualize the number of publications, authors, institutions, and countries of 3139 papers related to Brassicaceae plants from 2002 to 2022. The keywords and references were divided into two phases: Phase 1 (2002-2011) and Phase 2 (2012-2022) for quantitative and qualitative analysis. The results showed: An average annual publication volume of 149 articles, with an overall fluctuating upward trend; the research force was mainly led by Professor Ihsan A. Al-shehbaz from Missouri Botanical Garden; and the United States had the highest number of publications. In the first phase, research focused on the phylogeny of Brassicaceae plants, while the second phase delved into diverse research based on previous studies, research in areas such as polyploidy, molecular technique, physiology, and hyperaccumulator has been extended. Based on this research, we propounded some ideas for future studies on Brassicaceae plants and summarized the research gaps.

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.

Comparative chloroplast genome analyses provide insights into evolutionary history of Rhizophoraceae mangroves

Background

The Rhizophoraceae family comprises crucial mangrove plants that inhabit intertidal environments. In China, eight Rhizophoraceae mangrove species exist. Although complete chloroplast (Cp) genomes of four Rhizophoraceae mangrove plants have been reported, the Cp genomes of the remaining four species remain unclear, impeding a comprehensive understanding of the evolutionary history of this family.

Methods

Illumina high-throughput sequencing was employed to obtain the DNA sequences of Rhizophoraceae species. Cp genomes were assembled by NOVOPlasty and annotated using CpGAVAS software. Phylogenetic and divergence time analyses were conducted using MEGA and BEAST 2 software.

Results

Four novel Cp genomes of Rhizophoraceae mangrove species (Bruguiera sexangula, Bruguiera gymnorrhiza, Bruguiera × rhynchopetala and Rhizophora apiculata) were successfully assembled. The four Cp genomes ranged in length from 163,310 to 164,560 bp, with gene numbers varying from 124 to 128. The average nucleotide diversity (Pi) value of the eight Rhizophoraceae Cp genomes was 0.00596. Phylogenetic trees constructed based on the complete Cp genomes supported the monophyletic origin of Rhizophoraceae. Divergence time estimation based on the Cp genomes of representative species from Malpighiales showed that the origin of Rhizophoraceae occurred at approximately 58.54-50.02 million years ago (Mya). The divergence time within the genus Rhizophora (∼4.51 Mya) was much earlier than the divergence time within the genus Bruguiera (∼1.41 Mya), suggesting recent speciation processes in these genera. Our data provides new insights into phylogenetic relationship and evolutionary history of Rhizophoraceae mangrove plants.

Haplotype-resolved chromosomal-level assembly of wasabi (Eutrema japonicum) genome

In Japan, wasabi (Eutrema japonicum) is an important traditional condiment, and is recognized as an endemic species. In the present study, we generated a chromosome-level and haplotype-resolved reference genome for E. japonicum using PacBio CLR (continuous long reads), Illumina, and Hi-C sequencing data. The genome consists of 28 chromosomes that contain 1,512.1 Mb of sequence data, with a scaffold N50 length of 55.67 Mb. We also reported the subgenome and haplotype assignment of the 28 chromosomes by read-mapping and phylogenic analysis. Three validation methods (Benchmarking Universal Single-Copy Orthologs, Merqury, and Inspector) indicated that our obtained genome sequences were a high-quality and high-completeness genome assembly. Comparison of genome assemblies from previously published genomes showed that our obtained genome was of higher quality. Therefore, our genome will serve as a valuable genetic resource for both chemical ecology and evolution research of the genera Eutrema and Brassicaceae, as well as for wasabi breeding.

The evolution of the hypotetraploid Catolobus pendulus genome - the poorly known sister species of Capsella

The establishment of Arabidopsis as the most important plant model has also brought other crucifer species into the spotlight of comparative research. While the genus Capsella has become a prominent crucifer model system, its closest relative has been overlooked. The unispecific genus Catolobus is native to temperate Eurasian woodlands, from eastern Europe to the Russian Far East. Here, we analyzed chromosome number, genome structure, intraspecific genetic variation, and habitat suitability of Catolobus pendulus throughout its range. Unexpectedly, all analyzed populations were hypotetraploid (2n = 30, ~330 Mb). Comparative cytogenomic analysis revealed that the Catolobus genome arose by a whole-genome duplication in a diploid genome resembling Ancestral Crucifer Karyotype (ACK, n = 8). In contrast to the much younger Capsella allotetraploid genomes, the presumably autotetraploid Catolobus genome (2n = 32) arose early after the Catolobus/Capsella divergence. Since its origin, the tetraploid Catolobus genome has undergone chromosomal rediploidization, including a reduction in chromosome number from 2n = 32 to 2n = 30. Diploidization occurred through end-to-end chromosome fusion and other chromosomal rearrangements affecting a total of six of 16 ancestral chromosomes. The hypotetraploid Catolobus cytotype expanded toward its present range, accompanied by some longitudinal genetic differentiation. The sister relationship between Catolobus and Capsella allows comparative studies of tetraploid genomes of contrasting ages and different degrees of genome diploidization.

Re-engineering a lost trait: IPD3, a master regulator of arbuscular mycorrhizal symbiosis, affects genes for immunity and metabolism of non-host Arabidopsis when restored long after its evolutionary loss

Arbuscular mycorrhizal symbiosis (AM) is a beneficial trait originating with the first land plants, which has subsequently been lost by species scattered throughout the radiation of plant diversity to the present day, including the model Arabidopsis thaliana. To explore why an apparently beneficial trait would be repeatedly lost, we generated Arabidopsis plants expressing a constitutively active form of Interacting Protein of DMI3, a key transcription factor that enables AM within the Common Symbiosis Pathway, which was lost from Arabidopsis along with the AM host trait. We characterize the transcriptomic effect of expressing IPD3 in Arabidopsis with and without exposure to the AM fungus (AMF) Rhizophagus irregularis, and compare these results to the AM model Lotus japonicus and its ipd3 knockout mutant cyclops-4. Despite its long history as a non-AM species, restoring IPD3 in the form of its constitutively active DNA-binding domain to Arabidopsis altered expression of specific gene networks. Surprisingly, the effect of expressing IPD3 in Arabidopsis and knocking it out in Lotus was strongest in plants not exposed to AMF, which is revealed to be due to changes in IPD3 genotype causing a transcriptional state which partially mimics AMF exposure in non-inoculated plants. Our results indicate that despite the long interval since loss of AM and IPD3 in Arabidopsis, molecular connections to symbiosis machinery remain in place in this nonAM species, with implications for both basic science and the prospect of engineering this trait for agriculture.

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.

Ancestral self-compatibility facilitates the establishment of allopolyploids in Brassicaceae

Self-incompatibility systems based on self-recognition evolved in hermaphroditic plants to maintain genetic variation of offspring and mitigate inbreeding depression. Despite these benefits in diploid plants, for polyploids who often face a scarcity of mating partners, self-incompatibility can thwart reproduction. In contrast, self-compatibility provides an immediate advantage: a route to reproductive viability. Thus, diploid selfing lineages may facilitate the formation of new allopolyploid species. Here, we describe the mechanism of establishment of at least four allopolyploid species in Brassicaceae (Arabidopsis suecica, Arabidopsis kamchatica, Capsella bursa-pastoris, and Brassica napus), in a manner dependent on the prior loss of the self-incompatibility mechanism in one of the ancestors. In each case, the degraded S-locus from one parental lineage was dominant over the functional S-locus of the outcrossing parental lineage. Such dominant loss-of-function mutations promote an immediate transition to selfing in allopolyploids and may facilitate their establishment.

The genome wide analysis of Tryptophan Aminotransferase Related gene family, and their relationship with related agronomic traits in Brassica napus

Tryptophan Aminotransferase of Arabidopsis1/Tryptophan Aminotransferase-Related (TAA1/TAR) proteins are the enzymes that involved in auxin biosynthesis pathway. The TAA1/TAR gene family has been systematically characterized in several plants but has not been well reported in Brassica napus. In the present study, a total of 102 BnTAR genes with different number of introns were identified. It was revealed that these genes are distributed unevenly and occurred as clusters on different chromosomes except for A4, A5, A10 and C4 in B. napus. Most of the these BnTAR genes are conserved despite of existing of gene loss and gene gain. In addition, the segmental replication and whole-genome replication events were both play an important role in the BnTAR gene family formation. Expression profiles analysis indicated that the expression of BnTAR gene showed two patterns, part of them were mainly expressed in roots, stems and leaves of vegetative organs, and the others were mainly expressed in flowers and seeds of reproductive organs. Further analysis showed that many of BnTAR genes were located in QTL intervals of oil content or seed weight, for example BnAMI10 was located in cqOC-C5-4 and cqSW-A2-2, it indicated that some of the BnTAR genes might have relationship with these two characteristics. This study provides a multidimensional analysis of the TAA1/TAR gene family and a new insight into its biological function in B. napus.

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.

Environmental response in gene expression and DNA methylation reveals factors influencing the adaptive potential of Arabidopsis lyrata

Understanding what factors influence plastic and genetic variation is valuable for predicting how organisms respond to changes in the selective environment. Here, using gene expression and DNA methylation as molecular phenotypes, we study environmentally induced variation among Arabidopsis lyrata plants grown at lowland and alpine field sites. Our results show that gene expression is highly plastic, as many more genes are differentially expressed between the field sites than between populations. These environmentally responsive genes evolve under strong selective constraint – the strength of purifying selection on the coding sequence is high, while the rate of adaptive evolution is low. We find, however, that positive selection on cis-regulatory variants has likely contributed to the maintenance of genetically variable environmental responses, but such variants segregate only between distantly related populations. In contrast to gene expression, DNA methylation at genic regions is largely insensitive to the environment, and plastic methylation changes are not associated with differential gene expression. Besides genes, we detect environmental effects at transposable elements (TEs): TEs at the high-altitude field site have higher expression and methylation levels, suggestive of a broad-scale TE activation. Compared to the lowland population, plants native to the alpine environment harbor an excess of recent TE insertions, and we observe that specific TE families are enriched within environmentally responsive genes. Our findings provide insight into selective forces shaping plastic and genetic variation. We also highlight how plastic responses at TEs can rapidly create novel heritable variation in stressful conditions.

Diploid chromosome-level reference genome and population genomic analyses provide insights into Gypenoside biosynthesis and demographic evolution of Gynostemma pentaphyllum (Cucurbitaceae)

Gynostemma pentaphyllum (Thunb.) Makino is a perennial creeping herbaceous plant in the family Cucurbitaceae, which has great medicinal value and commercial potential, but urgent conservation efforts are needed due to the gradual decreases and fragmented distribution of its wild populations. Here, we report the high-quality diploid chromosome-level genome of G. pentaphyllum obtained using a combination of next-generation sequencing short reads, Nanopore long reads, and Hi-C sequencing technologies. The genome is anchored to 11 pseudo-chromosomes with a total size of 608.95 Mb and 26 588 predicted genes. Comparative genomic analyses indicate that G. pentaphyllum is estimated to have diverged from Momordica charantia 60.7 million years ago, with no recent whole-genome duplication event. Genomic population analyses based on genotyping-by-sequencing and ecological niche analyses indicated low genetic diversity but a strong population structure within the species, which could classify 32 G. pentaphyllum populations into three geographical groups shaped jointly by geographic and climate factors. Furthermore, comparative transcriptome analyses showed that the genes encoding enzyme involved in gypenoside biosynthesis had higher expression levels in the leaves and tendrils. Overall, the findings obtained in this study provide an effective molecular basis for further studies of demographic genetics, ecological adaption, and systematic evolution in Cucurbitaceae species, as well as contributing to molecular breeding, and the biosynthesis and biotransformation of gypenoside.

Genome-wide characterization of two Aubrieta taxa: Aubrieta canescens subsp. canescens and Au. macrostyla (Brassicaceae)

Aubrieta canescens complex is divided into two subspecies, Au. canescens subsp. canescens, Au. canescens subsp. cilicica and a distinct species, Au. macrostyla, based on molecular phylogeny. We generated a draft assembly of Au. canescens subsp. canescens and Au. macrostyla using paired-end shotgun sequencing. This is the first attempt at genome characterization for the genus. In the presented study, ~165 and ~157 Mbp of the genomes of Au. canescens subsp. canescens and Au. macrostyla were assembled, respectively, and a total of 32 425 and 31 372 gene models were predicted in the genomes of the target taxa, respectively. We corroborated the phylogenomic affinity of taxa with some core Brassicaceae species (Clades A and B) including Arabis alpina. The orthology-based tree suggested that Aubrieta species differentiated from A. alpina 1.3-2.0 mya (million years ago). The genome-wide syntenic comparison of two Aubrieta taxa revealed that Au. canescens subsp. canescens (46 %) and Au. macrostyla (45 %) have an almost identical syntenic gene pair ratio. These novel genome assemblies are the first steps towards the chromosome-level assembly of Au. canescens and understanding the genome diversity within the genus.

Comparative Analysis on the Evolution of Flowering Genes in Sugar Pathway in Brassicaceae

Sugar plays an important role in regulating the flowering of plants. However, studies of genes related to flowering regulation by the sugar pathway of Brassicaceae plants are scarce. In this study, we performed a comprehensive comparative genomics analysis of the flowering genes in the sugar pathway from seven members of the Brassicaceae, including: Arabidopsis thaliana, Arabidopsis lyrata, Astelia pumila, Camelina sativa, Brassica napus, Brassica oleracea, and Brassica rapa. We identified 105 flowering genes in the sugar pathway of these plants, and they were categorized into nine groups. Protein domain analysis demonstrated that the IDD8 showed striking structural variations in different Brassicaceae species. Selection pressure analysis revealed that sugar pathway genes related to flowering were subjected to strong purifying selection. Collinearity analysis showed that the identified flowering genes expanded to varying degrees, but SUS4 was absent from the genomes of Astelia pumila, Camelina sativa, Brassica napus, Brassica oleracea, and Brassica rapa. Tissue-specific expression of ApADG indicated functional differentiation. To sum up, genome-wide identification revealed the expansion, contraction, and diversity of flowering genes in the sugar pathway during Brassicaceae evolution. This study lays a foundation for further study on the evolutionary characteristics and potential biological functions of flowering genes in the sugar pathway of 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.

Dynamic enhancer transcription associates with reprogramming of immune genes during pattern triggered immunity in Arabidopsis

Background

Enhancers are cis-regulatory elements present in eukaryote genomes, which constitute indispensable determinants of gene regulation by governing the spatiotemporal and quantitative expression dynamics of target genes, and are involved in multiple life processes, for instance during development and disease states. The importance of enhancer activity has additionally been highlighted for immune responses in animals and plants; however, the dynamics of enhancer activities and molecular functions in plant innate immunity are largely unknown. Here, we investigated the involvement of distal enhancers in early innate immunity in Arabidopsis thaliana.

Results

A group of putative distal enhancers producing low-abundance transcripts either unidirectionally or bidirectionally are identified. We show that enhancer transcripts are dynamically modulated in plant immunity triggered by microbe-associated molecular patterns and are strongly correlated with open chromatin, low levels of methylated DNA, and increases in RNA polymerase II targeting and acetylated histone marks. Dynamic enhancer transcription is correlated with target early immune gene expression patterns. Cis motifs that are bound by immune-related transcription factors, such as WRKYs and SARD1, are highly enriched within upregulated enhancers. Moreover, a subset of core pattern-induced enhancers are upregulated by multiple patterns from diverse pathogens. The expression dynamics of putative immunity-related enhancers and the importance of WRKY binding motifs for enhancer function were also validated.

Conclusions

Our study demonstrates the general occurrence of enhancer transcription in plants and provides novel information on the distal regulatory landscape during early plant innate immunity, providing new insights into immune gene regulation and ultimately improving the mechanistic understanding of the plant immune system.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01362-8.

The Arabidopsis effector-triggered immunity landscape is conserved in oilseed crops

The bacterial phytopathogen Pseudomonas syringae causes disease on a wide array of plants, including the model plant Arabidopsis thaliana and its agronomically important relatives in the Brassicaceae family. To cause disease, P. syringae delivers effector proteins into plant cells through a type III secretion system. In response, plant nucleotide-binding leucine-rich repeat proteins recognize specific effectors and mount effector-triggered immunity (ETI). While ETI is pervasive across A. thaliana, with at least 19 families of P. syringae effectors recognized in this model species, the ETI landscapes of crop species have yet to be systematically studied. Here, we investigated the conservation of the A. thaliana ETI landscape in two closely related oilseed crops, Brassica napus (canola) and Camelina sativa (false flax). We show that the level of immune conservation is inversely related to the degree of evolutionary divergence from A. thaliana, with the more closely related C. sativa losing ETI responses to only one of the 19 P. syringae effectors tested, while the more distantly related B. napus loses ETI responses to four effectors. In contrast to the qualitative conservation of immune response, the quantitative rank order is not as well-maintained across the three species and diverges increasingly with evolutionary distance from A. thaliana. Overall, our results indicate that the A. thaliana ETI profile is qualitatively conserved in oilseed crops, but quantitatively distinct.

Chloroplast phylogenomics in Camelina (Brassicaceae) reveals multiple origins of polyploid species and the maternal lineage of C. sativa

The genus Camelina (Brassicaceae) comprises 7-8 diploid, tetraploid, and hexaploid species. Of particular agricultural interest is the biofuel crop, C. sativa (gold-of-pleasure or false flax), an allohexaploid domesticated from the widespread weed, C. microcarpa. Recent cytogenetics and genomics work has uncovered the identity of the parental diploid species involved in ancient polyploidization events in Camelina. However, little is known about the maternal subgenome ancestry of contemporary polyploid species. To determine the diploid maternal contributors of polyploid Camelina lineages, we sequenced and assembled 84 Camelina chloroplast genomes for phylogenetic analysis. Divergence time estimation was used to infer the timing of polyploidization events. Chromosome counts were also determined for 82 individuals to assess ploidy and cytotypic variation. Chloroplast genomes showed minimal divergence across the genus, with no observed gene-loss or structural variation. Phylogenetic analyses revealed C. hispida as a maternal diploid parent to the allotetraploid Camelina rumelica, and C. neglecta as the closest extant diploid contributor to the allohexaploids C. microcarpa and C. sativa. The tetraploid C. rumelica appears to have evolved through multiple independent hybridization events. Divergence times for polyploid lineages closely related to C. sativa were all inferred to be very recent, at only ~65 thousand years ago. Chromosome counts confirm that there are two distinct cytotypes within C. microcarpa (2n = 38 and 2n = 40). Based on these findings and other recent research, we propose a model of Camelina subgenome relationships representing our current understanding of the hybridization and polyploidization history of this recently-diverged genus.

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.

Transposition and duplication of MADS-domain transcription factor genes in annual and perennial Arabis species modulates flowering

Annual and perennial species differ in their timing and intensity of flowering, but the underlying mechanisms are poorly understood. We hybridized closely related annual and perennial plants and used genetics, transgenesis, and genomics to characterize differences in the activity and function of their flowering-time genes. We identify a gene encoding a transcription factor that moved between chromosomes and is retained in the annual but absent from the perennial. This gene strongly delays flowering, and we propose that it has been retained in the annual to compensate for reduced activity of closely related genes. This study highlights the value of using direct hybridization between closely related plant species to characterize functional differences in fast-evolving reproductive traits.

The timing of reproduction is an adaptive trait in many organisms. In plants, the timing, duration, and intensity of flowering differ between annual and perennial species. To identify interspecies variation in these traits, we studied introgression lines derived from hybridization of annual and perennial species, Arabis montbretiana and Arabis alpina, respectively. Recombination mapping identified two tandem A. montbretiana genes encoding MADS-domain transcription factors that confer extreme late flowering on A. alpina. These genes are related to the MADS AFFECTING FLOWERING (MAF) cluster of floral repressors of other Brassicaceae species and were named A. montbretiana (Am) MAF-RELATED (MAR) genes. AmMAR1 but not AmMAR2 prevented floral induction at the shoot apex of A. alpina, strongly enhancing the effect of the MAF cluster, and MAR1 is absent from the genomes of all A. alpina accessions analyzed. Exposure of plants to cold (vernalization) represses AmMAR1 transcription and overcomes its inhibition of flowering. Assembly of the tandem arrays of MAR and MAF genes of six A. alpina accessions and three related species using PacBio long-sequence reads demonstrated that the MARs arose within the Arabis genus by interchromosomal transposition of a MAF1-like gene followed by tandem duplication. Time-resolved comparative RNA-sequencing (RNA-seq) suggested that AmMAR1 may be retained in A. montbretiana to enhance the effect of the AmMAF cluster and extend the duration of vernalization required for flowering. Our results demonstrate that MAF genes transposed independently in different Brassicaceae lineages and suggest that they were retained to modulate adaptive flowering responses that differ even among closely related species.

Myxospermy Evolution in Brassicaceae: A Highly Complex and Diverse Trait with Arabidopsis as an Uncommon Model

The ability to extrude mucilage upon seed imbibition (myxospermy) occurs in several Angiosperm taxonomic groups, but its ancestral nature or evolutionary convergence origin remains misunderstood. We investigated seed mucilage evolution in the Brassicaceae family with comparison to the knowledge accumulated in Arabidopsis thaliana. The myxospermy occurrence was evaluated in 27 Brassicaceae species. Phenotyping included mucilage secretory cell morphology and topochemistry to highlight subtle myxospermy traits. In parallel, computational biology was driven on the one hundred genes constituting the so-called A. thaliana mucilage secretory cell toolbox to confront their sequence conservation to the observed phenotypes. Mucilage secretory cells show high morphology diversity; the three studied Arabidopsis species had a specific extrusion modality compared to the other studied Brassicaceae species. Orthologous genes from the A. thaliana mucilage secretory cell toolbox were mostly found in all studied species without correlation with the occurrence of myxospermy or even more sub-cellular traits. Seed mucilage may be an ancestral feature of the Brassicaceae family. It consists of highly diverse subtle traits, probably underlined by several genes not yet characterized in A. thaliana or by species-specific genes. Therefore, A. thaliana is probably not a sufficient reference for future myxospermy evo-devo studies.

The chloroplast genome of Amygdalus L. (Rosaceae) reveals the phylogenetic relationship and divergence time

Background

Limited access to genetic information has greatly hindered our understanding of the molecular evolution, phylogeny, and differentiation time of subg. Amygdalus. This study reported complete chloroplast (cp) genome sequences of subg. Amygdalus, which further enriched the available valuable resources of complete cp genomes of higher plants and deepened our understanding of the divergence time and phylogenetic relationships of subg. Amygdalus.

Results

The results showed that subg. Amygdalus species exhibited a tetrad structure with sizes ranging from 157,736 bp (P. kansuensis) to 158,971 bp (P. davidiana), a pair of inverted repeat regions (IRa/IRb) that ranged from 26,137–26,467 bp, a large single-copy region that ranged from 85,757–86,608 bp, and a small single-copy region that ranged from 19,020–19,133 bp. The average GC content of the complete cp genomes in the 12 species was 36.80%. We found that the structure of the subg. Amygdalus complete cp genomes was highly conserved, and the 12 subg. Amygdalus species had an rps19 pseudogene. There was not rearrangement of the complete cp genome in the 12 subg. Amygdalus species. All 12 subg. Amygdalus species clustered into one clade based on both Bayesian inference and maximum likelihood. The divergence time analyses based on the complete cp genome sequences showed that subg. Amygdalus species diverged approximately 15.65 Mya.

Conclusion

Our results provide data on the genomic structure of subg. Amygdalus and elucidates their phylogenetic relationships and divergence time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07968-6.

Genomic analysis of field pennycress (Thlaspi arvense) provides insights into mechanisms of adaptation to high elevation

Background

Understanding how organisms evolve and adapt to extreme habitats is of crucial importance in evolutionary ecology. Altitude gradients are an important determinant of the distribution pattern and range of organisms due to distinct climate conditions at different altitudes. High-altitude regions often provide extreme environments including low temperature and oxygen concentration, poor soil, and strong levels of ultraviolet radiation, leading to very few plant species being able to populate elevation ranges greater than 4000 m. Field pennycress (Thlaspi arvense) is a valuable oilseed crop and emerging model plant distributed across an elevation range of nearly 4500 m. Here, we generate an improved genome assembly to understand how this species adapts to such different environments.

Results

We sequenced and assembled de novo the chromosome-level pennycress genome of 527.3 Mb encoding 31,596 genes. Phylogenomic analyses based on 2495 single-copy genes revealed that pennycress is closely related to Eutrema salsugineum (estimated divergence 14.32–18.58 Mya), and both species form a sister clade to Schrenkiella parvula and genus Brassica. Field pennycress contains the highest percentage (70.19%) of transposable elements in all reported genomes of Brassicaceae, with the retrotransposon proliferation in the Middle Pleistocene being likely responsible for the expansion of genome size. Moreover, our analysis of 40 field pennycress samples in two high- and two low-elevation populations detected 1,256,971 high-quality single nucleotide polymorphisms. Using three complementary selection tests, we detected 130 candidate naturally selected genes in the Qinghai-Tibet Plateau (QTP) populations, some of which are involved in DNA repair and the ubiquitin system and potential candidates involved in high-altitude adaptation. Notably, we detected a single base mutation causing loss-of-function of the FLOWERING LOCUS C protein, responsible for the transition to early flowering in high-elevation populations.

Conclusions

Our results provide a genome-wide perspective of how plants adapt to distinct environmental conditions across extreme elevation differences and the potential for further follow-up research with extensive data from additional populations and species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01079-0.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5'-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

A prescient evolutionary model for genesis, duplication and differentiation of MIR160 homologs in Brassicaceae

MicroRNA160 is a class of nitrogen-starvation responsive genes which governs establishment of root system architecture by down-regulating AUXIN RESPONSE FACTOR genes (ARF10, ARF16 and ARF17) in plants. The high copy number of MIR160 variants discovered by us from land plants, especially polyploid crop Brassicas, posed questions regarding genesis, duplication, evolution and function. Absence of studies on impact of whole genome and segmental duplication on retention and evolution of MIR160 homologs in descendent plant lineages prompted us to undertake the current study. Herein, we describe ancestry and fate of MIR160 homologs in Brassicaceae in context of polyploidy driven genome re-organization, copy number and differentiation. Paralogy amongst Brassicaceae MIR160a, MIR160b and MIR160c was inferred using phylogenetic analysis of 468 MIR160 homologs from land plants. The evolutionarily distinct MIR160a was found to represent ancestral form and progenitor of MIR160b and MIR160c. Chronology of evolutionary events resulting in origin and diversification of genomic loci containing MIR160 homologs was delineated using derivatives of comparative synteny. A prescient model for causality of segmental duplications in establishment of paralogy in Brassicaceae MIR160, with whole genome duplication accentuating the copy number increase, is being posited in which post-segmental duplication events viz. differential gene fractionation, gene duplications and inversions are shown to drive divergence of chromosome segments. While mutations caused the diversification of MIR160a, MIR160b and MIR160c, duplicated segments containing these diversified genes suffered gene rearrangements via gene loss, duplications and inversions. Yet the topology of phylogenetic and phenetic trees were found congruent suggesting similar evolutionary trajectory. Over 80% of Brassicaceae genomes and subgenomes showed a preferential retention of single copy each of MIR160a, MIR160b and MIR160c suggesting functional relevance. Thus, our study provides a blue-print for reconstructing ancestry and phylogeny of MIRNA gene families at genomics level and analyzing the impact of polyploidy on organismal complexity. Such studies are critical for understanding the molecular basis of agronomic traits and deploying appropriate candidates for crop improvement.

Genomic basis of parallel adaptation varies with divergence in Arabidopsis and its relatives

Parallel adaptation provides valuable insight into the predictability of evolutionary change through replicated natural experiments. A steadily increasing number of studies have demonstrated genomic parallelism, yet the magnitude of this parallelism varies depending on whether populations, species, or genera are compared. This led us to hypothesize that the magnitude of genomic parallelism scales with genetic divergence between lineages, but whether this is the case and the underlying evolutionary processes remain unknown. Here, we resequenced seven parallel lineages of two Arabidopsis species, which repeatedly adapted to challenging alpine environments. By combining genome-wide divergence scans with model-based approaches, we detected a suite of 151 genes that show parallel signatures of positive selection associated with alpine colonization, involved in response to cold, high radiation, short season, herbivores, and pathogens. We complemented these parallel candidates with published gene lists from five additional alpine Brassicaceae and tested our hypothesis on a broad scale spanning ∼0.02 to 18 My of divergence. Indeed, we found quantitatively variable genomic parallelism whose extent significantly decreased with increasing divergence between the compared lineages. We further modeled parallel evolution over the Arabidopsis candidate genes and showed that a decreasing probability of repeated selection on the same standing or introgressed alleles drives the observed pattern of divergence-dependent parallelism. We therefore conclude that genetic divergence between populations, species, and genera, affecting the pool of shared variants, is an important factor in the predictability of genome evolution.

The tepary bean genome provides insight into evolution and domestication under heat stress

Tepary bean (Phaseolus acutifolis A. Gray), native to the Sonoran Desert, is highly adapted to heat and drought. It is a sister species of common bean (Phaseolus vulgaris L.), the most important legume protein source for direct human consumption, and whose production is threatened by climate change. Here, we report on the tepary genome including exploration of possible mechanisms for resilience to moderate heat stress and a reduced disease resistance gene repertoire, consistent with adaptation to arid and hot environments. Extensive collinearity and shared gene content among these Phaseolus species will facilitate engineering climate adaptation in common bean, a key food security crop, and accelerate tepary bean improvement.

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.

Assembly of the non-heading pak choi genome and comparison with the genomes of heading Chinese cabbage and the oilseed yellow sarson

Brassica rapa displays a wide range of morphological diversity which is exploited for a variety of food crops. Here we present a high-quality genome assembly for pak choi (Brassica rapa L. subsp. chinensis), an important non-heading leafy vegetable, and comparison with the genomes of heading type Chinese cabbage and the oilseed form, yellow sarson. Gene presence-absence variation (PAV) and genomic structural variations (SV) were identified, together with single nucleotide polymorphisms (SNPs). The structure and expression of genes for leaf morphology and flowering were compared between the three morphotypes revealing candidate genes for these traits in B. rapa. The pak choi genome assembly and its comparison with other B. rapa genome assemblies provides a valuable resource for the genetic improvement of this important vegetable crop and as a model to understand the diversity of morphological variation across Brassica species.

Kingdom-wide analysis of the evolution of the plant type III polyketide synthase superfamily

The emergence of type III polyketide synthases (PKSs) was a prerequisite for the conquest of land by the green lineage. Within the PKS superfamily, chalcone synthases (CHSs) provide the entry point reaction to the flavonoid pathway, while LESS ADHESIVE POLLEN 5 and 6 (LAP5/6) provide constituents of the outer exine pollen wall. To study the deep evolutionary history of this key family, we conducted phylogenomic synteny network and phylogenetic analyses of whole-genome data from 126 species spanning the green lineage including Arabidopsis thaliana, tomato (Solanum lycopersicum), and maize (Zea mays). This study thereby combined study of genomic location and context with changes in gene sequences. We found that the two major clades, CHS and LAP5/6 homologs, evolved early by a segmental duplication event prior to the divergence of Bryophytes and Tracheophytes. We propose that the macroevolution of the type III PKS superfamily is governed by whole-genome duplications and triplications. The combined phylogenetic and synteny analyses in this study provide insights into changes in the genomic location and context that are retained for a longer time scale with more recent functional divergence captured by gene sequence alterations.

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.

Gradual evolution of allopolyploidy in Arabidopsis suecica

Most diploid organisms have polyploid ancestors. The evolutionary process of polyploidization is poorly understood but has frequently been conjectured to involve some form of 'genome shock', such as genome reorganization and subgenome expression dominance. Here we study polyploidization in Arabidopsis suecica, a post-glacial allopolyploid species formed via hybridization of Arabidopsis thaliana and Arabidopsis arenosa. We generated a chromosome-level genome assembly of A. suecica and complemented it with polymorphism and transcriptome data from all species. Despite a divergence around 6 million years ago (Ma) between the ancestral species and differences in their genome composition, we see no evidence of a genome shock: the A. suecica genome is colinear with the ancestral genomes; there is no subgenome dominance in expression; and transposon dynamics appear stable. However, we find changes suggesting gradual adaptation to polyploidy. In particular, the A. thaliana subgenome shows upregulation of meiosis-related genes, possibly to prevent aneuploidy and undesirable homeologous exchanges that are observed in synthetic A. suecica, and the A. arenosa subgenome shows upregulation of cyto-nuclear processes, possibly in response to the new cytoplasmic environment of A. suecica, with plastids maternally inherited from A. thaliana. These changes are not seen in synthetic hybrids, and thus are likely to represent subsequent evolution.

The origin, evolution and diversification of multiple isoforms of light-dependent protochlorophyllide oxidoreductase (LPOR): focus on angiosperms

Light-dependent protochlorophyllide oxidoreductase (LPOR) catalyzes the reduction of protochlorophyllide to chlorophyllide, which is a key reaction for angiosperm development. Dark operative light-independent protochlorophyllide oxidoreductase (DPOR) is the other enzyme able to catalyze this reaction, however, it is not present in angiosperms. LPOR, which evolved later than DPOR, requires light to trigger the reaction. The ancestors of angiosperms lost DPOR genes and duplicated the LPORs, however, the LPOR evolution in angiosperms has not been yet investigated. In the present study, we built a phylogenetic tree using 557 nucleotide sequences of LPORs from both bacteria and plants to uncover the evolution of LPOR. The tree revealed that all modern sequences of LPOR diverged from a single sequence ∼1.36 billion years ago. The LPOR gene was then duplicated at least 10 times in angiosperms, leading to the formation of two or even more LPOR isoforms in multiple species. In the case of Arabidopsis thaliana, AtPORA and AtPORB originated in one duplication event, in contrary to the isoform AtPORC, which diverged first. We performed biochemical characterization of these isoforms in vitro, revealing differences in the lipid-driven properties. The results prone us to hypothesize that duplication events of LPOR gave rise to the isoforms having different lipid-driven activity, which may predispose them for functioning in different locations in plastids. Moreover, we showed that LPOR from Synechocystis operated in the lipid-independent manner, revealing differences between bacterial and plant LPORs. Based on the presented results, we propose a novel classification of LPOR enzymes based on their biochemical properties and phylogenetic relationships.

Genome sequencing and phylogenetic analysis of allotetraploid Salix matsudana Koidz

Polyploidy is a common phenomenon among willow species. In this study, genome sequencing was conducted for Salix matsudana Koidz (also named Chinese willow), an important greening and arbor tree species, and the genome of this species was compared with those of four other tree species in Salicaceae. The total genome sequence of S. matsudana was 655.72 Mb in size, with repeated sequences accounting for 45.97% of the total length. In total, 531.43 Mb of the genome sequence could be mapped onto 38 chromosomes using the published genetic map as a reference. The genome of S. matsudana could be divided into two groups, the A and B genomes, through homology analysis with the genome of Populus trichocarpa, and the A and B genomes contained 23,985 and 25,107 genes, respectively. 4DTv combined transposon analysis predicted that allotetraploidy in S. matsudana appeared ~4 million years ago. The results from this study will help reveal the evolutionary history of S. matsudana and lay a genetic basis for its breeding.

Genetic and genomic resources to study natural variation in Brassica rapa

The globally important crop Brassica rapa, a close relative of Arabidopsis, is an excellent system for modeling our current knowledge of plant growth on a morphologically diverse crop. The long history of B. rapa domestication across Asia and Europe provides a unique collection of locally adapted varieties that span large climatic regions with various abiotic and biotic stress-tolerance traits. This diverse gene pool provides a rich source of targets with the potential for manipulation toward the enhancement of productivity of crops both within and outside the Brassicaceae. To expand the genetic resources available to study natural variation in B. rapa, we constructed an Advanced Intercross Recombinant Inbred Line (AI-RIL) population using B. rapa subsp. trilocularis (Yellow Sarson) R500 and the B. rapa subsp. parachinensis (Cai Xin) variety L58. Our current understanding of genomic structure variation across crops suggests that a single reference genome is insufficient for capturing the genetic diversity within a species. To complement this AI-RIL population and current and future B. rapa genomic resources, we generated a de novo genome assembly of the B. rapa subsp. trilocularis (Yellow Sarson) variety R500, the maternal parent of the AI-RIL population. The genetic map for the R500 x L58 population generated using this de novo genome was used to map Quantitative Trait Loci (QTL) for seed coat color and revealed the improved mapping resolution afforded by this new assembly.

Formation and diversification of a paradigm biosynthetic gene cluster in plants

Numerous examples of biosynthetic gene clusters (BGCs), including for compounds of agricultural and medicinal importance, have now been discovered in plant genomes. However, little is known about how these complex traits are assembled and diversified. Here, we examine a large number of variants within and between species for a paradigm BGC (the thalianol cluster), which has evolved recently in a common ancestor of the Arabidopsis genus. Comparisons at the species level reveal differences in BGC organization and involvement of auxiliary genes, resulting in production of species-specific triterpenes. Within species, the thalianol cluster is primarily fixed, showing a low frequency of deleterious haplotypes. We further identify chromosomal inversion as a molecular mechanism that may shuffle more distant genes into the cluster, so enabling cluster compaction. Antagonistic natural selection pressures are likely involved in shaping the occurrence and maintenance of this BGC. Our work sheds light on the birth, life and death of complex genetic and metabolic traits in plants.

Genomic signatures of convergent adaptation to Alpine environments in three Brassicaceae species

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole-genome pooled population sequencing to study genome-wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over-representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species-specific.

Understanding Host-Pathogen Interactions in Brassica napus in the Omics Era

Brassica napus (canola/oilseed rape/rapeseed) is an economically important crop, mostly found in temperate and sub-tropical regions, that is cultivated widely for its edible oil. Major diseases of Brassica crops such as Blackleg, Clubroot, Sclerotinia Stem Rot, Downy Mildew, Alternaria Leaf Spot and White Rust have caused significant yield and economic losses in rapeseed-producing countries worldwide, exacerbated by global climate change, and, if not remedied effectively, will threaten global food security. To gain further insights into the host-pathogen interactions in relation to Brassica diseases, it is critical that we review current knowledge in this area and discuss how omics technologies can offer promising results and help to push boundaries in our understanding of the resistance mechanisms. Omics technologies, such as genomics, proteomics, transcriptomics and metabolomics approaches, allow us to understand the host and pathogen, as well as the interaction between the two species at a deeper level. With these integrated data in multi-omics and systems biology, we are able to breed high-quality disease-resistant Brassica crops in a more holistic, targeted and accurate way.

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.

Ancient orogenic and monsoon-driven assembly of the world's richest temperate alpine flora

Understanding how alpine biotas formed in response to historical environmental change may improve our ability to predict and mitigate the threats to alpine species posed by global warming. In the world's richest temperate alpine flora, that of the Tibet-Himalaya-Hengduan region, phylogenetic reconstructions of biome and geographic range evolution show that extant lineages emerged by the early Oligocene and diversified first in the Hengduan Mountains. By the early to middle Miocene, accelerated diversification and colonization of adjacent regions were likely driven jointly by mountain building and intensification of the Asian monsoon. The alpine flora of the Hengduan Mountains has continuously existed far longer than any other alpine flora on Earth and illustrates how modern biotas have been shaped by past geological and climatic events.

Polyploids increase overall diversity despite higher turnover than diploids in the Brassicaceae

Although polyploidy is widespread across the plant Tree of Life, its long-term evolutionary significance is still poorly understood. Here, we examine the effects of polyploidy in explaining the large-scale evolutionary patterns within angiosperms by focusing on a single family exhibiting extensive interspecific variation in chromosome numbers. We inferred ploidy from haploid chromosome numbers for 80% of species in the most comprehensive species-level chronogram for the Brassicaceae. After evaluating a total of 94 phylogenetic models of diversification, we found that ploidy influences diversification rates across the Brassicaceae. We also found that despite diversifying at a similar rate to diploids, polyploids have played a significant role in driving present-day differences in species richness among clades. Overall, in addition to highlighting the complexity in the evolutionary consequences of polyploidy, our results suggest that rare successful polyploids persist while significantly contributing to the long-term evolution of clades. Our findings further indicate that polyploidy has played a major role in driving the long-term evolution of the Brassicaceae and highlight the potential of polyploidy in shaping present-day diversity patterns across the plant Tree of Life.

Drivers of metabolic diversification: how dynamic genomic neighbourhoods generate new biosynthetic pathways in the Brassicaceae

Plants produce an array of specialized metabolites with important ecological functions. The mechanisms underpinning the evolution of new biosynthetic pathways are not well-understood. Here, we exploit available genome sequence resources to investigate triterpene biosynthesis across the Brassicaceae. Oxidosqualene cyclases (OSCs) catalyze the first committed step in triterpene biosynthesis. Systematic analysis of 13 sequenced Brassicaceae genomes was performed to identify all OSC genes. The genome neighbourhoods (GNs) around a total of 163 OSC genes were investigated to identify Pfam domains significantly enriched in these regions. All-vs-all comparisons of OSC neighbourhoods and phylogenomic analysis were used to investigate the sequence similarity and evolutionary relationships of the numerous candidate triterpene biosynthetic gene clusters (BGCs) observed. Functional analysis of three representative BGCs was carried out and their triterpene pathway products were elucidated. Our results indicate that plant genomes are remarkably plastic, and that dynamic GNs generate new biosynthetic pathways in different Brassicaceae lineages by shuffling the genes encoding a core palette of triterpene-diversifying enzymes, presumably in response to strong environmental selection pressure. These results illuminate a genomic basis for diversification of plant-specialized metabolism through natural combinatorics of enzyme families, which can be mimicked using synthetic biology to engineer diverse bioactive molecules.