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

Complementing model species with model clades

Model species continue to underpin groundbreaking plant science research. At the same time, the phylogenetic resolution of the land plant tree of life continues to improve. The intersection of these 2 research paths creates a unique opportunity to further extend the usefulness of model species across larger taxonomic groups. Here we promote the utility of the Arabidopsis thaliana model species, especially the ability to connect its genetic and functional resources, to species across the entire Brassicales order. We focus on the utility of using genomics and phylogenomics to bridge the evolution and diversification of several traits across the Brassicales to the resources in Arabidopsis, thereby extending scope from a model species by establishing a "model clade." These Brassicales-wide traits are discussed in the context of both the model species Arabidopsis and the family Brassicaceae. We promote the utility of such a "model clade" and make suggestions for building global networks to support future studies in the model order Brassicales.

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

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

The allotetraploid horseradish genome provides insights into subgenome diversification and formation of critical traits

Polyploidization can provide a wealth of genetic variation for adaptive evolution and speciation, but understanding the mechanisms of subgenome evolution as well as its dynamics and ultimate consequences remains elusive. Here, we report the telomere-to-telomere (T2T) gap-free reference genome of allotetraploid horseradish (Armoracia rusticana) sequenced using a comprehensive strategy. The (epi)genomic architecture and 3D chromatin structure of the A and B subgenomes differ significantly, suggesting that both the dynamics of the dominant long terminal repeat retrotransposons and DNA methylation have played critical roles in subgenome diversification. Investigation of the genetic basis of biosynthesis of glucosinolates (GSLs) and horseradish peroxidases reveals both the important role of polyploidization and subgenome differentiation in shaping the key traits. Continuous duplication and divergence of essential genes of GSL biosynthesis (e.g., FMO_GS-OX, IGMT, and GH1 gene family) contribute to the broad GSL profile in horseradish. Overall, the T2T assembly of the allotetraploid horseradish genome expands our understanding of polyploid genome evolution and provides a fundamental genetic resource for breeding and genetic improvement of horseradish.

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

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

Comparative genomics analysis of endangered wild Egyptian Moringa peregrina (Forssk.) Fiori plastome, with implications for the evolution of Brassicales order

Moringa is a mono-genus belonging to the Moringaceae family, which includes 13 species. Among them, Moringa peregrina is plant species native to the Arabian Peninsula, Southern Sinai in Egypt, and the Horn of Africa, and comprehensive studies on its nutritional, industrial, and medicinal values have been performed. Herein, we sequenced and analyzed the initial complete chloroplast genome of Moringa peregrina. Concurrently, we analyzed the new chloroplast genome along with 25 chloroplast genomes related to species representing eight families in the Brassicales order. The results indicate that the plastome sequence of M. peregrina consists of 131 genes, with an average GC content of 39.23%. There is a disparity in the IR regions of the 26 species ranging from 25,804 to 31,477 bp. Plastome structural variations generated 20 hotspot regions that could be considered prospective DNA barcode locations in the Brassicales order. Tandem repeats and SSR structures are reported as significant evidence of structural variations among the 26 tested specimens. Furthermore, selective pressure analysis was performed to estimate the substitution rate within the Moringaceae family, which revealing that the ndhA and accD genes are under positive selective pressure. The phylogenetic analysis of the Brassicales order produced an accurate monophyletic annotation cluster of the Moringaceae and Capparaceae species, offering unambiguous identification without overlapping groups between M. oleifera and M. peregrina, which are genetically strongly associated. Divergence time estimation suggests that the two Moringa species recently diversified, 0.467 Ma. Our findings highlight the first complete plastome of the Egyptian wild-type of M. peregrina, which can be used for determining plastome phylogenetic relationships and systematic evolution history within studies on the Moringaceae family.

Developing multifunctional crops by engineering Brassicaceae glucosinolate pathways

Glucosinolates (GSLs), found mainly in species of the Brassicaceae family, are one of the most well-studied classes of secondary metabolites. Produced by the action of myrosinase on GSLs, GSL-derived hydrolysis products (GHPs) primarily defend against biotic stress in planta. They also significantly affect the quality of crop products, with a subset of GHPs contributing unique food flavors and multiple therapeutic benefits or causing disagreeable food odors and health risks. Here, we explore the potential of these bioactive functions, which could be exploited for future sustainable agriculture. We first summarize our accumulated understanding of GSL diversity and distribution across representative Brassicaceae species. We then systematically discuss and evaluate the potential of exploited and unutilized genes involved in GSL biosynthesis, transport, and hydrolysis as candidate GSL engineering targets. Benefiting from available information on GSL and GHP functions, we explore options for multifunctional Brassicaceae crop ideotypes to meet future demand for food diversification and sustainable crop production. An integrated roadmap is subsequently proposed to guide ideotype development, in which maximization of beneficial effects and minimization of detrimental effects of GHPs could be combined and associated with various end uses. Based on several use-case examples, we discuss advantages and limitations of available biotechnological approaches that may contribute to effective deployment and could provide novel insights for optimization of future GSL engineering.

Phylogenomics and the flowering plant tree of life

The advances accelerated by next-generation sequencing and long-read sequencing technologies continue to provide an impetus for plant phylogenetic study. In the past decade, a large number of phylogenetic studies adopting hundreds to thousands of genes across a wealth of clades have emerged and ushered plant phylogenetics and evolution into a new era. In the meantime, a roadmap for researchers when making decisions across different approaches for their phylogenomic research design is imminent. This review focuses on the utility of genomic data (from organelle genomes, to both reduced representation sequencing and whole-genome sequencing) in phylogenetic and evolutionary investigations, describes the baseline methodology of experimental and analytical procedures, and summarizes recent progress in flowering plant phylogenomics at the ordinal, familial, tribal, and lower levels. We also discuss the challenges, such as the adverse impact on orthology inference and phylogenetic reconstruction raised from systematic errors, and underlying biological factors, such as whole-genome duplication, hybridization/introgression, and incomplete lineage sorting, together suggesting that a bifurcating tree may not be the best model for the tree of life. Finally, we discuss promising avenues for future plant phylogenomic studies.

