Cardamine hirsuta: a versatile genetic system for comparative studies

Topological analysis of 3D digital ovules identifies cellular patterns associated with ovule shape diversity

Tissue morphogenesis remains poorly understood. In plants, a central problem is how the 3D cellular architecture of a developing organ contributes to its final shape. We address this question through a comparative analysis of ovule morphogenesis, taking advantage of the diversity in ovule shape across angiosperms. Here, we provide a 3D digital atlas of Cardamine hirsuta ovule development at single cell resolution and compare it with an equivalent atlas of Arabidopsis thaliana. We introduce nerve-based topological analysis as a tool for unbiased detection of differences in cellular architectures and corroborate identified topological differences between two homologous tissues by comparative morphometrics and visual inspection. We find that differences in topology, cell volume variation and tissue growth patterns in the sheet-like integuments and the bulbous chalaza are associated with differences in ovule curvature. In contrast, the radialized conical ovule primordia and nucelli exhibit similar shapes, despite differences in internal cellular topology and tissue growth patterns. Our results support the notion that the structural organization of a tissue is associated with its susceptibility to shape changes during evolutionary shifts in 3D cellular architecture.

Shaping leaves through TALE homeodomain transcription factors

The first TALE homeodomain transcription factor gene to be described in plants was maize knotted1 (kn1). Dominant mutations in kn1 disrupt leaf development, with abnormal knots of tissue forming in the leaf blade. kn1 was found to be expressed in the shoot meristem but not in a peripheral region that gives rise to leaves. Furthermore, KN1 and closely related proteins were excluded from initiating and developing leaves. These findings were a prelude to a large body of work wherein TALE homeodomain proteins have been identified as vital regulators of meristem homeostasis and organ development in plants. KN1 homologues are widely represented across land plant taxa. Thus, studying the regulation and mechanistic action of this gene class has allowed investigations into the evolution of diverse plant morphologies. This review will focus on the function of TALE homeodomain transcription factors in leaf development in eudicots. Here, we discuss how TALE homeodomain proteins contribute to a spectrum of leaf forms, from the simple leaves of Arabidopsis thaliana to the compound leaves of Cardamine hirsuta and species beyond the Brassicaceae.

A chromosome-level genome assembly for the amphibious plant Rorippa aquatica reveals its allotetraploid origin and mechanisms of heterophylly upon submergence

The ability to respond to varying environments is crucial for sessile organisms such as plants. The amphibious plant Rorippa aquatica exhibits a striking type of phenotypic plasticity known as heterophylly, a phenomenon in which leaf form is altered in response to environmental factors. However, the underlying molecular mechanisms of heterophylly are yet to be fully understood. To uncover the genetic basis and analyze the evolutionary processes driving heterophylly in R. aquatica, we assembled the chromosome-level genome of the species. Comparative chromosome painting and chromosomal genomics revealed that allopolyploidization and subsequent post-polyploid descending dysploidy occurred during the speciation of R. aquatica. Based on the obtained genomic data, the transcriptome analyses revealed that ethylene signaling plays a central role in regulating heterophylly under submerged conditions, with blue light signaling acting as an attenuator of ethylene signal. The assembled R. aquatica reference genome provides insights into the molecular mechanisms and evolution of heterophylly.

Injection-based hairy root induction and plant regeneration techniques in Brassicaceae

BACKGROUND

Hairy roots constitute a valuable tissue culture system for species that are difficult to propagate through conventional seed-based methods. Moreover, the generation of transgenic plants derived from hairy roots can be facilitated by employing carefully designed hormone-containing media.

RESULTS

We initiated hairy root formation in the rare crucifer species Asperuginoides axillaris via an injection-based protocol using the Agrobacterium strain C58C1 harboring a hairy root-inducing (Ri) plasmid and successfully regenerated plants from established hairy root lines. Our study confirms the genetic stability of both hairy roots and their derived regenerants and highlights their utility as a permanent source of mitotic chromosomes for cytogenetic investigations. Additionally, we have developed an effective embryo rescue protocol to circumvent seed dormancy issues in A. axillaris seeds. By using inflorescence primary stems of Arabidopsis thaliana and Cardamine hirsuta as starting material, we also established hairy root lines that were subsequently used for regeneration studies.

CONCLUSION

We developed efficient hairy root transformation and regeneration protocols for various crucifers, namely A. axillaris, A. thaliana, and C. hirsuta. Hairy roots and derived regenerants can serve as a continuous source of plant material for molecular and cytogenetic analyses.

Interspersed expression of CUP-SHAPED COTYLEDON2 and REDUCED COMPLEXITY shapes Cardamine hirsuta complex leaf form

How genetically regulated growth shapes organ form is a key problem in developmental biology. Here, we investigate this problem using the leaflet-bearing complex leaves of Cardamine hirsuta as a model. Leaflet development requires the action of two growth-repressing transcription factors: REDUCED COMPLEXITY (RCO), a homeodomain protein, and CUP-SHAPED COTYLEDON2 (CUC2), a NAC-domain protein. However, how their respective growth-repressive actions are integrated in space and time to generate complex leaf forms remains unknown. By using live imaging, we show that CUC2 and RCO are expressed in an interspersed fashion along the leaf margin, creating a distinctive striped pattern. We find that this pattern is functionally important because forcing RCO expression in the CUC2 domain disrupts auxin-based marginal patterning and can abolish leaflet formation. By combining genetic perturbations with time-lapse imaging and cellular growth quantifications, we provide evidence that RCO-mediated growth repression occurs after auxin-based leaflet patterning and in association with the repression of cell proliferation. Additionally, through the use of genetic mosaics, we show that RCO is sufficient to repress both cellular growth and proliferation in a cell-autonomous manner. This mechanism of growth repression is different to that of CUC2, which occurs in proliferating cells. Our findings clarify how the two growth repressors RCO and CUC2 coordinate to subdivide developing leaf primordia into distinct leaflets and generate the complex leaf form. They also indicate different relationships between growth repression and cell proliferation in the patterning and post-patterning stages of organogenesis.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Development and diversity of lignin patterns

