Resolution of Brassicaceae Phylogeny Using Nuclear Genes Uncovers Nested Radiations and Supports Convergent Morphological Evolution

Early-Mid Pleistocene genetic differentiation and range expansions as exemplified by invasive Eurasian Bunias orientalis (Brassicaceae) indicates the Caucasus as key region

Turkish Warty cabbage, Bunias orientalis L. (Brassicaceae) is a perennial herb known for its 250 years of invasion history into Europe and worldwide temperate regions. Putative centers of origin were debated to be located in Turkey, the Caucasus or Eastern Europe. Based on the genetic variation from the nuclear and plastid genomes, we identified two major gene pools in the Caucasian-Irano-Turanian region and close to the Northern Caucasus, respectively. These gene pools are old and started to diverge and expand approximately 930 kya in the Caucasus. Pleistocene glaciation and deglaciation cycles favoured later expansion of a European gene pool 230 kya, which was effectively separated from the Caucasian-Irano-Turanian gene pool. Although the European gene pool is genetically less diverse, it has largely served as source for colonization of Western and Northern Europe in modern times with rare observations of genetic contributions from the Caucasian-Irano-Turanian gene pool such as in North-East America. This study largely utilized herbarium material to take advantage of a biodiversity treasure trove providing biological material and also giving access to detailed collection information.

Fruit fracture biomechanics and the release of Lepidium didymum pericarp-imposed mechanical dormancy by fungi

The biomechanical and ecophysiological properties of plant seed/fruit structures are fundamental to survival in distinct environments. Dispersal of fruits with hard pericarps (fruit coats) encasing seeds has evolved many times independently within taxa that have seed dispersal as their default strategy. The mechanisms by which the constraint of a hard pericarp determines germination timing in response to the environment are currently unknown. Here, we show that the hard pericarp of Lepidium didymum controls germination solely by a biomechanical mechanism. Mechanical dormancy is conferred by preventing full phase-II water uptake of the encased non-dormant seed. The lignified endocarp has biomechanically and morphologically distinct regions that serve as predetermined breaking zones. This pericarp-imposed mechanical dormancy is released by the activity of common fungi, which weaken these zones by degrading non-lignified pericarp cells. We propose that the hard pericarp with this biomechanical mechanism contributed to the global distribution of this species in distinct environments.

Fine-tuning timing: natural variation informs the mechanistic basis of the switch to flowering in Arabidopsis thaliana

The evolution of diverse life history strategies has allowed Arabidopsis thaliana to adapt to worldwide locations, spanning a range of latitudinal and environmental conditions. Arabidopsis thaliana accessions are either vernalization-requiring winter annuals or rapid cyclers, with extensive natural variation in vernalization requirement and response. Genetic and molecular analysis of this variation has enhanced our understanding of the mechanisms involved in life history determination, with translation to both natural and crop systems in the Brassicaceae and beyond.

Multispeed genome diploidization and diversification after an ancient allopolyploidization

Hybridization and genome doubling (allopolyploidy) have led to evolutionary novelties as well as to the origin of new clades and species. Despite the importance of allopolyploidization, the dynamics of postpolyploid diploidization (PPD) at the genome level has been only sparsely studied. The Microlepidieae (MICR) is a crucifer tribe of 17 genera and c. 56 species endemic to Australia and New Zealand. Our phylogenetic and cytogenomic analyses revealed that MICR originated via an intertribal hybridization between ancestors of Crucihimalayeae (n = 8; maternal genome) and Smelowskieae (n = 7; paternal genome), both native to the Northern Hemisphere. The reconstructed ancestral allopolyploid genome (n = 15) originated probably in northeastern Asia or western North America during the Late Miocene (c. 10.6-7 million years ago) and reached the Australian mainland via long-distance dispersal. In Australia, the allotetraploid genome diverged into at least three main subclades exhibiting different levels of PPD and diversity: 1.25-fold descending dysploidy (DD) of n = 15 → n = 12 (autopolyploidy → 24) in perennial Arabidella (3 species), 1.5-fold DD of n = 15 → n = 10 in the perennial Pachycladon (11 spp.) and 2.1-3.75-fold DD of n = 15 → n = 7-4 in the largely annual crown-group genera (42 spp. in 15 genera). These results are among the first to demonstrate multispeed genome evolution in taxa descending from a common allopolyploid ancestor. It is suggested that clade-specific PPD can operate at different rates and efficacies and can be tentatively linked to life histories and the extent of taxonomic diversity.