Genome-Wide Analysis of Lipoxygenase (LOX) Genes in Angiosperms

Lipoxygenases (LOXs) are enzymes that catalyze the addition of an oxygen molecule to unsaturated fatty acids, thus forming hydroperoxides. In plants, these enzymes are encoded by a multigene family found in several organs with varying activity patterns, by which they are classified as LOX9 or LOX13. They are involved in several physiological functions, such as growth, fruit development, and plant defense. Despite several studies on genes of the LOX family in plants, most studies are restricted to a single species or a few closely related species. This study aimed to analyze the diversity, evolution, and expression of LOX genes in angiosperm species. We identified 247 LOX genes among 23 species of angiosperms and basal plants. Phylogenetic analyses identified clades supporting LOX13 and two main clades for LOX9: LOX9_A and LOX9_B. Eudicot species such as Tarenaya hassleriana, Capsella rubella, and Arabidopsis thaliana did not present LOX9_B genes; however, LOX9_B was present in all monocots used in this study. We identified that there were potential new subcellular localization patterns and conserved residues of oxidation for LOX9 and LOX13 yet unexplored. In summary, our study provides a basis for the further functional and evolutionary study of lipoxygenases in angiosperms.

The miR393-Target Module Regulates Plant Development and Responses to Biotic and Abiotic Stresses

MicroRNAs (miRNAs), a class of endogenous small RNAs, are broadly involved in plant development, morphogenesis and responses to various environmental stresses, through manipulating the cleavage, translational expression, or DNA methylation of target mRNAs. miR393 is a conserved miRNA family present in many plants, which mainly targets genes encoding the transport inhibitor response1 (TIR1)/auxin signaling F-box (AFB) auxin receptors, and thus greatly affects the auxin signal perception, Aux/IAA degradation, and related gene expression. This review introduces the advances made on the miR393/target module regulating plant development and the plant’s responses to biotic and abiotic stresses. This module is valuable for genetic manipulation of optimized conditions for crop growth and development and would also be helpful in improving crop yield through molecular breeding.

Recurrent Plant-Specific Duplications of KNL2 and Its Conserved Function as a Kinetochore Assembly Factor

KINETOCHORE NULL2 (KNL2) plays key role in the recognition of centromeres and new CENH3 deposition. To gain insight into the origin and diversification of the KNL2 gene, we reconstructed its evolutionary history in the plant kingdom. Our results indicate that the KNL2 gene in plants underwent three independent ancient duplications in ferns, grasses and eudicots. Additionally, we demonstrated that previously unclassified KNL2 genes could be divided into two clades αKNL2 and βKNL2 in eudicots and γKNL2 and δKNL2 in grasses, respectively. KNL2s of all clades encode the conserved SANTA domain, but only the αKNL2 and γKNL2 groups additionally encode the CENPC-k motif. In the more numerous eudicot sequences, signatures of positive selection were found in both αKNL2 and βKNL2 clades, suggesting recent or ongoing adaptation. The confirmed centromeric localization of βKNL2 and mutant analysis suggests that it participates in loading of new CENH3, similarly to αKNL2. A high rate of seed abortion was found in heterozygous βKNL2 plants and the germinated homozygous mutants did not develop beyond the seedling stage. Taken together, our study provides a new understanding of the evolutionary diversification of the plant kinetochore assembly gene KNL2, and suggests that the plant-specific duplicated KNL2 genes are involved in centromere and/or kinetochore assembly for preserving genome stability.

Genomes, repeatomes and interphase chromosome organization in the meadowfoam family (Limnanthaceae, Brassicales)

The meadowfoam family (Limnanthaceae) is one of the smallest and genomically underexplored families of the Brassicales. The Limnanthaceae harbor about seven species in the genus Limnanthes (meadowfoam) and Floerkea proserpinacoides (false mermaidweed), all native to North America. Because all Limnanthes and Floerkea species have only five chromosome pairs, i.e., a chromosome number rare in Brassicales and shared with Arabidopsis thaliana (Arabidopsis), we examined the Limnanthaceae genomes as a potential model system. Using low-coverage whole-genome sequencing data, we reexamined phylogenetic relationships and characterized the repeatomes of Limnanthaceae genomes. Phylogenies based on complete chloroplast and 35S rDNA sequences corroborated the sister relationship between Floerkea and Limnanthes and two major clades in the latter genus. The genome size of Limnanthaceae species ranges from 1.5 to 2.1 Gb, apparently due to the large increase in DNA repeats, which constitute 60-70% of their genomes. Repeatomes are dominated by long terminal repeat retrotransposons, while tandem repeats represent only less than 0.5% of the genomes. The average chromosome size in Limnanthaceae species (340-420 Mb) is more than 10 times larger than in Arabidopsis (32 Mb). A three-dimensional fluorescence in situ hybridization analysis demonstrated that the five chromosome pairs in interphase nuclei of Limnanthes species adopt the Rabl-like configuration.

Bacterial-type ferroxidase tunes iron-dependent phosphate sensing during Arabidopsis root development