Different patterns of lignified cell walls are associated with diverse functions in a variety of plant tissues. These functions rely on the stiffness and hydrophobicity that lignin polymers impart to the cell wall. The precise pattern of subcellular lignin deposition is critical for the structure-function relationship in each lignified cell type. Here, we describe the role of xylem vessels as water pipes, Casparian strips as apoplastic barriers, and the role of asymmetrically lignified endocarp b cells in exploding seed pods. We highlight similarities and differences in the genetic mechanisms underpinning local lignin deposition in these diverse cell types. By bringing together examples from different developmental contexts and different plant species, we propose that comparative approaches can benefit our understanding of lignin patterning mechanisms.

Arabis alpina: A perennial model plant for ecological genomics and life-history evolution

Many model organisms were chosen and achieved prominence because of an advantageous combination of their life‐history characteristics, genetic properties and also practical considerations. Discoveries made in Arabidopsis thaliana, the most renowned noncrop plant model species, have markedly stimulated studies in other species with different biology. Within the family Brassicaceae, the arctic–alpine Arabis alpina has become a model complementary to Arabidopsis thaliana to study the evolution of life‐history traits, such as perenniality, and ecological genomics in harsh environments. In this review, we provide an overview of the properties that facilitated the rapid emergence of A. alpina as a plant model. We summarize the evolutionary history of A. alpina, including genomic aspects, the diversification of its mating system and demographic properties, and we discuss recent progress in the molecular dissection of developmental traits that are related to its perennial life history and environmental adaptation. From this published knowledge, we derive open questions that might inspire future research in A. alpina, other Brassicaceae species or more distantly related plant families.

Light signals generated by vegetation shade facilitate acclimation to low light in shade-avoider plants

When growing in search for light, plants can experience continuous or occasional shading by other plants. Plant proximity causes a decrease in the ratio of R to far-red light (low R:FR) due to the preferential absorbance of R light and reflection of FR light by photosynthetic tissues of neighboring plants. This signal is often perceived before actual shading causes a reduction in photosynthetically active radiation (low PAR). Here, we investigated how several Brassicaceae species from different habitats respond to low R:FR and low PAR in terms of elongation, photosynthesis, and photoacclimation. Shade-tolerant plants such as hairy bittercress (Cardamine hirsuta) displayed a good adaptation to low PAR but a poor or null response to low R:FR exposure. In contrast, shade-avoider species, such as Arabidopsis (Arabidopsis thaliana), showed a weak photosynthetic performance under low PAR but they strongly elongated when exposed to low R:FR. These responses could be genetically uncoupled. Most interestingly, exposure to low R:FR of shade-avoider (but not shade-tolerant) plants improved their photoacclimation to low PAR by triggering changes in photosynthesis-related gene expression, pigment accumulation, and chloroplast ultrastructure. These results indicate that low R:FR signaling unleashes molecular, metabolic, and developmental responses that allow shade-avoider plants (including most crops) to adjust their photosynthetic capacity in anticipation of eventual shading by nearby plants.

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.

A PHABULOSA-Controlled Genetic Pathway Regulates Ground Tissue Patterning in the Arabidopsis Root

In both animals and plants, development involves anatomical modifications. In the root of Arabidopsis thaliana, maturation of the ground tissue (GT)—a tissue comprising all cells between epidermal and vascular ones—is a paradigmatic example of these modifications, as it generates an additional tissue layer, the middle cortex (MC).1, 2, 3, 4 In early post-embryonic phases, the Arabidopsis root GT is composed of one layer of endodermis and one of cortex. A second cortex layer, the MC, is generated by asymmetric cell divisions in about 80% of Arabidopsis primary roots, in a time window spanning from 7 to 14 days post-germination (dpg). The cell cycle regulator CYCLIN D6;1 (CYCD6;1) plays a central role in this process, as its accumulation in the endodermis triggers the formation of MC.⁵ The phytohormone gibberellin (GA) is a key regulator of the timing of MC formation, as alterations in its signaling and homeostasis result in precocious endodermal asymmetric cell divisions.^3,6,7 However, little is known on how GAs are regulated during GT maturation. Here, we show that the HOMEODOMAIN LEUCINE ZIPPER III (HD-ZIPIII) transcription factor PHABULOSA (PHB) is a master regulator of MC formation, controlling the accumulation of CYCD6;1 in the endodermis in a cell non-autonomous manner. We show that PHB activates the GA catabolic gene GIBBERELLIN 2 OXIDASE 2 (GA2ox2) in the vascular tissue, thus regulating the stability of the DELLA protein GIBBERELLIN INSENSITIVE (GAI)—a GA signaling repressor—in the root and, hence, CYCD6;1 expression in the endodermis.

•

PHB regulates cell non-autonomously the timing of MC formation

•

A time-dependent rise of PHB expression controls the CYCD6;1 switch in the GT

•

PHB regulates GAI stability modulating GA levels

•

PHB regulates root GA levels activating GA2ox2 expression in the vasculature

Adjustment of the PIF7-HFR1 transcriptional module activity controls plant shade adaptation

