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

The water lily genome and the early evolution of flowering plants

Water lilies belong to the angiosperm order Nymphaeales. Amborellales, Nymphaeales and Austrobaileyales together form the so-called ANA-grade of angiosperms, which are extant representatives of lineages that diverged the earliest from the lineage leading to the extant mesangiosperms1-3. Here we report the 409-megabase genome sequence of the blue-petal water lily (Nymphaea colorata). Our phylogenomic analyses support Amborellales and Nymphaeales as successive sister lineages to all other extant angiosperms. The N. colorata genome and 19 other water lily transcriptomes reveal a Nymphaealean whole-genome duplication event, which is shared by Nymphaeaceae and possibly Cabombaceae. Among the genes retained from this whole-genome duplication are homologues of genes that regulate flowering transition and flower development. The broad expression of homologues of floral ABCE genes in N. colorata might support a similarly broadly active ancestral ABCE model of floral organ determination in early angiosperms. Water lilies have evolved attractive floral scents and colours, which are features shared with mesangiosperms, and we identified their putative biosynthetic genes in N. colorata. The chemical compounds and biosynthetic genes behind floral scents suggest that they have evolved in parallel to those in mesangiosperms. Because of its unique phylogenetic position, the N. colorata genome sheds light on the early evolution of angiosperms.

Maternal Inheritance of U's Triangle and Evolutionary Process of Brassica Mitochondrial Genomes

The sequences and genomic structures of plant mitochondrial (mt) genomes provide unique material for phylogenetic studies. The nature of uniparental inheritance renders an advantage when utilizing mt genomes for determining the parental sources of hybridized taxa. In this study, a concatenated matrix of mt genes was used to infer the phylogenetic relationships of six cultivated Brassica taxa and explore the maternal origins of three allotetraploids. The well-resolved sister relationships between two pairs of diploid and allotetraploid taxa suggest that Brassica carinata (car) possessed a maternal origin from Brassica nigra, while Brassica juncea (jun) was maternally derived from Brassica rapa (cam). Another allotetraploid taxon, Brassica napus (cv. Wester) may have been maternally derived from the common ancestor of B. rapa and Brassica oleracea (ole), and/or have undergone (an) extra hybridization event(s) along its evolutionary history. The characteristics of Brassica mt genomic structures also supported the phylogenetic results. Sinapis arvensis was nested inside the Brassica species, sister to the B. nigra-B. carinata lineage, and possessed an mt genome structure that mostly resembled B. nigra. Collectively, the evidence supported a systematic revision that placed S. arvensis within Brassica. Finally, ancestral mt genomes at each evolutionary node of Brassica were reconstructed, and the detailed and dynamic evolution of Brassica mt genomes was successfully reproduced. The mt genome of B. nigra structurally resembled that of the Brassica ancestor the most, with only one reversion of a block, and the Brassica oleracea underwent the most drastic changes. These findings suggested that repeat-mediated recombinations were largely responsible for the observed structural variations in the evolutionary history of Brassica mt genomes.

Field cress genome mapping: Integrating linkage and comparative maps with cytogenetic analysis for rDNA carrying chromosomes

Field cress (Lepidium campestre L.), despite its potential as a sustainable alternative oilseed plant, has been underutilized, and no prior attempts to characterize the genome at the genetic or molecular cytogenetic level have been conducted. Genetic maps are the foundation for anchoring and orienting annotated genome assemblies and positional cloning of candidate genes. Our principal goal was to construct a genetic map using integrated approaches of genetic, comparative and cytogenetic map analyses. In total, 503 F2 interspecific hybrid individuals were genotyped using 7,624 single nucleotide polymorphism markers. Comparative analysis demonstrated that ~57% of the sequenced loci in L. campestre were congruent with Arabidopsis thaliana (L.) genome and suggested a novel karyotype, which predates the ancestral crucifer karyotype. Aceto-orcein chromosome staining and fluorescence in situ hybridization (FISH) analyses confirmed that L. campestre, L. heterophyllum Benth. and their hybrids had a chromosome number of 2n = 2x = 16. Flow cytometric analysis revealed that both species possess 2C roughly 0.4 picogram DNA. Integrating linkage and comparative maps with cytogenetic map analyses assigned two linkage groups to their particular chromosomes. Future work could incorporate FISH utilizing A. thaliana mapped BAC clones to allow the chromosomes of field cress to be identified reliably.

Genome draft of the Arabidopsis relative Pachycladon cheesemanii reveals novel strategies to tolerate New Zealand's high ultraviolet B radiation environment

Background

Pachycladon cheesemanii is a close relative of Arabidopsis thaliana and is an allotetraploid perennial herb which is widespread in the South Island of New Zealand. It grows at altitudes of up to 1000 m where it is subject to relatively high levels of ultraviolet (UV)-B radiation. To gain first insights into how Pachycladon copes with UV-B stress, we sequenced its genome and compared the UV-B tolerance of two Pachycladon accessions with those of two A. thaliana accessions from different altitudes.

Results

A high-quality draft genome of P. cheesemanii was assembled with a high percentage of conserved single-copy plant orthologs. Synteny analysis with genomes from other species of the Brassicaceae family found a close phylogenetic relationship of P. cheesemanii with Boechera stricta from Brassicaceae lineage I. While UV-B radiation caused a greater growth reduction in the A. thaliana accessions than in the P. cheesemanii accessions, growth was not reduced in one P. cheesemanii accession. The homologues of A. thaliana UV-B radiation response genes were duplicated in P. cheesemanii, and an expression analysis of those genes indicated that the tolerance mechanism in P. cheesemanii appears to differ from that in A. thaliana.

Conclusion

Although the P. cheesemanii genome shows close similarity with that of A. thaliana, it appears to have evolved novel strategies allowing the plant to tolerate relatively high UV-B radiation.

A fossil-calibrated phylogeny reveals the biogeographic history of the Cladrastis clade, an amphi-Pacific early-branching group in papilionoid legumes