Access to inorganic phosphate (Pi), a principal intermediate of energy and nucleotide metabolism, profoundly affects cellular activities and plant performance. In most soils, antagonistic Pi-metal interactions restrict Pi bioavailability, which guides local root development to maximize Pi interception. Growing root tips scout the essential but immobile mineral nutrient; however, the mechanisms monitoring external Pi status are unknown. Here, we show that Arabidopsis LOW PHOSPHATE ROOT 1 (LPR1), one key determinant of Fe-dependent Pi sensing in root meristems, encodes a novel ferroxidase of high substrate specificity and affinity (apparent K_M ∼ 2 μM Fe²⁺). LPR1 typifies an ancient, Fe-oxidizing multicopper protein family that evolved early upon bacterial land colonization. The ancestor of streptophyte algae and embryophytes (land plants) acquired LPR1-type ferroxidase from soil bacteria via horizontal gene transfer, a hypothesis supported by phylogenomics, homology modeling, and biochemistry. Our molecular and kinetic data on LPR1 regulation indicate that Pi-dependent Fe substrate availability determines LPR1 activity and function. Guided by the metabolic lifestyle of extant sister bacterial genera, we propose that Arabidopsis LPR1 monitors subtle concentration differentials of external Fe availability as a Pi-dependent cue to adjust root meristem maintenance via Fe redox signaling and cell wall modification. We further hypothesize that the acquisition of bacterial LPR1-type ferroxidase by embryophyte progenitors facilitated the evolution of local Pi sensing and acquisition during plant terrestrialization.

•

Arabidopsis thaliana LPR1 multicopper oxidase typifies a novel ferroxidase cohort

•

Fe availability tunes LPR1-dependent root responses to phosphate (Pi) limitation

•

LPR1 specificity links Fe-Pi interactions to root Pi sensing via redox cycling

•

Streptophyte ancestors acquired LPR1-type ferroxidase from soil bacteria by HGT

Ancient Biosyntheses in an Oil Crop: Glucosinolate Profiles in Limnanthes alba and Its Relatives (Limnanthaceae, Brassicales)

The glucosinolate (GSL) profiles of four Limnanthaceae species, including the oil crop Limnanthes alba (meadowfoam), were investigated by an ultrahigh-performance liquid chromatography-quadrupole time-of-flight tandem mass spectrometry (UHPLC-QToF-MS/MS) analysis of desulfoGSLs after desulfation of native GSLs, supplemented by NMR of desulfated 2-hydroxy-2-methylpropylGSL and 3-methoxybenzylGSL. Leaves, roots, and seeds were investigated, providing an overview of biosynthetic capabilities in the genera Floerkea and Limnanthes. Methoxyl groups on benzylGSLs were in meta but not para positions; two 3,5-disubstituted benzylGSLs are tentatively proposed. 2-Hydroxy-2-methylpropylGSL was accompanied by an isomer that was not a previously reported GSL. The combined GSL profile of the family included GSLs derived from valine, leucine, isoleucine, phenylalanine, and tyrosine, and possibly methionine and tryptophan. Substituted indole GSLs and GSLs derived from chain-elongated amino acids or alanine were searched for but not detected. Hypothetic glycosides of GSLs were detected at low levels. Based on biochemical interpretation, we suggest biosynthetic schemes and gene families (CYP79C, GSOH) relevant for tailoring GSL profiles in Limnanthes crops.

Pif1 helicases and the evidence for a prokaryotic origin of Helitrons

Helitrons are the only group of rolling-circle transposons that encode a transposase with a helicase domain (Hel), which belongs to the Pif1 family. Because Pif1 helicases are important components of eukaryotic genomes, it has been suggested that Hel domains probably originated after a host eukaryotic Pif1 gene was captured by a Helitron ancestor. However, the few analyses exploring the evolution of Helitron transposases (RepHel) have focused on its Rep domain, which is also present in other mobile genetic elements. Here, we used phylogenetic and nonmetric multidimensional scaling analyses to investigate the relationship between Hel domains and Pif1-like helicases from a variety of organisms. Our results reveal that Hel domains are only distantly related to genomic helicases from eukaryotes and prokaryotes, and thus are unlikely to have originated from a captured Pif1 gene. Based on this evidence, and on recent studies indicating that Rep domains are more closely related to rolling-circle plasmids and phages, we suggest that Helitrons are descendants of a RepHel-encoding prokaryotic plasmid element that invaded eukaryotic genomes before the radiation of its major groups. We discuss how a Pif1-like helicase domain might have favored the transposition of Helitrons in eukaryotes beyond simply unwinding DNA intermediates. Finally, we demonstrate that some examples in the literature describing genomic helicases from eukaryotes actually consist of Hel domains from Helitrons, a finding that underscores how transposons can hamper the analysis of eukaryotic genes. This investigation also revealed that two groups of land plants appear to have lost genomic Pif1 helicases independently.

Ancestral Flowering Plant Chromosomes and Gene Orders Based on Generalized Adjacencies and Chromosomal Gene Co-Occurrences

Recurrent whole genome duplication and the ensuing loss of redundant genes-fractionation-complicate efforts to reconstruct the gene orders and chromosomes of the ancestors associated with the nodes of a phylogeny. Loss of genes disrupts the gene adjacencies key to current techniques. With our RACCROCHE pipeline, instead of starting with the inference of short ancestral segments, we suggest delaying the choice of gene adjacencies while we accumulate many more syntenically validated generalized (gapped) adjacencies. We obtain longer ancestral contigs using maximum weight matching (MWM). Similarly, we do not construct chromosomes by successively piecing together contigs into larger segments, but rather compile counts of pairwise contig co-occurrences on the set of extant genomes and use these to cluster the contigs. Chromosome-level contig assemblies for a monoploid genome emerge naturally at each node of the phylogeny and the contigs then can be ordered along the chromosome. Sampling alternative MWM solutions, visualizing heat maps, and applying gap statistics allow us to estimate the number of chromosomes in the reconstruction. We introduce several measures of quality: length of contigs, continuity of contig structure on successive ancestors, coverage of the extant genome by the reconstruction, and rearrangement relations among the inferred chromosomes. The reconstructed ancestors are visualized by painting the ancestral projections on the descendant genomes. We submit genomes drawn from a broad range of monocot orders to our pipeline, confirming the tetraploidization event "tau" in the stem lineage between the alismatids and the lilioids. We show additional applications to the Solanaceae and to four Brassica genomes, producing evidence about the monoploid ancestor in each case.