Shade caused by the proximity of neighboring vegetation triggers a set of acclimation responses to either avoid or tolerate shade. Comparative analyses between the shade-avoider Arabidopsis thaliana and the shade-tolerant Cardamine hirsuta revealed a role for the atypical basic-helix-loop-helix LONG HYPOCOTYL IN FR 1 (HFR1) in maintaining the shade tolerance in C. hirsuta, inhibiting hypocotyl elongation in shade and constraining expression profile of shade-induced genes. We showed that C. hirsuta HFR1 protein is more stable than its A. thaliana counterpart, likely due to its lower binding affinity to CONSTITUTIVE PHOTOMORPHOGENIC 1 (COP1), contributing to enhance its biological activity. The enhanced HFR1 total activity is accompanied by an attenuated PHYTOCHROME INTERACTING FACTOR (PIF) activity in C. hirsuta. As a result, the PIF-HFR1 module is differently balanced, causing a reduced PIF activity and attenuating other PIF-mediated responses such as warm temperature-induced hypocotyl elongation (thermomorphogenesis) and dark-induced senescence. By this mechanism and that of the already-known of phytochrome A photoreceptor, plants might ensure to properly adapt and thrive in habitats with disparate light amounts.

Genome-wide study of flowering-related MADS-box genes family in Cardamine hirsuta

MADS-box genes take part in diverse biological functions especially in development of reproductive structures and control of flowering time. Recently, Cardamine hirsuta has emerged as an exclusively powerful genetic system in comparative studies of development. Although the C. hirsuta genome sequence is available but a comprehensive analysis of its MADS-box family genes is still lacking. Here, we determined 50 Cardamine MADS-box genes through bioinformatics tools and classified them into 2 Mβ, 6 Mα and 2 Mγ and 40 MIKC-type (35 MIKCc and 5MIKC*) genes based on a phylogenetic analysis. The C. hirsuta MIKC subfamily could be further classified into 14 subgroups as Arabidopsis. However the number of MADS-box proteins was not equal among these subgroups. Based on the structural diversity among 50 MADS-box genes, 2 lineages were obtained, type I and type II. The lowest number of introns (0 or 1) was found in the Mα, Mβ, and Mγ groups of the type I genes. The most Cardamine MADS-box genes were randomly distributed on only three chromosomes. C. hirsuta had a relatively lower number of flowering MADS-box genes than A. thaliana and probably tandem duplication event resulted in the expansion of FLC, SQUA and TM3 family members in Arabidopsis. Moreover among the conserved motifs, ChMADS5 of SQUA, ChMADS34 of TM3 and ChMADS51 of AGL15 families had no K-domain. This study provides a basis for further functional investigation of MADS-box genes in C. hirsuta.

Outstanding questions in flower metabolism

The great diversity of flowers, their color, odor, taste, and shape, is mostly a result of the metabolic processes that occur in this reproductive organ when the flower and its tissues develop, grow, and finally die. Some of these metabolites serve to advertise flowers to animal pollinators, other confer protection towards abiotic stresses, and a large proportion of the molecules of the central metabolic pathways have bioenergetic and signaling functions that support growth and the transition to fruits and seeds. Although recent studies have advanced our general understanding of flower metabolism, several questions still await an answer. Here, we have compiled a list of open questions on flower metabolism encompassing molecular aspects, as well as topics of relevance for agriculture and the ecosystem. These questions include the study of flower metabolism through development, the biochemistry of nectar and its relevance to promoting plant-pollinator interaction, recycling of metabolic resources after flowers whiter and die, as well as the manipulation of flower metabolism by pathogens. We hope with this review to stimulate discussion on the topic of flower metabolism and set a reference point to return to in the future when assessing progress in the field.

Seed coat development in explosively dispersed seeds of Cardamine hirsuta

BACKGROUND AND AIMS

Seeds are dispersed by explosive coiling of the fruit valves in Cardamine hirsuta. This rapid coiling launches the small seeds on ballistic trajectories to spread over a 2 m radius around the parent plant. The seed surface interacts with both the coiling fruit valve during launch and subsequently with the air during flight. We aim to identify features of the seed surface that may contribute to these interactions by characterizing seed coat differentiation.

METHODS

Differentiation of the outermost seed coat layers from the outer integuments of the ovule involves dramatic cellular changes that we characterize in detail at the light and electron microscopical level including immunofluorescence and immunogold labelling.

KEY RESULTS

We found that the two outer integument (oi) layers of the seed coat contributed differently to the topography of the seed surface in the explosively dispersed seeds of C. hirsuta vs. the related species Arabidopsis thaliana where seed dispersal is non-explosive. The surface of A. thaliana seeds is shaped by the columella and the anticlinal cell walls of the epidermal oi2 layer. In contrast, the surface of C. hirsuta seeds is shaped by a network of prominent ridges formed by the anticlinal walls of asymmetrically thickened cells of the sub-epidermal oi1 layer, especially at the seed margin. Both the oi2 and oi1 cell layers in C. hirsuta seeds are characterized by specialized, pectin-rich cell walls that are deposited asymmetrically in the cell.

CONCLUSIONS

The two outermost seed coat layers in C. hirsuta have distinct properties: the sub-epidermal oi1 layer determines the topography of the seed surface, while the epidermal oi2 layer accumulates mucilage. These properties are influenced by polar deposition of distinct pectin polysaccharides in the cell wall. Although the ridged seed surface formed by oi1 cell walls is associated with ballistic dispersal in C. hirsuta, it is not restricted to explosively dispersed seeds in the Brassicaceae.

Plant science's next top models

BACKGROUND