An Arabidopsis introgression zone studied at high spatio-temporal resolution: interglacial and multiple genetic contact exemplified using whole nuclear and plastid genomes

Background

Gene flow between species, across ploidal levels, and even between evolutionary lineages is a common phenomenon in the genus Arabidopsis. However, apart from two genetically fully stabilized allotetraploid species that have been investigated in detail, the extent and temporal dynamics of hybridization are not well understood. An introgression zone, with tetraploid A. arenosa introgressing into A. lyrata subsp. petraea in the Eastern Austrian Forealps and subsequent expansion towards pannonical lowlands, was described previously based on morphological observations as well as molecular data using microsatellite and plastid DNA markers. Here we investigate the spatio-temporal context of this suture zone, making use of the potential of next-generation sequencing and whole-genome data. By utilizing a combination of nuclear and plastid genomic data, the extent, direction and temporal dynamics of gene flow are elucidated in detail and Late Pleistocene evolutionary processes are resolved.

Results

Analysis of nuclear genomic data significantly recognizes the clinal structure of the introgression zone, but also reveals that hybridization and introgression is more common and substantial than previously thought. Also tetraploid A. lyrata and A. arenosa subsp. borbasii from outside the previously defined suture zone show genomic signals of past introgression. A. lyrata is shown to serve usually as the maternal parent in these hybridizations, but one exception is identified from plastome-based phylogenetic reconstruction. Using plastid phylogenomics with secondary time calibration, the origin of A. lyrata and A. arenosa lineages is pre-dating the last three glaciation complexes (approx. 550,000 years ago). Hybridization and introgression followed during the last two glacial-interglacial periods (since approx. 300,000 years ago) with later secondary contact at the northern and southern border of the introgression zone during the Holocene.

Conclusions

Footprints of adaptive introgression in the Northeastern Forealps are older than expected and predate the Last Glaciation Maximum. This correlates well with high genetic diversity found within areas that served as refuge area multiple times. Our data also provide some first hints that early introgressed and presumably preadapted populations account for successful and rapid postglacial re-colonization and range expansion.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4220-6) contains supplementary material, which is available to authorized users.

Molecular Resources from Transcriptomes in the Brassicaceae Family

The rapidly falling costs and the increasing availability of large DNA sequence data sets facilitate the fast and affordable mining of large molecular markers data sets for comprehensive evolutionary studies. The Brassicaceae (mustards) are an important species-rich family in the plant kingdom with taxa distributed worldwide and a complex evolutionary history. We performed Simple Sequence Repeats (SSRs) mining using de novo assembled transcriptomes from 19 species across the Brassicaceae in order to study SSR evolution and provide comprehensive sets of molecular markers for genetic studies within the family. Moreover, we selected the genus Cochlearia to test the transferability and polymorphism of these markers among species. Additionally, we annotated Cochlearia pyrenaica transcriptome in order to identify the position of each of the mined SSRs. While we introduce a new set of tools that will further enable evolutionary studies across the Brassicaceae, we also discuss some broader aspects of SSR evolution. Overall, we developed 2012 ready-to-use SSR markers with their respective primers in 19 Brassicaceae species and a high quality annotated transcriptome for C. pyrenaica. As indicated by our transferability test with the genus Cochlearia these SSRs are transferable to species within the genus increasing exponentially the number of targeted species. Also, our polymorphism results showed substantial levels of variability for these markers. Finally, despite its complex evolutionary history, SSR evolution across the Brassicaceae family is highly conserved and we found no deviation from patterns reported in other Angiosperms.

Monophyletic Origin and Evolution of the Largest Crucifer Genomes

