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

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

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

Chromosomal Evolution and Apomixis in the Cruciferous Tribe Boechereae

The mustard family (Brassicaceae) comprises several dozen monophyletic clades usually ranked as tribes. The tribe Boechereae plays a prominent role in plant research due to the incidence of apomixis and its close relationship to Arabidopsis. This tribe, largely confined to western North America, harbors nine genera and c. 130 species, with >90% of species belonging to the genus Boechera. Hundreds of apomictic diploid and triploid Boechera hybrids have spurred interest in this genus, but the remaining Boechereae genomes remain virtually unstudied. Here we report on comparative genome structure of six genera (Borodinia, Cusickiella, Phoenicaulis, Polyctenium, Nevada, and Sandbergia) and three Boechera species as revealed by comparative chromosome painting (CCP). All analyzed taxa shared the same seven-chromosome genome structure. Comparisons with the sister Halimolobeae tribe (n = 8) showed that the ancestral Boechereae genome (n = 7) was derived from an older n = 8 genome by descending dysploidy followed by the divergence of extant Boechereae taxa. As tribal divergence post-dated the origin of four tribe-specific chromosomes, it is proposed that these chromosomal rearrangements were a key evolutionary innovation underlaying the origin and diversification of the Boechereae in North America. Although most Boechereae genera exhibit genomic conservatism, intra-tribal cladogenesis has occasionally been accompanied by chromosomal rearrangements (particularly inversions). Recently, apomixis was reported in the Boechereae genera Borodinia and Phoenicaulis. Here, we report sexual reproduction in diploid Nevada, diploid Sandbergia, and tetraploid Cusickiella and aposporous apomixis in tetraploids of Polyctenium and Sandbergia. In sum, apomixis is now known to occur in five of the nine Boechereae genera.

Developmental Mechanisms of Fleshy Fruit Diversity in Rosaceae

Rosaceae (the rose family) is an economically important family that includes species prized for high-value fruits and ornamentals. The family also exhibits diverse fruit types, including drupe (peach), pome (apple), drupetum (raspberry), and achenetum (strawberry). Phylogenetic analysis and ancestral fruit-type reconstruction suggest independent evolutionary paths of multiple fleshy fruit types from dry fruits. A recent whole genome duplication in the Maleae/Pyreae tribe (with apple, pear, hawthorn, and close relatives; referred to as Maleae here) may have contributed to the evolution of pome fruit. MADS-box genes, known to regulate floral organ identity, are emerging as important regulators of fruit development. The differential competence of floral organs to respond to fertilization signals may explain the different abilities of floral organs to form fleshy fruit. Future comparative genomics and functional studies in closely related Rosaceae species with distinct fruit types will test hypotheses and provide insights into mechanisms of fleshy fruit diversity. These efforts will be facilitated by the wealth of genome data and resources in Rosaceae.

Independent evolution of ancestral and novel defenses in a genus of toxic plants (Erysimum, Brassicaceae)

Phytochemical diversity is thought to result from coevolutionary cycles as specialization in herbivores imposes diversifying selection on plant chemical defenses. Plants in the speciose genus Erysimum (Brassicaceae) produce both ancestral glucosinolates and evolutionarily novel cardenolides as defenses. Here we test macroevolutionary hypotheses on co-expression, co-regulation, and diversification of these potentially redundant defenses across this genus. We sequenced and assembled the genome of E. cheiranthoides and foliar transcriptomes of 47 additional Erysimum species to construct a phylogeny from 9868 orthologous genes, revealing several geographic clades but also high levels of gene discordance. Concentrations, inducibility, and diversity of the two defenses varied independently among species, with no evidence for trade-offs. Closely related, geographically co-occurring species shared similar cardenolide traits, but not glucosinolate traits, likely as a result of specific selective pressures acting on each defense. Ancestral and novel chemical defenses in Erysimum thus appear to provide complementary rather than redundant functions.

Brassicaceae-Derived Anticancer Agents: Towards a Green Approach to Beat Cancer

Cancer is the main cause of mortality and morbidity worldwide. Although a large variety of therapeutic approaches have been developed and translated into clinical protocols, the toxic side effects of cancer treatments negatively impact patients, allowing cancer to grow. Brassica metabolites are emerging as new weapons for anti-cancer therapeutics. The beneficial role of the consumption of brassica vegetables, the most-used vegetables in the Mediterranean diet, particularly broccoli, in the prevention of chronic diseases, including cardiovascular diseases, diabetes, and obesity, has been well-documented. In this review, we discuss the anti-tumor effects of the bioactive compounds from Brassica vegetables with regard to the compounds and types of cancer against which they show activity, providing current knowledge on the anti-cancer effects of Brassica metabolites against major types of tumors. In addition, we discuss the impacts of industrial and domestic processing on the compounds’ functional properties before their consumption as well as the main strategies used to increase the content of health-promoting metabolites in Brassica plants through biofortification. Finally, the impacts of microbiota on the compounds’ bioactivity are considered. This information will be helpful for the further development of efficacious anti-cancer drugs.