Model organisms are at the core of life science research. Notable examples include the mouse as a model for humans, baker's yeast for eukaryotic unicellular life and simple genetics, or the enterobacteria phage λ in virology. Plant research was an exception to this rule, with researchers relying on a variety of non-model plants until the eventual adoption of Arabidopsis thaliana as primary plant model in the 1980s. This proved to be an unprecedented success, and several secondary plant models have since been established. Currently, we are experiencing another wave of expansion in the set of plant models.

SCOPE

Since the 2000s, new model plants have been established to study numerous aspects of plant biology, such as the evolution of land plants, grasses, invasive and parasitic plant life, adaptation to environmental challenges, and the development of morphological diversity. Concurrent with the establishment of new plant models, the advent of the 'omics' era in biology has led to a resurgence of the more complex non-model plants. With this review, we introduce some of the new and fascinating plant models, outline why they are interesting subjects to study, the questions they will help to answer, and the molecular tools that have been established and are available to researchers.

CONCLUSIONS

Understanding the molecular mechanisms underlying all aspects of plant biology can only be achieved with the adoption of a comprehensive set of models, each of which allows the assessment of at least one aspect of plant life. The model plants described here represent a step forward towards our goal to explore and comprehend the diversity of plant form and function. Still, several questions remain unanswered, but the constant development of novel technologies in molecular biology and bioinformatics is already paving the way for the next generation of plant models.

Variation in Below-to Aboveground Systemic Induction of Glucosinolates Mediates Plant Fitness Consequences under Herbivore Attack

Plants defend themselves against herbivore attack by constitutively producing toxic secondary metabolites, as well as by inducing them in response to herbivore feeding. Induction of secondary metabolites can cross plant tissue boundaries, such as from root to shoot. However, whether the potential for plants to systemically induce secondary metabolites from roots to shoots shows genetic variability, and thus, potentially, is under selection conferring fitness benefits to the plants is an open question. To address this question, we induced 26 maternal plant families of the wild species Cardamine hirsuta belowground (BG) using the wound-mimicking phytohormone jasmonic acid (JA). We measured resistance against a generalist (Spodoptera littoralis) and a specialist (Pieris brassicae) herbivore species, as well as the production of glucosinolates (GSLs) in plants. We showed that BG induction increased AG resistance against the generalist but not against the specialist, and found substantial plant family-level variation for resistance and GSL induction. We further found that the systemic induction of several GSLs tempered the negative effects of herbivory on total seed set production. Using a widespread natural system, we thus confirm that BG to AG induction has a strong genetic component, and can be under positive selection by increasing plant fitness. We suggest that natural variation in systemic induction is in part dictated by allocation trade-offs between constitutive and inducible GSL production, as well as natural variation in AG and BG herbivore attack in nature.

Trait-dependent resemblance of the flowering phenology and floral morphology of the allopolyploid Cardamine flexuosa to those of the parental diploids in natural habitats

Allopolyploids possess complete sets of genomes derived from different parental species and exhibit a range of variation in various traits. Reproductive traits may play a key role in the reproductive isolation between allopolyploids and their parental species, thus affecting the thriving of allopolyploids. However, empirical data, especially in natural habitats, comparing reproductive trait variation between allopolyploids and their parental species remain rare. Here, we documented the flowering phenology and floral morphology of the allopolyploid wild plant Cardamine flexuosa and its diploid parents C. amara and C. hirsuta in their native range in Switzerland. The flowering of C. flexuosa started at an intermediate time compared with those of the parents and the flowering period of C. flexuosa overlapped with those of the parents. Cardamine flexuosa resembled C. hirsuta in the size of flowers and petals and the length/width ratio of petals, while it resembled C. amara in the length/width ratio of flowers. These results provide empirical evidence of the trait-dependent variation of allopolyploid phenotypes in natural habitats at the local scale. They also suggest that the variation in some reproductive traits in C. flexuosa is associated with self-fertilization. Therefore, it is helpful to consider the mating system in furthering the understanding of the processes that may have shaped trait variation in polyploids in nature.

CRISPR/Cas9-Mediated Mutagenesis of RCO in Cardamine hirsuta

The small crucifer Cardamine hirsuta bears complex leaves divided into leaflets. This is in contrast to its relative, the reference plant Arabidopsis thaliana, which has simple leaves. Comparative studies between these species provide attractive opportunities to study the diversification of form. Here, we report on the implementation of the CRISPR/Cas9 genome editing methodology in C. hirsuta and with it the generation of novel alleles in the RCO gene, which was previously shown to play a major role in the diversification of form between the two species. Thus, genome editing can now be deployed in C. hirsuta, thereby increasing its versatility as a model system to study gene function and evolution.

The story of promiscuous crucifers: origin and genome evolution of an invasive species, Cardamine occulta (Brassicaceae), and its relatives

BACKGROUND AND AIMS

Cardamine occulta (Brassicaceae) is an octoploid weedy species (2n = 8x = 64) originated in Eastern Asia. It has been introduced to other continents including Europe and considered to be an invasive species. Despite its wide distribution, the polyploid origin of C. occulta remained unexplored. The feasibility of comparative chromosome painting (CCP) in crucifers allowed us to elucidate the origin and genome evolution in Cardamine species. We aimed to investigate the genome structure of C. occulta in comparison with its tetraploid (2n = 4x = 32, C. kokaiensis and C. scutata) and octoploid (2n = 8x = 64, C. dentipetala) relatives.

METHODS

Genomic in situ hybridization (GISH) and large-scale CCP were applied to uncover the parental genomes and chromosome composition of the investigated Cardamine species.

KEY RESULTS