Clade E, or the Hesperis clade, is one of the major Brassicaceae (Crucifereae) clades, comprising some 48 genera and 351 species classified into seven tribes and is distributed predominantly across arid and montane regions of Asia. Several taxa have socioeconomic significance, being important ornamental but also weedy and invasive species. From the comparative genomic perspective, the clade is noteworthy as it harbors species with the largest crucifer genomes but low numbers of chromosomes (n = 5-7). By applying comparative cytogenetic analysis and whole-chloroplast phylogenetics, we constructed, to our knowledge, the first partial and complete cytogenetic maps for selected representatives of clade E tribes and investigated their relationships in a family-wide context. The Hesperis clade is a well-supported monophyletic lineage comprising seven tribes: Anchonieae, Buniadeae, Chorisporeae, Dontostemoneae, Euclidieae, Hesperideae, and Shehbazieae. The clade diverged from other Brassicaceae crown-group clades during the Oligocene, followed by subsequent Miocene tribal diversifications in central/southwestern Asia. The inferred ancestral karyotype of clade E (CEK; n = 7) originated from an older n = 8 genome, which also was the purported progenitor of tribe Arabideae (KAA genome). In most taxa of clade E, the seven linkage groups of CEK either remained conserved (Chorisporeae) or were reshuffled by chromosomal translocations (Euclidieae). In 50% of Anchonieae and Hesperideae species, the CEK genome has undergone descending dysploidy toward n = 6 (-5). These genomic data elucidate early genome evolution in Brassicaceae and pave the way for future whole-genome sequencing and assembly efforts in this as yet genomically neglected group of crucifer plants.

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

Hybridization and polyploidy followed by genome-wide diploidization had a significant impact on the diversification of land plants. The ancient At-α whole-genome duplication (WGD) preceded the diversification of crucifers (Brassicaceae). Some genera and tribes also experienced younger, mesopolyploid WGDs concealed by subsequent genome diploidization. Here we tested if multiple base chromosome numbers originated due to genome diploidization after independent mesopolyploid WGDs and how diploidization affected post-polyploid gene retention. Sixteen species representing 10 Brassicaceae tribes were analyzed by comparative chromosome painting and/or whole-transcriptome analysis of gene age distributions and phylogenetic analyses of gene duplications. Overall, we found evidence for at least 13 independent mesopolyploidies followed by different degrees of diploidization across the Brassicaceae. New mesotetraploid events were uncovered for the tribes Anastaticeae, Iberideae and Schizopetaleae, and mesohexaploid WGDs for Cochlearieae and Physarieae. In contrast, we found convergent patterns of gene retention and loss among these independent WGDs. Our combined analyses of genomic data for Brassicaceae indicate that extant chromosome number variation in many plant groups, and especially monophyletic taxa with multiple base chromosome numbers, can result from clade-specific genome duplications followed by diploidization. Our observation of parallel gene retention and loss across multiple independent WGDs provides one of the first multi-species tests of the predictability of patterns of post-polyploid genome evolution.

Evolution of Rosaceae Fruit Types Based on Nuclear Phylogeny in the Context of Geological Times and Genome Duplication

Fruits are the defining feature of angiosperms, likely have contributed to angiosperm successes by protecting and dispersing seeds, and provide foods to humans and other animals, with many morphological types and important ecological and agricultural implications. Rosaceae is a family with ∼3000 species and an extraordinary spectrum of distinct fruits, including fleshy peach, apple, and strawberry prized by their consumers, as well as dry achenetum and follicetum with features facilitating seed dispersal, excellent for studying fruit evolution. To address Rosaceae fruit evolution and other questions, we generated 125 new transcriptomic and genomic datasets and identified hundreds of nuclear genes to reconstruct a well-resolved Rosaceae phylogeny with highly supported monophyly of all subfamilies and tribes. Molecular clock analysis revealed an estimated age of ∼101.6 Ma for crown Rosaceae and divergence times of tribes and genera, providing a geological and climate context for fruit evolution. Phylogenomic analysis yielded strong evidence for numerous whole genome duplications (WGDs), supporting the hypothesis that the apple tribe had a WGD and revealing another one shared by fleshy fruit-bearing members of this tribe, with moderate support for WGDs in the peach tribe and other groups. Ancestral character reconstruction for fruit types supports independent origins of fleshy fruits from dry-fruit ancestors, including the evolution of drupes (e.g., peach) and pomes (e.g., apple) from follicetum, and drupetum (raspberry and blackberry) from achenetum. We propose that WGDs and environmental factors, including animals, contributed to the evolution of the many fruits in Rosaceae, which provide a foundation for understanding fruit evolution.

Resolution of deep eudicot phylogeny and their temporal diversification using nuclear genes from transcriptomic and genomic datasets