Plastid phylogenomic insights into relationships of all flowering plant families

BACKGROUND

Flowering plants (angiosperms) are dominant components of global terrestrial ecosystems, but phylogenetic relationships at the familial level and above remain only partially resolved, greatly impeding our full understanding of their evolution and early diversification. The plastome, typically mapped as a circular genome, has been the most important molecular data source for plant phylogeny reconstruction for decades.

RESULTS

Here, we assembled by far the largest plastid dataset of angiosperms, composed of 80 genes from 4792 plastomes of 4660 species in 2024 genera representing all currently recognized families. Our phylogenetic tree (PPA II) is essentially congruent with those of previous plastid phylogenomic analyses but generally provides greater clade support. In the PPA II tree, 75% of nodes at or above the ordinal level and 78% at or above the familial level were resolved with high bootstrap support (BP ≥ 90). We obtained strong support for many interordinal and interfamilial relationships that were poorly resolved previously within the core eudicots, such as Dilleniales, Saxifragales, and Vitales being resolved as successive sisters to the remaining rosids, and Santalales, Berberidopsidales, and Caryophyllales as successive sisters to the asterids. However, the placement of magnoliids, although resolved as sister to all other Mesangiospermae, is not well supported and disagrees with topologies inferred from nuclear data. Relationships among the five major clades of Mesangiospermae remain intractable despite increased sampling, probably due to an ancient rapid radiation.

CONCLUSIONS

We provide the most comprehensive dataset of plastomes to date and a well-resolved phylogenetic tree, which together provide a strong foundation for future evolutionary studies of flowering plants.

Nuclear organization in crucifer genomes: nucleolus-associated telomere clustering is not a universal interphase configuration in Brassicaceae

Arabidopsis thaliana has become a major plant research model, where interphase nuclear organization exhibits unique features, including nucleolus-associated telomere clustering. The chromocenter (CC)-loop model, or rosette-like configuration, describes intranuclear chromatin organization in Arabidopsis as megabase-long loops anchored in, and emanating from, peripherally positioned CCs, with those containing telomeres associating with the nucleolus. To investigate whether the CC-loop organization is universal across the mustard family (crucifers), the nuclear distributions of centromeres, telomeres and nucleoli were analyzed by fluorescence in situ hybridization in seven diploid species (2n = 10-16) representing major crucifer clades with an up to 26-fold variation in genome size (160-4260 Mb). Nucleolus-associated telomere clustering was confirmed in Arabidopsis (157 Mb) and was newly identified as the major nuclear phenotype in other species with a small genome (215-381 Mb). In large-genome species (2611-4264 Mb), centromeres and telomeres adopted a Rabl-like configuration or dispersed distribution in the nuclear interior; telomeres only rarely associated with the nucleolus. In Arabis cypria (381 Mb) and Bunias orientalis (2611 Mb), tissue-specific patterns deviating from the major nuclear phenotypes were observed in anther and stem tissues, respectively. The rosette-like configuration, including nucleolus-associated telomere clustering in small-genome species from different infrafamiliar clades, suggests that genomic properties rather than phylogenetic position determine the interphase nuclear organization. Our data suggest that nuclear genome size, average chromosome size and degree of longitudinal chromosome compartmentalization affect interphase chromosome organization in crucifer genomes.

Comparative phylogenetics of repetitive elements in a diverse order of flowering plants (Brassicales)

Genome sizes of plants have long piqued the interest of researchers due to the vast differences among organisms. However, the mechanisms that drive size differences have yet to be fully understood. Two important contributing factors to genome size are expansions of repetitive elements, such as transposable elements (TEs), and whole-genome duplications (WGD). Although studies have found correlations between genome size and both TE abundance and polyploidy, these studies typically test for these patterns within a genus or species. The plant order Brassicales provides an excellent system to further test if genome size evolution patterns are consistent across larger time scales, as there are numerous WGDs. This order is also home to one of the smallest plant genomes, Arabidopsis thaliana-chosen as the model plant system for this reason-as well as to species with very large genomes. With new methods that allow for TE characterization from low-coverage genome shotgun data and 71 taxa across the Brassicales, we confirm the correlation between genome size and TE content, however, we are unable to reconstruct phylogenetic relationships and do not detect any shift in TE abundance associated with WGD.

Hundreds of nuclear and plastid loci yield novel insights into orchid relationships

PREMISE

The inference of evolutionary relationships in the species-rich family Orchidaceae has hitherto relied heavily on plastid DNA sequences and limited taxon sampling. Previous studies have provided a robust plastid phylogenetic framework, which was used to classify orchids and investigate the drivers of orchid diversification. However, the extent to which phylogenetic inference based on the plastid genome is congruent with the nuclear genome has been only poorly assessed.

METHODS

We inferred higher-level phylogenetic relationships of orchids based on likelihood and ASTRAL analyses of 294 low-copy nuclear genes sequenced using the Angiosperms353 universal probe set for 75 species (representing 69 genera, 16 tribes, 24 subtribes) and a concatenated analysis of 78 plastid genes for 264 species (117 genera, 18 tribes, 28 subtribes). We compared phylogenetic informativeness and support for the nuclear and plastid phylogenetic hypotheses.

RESULTS

Phylogenetic inference using nuclear data sets provides well-supported orchid relationships that are highly congruent between analyses. Comparisons of nuclear gene trees and a plastid supermatrix tree showed that the trees are mostly congruent, but revealed instances of strongly supported phylogenetic incongruence in both shallow and deep time. The phylogenetic informativeness of individual Angiosperms353 genes is in general better than that of most plastid genes.

CONCLUSIONS