All investigated species descended from a common ancestral Cardamine genome (n = 8), structurally resembling the Ancestral Crucifer Karyotype (n = 8), but differentiated by a translocation between chromosomes AK6 and AK8. Allotetraploid C. scutata originated by hybridization between two diploid species, C. parviflora and C. amara (2n = 2x = 16). By contrast, C. kokaiensis has an autotetraploid origin from a parental genome related to C. parviflora. Interestingly, octoploid C. occulta probably originated through hybridization between the tetraploids C. scutata and C. kokaiensis. The octoploid genome of C. dentipetala probably originated from C. scutata via autopolyploidization. Except for five species-specific centromere repositionings and one pericentric inversion post-dating the polyploidization events, the parental subgenomes remained stable in the tetra- and octoploids.

CONCLUSIONS

Comparative genome structure, origin and evolutionary history was reconstructed in C. occulta and related species. For the first time, whole-genome cytogenomic maps were established for octoploid plants. Post-polyploid evolution in Asian Cardamine polyploids has not been associated with descending dysploidy and intergenomic rearrangements. The combination of different parental (sub)genomes adapted to distinct habitats provides an evolutionary advantage to newly formed polyploids by occupying new ecological niches.

Genome invasion by a hypomethylated satellite repeat in Australian crucifer Ballantinia antipoda

Repetitive sequences are ubiquitous components of all eukaryotic genomes. They contribute to genome evolution and the regulation of gene transcription. However, the uncontrolled activity of repetitive sequences can negatively affect genome functions and stability. Therefore, repetitive DNAs are embedded in a highly repressive heterochromatic environment in plant cell nuclei. Here, we analyzed the sequence, composition and the epigenetic makeup of peculiar non-pericentromeric heterochromatic segments in the genome of the Australian crucifer Ballantinia antipoda. By the combination of high throughput sequencing, graph-based clustering and cytogenetics, we found that the heterochromatic segments consist of a mixture of unique sequences and an A-T-rich 174 bp satellite repeat (BaSAT1). BaSAT1 occupies about 10% of the B. antipoda nuclear genome in >250 000 copies. Unlike many other highly repetitive sequences, BaSAT1 repeats are hypomethylated; this contrasts with the normal patterns of DNA methylation in the B. antipoda genome. Detailed analysis of several copies revealed that these non-methylated BaSAT1 repeats were also devoid of heterochromatic histone H3K9me2 methylation. However, the factors decisive for the methylation status of BaSAT1 repeats remain currently unknown. In summary, we show that even highly repetitive sequences can exist as hypomethylated in the plant nuclear genome.

A high-density genetic map and QTL mapping of leaf traits and glucosinolates in Barbarea vulgaris

Background

Barbarea vulgaris is a wild cruciferous plant and include two distinct types: the G- and P-types named after their glabrous and pubescent leaves, respectively. The types differ significantly in resistance to a range of insects and diseases as well as glucosinolates and other chemical defenses. A high-density linkage map was needed for further progress to be made in the molecular research of this plant.

Results

We performed restriction site-associated DNA sequencing (RAD-Seq) on an F₂ population generated from G- and P-type B. vulgaris. A total of 1545 SNP markers were mapped and ordered in eight linkage groups, which represents the highest density linkage map to date for the crucifer tribe Cardamineae. A total of 722 previously published genome contigs (50.2 Mb, 30% of the total length) can be anchored to this high density genetic map, an improvement compared to a previously published map (431 anchored contigs, 38.7 Mb, 23% of the assembly genome). Most of these (572 contigs, 31.2 Mb) were newly anchored to the map, representing a significant improvement. On the basis of the present high-density genetic map, 37 QTL were detected for eleven traits, each QTL explaining 2.9–71.3% of the phenotype variation. QTL of glucosinolates, leaf size and color traits were in most cases overlapping, possibly implying a functional connection.

Conclusions

This high-density linkage map and the QTL obtained in this study will be useful for further understanding of the genetic of the B. vulgaris and molecular basis of these traits, many of which are shared in the related crop watercress.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5769-z) contains supplementary material, which is available to authorized users.

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.

Fruit Development: Turning Sticks into Hearts

Fruit morphological diversity reflects the versatility of these angiosperm-specific structures, which facilitate plant progeny dispersal from their sessile parents. A recent study links regulatory changes in a key genetic network for fruit patterning with the origin of heart-shaped pods in Brassicaceae.

Genetic basis of plasticity in plants

The ability of an organism to change its phenotype in response to different environments, termed plasticity, is a particularly important characteristic to enable sessile plants to adapt to rapid changes in their surroundings. Plasticity is a quantitative trait that can provide a fitness advantage and mitigate negative effects due to environmental perturbations. Yet, its genetic basis is not fully understood. Alongside technological limitations, the main challenge in studying plasticity has been the selection of suitable approaches for quantification of phenotypic plasticity. Here, we propose a categorization of the existing quantitative measures of phenotypic plasticity into nominal and relative approaches. Moreover, we highlight the recent advances in the understanding of the genetic architecture underlying phenotypic plasticity in plants. We identify four pillars for future research to uncover the genetic basis of phenotypic plasticity, with emphasis on development of computational approaches and theories. These developments will allow us to perform specific experiments to validate the causal genes for plasticity and to discover their role in plant fitness and evolution.

Shade Avoidance and Neighbor Detection

Plants detect neighboring vegetation as potential competitors for resources. Vegetation proximity is perceived by changes in the red (R) to far-red (FR) ratio (R:FR) through the phytochrome photoreceptors. To face this challenge, many plants have evolved the strategy to avoid shade, displaying a series of responses known as the shade avoidance syndrome (SAS). The SAS responses have been mostly studied at the seedling stage, and cover hypocotyl elongation as well as cotyledon and primary leaf expansion. In adult stages, SAS responses include an increase in petiole elongation and a decrease in leaf expansion, and an increase in plant height. Thus, the analysis of these responses provides a valuable and simple way to study how vegetation proximity affects plant development in both seedlings and adult plants. Here we describe a simple protocol to simulate shade in the laboratory and to evaluate these responses. Overall, our protocol can be easily used to expand the set of SAS responses of plants at different stages of development.