Explosive diversification is widespread in eukaryotes, making it difficult to resolve phylogenetic relationships. Eudicots contain c. 75% of extant flowering plants, are important for human livelihood and terrestrial ecosystems, and have probably experienced explosive diversifications. The eudicot phylogenetic relationships, especially among those of the Pentapetalae, remain unresolved. Here, we present a highly supported eudicot phylogeny and diversification rate shifts using 31 newly generated transcriptomes and 88 other datasets covering 70% of eudicot orders. A highly supported eudicot phylogeny divided Pentapetalae into two groups: one with rosids, Saxifragales, Vitales and Santalales; the other containing asterids, Caryophyllales and Dilleniaceae, with uncertainty for Berberidopsidales. Molecular clock analysis estimated that crown eudicots originated c. 146 Ma, considerably earlier than earliest tricolpate pollen fossils and most other molecular clock estimates, and Pentapetalae sequentially diverged into eight major lineages within c. 15 Myr. Two identified increases of diversification rate are located in the stems leading to Pentapetalae and asterids, and lagged behind the gamma hexaploidization. The nuclear genes from newly generated transcriptomes revealed a well-resolved eudicot phylogeny, sequential separation of major core eudicot lineages and temporal mode of diversifications, providing new insights into the evolutionary trend of morphologies and contributions to the diversification of eudicots.

Young inversion with multiple linked QTLs under selection in a hybrid zone

Fixed chromosomal inversions can reduce gene flow and promote speciation in two ways: by suppressing recombination and by carrying locally favoured alleles at multiple loci. However, it is unknown whether favoured mutations slowly accumulate on older inversions or if young inversions spread because they capture pre-existing adaptive quantitative trait loci (QTLs). By genetic mapping, chromosome painting and genome sequencing, we have identified a major inversion controlling ecologically important traits in Boechera stricta. The inversion arose since the last glaciation and subsequently reached local high frequency in a hybrid speciation zone. Furthermore, the inversion shows signs of positive directional selection. To test whether the inversion could have captured existing, linked QTLs, we crossed standard, collinear haplotypes from the hybrid zone and found multiple linked phenology QTLs within the inversion region. These findings provide the first direct evidence that linked, locally adapted QTLs may be captured by young inversions during incipient speciation.

Evolution of the 3R-MYB Gene Family in Plants

Plant 3R-MYB transcription factors are an important subgroup of the MYB super family in plants; however, their evolutionary history and functions remain poorly understood. We identified 225 3R-MYB proteins from 65 plant species, including algae and all major lineages of land plants. Two segmental duplication events preceding the common ancestor of angiosperms have given rise to three subgroups of the 3R-MYB proteins. Five conserved introns in the domain region of the 3R-MYB genes were identified, which arose through a step-wise pattern of intron gain during plant evolution. Alternative splicing (AS) analysis of selected species revealed that transcripts from more than 60% of 3R-MYB genes undergo AS. AS could regulate transcriptional activity for some of the plant 3R-MYBs by generating different regulatory motifs. The 3R-MYB genes of all subgroups appear to be enriched for Mitosis-Specific Activator element core sequences within their upstream promoter region, which suggests a functional involvement in cell cycle. Notably, expression of 3R-MYB genes from different species exhibits differential regulation under various abiotic stresses. These data suggest that the plant 3R-MYBs function in both cell cycle regulation and abiotic stress response, which may contribute to the adaptation of plants to a sessile lifestyle.

Frequency and maintenance of unreduced gametes in natural plant populations: associations with reproductive mode, life history and genome size

Fertilization involving unreduced (2n) gametes is considered the dominant mechanism of polyploid formation in angiosperms; however, our knowledge of the prevalence of and evolutionary mechanisms maintaining 2n gametes in natural populations is limited. We hypothesize that 2n gametes are deleterious consequences of meiotic errors maintained by mutation-selection balance and should increase in species with relaxed opportunities for selection on sexual processes (asexuality), reduced efficacy of selection (asexuality, selfing) and increased genome instability (high chromosome number). We used flow cytometry to estimate male 2n gamete production in 60 populations from 24 species of Brassicaceae. We quantified variation in 2n gamete production within and among species, and examined associations with life history, reproductive mode, genome size and chromosomal number while accounting for phylogeny. Most individuals produced < 2% 2n male gametes, whereas a small number had > 5% (up to 85%) production. Variation in 2n gamete production was significant among species and related to reproductive system; asexual species produced significantly more 2n gametes than mixed-mating and outcrossing species. Our results, unique in their multi-species perspective, are consistent with 2n gametes being deleterious but maintained when opportunities for selection are limited. Rare individuals with elevated 2n gamete production may be key contributors to polyploid formation.