Our study provides the first robust nuclear phylogenomic framework for Orchidaceae and an assessment of intragenomic nuclear discordance, plastid-nuclear tree incongruence, and phylogenetic informativeness across the family. Our results also demonstrate what has long been known but rarely thoroughly documented: nuclear and plastid phylogenetic trees can contain strongly supported discordances, and this incongruence must be reconciled prior to interpretation in evolutionary studies, such as taxonomy, biogeography, and character evolution.

Glucosinolate profiles and phylogeny in Barbarea compared to other tribe Cardamineae (Brassicaceae) and Reseda (Resedaceae), based on a library of ion trap HPLC-MS/MS data of reference desulfoglucosinolates

A library of ion trap MS2 spectra and HPLC retention times reported here allowed distinction in plants of at least 70 known glucosinolates (GSLs) and some additional proposed GSLs. We determined GSL profiles of selected members of the tribe Cardamineae (Brassicaceae) as well as Reseda (Resedaceae) used as outgroup in evolutionary studies. We included several accessions of each species and a range of organs, and paid attention to minor peaks and GSLs not detected. In this way, we obtained GSL profiles of Barbarea australis, Barbarea grayi, Planodes virginica selected for its apparent intermediacy between Barbarea and the remaining tribe and family, and Rorippa sylvestris and Nasturtium officinale, for which the presence of acyl derivatives of GSLs was previously untested. We also screened Armoracia rusticana, with a remarkably diverse GSL profile, the emerging model species Cardamine hirsuta, for which we discovered a GSL polymorphism, and Reseda luteola and Reseda odorata. The potential for aliphatic GSL biosynthesis in Barbarea vulgaris was of interest, and we subjected P-type and G-type B. vulgaris to several induction regimes in an attempt to induce aliphatic GSL. However, aliphatic GSLs were not detected in any of the B. vulgaris types. We characterized the investigated chemotypes phylogenetically, based on nuclear rDNA internal transcribed spacer (ITS) sequences, in order to understand their relation to the species B. vulgaris in general, and found them to be representative of the species as it occurs in Europe, as far as documented in available ITS-sequence repositories. In short, we provide GSL profiles of a wide variety of tribe Cardamineae plants and conclude aliphatic GSLs to be absent or below our limit of detection in two major evolutionary lines of B. vulgaris. Concerning analytical chemistry, we conclude that availability of authentic reference compounds or reference materials is critical for reliable GSL analysis and characterize two publicly available reference materials: seeds of P. virginica and N. officinale.

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

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

Mix and match: Patchwork domain evolution of the land plant-specific Ca2+-permeable mechanosensitive channel MCA

Multidomain proteins can have a complex evolutionary history that may involve de novo domain evolution, recruitment and / or recombination of existing domains and domain losses. Here, the domain evolution of the plant-specific Ca²⁺-permeable mechanosensitive channel protein, MID1-COMPLEMENTING ACTIVITY (MCA), was investigated. MCA, a multidomain protein, possesses a Ca²⁺-influx-MCA^func domain and a PLAC8 domain. Profile Hidden Markov Models (HMMs) of domains were assessed in 25 viridiplantae proteomes. While PLAC8 was detected in plants, animals, and fungi, MCA^func was found in streptophytes but not in chlorophytes. Full MCA proteins were only found in embryophytes. We identified the MCA^func domain in all streptophytes including charophytes where it appeared in E3 ubiquitin ligase-like proteins. Our Maximum Likelihood (ML) analyses suggested that the MCA^func domain evolved early in the history of streptophytes. The PLAC8 domain showed similarity to Plant Cadmium Resistance (PCR) genes, and the coupling of MCA^func and PLAC8 seemed to represent a single evolutionary event. This combination is unique in MCA, and does not exist in other plant mechanosensitive channels. Within angiosperms, gene duplications increased the number of MCAs. Considering their role in mechanosensing in roots, MCA might be instrumental for the rise of land plants. This study provides a textbook example of de novo domain emergence, recombination, duplication, and losses, leading to the convergence of function of proteins in plants.

Evolution of Conserved Noncoding Sequences in Arabidopsis thaliana

Recent pangenome studies have revealed a large fraction of the gene content within a species exhibits presence-absence variation (PAV). However, coding regions alone provide an incomplete assessment of functional genomic sequence variation at the species level. Little to no attention has been paid to noncoding regulatory regions in pangenome studies, though these sequences directly modulate gene expression and phenotype. To uncover regulatory genetic variation, we generated chromosome-scale genome assemblies for thirty Arabidopsis thaliana accessions from multiple distinct habitats and characterized species level variation in Conserved Noncoding Sequences (CNS). Our analyses uncovered not only PAV and positional variation (PosV) but that diversity in CNS is nonrandom, with variants shared across different accessions. Using evolutionary analyses and chromatin accessibility data, we provide further evidence supporting roles for conserved and variable CNS in gene regulation. Additionally, our data suggests that transposable elements contribute to CNS variation. Characterizing species-level diversity in all functional genomic sequences may later uncover previously unknown mechanistic links between genotype and phenotype.

Glucosinolate Biosynthesis and the Glucosinolate–Myrosinase System in Plant Defense

Insect pests represent a major global challenge to important agricultural crops. Insecticides are often applied to combat such pests, but their use has caused additional challenges such as environmental contamination and human health issues. Over millions of years, plants have evolved natural defense mechanisms to overcome insect pests and pathogens. One such mechanism is the production of natural repellents or specialized metabolites like glucosinolates. There are three types of glucosinolates produced in the order Brassicales: aliphatic, indole, and benzenic glucosinolates. Upon insect herbivory, a “mustard oil bomb” consisting of glucosinolates and their hydrolyzing enzymes (myrosinases) is triggered to release toxic degradation products that act as insect deterrents. This review aims to provide a comprehensive summary of glucosinolate biosynthesis, the “mustard oil bomb”, and how these metabolites function in plant defense against pathogens and insects. Understanding these defense mechanisms will not only allow us to harness the benefits of this group of natural metabolites for enhancing pest control in Brassicales crops but also to transfer the “mustard oil bomb” to non-glucosinolate producing crops to boost their defense and thereby reduce the use of chemical pesticides.