Leaf development and evolution

Plant leaves are differentiated organs that arise sequentially from a population of pluripotent stem cells at the shoot apical meristem (SAM). There is substantial diversity in leaf shape, much of which depends on the size and arrangement of outgrowths at the leaf margin. These outgrowths are generated by a patterning mechanism similar to the phyllotactic processes producing organs at the SAM, which involves the transcription factors CUP-SHAPED COTYLEDON and the phytohormone auxin. In the leaf, this patterning mechanism creates sequential protrusions and indentations along the margin. The size, shape, and distribution of these protrusions also depend on the overall growth of the leaf lamina. Globally, growth is regulated by a complex genetic network controlling the distribution of cell proliferation and the timing of differentiation. Evolutionary changes in margin form arise from changes in two different classes of homeobox genes that modify the outcome of marginal patterning in diverse ways, and are under intense investigation.

Building the differences: a case for the ground tissue patterning in plants

A key question in biology is to understand how interspecies morphological diversities originate. Plant roots present a huge interspecific phenotypical variability, mostly because roots largely contribute to adaptation to different kinds of soils. One example is the interspecific cortex layer number variability, spanning from one to several. Here, we review the latest advances in the understanding of the mechanisms expanding and/or restricting cortical layer number in Arabidopsis thaliana and their involvement in cortex pattern variability among multi-cortical layered species such as Cardamine hirsuta or Oryza sativa.

LMI1 homeodomain protein regulates organ proportions by spatial modulation of endoreduplication

Here, Vuolo et al. investigated the mechanisms controlling the relative size of leaves compared with their lateral appendages (stipules). Using genetics, live imaging, and modeling, they demonstrate that the LATE MERISTEM IDENTITY1 (LMI1) homeodomain protein regulates stipule proportions via an endoreduplication-dependent trade-off that limits tissue size despite increasing cell growth.

How the interplay between cell- and tissue-level processes produces correctly proportioned organs is a key problem in biology. In plants, the relative size of leaves compared with their lateral appendages, called stipules, varies tremendously throughout development and evolution, yet relevant mechanisms remain unknown. Here we use genetics, live imaging, and modeling to show that in Arabidopsis leaves, the LATE MERISTEM IDENTITY1 (LMI1) homeodomain protein regulates stipule proportions via an endoreduplication-dependent trade-off that limits tissue size despite increasing cell growth. LM1 acts through directly activating the conserved mitosis blocker WEE1, which is sufficient to bypass the LMI1 requirement for leaf proportionality.

The role of APETALA1 in petal number robustness

Invariant floral forms are important for reproductive success and robust to natural perturbations. Petal number, for example, is invariant in Arabidopsis thaliana flowers. However, petal number varies in the closely related species Cardamine hirsuta, and the genetic basis for this difference between species is unknown. Here we show that divergence in the pleiotropic floral regulator APETALA1 (AP1) can account for the species-specific difference in petal number robustness. This large effect of AP1 is explained by epistatic interactions: A. thaliana AP1 confers robustness by masking the phenotypic expression of quantitative trait loci controlling petal number in C. hirsuta. We show that C. hirsuta AP1 fails to complement this function of A. thaliana AP1, conferring variable petal number, and that upstream regulatory regions of AP1 contribute to this divergence. Moreover, variable petal number is maintained in C. hirsuta despite sufficient standing genetic variation in natural accessions to produce plants with four-petalled flowers.

Transcription factor dosage: more or less sufficient for growth

Recent findings highlight three instances in which major aspects of plant development are controlled by dosage-dependent protein levels. In the shoot apical meristem the mobile transcription factor WUS displays an intricate function with respect to target regulation that involves WUS dosage, binding site affinity and protein dimerization. The size of the root meristem is controlled by dosage-dependent PLT protein activity. Recent identification of targets and feedbacks provide new insights and entry into possible mechanisms of dosage read-out. Finally, HD-ZIPIII dosage, enforced by a gradient of mobile miRNAs, presents a relatively unexplored case in the radial patterning of vasculature and ground tissue. We evaluate our current knowledge of these three examples and address molecular mechanisms of dosage translation where possible.

Agrobacterium-mediated floral dip transformation of the model polyploid species Arabidopsis kamchatica

Polyploidization has played an important role in the speciation and diversification of plant species. However, genetic analyses of polyploids are challenging because the vast majority of the model species are diploids. The allotetraploid Arabidopsis kamchatica, which originated through the hybridization of the diploid Arabidopsis halleri and Arabidopsis lyrata, is an emerging model system for studying various aspects of polyploidy. However, a transgenic method that allows the insertion of a gene of interest into A. kamchatica is still lacking. In this study, we investigated the early development of pistils in A. kamchatica and confirmed the formation of open pistils in young flower buds (stages 8-9), which is important for allowing Agrobacterium to access female reproductive tissues. We established a simple Agrobacterium-mediated floral dip transformation method to transform a gene of interest into A. kamchatica by dipping A. kamchatica inflorescences bearing many young flower buds into a 5% sucrose solution containing 0.05% Silwet L-77 and Agrobacterium harboring the gene of interest. We showed that a screenable marker comprising fluorescence-accumulating seed technology with green fluorescent protein was useful for screening the transgenic seeds of two accessions of A. kamchatica subsp. kamchatica and an accession of A. kamchatica subsp. kawasakiana.