Using mustard genomes to explore the genetic basis of evolutionary change

Recent advances in sequencing technologies and gene manipulation tools have driven mustard species into the spotlight of comparative research and have offered powerful insight how phenotypic space is explored during evolution. Evidence emerged for genome-wide signal of transcription factors and gene duplication contributing to trait divergence, e.g., PLETHORA5/7 in leaf complexity. Trait divergence is often manifested in differential expression due to cis-regulatory divergence, as in KNOX genes and REDUCED COMPLEXITY, and can be coupled with protein divergence. Fruit shape in Capsella rubella results from anisotropic growth during three distinct phases. Brassicaceae exhibit novel fruit dispersal strategy, explosive pod shatter, where the rapid movement depends on slow build-up of tension and its rapid release facilitated by asymmetric cell wall thickenings.

Plastome phylogeny and early diversification of Brassicaceae

Background

The family Brassicaceae encompasses diverse species, many of which have high scientific and economic importance. Early diversifications and phylogenetic relationships between major lineages or clades remain unclear. Here we re-investigate Brassicaceae phylogeny with complete plastomes from 51 species representing all four lineages or 5 of 6 major clades (A, B, C, E and F) as identified in earlier studies.

Results

Bayesian and maximum likelihood phylogenetic analyses using a partitioned supermatrix of 77 protein coding genes resulted in nearly identical tree topologies exemplified by highly supported relationships between clades. All four lineages were well identified and interrelationships between them were resolved. The previously defined Clade C was found to be paraphyletic (the genus Megadenia formed a separate lineage), while the remaining clades were monophyletic. Clade E (lineage III) was sister to clades B + C rather than to all core Brassicaceae (clades A + B + C or lineages I + II), as suggested by a previous transcriptome study. Molecular dating based on plastome phylogeny supported the origin of major lineages or clades between late Oligocene and early Miocene, and the following radiative diversification across the family took place within a short timescale. In addition, gene losses in the plastomes occurred multiple times during the evolutionary diversification of the family.

Conclusions

Plastome phylogeny illustrates the early diversification of cruciferous species. This phylogeny will facilitate our further understanding of evolution and adaptation of numerous species in the model family Brassicaceae.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3555-3) contains supplementary material, which is available to authorized users.

Plastome phylogeny and early diversification of Brassicaceae

BACKGROUND

The family Brassicaceae encompasses diverse species, many of which have high scientific and economic importance. Early diversifications and phylogenetic relationships between major lineages or clades remain unclear. Here we re-investigate Brassicaceae phylogeny with complete plastomes from 51 species representing all four lineages or 5 of 6 major clades (A, B, C, E and F) as identified in earlier studies.

RESULTS

Bayesian and maximum likelihood phylogenetic analyses using a partitioned supermatrix of 77 protein coding genes resulted in nearly identical tree topologies exemplified by highly supported relationships between clades. All four lineages were well identified and interrelationships between them were resolved. The previously defined Clade C was found to be paraphyletic (the genus Megadenia formed a separate lineage), while the remaining clades were monophyletic. Clade E (lineage III) was sister to clades B + C rather than to all core Brassicaceae (clades A + B + C or lineages I + II), as suggested by a previous transcriptome study. Molecular dating based on plastome phylogeny supported the origin of major lineages or clades between late Oligocene and early Miocene, and the following radiative diversification across the family took place within a short timescale. In addition, gene losses in the plastomes occurred multiple times during the evolutionary diversification of the family.

CONCLUSIONS

Plastome phylogeny illustrates the early diversification of cruciferous species. This phylogeny will facilitate our further understanding of evolution and adaptation of numerous species in the model family Brassicaceae.

Anastatica hierochuntica, an Arabidopsis Desert Relative, Is Tolerant to Multiple Abiotic Stresses and Exhibits Species-Specific and Common Stress Tolerance Strategies with Its Halophytic Relative, Eutrema (Thellungiella) salsugineum