Expression Partitioning of Duplicate Genes at Single Cell Resolution in Arabidopsis Roots

Gene duplication is a key evolutionary phenomenon, prevalent in all organisms but particularly so in plants, where whole genome duplication (WGD; polyploidy) is a major force in genome evolution. Much effort has been expended in attempting to understand the evolution of duplicate genes, addressing such questions as why some paralog pairs rapidly return to single copy status whereas, in other pairs, both paralogs are retained and may diverge in expression pattern or function. The effect of a gene - its site of expression and thus the initial locus of its function - occurs at the level of a cell comprising a single cell type at a given state of the cell's development. Using Arabidopsis thaliana single cell transcriptomic data we categorized patterns of expression for 11,470 duplicate gene pairs across 36 cell clusters comprising nine cell types and their developmental states. Among these 11,470 pairs, 10,187 (88.8%) had at least one copy expressed in at least one of the 36 cell clusters. Pairs produced by WGD more often had both paralogs expressed in root cells than did pairs produced by small scale duplications. Three quarters of gene pairs expressed in the 36 cell clusters (7,608/10,187) showed extreme expression bias in at least one cluster, including 352 cases of reciprocal bias, a pattern consistent with expression subfunctionalization. More than twice as many pairs showed reciprocal expression bias between cell states than between cell types or between roots and leaves. A group of 33 gene pairs with reciprocal expression bias showed evidence of concerted divergence of gene networks in stele vs. epidermis. Pairs with both paralogs expressed without bias were less likely to have paralogs with divergent mutant phenotypes; such bias-free pairs showed evidence of preservation by maintenance of dosage balance. Overall, we found considerable evidence of shifts in gene expression following duplication, including in >80% of pairs encoding 7,653 genes expressed ubiquitously in all root cell types and states for which we inferred the polarity of change.

Phylogeny and multiple independent whole-genome duplication events in the Brassicales

PREMISE

Whole-genome duplications (WGDs) are prevalent throughout the evolutionary history of plants. For example, dozens of WGDs have been phylogenetically localized across the order Brassicales, specifically, within the family Brassicaceae. A WGD event has also been identified in the Cleomaceae, the sister family to Brassicaceae, yet its placement, as well as that of WGDs in other families in the order, remains unclear.

METHODS

Phylo-transcriptomic data were generated and used to infer a nuclear phylogeny for 74 Brassicales taxa. Genome survey sequencing was also performed on 66 of those taxa to infer a chloroplast phylogeny. These phylogenies were used to assess and confirm relationships among the major families of the Brassicales and within Brassicaceae. Multiple WGD inference methods were then used to assess the placement of WGDs on the nuclear phylogeny.

RESULTS

Well-supported chloroplast and nuclear phylogenies for the Brassicales and the putative placement of the Cleomaceae-specific WGD event Th-ɑ are presented. This work also provides evidence for previously hypothesized WGDs, including a well-supported event shared by at least two members of the Resedaceae family, and a possible event within the Capparaceae.

CONCLUSIONS

Phylogenetics and the placement of WGDs within highly polyploid lineages continues to be a major challenge. This study adds to the conversation on WGD inference difficulties by demonstrating that sampling is especially important for WGD identification and phylogenetic placement. Given its economic importance and genomic resources, the Brassicales continues to be an ideal group for assessing WGD inference methods.

Phylogenetic and evolutionary features of the plastome of Tropaeolum pentaphyllum Lam. (Tropaeolaceae)

Complete plastome sequence of Tropaeolum pentaphyllum revealed molecular markers, hotspots of nucleotide polymorphism, RNA editing sites and phylogenetic aspects Tropaeolaceae Juss. ex DC. comprises approximately 95 species across North and South Americas. Tropaeolum pentaphyllum Lam. is an unconventional and endangered species with occurrence in some countries of South America. Although this species presents nutritional, medicinal and ornamental uses, genetic studies involving natural populations or promising genotypes are practically non-existent. Here, we report the nucleotide sequence of T. pentaphyllum plastome. It represents the first complete plastome sequence of the family Tropaeolaceae to be fully sequenced and analyzed in detail. The sequencing data revealed that the T. pentaphyllum plastome is highly similar to the plastomes of other Brassicales. Notwithstanding, our analyses detected some specific features concerning events of IR expansion and structural changes in some genes such as matK, rpoA, and rpoC2. We also detected 251 SSR loci, nine hotspots of nucleotide polymorphism, and two specific RNA editing sites in the plastome of T. pentaphyllum. Moreover, plastid phylogenomic inference indicated a closed relationship between the families Tropaeolaceae and Akaniaceae, which formed a sister group to Moringaceae-Caricaceae. Finally, our data bring new molecular markers and evolutionary features to be applied in the natural population, germplasm collection, and genotype selection aiming conservation, genetic diversity evaluation, and exploitation of this endangered species.

Gene Balance Predicts Transcriptional Responses Immediately Following Ploidy Change in Arabidopsis thaliana

The gene balance hypothesis postulates that there is selection on gene copy number (gene dosage) to preserve the stoichiometric balance among interacting proteins. This presupposes that gene product abundance is governed by gene dosage and that gene dosage responses are consistent for interacting genes in a dosage-balance-sensitive network or complex. Gene dosage responses, however, have rarely been quantified, and the available data suggest that they are highly variable. We sequenced the transcriptomes of two synthetic autopolyploid accessions of Arabidopsis (Arabidopsis thaliana) and their diploid progenitors, as well as one natural tetraploid and its synthetic diploid produced via haploid induction, to estimate transcriptome size and dosage responses immediately following ploidy change. Similar to what has been observed in previous studies, overall transcriptome size does not exhibit a simple doubling in response to genome doubling, and individual gene dosage responses are highly variable in all three accessions, indicating that expression is not strictly coupled with gene dosage. Nonetheless, putatively dosage balance-sensitive gene groups (Gene Ontology terms, metabolic networks, gene families, and predicted interacting proteins) exhibit smaller and more coordinated dosage responses than do putatively dosage-insensitive gene groups, suggesting that constraints on dosage balance operate immediately following whole-genome duplication and that duplicate gene retention patterns are shaped by selection to preserve dosage balance.