Differential spatial distribution of miR165/6 determines variability in plant root anatomy

A clear example of interspecific variation is the number of root cortical layers in plants. The genetic mechanisms underlying this variability are poorly understood, partly because of the lack of a convenient model. Here, we demonstrate that Cardamine hirsuta, unlike Arabidopsis thaliana, has two cortical layers that are patterned during late embryogenesis. We show that a miR165/6-dependent distribution of the HOMEODOMAIN LEUCINE ZIPPER III (HD-ZIPIII) transcription factor PHABULOSA (PHB) controls this pattern. Our findings reveal that interspecies variation in miRNA distribution can determine differences in anatomy in plants.

MowJoe: a method for automated-high throughput dissected leaf phenotyping

BACKGROUND

Accurate and automated phenotyping of leaf images is necessary for high throughput studies of leaf form like genome-wide association analysis and other forms of quantitative trait locus mapping. Dissected leaves (also referred to as compound) that are subdivided into individual units are an attractive system to study diversification of form. However, there are only few software tools for their automated analysis. Thus, high-throughput image processing algorithms are needed that can partition these leaves in their phenotypically relevant units and calculate morphological features based on these units.

RESULTS

We have developed MowJoe, an image processing algorithm that dissects a dissected leaf into leaflets, petiolule, rachis and petioles. It employs image skeletonization to convert leaves into graphs, and thereafter applies algorithms operating on graph structures. This partitioning of a leaf allows the derivation of morphological features such as leaf size, or eccentricity of leaflets. Furthermore, MowJoe automatically places landmarks onto the terminal leaflet that can be used for further leaf shape analysis. It generates specific output files that can directly be imported into downstream shape analysis tools. We applied the algorithm to two accessions of Cardamine hirsuta and show that our features are able to robustly discriminate between these accessions.

CONCLUSION

MowJoe is a tool for the semi-automated, quantitative high throughput shape analysis of dissected leaf images. It provides the statistical power for the detection of the genetic basis of quantitative morphological variations.

Conservation vs divergence in LEAFY and APETALA1 functions between Arabidopsis thaliana and Cardamine hirsuta

A conserved genetic toolkit underlies the development of diverse floral forms among angiosperms. However, the degree of conservation vs divergence in the configuration of these gene regulatory networks is less clear. We addressed this question in a parallel genetic study between the closely related species Arabidopsis thaliana and Cardamine hirsuta. We identified leafy (lfy) and apetala1 (ap1) alleles in a mutant screen for floral regulators in C. hirsuta. C. hirsuta lfy mutants showed a complete homeotic conversion of flowers to leafy shoots, mimicking lfy ap1 double mutants in A. thaliana. Through genetic and molecular experiments, we showed that AP1 activation is fully dependent on LFY in C. hirsuta, by contrast to A. thaliana. Additionally, we found that LFY influences heteroblasty in C. hirsuta, such that loss or gain of LFY function affects its progression. Overexpression of UNUSUAL FLORAL ORGANS also alters C. hirsuta leaf shape in an LFY-dependent manner. We found that LFY and AP1 are conserved floral regulators that act nonredundantly in C. hirsuta, such that LFY has more obvious roles in floral and leaf development in C. hirsuta than in A. thaliana.

The Plant-Dependent Life Cycle of Thecaphora thlaspeos: A Smut Fungus Adapted to Brassicaceae

Smut fungi are globally distributed plant pathogens that infect agriculturally important crop plants such as maize or potato. To date, molecular studies on plant responses to smut fungi are challenging due to the genetic complexity of their host plants. Therefore, we set out to investigate the known smut fungus of Brassicaceae hosts, Thecaphora thlaspeos. T. thlaspeos infects different Brassicaceae plant species throughout Europe, including the perennial model plant Arabis alpina. In contrast to characterized smut fungi, mature and dry T. thlaspeos teliospores germinated only in the presence of a plant signal. An infectious filament emerges from the teliospore, which can proliferate as haploid filamentous cultures. Haploid filaments from opposite mating types mate, similar to sporidia of the model smut fungus Ustilago maydis. Consistently, the a and b mating locus genes are conserved. Infectious filaments can penetrate roots and aerial tissues of host plants, causing systemic colonization along the vasculature. Notably, we could show that T. thlaspeos also infects Arabidopsis thaliana. Exploiting the genetic resources of A. thaliana and Arabis alpina will allow us to characterize plant responses to smut infection in a comparative manner and, thereby, characterize factors for endophytic growth as well as smut fungi virulence in dicot plants.

Using mustard genomes to explore the genetic basis of evolutionary change

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

Field Guide to Plant Model Systems

For the past several decades, advances in plant development, physiology, cell biology, and genetics have relied heavily on the model (or reference) plant Arabidopsis thaliana. Arabidopsis resembles other plants, including crop plants, in many but by no means all respects. Study of Arabidopsis alone provides little information on the evolutionary history of plants, evolutionary differences between species, plants that survive in different environments, or plants that access nutrients and photosynthesize differently. Empowered by the availability of large-scale sequencing and new technologies for investigating gene function, many new plant models are being proposed and studied.

Cells, walls, and endless forms