The search for novel stress tolerance determinants has led to increasing interest in plants native to extreme environments - so called "extremophytes." One successful strategy has been comparative studies between Arabidopsis thaliana and extremophyte Brassicaceae relatives such as the halophyte Eutrema salsugineum located in areas including cold, salty coastal regions of China. Here, we investigate stress tolerance in the desert species, Anastatica hierochuntica (True Rose of Jericho), a member of the poorly investigated lineage III Brassicaceae. We show that A. hierochuntica has a genome approximately 4.5-fold larger than Arabidopsis, divided into 22 diploid chromosomes, and demonstrate that A. hierochuntica exhibits tolerance to heat, low N and salt stresses that are characteristic of its habitat. Taking salt tolerance as a case study, we show that A. hierochuntica shares common salt tolerance mechanisms with E. salsugineum such as tight control of shoot Na+ accumulation and resilient photochemistry features. Furthermore, metabolic profiling of E. salsugineum and A. hierochuntica shoots demonstrates that the extremophytes exhibit both species-specific and common metabolic strategies to cope with salt stress including constitutive up-regulation (under control and salt stress conditions) of ascorbate and dehydroascorbate, two metabolites involved in ROS scavenging. Accordingly, A. hierochuntica displays tolerance to methyl viologen-induced oxidative stress suggesting that a highly active antioxidant system is essential to cope with multiple abiotic stresses. We suggest that A. hierochuntica presents an excellent extremophyte Arabidopsis relative model system for understanding plant survival in harsh desert conditions.

Glucosinolate diversity within a phylogenetic framework of the tribe Cardamineae (Brassicaceae) unraveled with HPLC-MS/MS and NMR-based analytical distinction of 70 desulfoglucosinolates

As a basis for future investigations of evolutionary trajectories and biosynthetic mechanisms underlying variations in glucosinolate structures, we screened members of the crucifer tribe Cardamineae by HPLC-MS/MS, isolated and identified glucosinolates by NMR, searched the literature for previous data for the tribe, and collected HPLC-MS/MS data for nearly all glucosinolates known from the tribe as well as some related structures (70 in total). This is a considerable proportion of the approximately 142 currently documented natural glucosinolates. Calibration with authentic references allowed distinction (or elucidation) of isomers in many cases, such as distinction of β-hydroxyls, methylthios, methylsulfinyls and methylsulfonyls. A mechanism for fragmentation of secondary β-hydroxyls in MS was elucidated, and two novel glucosinolates were discovered: 2-hydroxy-3-methylpentylglucosinolate in roots of Cardamine pratensis and 2-hydroxy-8-(methylsulfinyl)octylglucosinolate in seeds of Rorippa amphibia. A large number of glucosinolates (ca. 54 with high structural certainty and a further 28 or more suggested from tandem MS), representing a wide structural variation, is documented from the tribe. This included glucosinolates apparently derived from Met, Phe, Trp, Val/Leu, Ile and higher homologues. Normal side chain elongation and side chain decoration by oxidation or methylation was observed, as well as rare abnormal side chain decoration (hydroxylation of aliphatics at the δ rather than β-position). Some species had diverse profiles, e.g. R. amphibia and C. pratensis (19 and 16 individual glucosinolates, respectively), comparable to total diversity in literature reports of Armoracia rusticana (17?), Barbarea vulgaris (20-24), and Rorippa indica (>20?). The ancestor or the tribe would appear to have used Trp, Met, and homoPhe as glucosinolate precursor amino acids, and to exhibit oxidation of thio to sulfinyl, formation of alkenyls, β-hydroxylation of aliphatic chains and hydroxylation and methylation of indole glucosinolates. Two hotspots of apparent biochemical innovation and loss were identified: C. pratensis and the genus Barbarea. Diversity in other species mainly included structures also known from other crucifers. In addition to a role of gene duplication, two contrasting genetic/biochemical mechanisms for evolution of such combined diversity and redundancy are discussed: (i) involvement of widespread genes with expression varying during evolution, and (ii) mutational changes in substrate specificities of CYP79F and GS-OH enzymes.

Multiple Polyploidization Events across Asteraceae with Two Nested Events in the Early History Revealed by Nuclear Phylogenomics