Genetic analysis of the Arabidopsis TIR1/AFB auxin receptors reveals both overlapping and specialized functions

The TIR1/AFB auxin co-receptors mediate diverse responses to the plant hormone auxin. The Arabidopsis genome encodes six TIR1/AFB proteins representing three of the four clades that were established prior to angiosperm radiation. To determine the role of these proteins in plant development we performed an extensive genetic analysis involving the generation and characterization of all possible multiply-mutant lines. We find that loss of all six TIR1/AFB proteins results in early embryo defects and eventually seed abortion, and yet a single wild-type allele of TIR1 or AFB2 is sufficient to support growth throughout development. Our analysis reveals extensive functional overlap between even the most distantly related TIR1/AFB genes except for AFB1. Surprisingly, AFB1 has a specialized function in rapid auxin-dependent inhibition of root growth and early phase of root gravitropism. This activity may be related to a difference in subcellular localization compared to the other members of the family.

Glucosinolate structural diversity, identification, chemical synthesis and metabolism in plants

The glucosinolates (GSLs) is a well-defined group of plant metabolites characterized by having an S-β-d-glucopyrano unit anomerically connected to an O-sulfated (Z)-thiohydroximate function. After enzymatic hydrolysis, the sulfated aglucone can undergo rearrangement to an isothiocyanate, or form a nitrile or other products. The number of GSLs known from plants, satisfactorily characterized by modern spectroscopic methods (NMR and MS) by mid-2018, is 88. In addition, a group of partially characterized structures with highly variable evidence counts for approximately a further 49. This means that the total number of characterized GSLs from plants is somewhere between 88 and 137. The diversity of GSLs in plants is critically reviewed here, resulting in significant discrepancies with previous reviews. In general, the well-characterized GSLs show resemblance to C-skeletons of the amino acids Ala, Val, Leu, Trp, Ile, Phe/Tyr and Met, or to homologs of Ile, Phe/Tyr or Met. Insufficiently characterized, still hypothetic GSLs include straight-chain alkyl GSLs and chain-elongated GSLs derived from Leu. Additional reports (since 2011) of insufficiently characterized GSLs are reviewed. Usually the crucial missing information is correctly interpreted NMR, which is the most effective tool for GSL identification. Hence, modern use of NMR for GSL identification is also reviewed and exemplified. Apart from isolation, GSLs may be obtained by organic synthesis, allowing isotopically labeled GSLs and any kind of side chain. Enzymatic turnover of GSLs in plants depends on a considerable number of enzymes and other protein factors and furthermore depends on GSL structure. Identification of GSLs must be presented transparently and live up to standard requirements in natural product chemistry. Unfortunately, many recent reports fail in these respects, including reports based on chromatography hyphenated to MS. In particular, the possibility of isomers and isobaric structures is frequently ignored. Recent reports are re-evaluated and interpreted as evidence of the existence of "isoGSLs", i.e. non-GSL isomers of GSLs in plants. For GSL analysis, also with MS-detection, we stress the importance of using authentic standards.

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.

Distinguishing successive ancient polyploidy levels based on genome-internal syntenic alignment

BACKGROUND

A basic tool for studying the polyploidization history of a genome, especially in plants, is the distribution of duplicate gene similarities in syntenically aligned regions of a genome. This distribution can usually be decomposed into two or more components identifiable by peaks, or local maxima, each representing a different polyploidization event. The distributions may be generated by means of a discrete time branching process, followed by a sequence divergence model. The branching process, as well as the inference of fractionation rates based on it, requires knowledge of the ploidy level of each event, which cannot be directly inferred from the pair similarity distribution.

RESULTS

For a sequence of two events of unknown ploidy, either tetraploid, giving rise to whole genome doubling (WGD), or hexaploid, giving rise to whole genome tripling (WGT), we base our analysis on triples of similar genes. We calculate the probability of the four triplet types with origins in one or the other event, or both, and impose a mutational model so that the distribution resembles the original data. Using a ML transition point in the similarities between the two events as a discriminator for the hypothesized origin of each similarity, we calculate the predicted number of triplets of each type for each model combining WGT and/or WGD. This yields a predicted profile of triplet types for each model. We compare the observed and predicted triplet profiles for each model to confirm the polyploidization history of durian, poplar and cabbage.

CONCLUSIONS

We have developed a way of inferring the ploidy of up to three successive WGD and/or WGT events by estimating the time of origin of each of the similarities in triples of genes. This may be generalized to a larger number of events and to higher ploidies.

How to build a fruit: Transcriptomics of a novel fruit type in the Brassiceae