The early-branching Cladrastis clade of papilionoid legumes (Leguminosae, Papilionoideae) has an intriguing amphi-Pacific disjunct distribution in eastern Asia and temperate-tropical Americas. Here we used nuclear and three plastid regions to reconstruct the phylogenetic relationships and divergence times in the Cladrastis clade, as well as the evolution of morphological characters that might have been key in its biogeographic history. The ancestral character state estimation revealed that the most recent common ancestor of the Cladrastis clade was deciduous trees possessing compressed, winged fruits. The Cladrastis clade was inferred to have originated in the mid-latitude thermophilic forests of North America in the early Eocene, followed by the split between ancestors of wing-fruited Platyosprion and the non-wing-fruited group, and later the divergence of Cladrastis s.s. from the non-wing-fruited group in middle Eocene. Platyosprion and Cladrastis s.s. display an "out-of-North-America" biogeographic pattern and might have migrated to Asia via the Bering land bridge (BLB) or the North Atlantic land bridges (NALB) during middle to late Eocene. Our results, coupled with the relatively well documented fossil record for the clade, suggest that Platyosprion experienced an extinction event in North America caused by climatic cooling around the Eocene-Oligocene transition, which drove a major vegetation shift in western North America, in turn serving as a barrier for the vicariance of Pickeringia and Styphnolobium. The evolution of shrubby habit and sclerophyllous leaves in the former might be adaption to the chaparral vegetation in southwestern North America; the latter gained the trait of moniliform, succulent fruit. Styphnolobium further dispersed southward to tropical North America in the Oligocene, and eastward to Asia through BLB during middle Miocene. Subsequent sundering of BLB facilitated the vicariance of St. affine and St. japonicum.

Genome Improvement and Genetic Map Construction for Aethionema arabicum, the First Divergent Branch in the Brassicaceae Family