Biodiversity results from multiple evolutionary mechanisms, including genetic variation and natural selection. Whole-genome duplications (WGDs), or polyploidizations, provide opportunities for large-scale genetic modifications. Many evolutionarily successful lineages, including angiosperms and vertebrates, are ancient polyploids, suggesting that WGDs are a driving force in evolution. However, this hypothesis is challenged by the observed lower speciation and higher extinction rates of recently formed polyploids than diploids. Asteraceae includes about 10% of angiosperm species, is thus undoubtedly one of the most successful lineages and paleopolyploidization was suggested early in this family using a small number of datasets. Here, we used genes from 64 new transcriptome datasets and others to reconstruct a robust Asteraceae phylogeny, covering 73 species from 18 tribes in six subfamilies. We estimated their divergence times and further identified multiple potential ancient WGDs within several tribes and shared by the Heliantheae alliance, core Asteraceae (Asteroideae-Mutisioideae), and also with the sister family Calyceraceae. For two of the WGD events, there were subsequent great increases in biodiversity; the older one proceeded the divergence of at least 10 subfamilies within 10 My, with great variation in morphology and physiology, whereas the other was followed by extremely high species richness in the Heliantheae alliance clade. Our results provide different evidence for several WGDs in Asteraceae and reveal distinct association among WGD events, dramatic changes in environment and species radiations, providing a possible scenario for polyploids to overcome the disadvantages of WGDs and to evolve into lineages with high biodiversity.

The Cardamine hirsuta genome offers insight into the evolution of morphological diversity

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

Patterns of Substitution Rate Variation at Many Nuclear Loci in Two Species Trios in the Brassicaceae Partitioned with ANOVA

There are marked variations among loci and among lineages in rates of nucleotide substitution. The generation time hypothesis (GTH) is a neutral explanation for substitution rate heterogeneity that has genomewide application, predicting that species with shorter generation times accumulate DNA sequence substitutions faster than species with longer generation times do since faster genome replication provides more opportunities for mutations to occur and reach fixation by genetic drift. Relatively few studies have rigorously evaluated the GTH in plants, and there are numerous alternative hypotheses for plant substitution rate variation. One major challenge has been finding pairs of closely related plant species with contrasting generation times and appropriate outgroup taxa that all also have DNA sequence data for numerous loci. To test for causes of rate variation, we obtained sequence data for 256 genes for Arabidopsis thaliana, normally reproducing every year, and the biennial Arabidopsis lyrata with three closely related outgroup taxa (Brassica rapa, Capsella grandiflora, and Neslia paniculata) as well as the biennial Brassica oleracea and the annual B. rapa lineage with the outgroup N. paniculata. A sign test indicated that more loci than expected by chance have faster rates of substitution on the branch leading to the annual than to the perennial for one three-species trio but not another. Tajima's 1D and 2D tests, and a likelihood ratio test that incorporated saturation correction, rejected rate homogeneity for up to 26 genes (up to 14 genes when correcting for multiple tests), consistently showing faster rates for the annual lineage in the Arabidopsis species trio. ANOVA showed significant rate heterogeneity between the Arabidopsis and Brassica species trios (about 6 % of rate variation) and among loci (about 26-32 % of rate variation). The lineage-by-locus interaction which would be caused by locus- and lineage-specific natural selection explained about 13 % of substitution rate variation in one ANOVA model using substitution rates from genes partitioned into odd and even codons but was not a significant effect without partitioned genes. Annual/perennial lineage and species trio by annual/perennial lineage each explained about 1 % of substitution rate variation.

Molecular phylogeny reveals the non-monophyly of tribe Yinshanieae (Brassicaceae) and description of a new tribe, Hillielleae

The taxonomic treatment within the unigeneric tribe Yinshanieae (Brassicaceae) is controversial, owing to differences in generic delimitation applied to its species. In this study, sequences from nuclear ITS and chloroplast trnL-F regions were used to test the monophyly of Yinshanieae, while two nuclear markers (ITS, ETS) and four chloroplast markers (trnL-F, trnH-psbA, rps16, rpL32-trnL) were used to elucidate the phylogenetic relationships within the tribe. Using maximum parsimony, maximum likelihood, and Bayesian inference methods, we reconstructed the phylogeny of Brassicaceae and Yinshanieae. The results show that Yinshanieae is not a monophyletic group, with the taxa splitting into two distantly related clades: one clade contains four taxa and falls in Lineage I, whereas the other includes all species previously placed in Hilliella and is embedded in the Expanded Lineage II. The tribe Yinshanieae is redefined, and a new tribe, Hillielleae, is proposed based on combined evidence from molecular phylogeny, morphology, and cytology.