A key question in biology is how the endless diversity of forms found in nature evolved. Understanding the cellular basis of this diversity has been aided by advances in non-model experimental systems, quantitative image analysis tools, and modeling approaches. Recent work in plants highlights the importance of cell wall and cuticle modifications for the emergence of diverse forms and functions. For example, explosive seed dispersal in Cardamine hirsuta depends on the asymmetric localization of lignified cell wall thickenings in the fruit valve. Similarly, the iridescence of Hibiscus trionum petals relies on regular striations formed by cuticular folds. Moreover, NAC transcription factors regulate the differentiation of lignified xylem vessels but also the water-conducting cells of moss that lack a lignified secondary cell wall, pointing to the origin of vascular systems. Other novel forms are associated with modified cell growth patterns, including oriented cell expansion or division, found in the long petal spurs of Aquilegia flowers, and the Sarracenia purpurea pitcher leaf, respectively. Another good example is the regulation of dissected leaf shape in C. hirsuta via local growth repression, controlled by the REDUCED COMPLEXITY HD-ZIP class I transcription factor. These studies in non-model species often reveal as much about fundamental processes of development as they do about the evolution of form.

Coupled enhancer and coding sequence evolution of a homeobox gene shaped leaf diversity

In this study, Vuolo et al. investigate the mechanisms underlying the genetic basis for morphological diversity in leaf shape. They show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base.

Here we investigate mechanisms underlying the diversification of biological forms using crucifer leaf shape as an example. We show that evolution of an enhancer element in the homeobox gene REDUCED COMPLEXITY (RCO) altered leaf shape by changing gene expression from the distal leaf blade to its base. A single amino acid substitution evolved together with this regulatory change, which reduced RCO protein stability, preventing pleiotropic effects caused by its altered gene expression. We detected hallmarks of positive selection in these evolved regulatory and coding sequence variants and showed that modulating RCO activity can improve plant physiological performance. Therefore, interplay between enhancer and coding sequence evolution created a potentially adaptive path for morphological evolution.

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

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

How diploidization turned a tetraploid into a pseudotriploid

PREMISE OF THE STUDY

Despite being highly fertile and occupying a large geographic region, the North American heartleaf bittercress (Cardamine cordifolia; Brassicaceae) has a puzzling triploid-like chromosome number (2n = 3x = 24). As most triploids are sterile, we embarked on a detailed analysis of the C. cordifolia genome to elucidate its origin and structure.

METHODS

Mitotic and meiotic chromosome complement of C. cordifolia was analyzed by comparative chromosome painting using chromosome-specific BAC contigs of Arabidopsis thaliana. Resulting chromosome patterns were documented by multicolor fluorescence microscopy and compared with known ancestral and extant Brassicaceae genomes.

KEY RESULTS

We discovered that C. cordifolia is not a triploid hybrid but a diploidized tetraploid with the prevalence of regular, diploid-like meiotic pairing. The ancestral tetraploid chromosome number (2n = 32) was reduced to a triploid-like number (2n = 24) through four terminal chromosome translocations.

CONCLUSIONS

The structure of the pseudotriploid C. cordifolia genome results from a stepwise diploidization process after whole-genome duplication. We showed that translocation-based descending dysploidy (from n = 16 to n = 12) was mediated by the formation of five new chromosomes. The genome of C. cordifolia represents the diploidization process in statu nascendi and provides valuable insights into mechanisms of postpolyploidy rediploidization in land plants. Our data further suggest that chromosome number alone does not need to be a reliable proxy of species' evolutionary past and that the same chromosome number may originate either by polyploidization (hybridization) or due to descending dysploidy.

Morphomechanical Innovation Drives Explosive Seed Dispersal

How mechanical and biological processes are coordinated across cells, tissues, and organs to produce complex traits is a key question in biology. Cardamine hirsuta, a relative of Arabidopsis thaliana, uses an explosive mechanism to disperse its seeds. We show that this trait evolved through morphomechanical innovations at different spatial scales. At the organ scale, tension within the fruit wall generates the elastic energy required for explosion. This tension is produced by differential contraction of fruit wall tissues through an active mechanism involving turgor pressure, cell geometry, and wall properties of the epidermis. Explosive release of this tension is controlled at the cellular scale by asymmetric lignin deposition within endocarp b cells-a striking pattern that is strictly associated with explosive pod shatter across the Brassicaceae plant family. By bridging these different scales, we present an integrated mechanism for explosive seed dispersal that links evolutionary novelty with complex trait innovation. VIDEO ABSTRACT.

Breaking the mold: understanding the evolution and development of lateral organs in diverse plant models

The formation of complex three-dimensional shape differs significantly between plants and animals due to the presence of the cell wall in the former, which prevents all cell migration. Instead, in lateral plant organs such as leaves or petals, shape is controlled by a series of developmental phases in which the organ acquires polarity, cells undergo proliferation, and, lastly, cells expand to their final shape and size. Although these processes were first described based on mutagenesis approaches in major model systems like Arabidopsis thaliana, further insight into their complexity is best provided by studies of natural variation in organ shape in alternative model systems that sample a broader range of plant form. Weaving together work from both forward and evolutionary genetics, this review focuses on how modification in polarity establishment, cell proliferation and cell expansion drives modifications in the fundamental lateral organ developmental program to create diversity in shape.

Cardamine hirsuta: a comparative view

Current advances in developmental genetics are increasingly underpinned by comparative approaches as more powerful experimental tools become available in non-model organisms. Cardamine hirsuta is related to the model plant Arabidopsis thaliana and comparisons between these two experimentally tractable species have advanced our understanding of development and diversity. The power of forward genetics to uncover new biology was evident in the isolation of REDUCED COMPLEXITY, a gene which is present in C. hirsuta but lost in A. thaliana, and shapes crucifer leaf diversity. Transferring two Knotted1-like homeobox genes between C. hirsuta and A. thaliana revealed a constraint imposed by pleiotropy on the evolutionary potential of cis regulatory change to modify leaf shape. FLOWERING LOCUS C was identified as a heterochronic gene that underlies natural leaf shape variation in C. hirsuta.