The genus Aethionema is a sister-group to the core-group of the Brassicaceae family that includes Arabidopsis thaliana and the Brassica crops. Thus, Aethionema is phylogenetically well-placed for the investigation and understanding of genome and trait evolution across the family. We aimed to improve the quality of the reference genome draft version of the annual species Aethionema arabicum. Second, we constructed the first Ae. arabicum genetic map. The improved reference genome and genetic map enabled the development of each other. We started with the initially published genome (version 2.5). PacBio and MinION sequencing together with genetic map v2.5 were incorporated to produce the new reference genome v3.0. The improved genome contains 203 MB of sequence, with approximately 94% of the assembly made up of called (non-gap) bases, assembled into 2,883 scaffolds (with only 6% of the genome made up of non-called bases (Ns)). The N₅₀ (10.3 MB) represents an 80-fold increase over the initial genome release. We generated a Recombinant Inbred Line (RIL) population that was derived from two ecotypes: Cyprus and Turkey (the reference genotype. Using a Genotyping by Sequencing (GBS) approach, we generated a high-density genetic map with 749 (v2.5) and then 632 SNPs (v3.0) was generated. The genetic map and reference genome were integrated, thus greatly improving the scaffolding of the reference genome into 11 linkage groups. We show that long-read sequencing data and genetics are complementary, resulting in an improved genome assembly in Ae. arabicum. They will facilitate comparative genetic mapping work for the Brassicaceae family and are also valuable resources to investigate wide range of life history traits in Aethionema.

Evolution and functional differentiation of recently diverged phytochelatin synthase genes from Arundo donax L

Phytochelatin synthases (PCSs) play pivotal roles in the detoxification of heavy metals and metalloids in plants; however, little information on the evolution of recently duplicated PCS genes in plant species is available. Here we characterize the evolution and functional differentiation of three PCS genes from the giant reed (Arundo donax L.), a biomass/bioenergy crop with remarkable resistance to cadmium and other heavy metals. Phylogenetic reconstruction with PCS genes from fully sequenced monocotyledonous genomes indicated that the three A. donax PCSs, namely AdPCS1-3, form a monophyletic clade. The AdPCS1-3 genes were expressed at low levels in many A. donax organs and displayed different levels of cadmium-responsive expression in roots. Overexpression of AdPCS1-3 in Arabidopsis thaliana and yeast reproduced the phenotype of functional PCS genes. Mass spectrometry analyses confirmed that AdPCS1-3 are all functional enzymes, but with significant differences in the amount of the phytochelatins synthesized. Moreover, heterogeneous evolutionary rates characterized the AdPCS1-3 genes, indicative of relaxed natural selection. These results highlight the elevated functional differentiation of A. donax PCS genes from both a transcriptional and an enzymatic point of view, providing evidence of the high evolvability of PCS genes and of plant responsiveness to heavy metal stress.

Trichoderma harzianum favours the access of arbuscular mycorrhizal fungi to non-host Brassicaceae roots and increases plant productivity

The family Brassicaceae includes plants that are non-host for arbuscular mycorrhizal fungi (AMF) such as the model plant Arabidopsis thaliana (arabidopsis) and the economically important crop plant Brassica napus (rapeseed). It is well known that Trichoderma species have the ability to colonize the rhizosphere of Brassicaceae plants, promoting growth and development as well as stimulating systemic defenses. The aim of the present work is to ascertain that Brassicaceae plants increase productivity when AMF and Trichoderma are combinedly applied, and how such an effect can be ruled. This simultaneous application of a Trichoderma harzianum biocontrol strain and an AMF formulation produces a significant increase in the colonization by Trichoderma and the presence of AMF in arabidopsis and rapeseed roots, such colonization accompanied by improved productivity in both Brassicaceae species. Expression profiling of defense-related marker genes suggests that the phytohormone salicylic acid plays a key role in the modulation of the root colonization process when both fungi are jointly applied.

Intergeneric Relationships within the Family Salicaceae s.l. based on Plastid Phylogenomics

Many Salicaceae s.l. plants are recognized for their important role in the production of products such as wood, oils, and medicines, and as a model organism in life studies. However, the difference in plastid sequence, phylogenetic relationships, and lineage diversification of the family Salicaceae s.l. remain poorly understood. In this study, we compare 24 species representing 18 genera of the family. Simple sequence repeats (SSRs) are considered effective molecular markers for plant species identification and population genetics. Among them, a total of 1798 SSRs were identified, among which mononucleotide repeat was the most common with 1455 accounts representing 80.92% of the total. Most of the SSRs are located in the non-coding region. We also identified five other types of repeats, including 1750 tandems, 434 forward, 407 palindromic, 86 reverse, and 30 complementary repeats. The species in Salicaceae s.l. have a conserved plastid genome. Each plastome presented a typical quadripartite structure and varied in size due to the expansion and contraction of the inverted repeat (IR) boundary, lacking major structural variations, but we identified six divergence hotspot regions. We obtained phylogenetic relationships of 18 genera in Salicaceae s.l. and the 24 species formed a highly supported lineage. Casearia was identified as the basal clade. The divergence time between Salicaceae s.l. and the outgroup was estimated as ~93 Mya; Salix, and Populus diverged around 34 Mya, consistent with the previously reported time. Our research will contribute to a better understanding of the phylogenetic relationships among the members of the Salicaceae s.l.

The large genome size variation in the Hesperis clade was shaped by the prevalent proliferation of DNA repeats and rarer genome downsizing

BACKGROUND AND AIMS

Most crucifer species (Brassicaceae) have small nuclear genomes (mean 1C-value 617 Mb). The species with the largest genomes occur within the monophyletic Hesperis clade (Mandáková et al., Plant Physiology174: 2062-2071; also known as Clade E or Lineage III). Whereas most chromosome numbers in the clade are 6 or 7, monoploid genome sizes vary 16-fold (256-4264 Mb). To get an insight into genome size evolution in the Hesperis clade (~350 species in ~48 genera), we aimed to identify, quantify and localize in situ the repeats from which these genomes are built. We analysed nuclear repeatomes in seven species, covering the phylogenetic and genome size breadth of the clade, by low-pass whole-genome sequencing.

METHODS

Genome size was estimated by flow cytometry. Genomic DNA was sequenced on an Illumina sequencer and DNA repeats were identified and quantified using RepeatExplorer; the most abundant repeats were localized on chromosomes by fluorescence in situ hybridization. To evaluate the feasibility of bacterial artificial chromosome (BAC)-based comparative chromosome painting in Hesperis-clade species, BACs of arabidopsis were used as painting probes.

KEY RESULTS

Most biennial and perennial species of the Hesperis clade possess unusually large nuclear genomes due to the proliferation of long terminal repeat retrotransposons. The prevalent genome expansion was rarely, but repeatedly, counteracted by purging of transposable elements in ephemeral and annual species.

CONCLUSIONS

The most common ancestor of the Hesperis clade has experienced genome upsizing due to transposable element amplification. Further genome size increases, dominating diversification of all Hesperis-clade tribes, contrast with the overall stability of chromosome numbers. In some subclades and species genome downsizing occurred, presumably as an adaptive transition to an annual life cycle. The amplification versus purging of transposable elements and tandem repeats impacted the chromosomal architecture of the Hesperis-clade species.

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

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

Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta

BACKGROUND

Genomic variation is widespread, and both neutral and selective processes can generate similar patterns in the genome. These processes are not mutually exclusive, so it is difficult to infer the evolutionary mechanisms that govern population and species divergence. Boechera stricta is a perennial relative of Arabidopsis thaliana native to largely undisturbed habitats with two geographic and ecologically divergent subspecies. Here, we delineate the evolutionary processes driving the genetic diversity and population differentiation in this species.

RESULTS

Using whole-genome re-sequencing data from 517 B. stricta accessions, we identify four genetic groups that diverged around 30-180 thousand years ago, with long-term small effective population sizes and recent population expansion after the Last Glacial Maximum. We find three genomic regions with elevated nucleotide diversity, totaling about 10% of the genome. These three regions of elevated nucleotide diversity show excess of intermediate-frequency alleles, higher absolute divergence (dXY), and lower relative divergence (FST) than genomic background, and significant enrichment in immune-related genes, reflecting long-term balancing selection. Scattered across the genome, we also find regions with both high FST and dXY among the groups, termed FST-islands. Population genetic signatures indicate that FST-islands with elevated divergence, which have experienced directional selection, are derived from divergent sorting of ancient polymorphisms.

CONCLUSIONS

Our results suggest that long-term balancing selection on disease resistance genes may have maintained ancestral haplotypes across different geographical lineages, and unequal sorting of balanced polymorphisms may have generated genomic regions with elevated divergence. This study highlights the importance of ancestral balanced polymorphisms as crucial components of genome-wide variation.

A Dead Gene Walking: Convergent Degeneration of a Clade of MADS-Box Genes in Crucifers

Genes are "born," and eventually they "die." These processes shape the phenotypic evolution of organisms and are hence of great biological interest. If genes die in plants, they generally do so quite rapidly. Here, we describe the fate of GOA-like genes that evolve in a dramatically different manner. GOA-like genes belong to the subfamily of Bsister genes of MIKC-type MADS-box genes. Typical MIKC-type genes encode conserved transcription factors controlling plant development. We show that ABS-like genes, a clade of Bsister genes, are indeed highly conserved in crucifers (Brassicaceae) maintaining the ancestral function of Bsister genes in ovule and seed development. In contrast, their closest paralogs, the GOA-like genes, have been undergoing convergent gene death in Brassicaceae. Intriguingly, erosion of GOA-like genes occurred after millions of years of coexistence with ABS-like genes. We thus describe Delayed Convergent Asymmetric Degeneration, a so far neglected but possibly frequent pattern of duplicate gene evolution that does not fit classical scenarios. Delayed Convergent Asymmetric Degeneration of GOA-like genes may have been initiated by a reduction in the expression of an ancestral GOA-like gene in the stem group of Brassicaceae and driven by dosage subfunctionalization. Our findings have profound implications for gene annotations in genomics, interpreting patterns of gene evolution and using genes in phylogeny reconstructions of species.

Genetic Interactions and Molecular Evolution of the Duplicated Genes ICARUS2 and ICARUS1 Help Arabidopsis Plants Adapt to Different Ambient Temperatures

Understanding how plants adapt to ambient temperatures has become a major challenge prompted by global climate change. This has led to the identification of several genes regulating the thermal plasticity of plant growth and flowering time. However, the mechanisms accounting for the natural variation and evolution of such developmental plasticity remain mostly unknown. In this study, we determined that natural variation at ICARUS2 (ICA2), which interacts genetically with its homolog ICA1, alters growth and flowering time plasticity in relation to temperature in Arabidopsis (Arabidopsis thaliana). Transgenic analyses demonstrated multiple functional effects for ICA2 and supported the notion that structural polymorphisms in ICA2 likely underlie its natural variation. Two major ICA2 haplogroups carrying distinct functionally active alleles showed high frequency, strong geographic structure, and significant associations with climatic variables related to annual and daily fluctuations in temperature. Genome analyses across the plant phylogeny indicated that the prevalent plant ICA genes encoding two tRNAHis guanylyl transferase 1 units evolved ∼120 million years ago during the early divergence of mono- and dicotyledonous clades. In addition, ICA1/ICA2 duplication occurred specifically in the Camelineae tribe (Brassicaceae). Thus, ICA2 appears to be ubiquitous across plant evolution and likely contributes to climate adaptation through modifications of thermal developmental plasticity in Arabidopsis.

Expansion and Functional Divergence of Inositol Polyphosphate 5-Phosphatases in Angiosperms

Inositol polyphosphate 5-phosphatase (5PTase), a key enzyme that hydrolyzes the 5` position of the inositol ring, has essential functions in growth, development, and stress responses in plants, yeasts, and animals. However, the evolutionary history and patterns of 5PTases have not been examined systematically. Here, we report a comprehensive molecular evolutionary analysis of the 5PTase gene family and define four groups. These four groups are different from former classifications, which were based on in vitro substrate specificity. Most orthologous groups appear to be conserved as single or low-copy genes in all lineages in Groups II-IV, whereas 5PTase genes in Group I underwent several duplication events in angiosperm, resulting in multiple gene copies. Whole-genome duplication (WGD) was the main mechanism for 5PTase duplications in angiosperm. Plant 5PTases have more members than that of animals, and most plant 5PTase genes appear to have evolved under strong purifying selection. The paralogs have diverged in substrate specificity and expression pattern, showing evidence of selection pressure. Meanwhile, the increase in 5PTases and divergences in sequence, expression, and substrate might have contributed to the divergent functions of 5PTase genes, allowing the angiosperms to successfully adapt to a great number of ecological niches.

Brassicaceae flowers: diversity amid uniformity

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

Origin of angiosperms and the puzzle of the Jurassic gap

Angiosperms are by far the most species-rich clade of land plants, but their origin and early evolutionary history remain poorly understood. We reconstructed angiosperm phylogeny based on 80 genes from 2,881 plastid genomes representing 85% of extant families and all orders. With a well-resolved plastid tree and 62 fossil calibrations, we dated the origin of the crown angiosperms to the Upper Triassic, with major angiosperm radiations occurring in the Jurassic and Lower Cretaceous. This estimated crown age is substantially earlier than that of unequivocal angiosperm fossils, and the difference is here termed the 'Jurassic angiosperm gap'. Our time-calibrated plastid phylogenomic tree provides a highly relevant framework for future comparative studies of flowering plant evolution.

Resolving the backbone of the Brassicaceae phylogeny for investigating trait diversity

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

The Boechera Genus as a Resource for Apomixis Research

The genera Boechera (A. Löve et D. Löve) and Arabidopsis, the latter containing the model plant Arabidopsis thaliana, belong to the same clade within the Brassicaceae family. Boechera is the only among the more than 370 genera in the Brassicaceae where apomixis is well documented. Apomixis refers to the asexual reproduction through seed, and a better understanding of the underlying mechanisms has great potential for applications in agriculture. The Boechera genus currently includes 110 species (of which 38 are reported to be triploid and thus apomictic), which are distributed mostly in the North America. The apomictic lineages of Boechera occur at both the diploid and triploid level and show signs of a hybridogenic origin, resulting in a modification of their chromosome structure, as reflected by alloploidy, aneuploidy, substitutions of homeologous chromosomes, and the presence of aberrant chromosomes. In this review, we discuss the advantages of the Boechera genus to study apomixis, consider its modes of reproduction as well as the inheritance and possible mechanisms controlling apomixis. We also consider population genetic aspects and a possible role of hybridization at the origin of apomixis in Boechera. The molecular tools available to study Boechera, such as transformation techniques, laser capture microdissection, analysis of transcriptomes etc. are also discussed. We survey available genome assemblies of Boechera spp. and point out the challenges to assemble the highly heterozygous genomes of apomictic species. Due to these challenges, we argue for the application of an alternative reference-free method for the comparative analysis of such genomes, provide an overview of genomic sequencing data in the genus Boechera suitable for such analysis, and provide examples of its application.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Angiosperm-Wide and Family-Level Analyses of AP2/ERF Genes Reveal Differential Retention and Sequence Divergence After Whole-Genome Duplication

Plants are immobile and often face stressful environmental conditions, prompting the evolution of genes regulating environmental responses. Such evolution is achieved largely through gene duplication and subsequent divergence. One of the most important gene families involved in regulating plant environmental responses and development is the AP2/ERF superfamily; however, the evolutionary history of these genes is unclear across angiosperms and in major angiosperm families adapted to various ecological niches. Specifically, the impact on gene copy number of whole-genome duplication events occurring around the time of the origins of several plant families is unknown. Here, we present the first angiosperm-wide comparative study of AP2/ERF genes, identifying 75 Angiosperm OrthoGroups (AOGs), each derived from an ancestral angiosperm gene copy. Among these AOGs, 21 retain duplicates with increased copy number in many angiosperm lineages, while the remaining 54 AOGs tend to maintain low copy number. Further analyses of multiple species in the Brassicaceae family indicated that family-specific duplicates experienced differential selective pressures in coding regions, with some paralogs showing signs of positive selection. Further, cis regulatory elements also exhibit extensive divergence between duplicates in Arabidopsis. Moreover, comparison of expression levels suggested that AP2/ERF genes with frequently retained duplicates are enriched for broad expression patterns, offering increased opportunities for functional diversification via changes in expression patterns, and providing a mechanism for repeated duplicate retention in some AOGs. Our results represent the most comprehensive evolutionary history of the AP2/ERF gene family, and support the hypothesis that AP2/ERF genes with broader expression patterns are more likely to be retained as duplicates than those with narrower expression profiles, which could lead to a higher chance of duplicate gene subfunctionalization. The greater tendency of some AOGs to retain duplicates, allowing expression and functional divergence, may facilitate the evolution of complex signaling networks in response to new environmental conditions.

Population history provides foundational knowledge for utilizing and developing native plant restoration materials

A species' population structure and history are critical pieces of information that can help guide the use of available native plant materials in restoration treatments and decide what new native plant materials should be developed to meet future restoration needs. In the western United States, Pseudoroegneria spicata (bluebunch wheatgrass; Poaceae) is an important component of grassland and shrubland plant communities and commonly used for restoration due to its drought resistance and competitiveness with exotic weeds. We used next-generation sequencing data to investigate the processes that shaped P. spicata's geographic pattern of genetic variation across the Intermountain West. Pseudoroegneria spicata's genetic diversity is partitioned into populations that likely differentiated since the Last Glacial Maximum. Adjacent populations display varying magnitudes of historical gene flow, with migration rates ranging from multiple migrants per generation to multiple generations per migrant. When considering the commercial germplasm sources available for restoration, genetic identities remain representative of the wildland localities from which germplasm sources were originally developed, and they maintain high levels of heterozygosity and nucleotide diversity. However, the commercial germplasm sources represent a small fraction of the overall genetic diversity of P. spicata in the Intermountain West. Given the low migration rates and long divergence times between some pairs of P. spicata populations, using commercial germplasm sources could facilitate undesirable restoration outcomes when used in certain geographic areas, even if the environment in which the commercial materials thrive is similar to that of the restoration site. As such, population structure and history can be used to provide guidance on what geographic areas may need additional native plant materials so that restoration efforts support species and community resilience and improve outcomes.

Salt tolerance mechanisms in three Irano-Turanian Brassicaceae halophytes relatives of Arabidopsis thaliana

Salt tolerance mechanisms were studied in three Irano-Turanian halophytic species from the Brassicaceae ‎‎(Lepidium latifolium, L. perfoliatum and Schrenkiella parvula) and compared with the glycophyte Arabidopsis thaliana. According to seed germination under salt stress, L. perfoliatum was the most tolerant species, while L. latifolium and S. parvula were rather susceptible. Contrastingly, based on biomass production L. perfoliatum was more salt sensitive than the other two species. In S. parvula biomass was increased up to 2.8-fold by 100 mM NaCl; no significant growth reduction was observed even when exposed to 400 mM NaCl. Stable activities of antioxidative defense enzymes, nil or negligible accumulation of superoxide anion and hydrogen peroxide, as well as stable membrane integrity in the three halophytes revealed that no oxidative stress occurred in these tolerant species under salt stress. Proline levels increased in response to salt treatment. However, it contributed only by 0.3‒2.0% to the total osmolyte concentration in the three halophytes (at 400 mM NaCl) and even less (0.04%) in the glycophyte, A. thaliana (at 100 mM NaCl). Soluble sugars in all three halophytes and free amino acids pool in S. parvula decreased under salt treatment in contrast to the glycophyte, A. thaliana. The contribution of organic osmolytes to the total osmolyte pool increased by salt treatment in the roots, while decreased in halophyte and glycophyte, A. thaliana leaves. Interestingly, this reduction was compensated by a higher relative contribution of K in the leaves of the halophytes, but of Na in A. thaliana. Taken together, biomass data and biochemical indicators show that S. parvula is more salt tolerant than the two Lepidium species. Our data indicate that L. latifolium, as a perennial halophyte with a large biomass, is highly suitable for both restoration of saline habitats and saline agriculture.

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

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

Molecular genetic and biochemical evidence for adaptive evolution of leaf abaxial epicuticular wax crystals in the genus Lithocarpus (Fagaceae)

BACKGROUND

Leaf epicuticular wax is an important functional trait for physiological regulation and pathogen defense. This study tests how selective pressure may have forced the trait of leaf abaxial epicuticular wax crystals (LAEWC) and whether the presence/absence of LAEWC is associated with other ecophysiological traits. Scanning Electron Microscopy was conducted to check for LAEWC in different Lithocarpus species. Four wax biosynthesis related genes, including two wax backbone genes ECERIFERUM 1 (CER1) and CER3, one regulatory gene CER7 and one transport gene CER5, were cloned and sequenced. Ecophysiological measurements of secondary metabolites, photosynthesis, water usage efficiency, and nutrition indices were also determined. Evolutionary hypotheses of leaf wax character transition associated with the evolution of those ecophysiological traits as well as species evolution were tested by maximum likelihood.

RESULTS

Eight of 14 studied Lithocarpus species have obvious LAEWC appearing with various types of trichomes. Measurements of ecophysiological traits show no direct correlations with the presence/absence of LAEWC. However, the content of phenolic acids is significantly associated with the gene evolution of the wax biosynthetic backbone gene CER1, which was detected to be positively selected when LAEWC was gained during the late-Miocene-to-Pliocene period.

CONCLUSIONS

Changes of landmass and vegetation type accelerated the diversification of tropical and subtropical forest trees and certain herbivores during the late Miocene. As phenolic acids were long thought to be associated with defense against herbivories, co-occurrence of LAEWC and phenolic acids may suggest that LAEWC might be an adaptive defensive mechanism in Lithocarpus.

Species Tree Estimation Using ASTRAL: How Many Genes Are Enough?

Species tree reconstruction from genomic data is increasingly performed using methods that account for sources of gene tree discordance such as incomplete lineage sorting. One popular method for reconstructing species trees from unrooted gene tree topologies is ASTRAL. In this paper, we derive theoretical sample complexity results for the number of genes required by ASTRAL to guarantee reconstruction of the correct species tree with high probability. We also validate those theoretical bounds in a simulation study. Our results indicate that ASTRAL requires gene trees to reconstruct the species tree correctly with high probability where is the number of species and is the length of the shortest branch in the species tree. Our simulations, some under the anomaly zone, show trends consistent with the theoretical bounds and also provide some practical insights on the conditions where ASTRAL works well.

Climate and host-plant associations shaped the evolution of ceutorhynch weevils throughout the Cenozoic

Using molecular phylogenetic data and methods we inferred divergence times and diversification patterns for the weevil subfamily Ceutorhynchinae in the context of host‐plant associations and global climate over evolutionary time. We detected four major diversification shifts that correlate with both host shifts and major climate events. Ceutorhynchinae experienced an increase in diversification rate at ∼53 Ma, during the Early Eocene Climate Optimum, coincident with a host shift to Lamiaceae. A second major diversification phase occurred at the end of the Eocene (∼34 Ma). This contrasts with the overall deterioration in climate equability at the Eocene‐Oligocene boundary, but tracks the diversification of important host plant clades in temperate (higher) latitudes, leading to increased diversification rates in the weevil clades infesting temperate hosts. A third major phase of diversification is correlated with the rising temperatures of the Late Oligocene Warming Event (∼26.5 Ma); diversification rates then declined shortly after the Middle Miocene Climate Transition (∼14.9 Ma). Our results indicate that biotic and abiotic factors together explain the evolution of Ceutorhynchinae better than each of these drivers viewed in isolation.

The Aphid-Transmitted Turnip yellows virus Differentially Affects Volatiles Emission and Subsequent Vector Behavior in Two Brassicaceae Plants

Aphids are important pests which cause direct damage by feeding or indirect prejudice by transmitting plant viruses. Viruses are known to induce modifications of plant cues in ways that can alter vector behavior and virus transmission. In this work, we addressed whether the modifications induced by the aphid-transmitted Turnip yellows virus (TuYV) in the model plant Arabidopsis thaliana also apply to the cultivated plant Camelina sativa, both belonging to the Brassicaceae family. In most experiments, we observed a significant increase in the relative emission of volatiles from TuYV-infected plants. Moreover, due to plant size, the global amounts of volatiles emitted by C. sativa were higher than those released by A. thaliana. In addition, the volatiles released by TuYV-infected C. sativa attracted the TuYV vector Myzus persicae more efficiently than those emitted by non-infected plants. In contrast, no such preference was observed for A. thaliana. We propose that high amounts of volatiles rather than specific metabolites are responsible for aphid attraction to infected C. sativa. This study points out that the data obtained from the model pathosystem A. thaliana/TuYV cannot be straightforwardly extrapolated to a related plant species infected with the same virus.

Plant phylogenomics based on genome-partitioning strategies: Progress and prospects

The rapid expansion of next-generation sequencing (NGS) has generated a powerful array of approaches to address fundamental questions in biology. Several genome-partitioning strategies to sequence selected subsets of the genome have emerged in the fields of phylogenomics and evolutionary genomics. In this review, we summarize the applications, advantages and limitations of four NGS-based genome-partitioning approaches in plant phylogenomics: genome skimming, transcriptome sequencing (RNA-seq), restriction site associated DNA sequencing (RAD-Seq), and targeted capture (Hyb-seq). Of these four genome-partitioning approaches, targeted capture (especially Hyb-seq) shows the greatest promise for plant phylogenetics over the next few years. This review will aid researchers in their selection of appropriate genome-partitioning approaches to address questions of evolutionary scale, where we anticipate continued development and expansion of whole-genome sequencing strategies in the fields of plant phylogenomics and evolutionary biology research.

Fragmentary Gene Sequences Negatively Impact Gene Tree and Species Tree Reconstruction

Species tree reconstruction from genome-wide data is increasingly being attempted, in most cases using a two-step approach of first estimating individual gene trees and then summarizing them to obtain a species tree. The accuracy of this approach, which promises to account for gene tree discordance, depends on the quality of the inferred gene trees. At the same time, phylogenomic and phylotranscriptomic analyses typically use involved bioinformatics pipelines for data preparation. Errors and shortcomings resulting from these preprocessing steps may impact the species tree analyses at the other end of the pipeline. In this article, we first show that the presence of fragmentary data for some species in a gene alignment, as often seen on real data, can result in substantial deterioration of gene trees, and as a result, the species tree. We then investigate a simple filtering strategy where individual fragmentary sequences are removed from individual genes but the rest of the gene is retained. Both in simulations and by reanalyzing a large insect phylotranscriptomic data set, we show the effectiveness of this simple filtering strategy.

Aberrant imprinting may underlie evolution of parthenogenesis

Genomic imprinting confers parent-of-origin-specific gene expression, thus non-equivalent and complementary function of parental genomes. As a consequence, genomic imprinting poses an epigenetic barrier to parthenogenesis in sexual organisms. We report aberrant imprinting in Boechera, a genus in which apomicts evolved from sexuals multiple times. Maternal activation of a MADS-box gene, a homolog of which is imprinted and paternally expressed in the sexual relative Arabidopsis, is accompanied by locus-specific DNA methylation changes in apomicts where parental imprinting seems to be relaxed.

Aberrant imprinting may underlie evolution of parthenogenesis

Genomic imprinting confers parent-of-origin-specific gene expression, thus non-equivalent and complementary function of parental genomes. As a consequence, genomic imprinting poses an epigenetic barrier to parthenogenesis in sexual organisms. We report aberrant imprinting in Boechera, a genus in which apomicts evolved from sexuals multiple times. Maternal activation of a MADS-box gene, a homolog of which is imprinted and paternally expressed in the sexual relative Arabidopsis, is accompanied by locus-specific DNA methylation changes in apomicts where parental imprinting seems to be relaxed.

Panax ginseng genome examination for ginsenoside biosynthesis

Ginseng, which contains ginsenosides as bioactive compounds, has been regarded as an important traditional medicine for several millennia. However, the genetic background of ginseng remains poorly understood, partly because of the plant's large and complex genome composition. We report the entire genome sequence of Panax ginseng using next-generation sequencing. The 3.5-Gb nucleotide sequence contains more than 60% repeats and encodes 42 006 predicted genes. Twenty-two transcriptome datasets and mass spectrometry images of ginseng roots were adopted to precisely quantify the functional genes. Thirty-one genes were identified to be involved in the mevalonic acid pathway. Eight of these genes were annotated as 3-hydroxy-3-methylglutaryl-CoA reductases, which displayed diverse structures and expression characteristics. A total of 225 UDP-glycosyltransferases (UGTs) were identified, and these UGTs accounted for one of the largest gene families of ginseng. Tandem repeats contributed to the duplication and divergence of UGTs. Molecular modeling of UGTs in the 71st, 74th, and 94th families revealed a regiospecific conserved motif located at the N-terminus. Molecular docking predicted that this motif captures ginsenoside precursors. The ginseng genome represents a valuable resource for understanding and improving the breeding, cultivation, and synthesis biology of this key herb.

A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families

Ferns account for 80% of nonflowering vascular plant species and are the sister lineage of seed plants. Recent molecular phylogenetics have greatly advanced understanding of fern tree of life, but relationships among some major lineages remain unclear. To better resolve the phylogenetic relationships of ferns, we generated transcriptomes from 125 ferns and two lycophytes, with three additional public datasets, to represent all 11 orders and 85% of families of ferns. Our nuclear phylogeny provides strong supports for the monophyly of all four subclasses and nearly all orders and families, and for relationships among these lineages. The only exception is Gleicheniales, which was highly supported as being paraphyletic with Dipteridaceae sister to a clade with Gleicheniaceae + Hymenophyllales. In addition, new and strongly supported phylogenetic relationships are found for suborders and families in Polypodiales. We provide the first dated fern phylogenomic tree using many nuclear genes from a large majority of families, with an estimate for separation of the ancestors of ferns and seed plants in early Devonian at ∼400 Mya and subsequent gradual divergences of fern orders from ∼380 to 200 Mya. Moreover, the newly obtained fern phylogeny provides a framework for gene family analyses, which indicate that the vast majority of transcription factor families found in seed plants were already present in the common ancestor of extant vascular plants. In addition, fern transcription factor genes show similar duplication patterns to those in seed plants, with some showing stable copy number and others displaying independent expansions in both ferns and seed plants. This study provides a robust phylogenetic and gene family evolution framework, as well as rich molecular resources for understanding the morphological and functional evolution in ferns.

New insights into the taxonomy of tribe Euclidieae (Brassicaceae), evidence from nrITS sequence data

As currently delimitated, the species-rich mustard tribe Euclidieae DC. (Brassicaceae) comprises 28 genera and 152 species distributed primarily in Asia. To date, no tribe-wide comprehensive phylogenetic analysis has been conducted. In this study, sequence data from the nuclear ribosomal internal transcribed spacer (nrITS) region of 82 species in all 28 genera of Euclidieae were used to test its monophyly and infer inter- and intra-generic relationships within. Phylogenetic analyses revealed that Rhammatophyllum and Sisymbriopsis are embedded within Solms-laubachia s.l., and Solms-laubachia lanuginosa (Eurycarpus lanuginosus) fell outside the tribe. Therefore, Solms-laubachia s.l. as currently recognized is not monophyletic and its generic delimitation needed further study. Besides, our results suggest that the genera Lepidostemon, Neotorularia, and Tetracme are polyphyletic.

Identification of genes regulating traits targeted for domestication of field cress (Lepidium campestre) as a biennial and perennial oilseed crop

BACKGROUND

The changing climate and the desire to use renewable oil sources necessitate the development of new oilseed crops. Field cress (Lepidium campestre) is a species in the Brassicaceae family that has been targeted for domestication not only as an oilseed crop that produces seeds with a desirable industrial oil quality but also as a cover/catch crop that provides valuable ecosystem services. Lepidium is closely related to Arabidopsis and display significant proportions of syntenic regions in their genomes. Arabidopsis genes are among the most characterized genes in the plant kingdom and, hence, comparative genomics of Lepidium-Arabidopsis would facilitate the identification of Lepidium candidate genes regulating various desirable traits.

RESULTS

Homologues of 30 genes known to regulate vernalization, flowering time, pod shattering, oil content and quality in Arabidopsis were identified and partially characterized in Lepidium. Alignments of sequences representing field cress and two of its closely related perennial relatives: L. heterophyllum and L. hirtum revealed 243 polymorphic sites across the partial sequences of the 30 genes, of which 95 were within the predicted coding regions and 40 led to a change in amino acids of the target proteins. Within field cress, 34 polymorphic sites including nine non-synonymous substitutions were identified. The phylogenetic analysis of the data revealed that field cress is more closely related to L. heterophyllum than to L. hirtum.

CONCLUSIONS

There is significant variation within and among Lepidium species within partial sequences of the 30 genes known to regulate traits targeted in the present study. The variation within these genes are potentially useful to speed-up the process of domesticating field cress as future oil crop. The phylogenetic relationship between the Lepidium species revealed in this study does not only shed some light on Lepidium genome evolution but also provides important information to develop efficient schemes for interspecific hybridization between different Lepidium species as part of the domestication efforts.

Phytochrome B regulates resource allocation in Brassica rapa

Crop biomass and yield are tightly linked to how the light signaling network translates information about the environment into allocation of resources, including photosynthates. Once activated, the phytochrome (phy) class of photoreceptors signal and re-deploy carbon resources to alter growth, plant architecture, and reproductive timing. Most of the previous characterization of the light-modulated growth program has been performed in the reference plant Arabidopsis thaliana. Here, we use Brassica rapa as a crop model to test for conservation of the phytochrome-carbon network. In response to elevated levels of CO2, B. rapa seedlings showed increases in hypocotyl length, shoot and root fresh weight, and the number of lateral roots. All of these responses were dependent on nitrogen and polar auxin transport. In addition, we identified putative B. rapa orthologs of PhyB and isolated two nonsense alleles. BrphyB mutants had significantly decreased or absent CO2-stimulated growth responses. Mutant seedlings also showed misregulation of auxin-dependent genes and genes involved in chloroplast development. Adult mutant plants had reduced chlorophyll levels, photosynthetic rate, stomatal index, and seed yield. These findings support a recently proposed holistic role for phytochromes in regulating resource allocation, biomass production, and metabolic state in the developing plant.

Kale (Brassica oleracea var. acephala) as a superfood: Review of the scientific evidence behind the statement

Kale (Brassica oleracea var. acephala) is a cruciferous vegetable, characterized by leaves along the stem, which, in recent years, have gained a great popularity as a ´superfood´. Consequently, in a popular culture it is listed in many ´lists of the healthiest vegetables´. Without the doubt, a scientific evidences support the fact that cruciferous vegetables included in human diet can positively affect health and well-being, but remains unclear why kale is declared superior in comparison with other cruciferous. It is questionable if this statement about kale is triggered by scientific evidence or by some other factors. Our review aims to bring an overview of kale's botanical characteristics, agronomic requirements, contemporary and traditional use, macronutrient and phytochemical content and biological activity, in order to point out the reasons for tremendous kale popularity.

A segmental duplication in the common ancestor of Brassicaceae is responsible for the origin of the paralogs KCS6-KCS5, which are not shared with other angiosperms

Novel morphological structures allowed adaptation to dry conditions in early land plants. The cuticle, one such novelty, plays diverse roles in tolerance to abiotic and biotic stresses and plant development. Cuticular waxes represent a major constituent of the cuticle and are comprised of an assortment of chemicals that include, among others, very long chain fatty acids (VLCFAs). Members of the β-ketoacyl coenzyme A synthases (KCS) gene family code for enzymes that are essential for fatty acid biosynthesis. The gene KCS6 (CUT1) is known to be a key player in the production of VLCFA precursors essential for the synthesis of cuticular waxes in the model plant Arabidopsis thaliana (Brassicaceae). Despite its functional importance, relatively little is known about the evolutionary history of KCS6 or its paralog KCS5 in Brassicaceae or beyond. This lacuna becomes important when we extrapolate understanding of mechanisms gained from the model plant to its containing clades Brassicaceae, flowering plants, or beyond. The Brassicaceae, with several sequenced genomes and a known history of paleoploidy, mesopolyploidy and neopolyploidy, offer a system in which to study the evolution and diversification of the KCS6-KCS5 paralogy. Our phylogenetic analyses across green plants, combined with comparative genomic, microsynteny and evolutionary rates analyses across nine genomes of Brassicaceae, reveal that (1) the KCS6-KCS5 paralogy arose as the result of a large segmental duplication in the ancestral Brassicaceae, (2) the KCS6-KCS5 lineage is represented by a single copy in other flowering plant lineages, (3) the duplicated segments undergo different degrees of retention and loss, and (4) most of the genes in the KCS6 and KCS5 gene blocks (including KCS6 and KCS5 themselves) are under purifying selection. The last also true for most members of the KCS gene family in Brassicaceae, except for KCS8, KCS9 and KCS17, which are under positive selection and may be undergoing functional evolution, meriting further investigation. Overall, our results clearly establish that the ancestral KCS6/5 gene duplicated in the Brassicaceae lineage. It is possible that any specialized functions of KCS5 found in Brassicaceae are either part of a set of KCS6/5 gene functions in the rest of the flowering plants, or unique to Brassicaceae.

The complete chloroplast genome of Sisymbrium irio

The complete chloroplast genome of Sisymbrium irio was determined. The length of the complete chloroplast genome is 154,001 bp. The whole chloroplast genome consists of 83,891 bp long single copy (LSC) and 17,630 bp small single copy (SSC) regions, separated by a pair of 26,240 bp inverted repeat (IR) regions. The S. irio chloroplast genome encodes 112 annotated known unique genes including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. The phylogenetic position of S. irio is sister to Brassiceae and Thlaspideae.

Preferential retention of genes from one parental genome after polyploidy illustrates the nature and scope of the genomic conflicts induced by hybridization

Polyploidy is increasingly seen as a driver of both evolutionary innovation and ecological success. One source of polyploid organisms' successes may be their origins in the merging and mixing of genomes from two different species (e.g., allopolyploidy). Using POInT (the Polyploid Orthology Inference Tool), we model the resolution of three allopolyploidy events, one from the bakers' yeast (Saccharomyces cerevisiae), one from the thale cress (Arabidopsis thaliana) and one from grasses including Sorghum bicolor. Analyzing a total of 21 genomes, we assign to every gene a probability for having come from each parental subgenome (i.e., derived from the diploid progenitor species), yielding orthologous segments across all genomes. Our model detects statistically robust evidence for the existence of biased fractionation in all three lineages, whereby genes from one of the two subgenomes were more likely to be lost than those from the other subgenome. We further find that a driver of this pattern of biased losses is the co-retention of genes from the same parental genome that share functional interactions. The pattern of biased fractionation after the Arabidopsis and grass allopolyploid events was surprisingly constant in time, with the same parental genome favored throughout the lineages' history. In strong contrast, the yeast allopolyploid event shows evidence of biased fractionation only immediately after the event, with balanced gene losses more recently. The rapid loss of functionally associated genes from a single subgenome is difficult to reconcile with the action of genetic drift and suggests that selection may favor the removal of specific duplicates. Coupled to the evidence for continuing, functionally-associated biased fractionation after the A. thaliana At-α event, we suggest that, after allopolyploidy, there are functional conflicts between interacting genes encoded in different subgenomes that are ultimately resolved through preferential duplicate loss.

Widespread Whole Genome Duplications Contribute to Genome Complexity and Species Diversity in Angiosperms

Gene duplications provide evolutionary potentials for generating novel functions, while polyploidization or whole genome duplication (WGD) doubles the chromosomes initially and results in hundreds to thousands of retained duplicates. WGDs are strongly supported by evidence commonly found in many species-rich lineages of eukaryotes, and thus are considered as a major driving force in species diversification. We performed comparative genomic and phylogenomic analyses of 59 public genomes/transcriptomes and 46 newly sequenced transcriptomes covering major lineages of angiosperms to detect large-scale gene duplication events by surveying tens of thousands of gene family trees. These analyses confirmed most of the previously reported WGDs and provided strong evidence for novel ones in many lineages. The detected WGDs supported a model of exponential gene loss during evolution with an estimated half-life of approximately 21.6 million years, and were correlated with both the emergence of lineages with high degrees of diversification and periods of global climate changes. The new datasets and analyses detected many novel WGDs widely spread during angiosperm evolution, uncovered preferential retention of gene functions in essential cellular metabolisms, and provided clues for the roles of WGD in promoting angiosperm radiation and enhancing their adaptation to environmental changes.

Brassicales phylogeny inferred from 72 plastid genes: A reanalysis of the phylogenetic localization of two paleopolyploid events and origin of novel chemical defenses

PREMISE OF THE STUDY

Previous phylogenetic studies employing molecular markers have yielded various insights into the evolutionary history across Brassicales, but many relationships between families remain poorly supported or unresolved. A recent phylotranscriptomic approach utilizing 1155 nuclear markers obtained robust estimates for relationships among 14 of 17 families. Here we report a complete family-level phylogeny estimated using the plastid genome.

METHODS

We conducted phylogenetic analyses on a concatenated data set comprising 44,926 bp from 72 plastid genes for species distributed across all 17 families. Our analysis includes three additional families, Tovariaceae, Salvadoraceae, and Setchellanthaceae, that were omitted in the previous phylotranscriptomic study.

KEY RESULTS

Our phylogenetic analyses obtained fully resolved and strongly supported estimates for all nodes across Brassicales. Importantly, these findings are congruent with the topology reported in the phylotranscriptomic study. This consistency suggests that future studies could utilize plastid genomes as markers for resolving relationships within some notoriously difficult clades across Brassicales. We used this new phylogenetic framework to verify the placement of the At-α event near the origin of Brassicaceae, with median date estimates of 31.8 to 42.8 million years ago and restrict the At-β event to one of two nodes with median date estimates between 85 to 92.2 million years ago. These events ultimately gave rise to novel chemical defenses and are associated with subsequent shifts in net diversification rates.

CONCLUSIONS

We anticipate that these findings will aid future comparative evolutionary studies across Brassicales, including selecting candidates for whole-genome sequencing projects.

Silique valves as sails in anemochory of Lunaria (Brassicaceae)

The generally held opinion that seeds of Lunaria remain at the replum after detachment of the two valves and then wind causes a shaking or rattling of the replum with its diaphragm, thus launching the seeds, is challenged. In a sparse forest in the Swabian Alb, the first author noticed flying valves of Lunaria rediviva to which the narrow-winged flat seeds are attached. Investigations with SEM and histology have shown that the valves secrete a glue only at those sites where the seeds rest on the valves before valve tissues die. Further analysis has shown (using the periodic acid-Schiff reaction) that the glue consists of polysaccharides. After detachment and dispersal of the valves, the adhesive strength continuously decreases. This is the first report for a sticky valve exudate in the Brassicaceae. Because of the adhesion of Lunaria seeds to their valves for some time, the 1st order diaspore is a mericarp, in a broad sense, and can be interpreted as an adaptation to long-distance dispersal by stronger winds. In this context, the 'flying carpets' of Lunaria are more effective and transport more than one seed. Molecular studies assigned Lunaria to the tribe Biscutelleae, which now contains the angustiseptate genera Biscutella and Megadenia as well as the latiseptate genera Lunaria and Ricotia. The valves in Ricotia can easily be detached (studied in herbarium material and a living plant), but, in contrast to Lunaria, the ripe seeds remain at the replum and its diaphragm, respectively.