Comparative paleogenomics of crucifers: ancestral genomic blocks revisited

A decade ago the concept of the Ancestral Crucifer Karyotype (ACK) and the definition of 24 conserved genomic blocks was presented. Subsequently, 35 cytogenetic reconstructions and/or draft genome sequences of crucifer species (members of the Brassicaceae family) have been analyzed in the context of this system; placing crucifers at the forefront of plant phylogenomics. In this review, we highlight how the ACK and genomic blocks have facilitated and guided genomic analysis of crucifers in the last 10 years and provide an update of this robust model.

Biogeography and diversification of Brassicales: A 103million year tale

Brassicales is a diverse order perhaps most famous because it houses Brassicaceae and, its premier member, Arabidopsis thaliana. This widely distributed and species-rich lineage has been overlooked as a promising system to investigate patterns of disjunct distributions and diversification rates. We analyzed plastid and mitochondrial sequence data from five gene regions (>8000bp) across 151 taxa to: (1) produce a chronogram for major lineages in Brassicales, including Brassicaceae and Arabidopsis, based on greater taxon sampling across the order and previously overlooked fossil evidence, (2) examine biogeographical ancestral range estimations and disjunct distributions in BioGeoBEARS, and (3) determine where shifts in species diversification occur using BAMM. The evolution and radiation of the Brassicales began 103Mya and was linked to a series of inter-continental vicariant, long-distance dispersal, and land bridge migration events. North America appears to be a significant area for early stem lineages in the order. Shifts to Australia then African are evident at nodes near the core Brassicales, which diverged 68.5Mya (HPD=75.6-62.0). This estimated age combined with fossil evidence, indicates that some New World clades embedded amongst Old World relatives (e.g., New World capparoids) are the result of different long distance dispersal events, whereas others may be best explained by land bridge migration (e.g., Forchhammeria). Based on these analyses, the Brassicaceae crown group diverged in Europe/Northern Africa in the Eocene, circa 43.4Mya (HPD=46.6-40.3) and Arabidopsis separated from close congeners circa 10.4Mya. These ages fall between divergent dates that were previously published, suggesting we are slowly converging on a robust age estimate for the family. Three significant shifts in species diversification are observed in the order: (1) 58Mya at the crown of Capparaceae, Cleomaceae and Brassicaceae, (2) 38Mya at the crown of Resedaceae+Stixis clade, and (3) 21Mya at the crown of the tribes Brassiceae and Sisymbrieae within Brassicaceae.

Arabidopsis Cell Division Cycle 20.1 Is Required for Normal Meiotic Spindle Assembly and Chromosome Segregation

Cell division requires proper spindle assembly; a surveillance pathway, the spindle assembly checkpoint (SAC), monitors whether the spindle is normal and correctly attached to kinetochores. The SAC proteins regulate mitotic chromosome segregation by affecting CDC20 (Cell Division Cycle 20) function. However, it is unclear whether CDC20 regulates meiotic spindle assembly and proper homolog segregation. Here, we show that the Arabidopsis thaliana CDC20.1 gene is indispensable for meiosis and male fertility. We demonstrate that cdc20.1 meiotic chromosomes align asynchronously and segregate unequally and the metaphase I spindle has aberrant morphology. Comparison of the distribution of meiotic stages at different time points between the wild type and cdc20.1 reveals a delay of meiotic progression from diakinesis to anaphase I. Furthermore, cdc20.1 meiocytes exhibit an abnormal distribution of a histone H3 phosphorylation mark mediated by the Aurora kinase, providing evidence that CDC20.1 regulates Aurora localization for meiotic chromosome segregation. Further evidence that CDC20.1 and Aurora are functionally related was provided by meiosis-specific knockdown of At-Aurora1 expression, resulting in meiotic chromosome segregation defects similar to those of cdc20.1. Taken together, these results suggest a critical role for CDC20.1 in SAC-dependent meiotic chromosome segregation.