Comparative gene expression studies are invaluable for predicting how existing genetic pathways may be modified or redeployed to produce novel and variable phenotypes. Fruits are ecologically important organs because of their impact on plant fitness and seed dispersal, modifications in which results in morphological variation across species. A novel fruit type in the Brassicaceae known as heteroarthrocarpy enables distinct dispersal methods in a single fruit through segmentation via a lateral joint and variable dehiscence at maturity. Given the close relationship to Arabidopsis, species that exhibit heteroarthrocarpy are powerful models to elucidate how differences in gene expression of a fruit patterning pathway may result in novel fruit types. Transcriptomes of distal, joint, and proximal regions from Erucaria erucarioides and Cakile lanceolata were analyzed to elucidate within fruit and between species differences in whole transcriptome, gene ontology, and fruit patterning expression profiles. Whole transcriptome expression profiles vary between fruit regions in patterns that are consistent with fruit anatomy. These transcriptomic variances do not correlate with changes in gene ontology, as they remain generally stable within and between both species. Upstream regulators in the fruit patterning pathway, FILAMENTOUS FLOWER and YABBY3, are expressed in the distal and proximal regions of E. erucarioides, but not in the joint, implicating alterations in the pathway in heteroarthrocarpic fruits. Downstream gene, INDEHISCENT, is significantly upregulated in the abscissing joint region of C. lanceolata, which suggests repurposing of valve margin genes for novel joint disarticulation in an otherwise indehiscent fruit. In summary, these data are consistent with modifications in fruit patterning genes producing heteroarthrocarpic fruits through different components of the pathway relative to other indehiscent, non-heteroarthrocarpic, species within the family. Our understanding of fruit development in Arabidopsis is now extended to atypical siliques within the Brassicaceae, facilitating future studies on seed shattering in important Brassicaceous crops and pernicious weeds.

The plastid genome of Pentadiplandra brazzeana Baillon (Pentadiplandraceae)

Pentadiplandra brazzeana is the only species of Pentadiplandraceae. This species is shrubs or tuberous liana with a natural distribution in tropical West Africa. This study firstly determined its complete plastome. The plastome is totally 156,625 bp in length, which contains a pair of 26,751-bp-long inverted repeat regions (IRs), a large single copy region of 85,318 bp, and a small single copy region of 17,805 bp. A total of 112 unique genes were identified in this plastome, of which 78 are protein-coding genes, 30 are tRNA genes, and 4 are rRNA genes. Phylogenetic analysis based on 82 genes fully resolved relationships among sampled families of Brassicales, which are consistent with previous studies.

[Biosynthesis of indole glucosinolates and ecological role of secondary modification pathways]

Indole glucosinolates are plant secondary metabolites derived from the amino acid tryptophan. They are part of a large group of sulfur-containing molecules almost exclusively found among Brassicales, which include the mustard family (Brassicaceae) with many edible plant species of major nutritional importance. These compounds mediate numerous interactions between these plants and their natural enemies and are therefore of major biological and economical interest. This literature review aims at taking stock of recent advances of our knowledge about the biosynthetic pathways of indole glucosinolates, but also about the defense strategies and ecological processes involving these metabolites.

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.

Evolution of Glucosinolate Diversity via Whole-Genome Duplications, Gene Rearrangements, and Substrate Promiscuity

Over several decades, glucosinolates have become a model system for the study of specialized metabolic diversity in plants. The near-complete identification of biosynthetic enzymes, regulators, and transporters has provided support for the role of gene duplication and subsequent changes in gene expression, protein function, and substrate specificity as the evolutionary bases of glucosinolate diversity. Here, we provide examples of how whole-genome duplications, gene rearrangements, and substrate promiscuity potentiated the evolution of glucosinolate biosynthetic enzymes, regulators, and transporters by natural selection. This in turn may have led to the repeated evolution of glucosinolate metabolism and diversity in higher plants.

Whole-Genome Duplication and Plant Macroevolution

Whole-genome duplication (WGD) is characteristic of almost all fundamental lineages of land plants. Unfortunately, the timings of WGD events are loosely constrained and hypotheses of evolutionary consequence are poorly formulated, making them difficult to test. Using examples from across the plant kingdom, we show that estimates of timing can be improved through the application of molecular clock methodology to multigene datasets. Further, we show that phenotypic change can be quantified in morphospaces and that relative phenotypic disparity can be compared in the light of WGD. Together, these approaches facilitate tests of hypotheses on the role of WGD in plant evolution, underscoring the potential of plants as a model system for investigating the role WGD in macroevolution.

Lessons from Cleomaceae, the Sister of Crucifers

Cleomaceae is a diverse group well-suited to addressing fundamental genomic and evolutionary questions as the sister group to Brassicaceae, facilitating transfer of knowledge from the model Arabidopsis thaliana. Phylogenetic and taxonomic revisions provide a framework for examining the evolution of substantive morphological and physiology diversity in Cleomaceae, but not necessarily in Brassicaceae. The investigation of both nested and contrasting whole-genome duplications (WGDs) between Cleomaceae and Brassicaceae allows comparisons of independently duplicated genes and investigation of whether they may be drivers of the observed innovations. Further, a wealth of outstanding genetic research has provided insight into how the important alternative carbon fixation pathway, C₄ photosynthesis, has evolved via differential expression of a suite of genes, of which the underlying mechanisms are being elucidated.

Comparative Analyses of Chloroplast Genomes of Cucurbitaceae Species: Lights into Selective Pressures and Phylogenetic Relationships

Cucurbitaceae is the fourth most important economic plant family with creeping herbaceous species mainly distributed in tropical and subtropical regions. Here, we described and compared the complete chloroplast genome sequences of ten representative species from Cucurbitaceae. The lengths of the ten complete chloroplast genomes ranged from 155,293 bp (C. sativus) to 158,844 bp (M. charantia), and they shared the most common genomic features. 618 repeats of three categories and 813 microsatellites were found. Sequence divergence analysis showed that the coding and IR regions were highly conserved. Three protein-coding genes (accD, clpP, and matK) were under selection and their coding proteins often have functions in chloroplast protein synthesis, gene transcription, energy transformation, and plant development. An unconventional translation initiation codon of psbL gene was found and provided evidence for RNA editing. Applying BI and ML methods, phylogenetic analysis strongly supported the position of Gomphogyne, Hemsleya, and Gynostemma as the relatively original lineage in Cucurbitaceae. This study suggested that the complete chloroplast genome sequences were useful for phylogenetic studies. It would also determine potential molecular markers and candidate DNA barcodes for coming studies and enrich the valuable complete chloroplast genome resources of Cucurbitaceae.