Comparative Mitogenome Analysis of the Genus Trifolium Reveals Independent Gene Fission of ccmFn and Intracellular Gene Transfers in Fabaceae

The genus Trifolium is the largest of the tribe Trifolieae in the subfamily Papilionoideae (Fabaceae). The paucity of mitochondrial genome (mitogenome) sequences has hindered comparative analyses among the three genomic compartments of the plant cell (nucleus, mitochondrion and plastid). We assembled four mitogenomes from the two subgenera (Chronosemium and Trifolium) of the genus. The four Trifolium mitogenomes were compact (294,911-348,724 bp in length) and contained limited repetitive (6.6-8.6%) DNA. Comparison of organelle repeat content highlighted the distinct evolutionary trajectory of plastid genomes in a subset of Trifolium species. Intracellular gene transfer (IGT) was analyzed among the three genomic compartments revealing functional transfer of mitochondrial rps1 to nuclear genome along with other IGT events. Phylogenetic analysis based on mitochondrial and nuclear rps1 sequences revealed that the functional transfer in Trifolieae was independent from the event that occurred in robinioid clade that includes genus Lotus. A novel, independent fission event of ccmFn in Trifolium was identified, caused by a 59 bp deletion. Fissions of this gene reported previously in land plants were reassessed and compared with Trifolium.

Phylogenomic Insights into Deep Phylogeny of Angiosperms Based on Broad Nuclear Gene Sampling

Angiosperms (flowering plants) are the most diverse and species-rich group of plants. The vast majority (∼99.95%) of angiosperms form a clade called Mesangiospermae, which is subdivided into five major groups: eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales. The relationships among these Mesangiospermae groups have been the subject of long debate. In this study, we assembled a phylogenomic dataset of 1594 genes from 151 angiosperm taxa, including representatives of all five lineages, to investigate the phylogeny of major angiosperm lineages under both coalescent- and concatenation-based methods. We dissected the phylogenetic signal and found that more than half of the genes lack phylogenetic information for the backbone of angiosperm phylogeny. We further removed the genes with weak phylogenetic signal and showed that eudicots, Ceratophyllales, and Chloranthales form a clade, with magnoliids and monocots being the next successive sister lineages. Similar frequencies of gene tree conflict are suggestive of incomplete lineage sorting along the backbone of the angiosperm phylogeny. Our analyses suggest that a fully bifurcating species tree may not be the best way to represent the early radiation of angiosperms. Meanwhile, we inferred that the crown-group angiosperms originated approximately between 255.1 and 222.2 million years ago, and Mesangiospermae diversified into the five extant groups in a short time span (∼27 million years) at the Early to Late Jurassic.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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

The complete chloroplast genome sequence of yellow mustard (Sinapis alba L.) and its phylogenetic relationship to other Brassicaceae species

As a member of the large Brassicaceae family, yellow mustard (Sinapis alba L.) has been used as an important gene pool for the genetic improvement of cash crops in Brassicaceae. Understanding the phylogenetic relationship between Sinapis alba (S. alba) and other Brassicaceae crops can provide guidance on the introgression of its favorable alleles into related species. The chloroplast (cp) genome is an ideal model for assessing genome evolution and the phylogenetic relationships of complex angiosperm families. Herein, we de novo assembled the complete cp genome of S. alba by integrating the PacBio and Illumina sequencing platforms. A 153,760 bp quadripartite cycle without any gap was obtained, including a pair of inverted repeats (IRa and IRb) of 26,221 bp, separated by a large single copy (LSC) region of 83,506 bp and a small single copy (SSC) region of 17,821 bp. A total of 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified in this cp genome, as were 89 simple sequence repeat (SSR) loci of 18 types. The codon usage analysis revealed a preferential use of the Leu codon with the A/U ending. The phylogenetic analysis using 82 Brassicaceae species demonstrated that S. alba had a close relationship with important Brassica and Raphanus species; moreover, it likely originated from a separate evolutionary pathway compared with the congeneric Sinapis arvensis. The synonymous (Ks) and non-synonymous (Ks) substitution rate analysis showed that genes encoding "Subunits of cytochrome b/f complex" were under the lowest purifying selection pressure, whereas those associated with "Maturase", "Subunit of acetyl-CoA", and "Subunits of NADH-dehydrogenase" underwent relatively higher purifying selection pressures. Our results provide valuable information for fully utilizing the S. alba cp genome as a potential genetic resource for the genetic improvement of Brassica and Raphanus species.

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.