Alternate wiring of a KNOXI genetic network underlies differences in leaf development of A. thaliana and C. hirsuta

Leaf form diversity and evolution: a never-ending story in plant biology

Leaf form can vary at different levels, such as inter/intraspecies, and diverse leaf shapes reflect their remarkable ability to adapt to various environmental conditions. Over the past two decades, considerable progress has been made in unraveling the molecular mechanisms underlying leaf form diversity, particularly the regulatory mechanisms of leaf complexity. However, the mechanisms identified thus far are only part of the entire process, and numerous questions remain unanswered. This review aims to provide an overview of the current understanding of the molecular mechanisms driving leaf form diversity while highlighting the existing gaps in our knowledge. By focusing on the unanswered questions, this review aims to shed light on areas that require further research, ultimately fostering a more comprehensive understanding of leaf form diversity.

A CUC1/auxin genetic module links cell polarity to patterned tissue growth and leaf shape diversity in crucifer plants

How tissue-level information encoded by fields of regulatory gene activity is translated into the patterns of cell polarity and growth that generate the diverse shapes of different species remains poorly understood. Here, we investigate this problem in the case of leaf shape differences between Arabidopsis thaliana, which has simple leaves, and its relative Cardamine hirsuta that has complex leaves divided into leaflets. We show that patterned expression of the transcription factor CUP-SHAPED COTYLEDON1 in C. hirsuta (ChCUC1) is a key determinant of leaf shape differences between the two species. Through inducible genetic perturbations, time-lapse imaging of growth, and computational modeling, we find that ChCUC1 provides instructive input into auxin-based leaf margin patterning. This input arises via transcriptional regulation of multiple auxin homeostasis components, including direct activation of WAG kinases that are known to regulate the polarity of PIN-FORMED auxin transporters. Thus, we have uncovered a mechanism that bridges biological scales by linking spatially distributed and species-specific transcription factor expression to cell-level polarity and growth, to shape diverse leaf forms.

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.

GRAS transcription factor PINNATE-LIKE PENTAFOLIATA2 controls compound leaf morphogenesis in Medicago truncatula

The milestone of compound leaf development is the generation of separate leaflet primordia during the early stages, which involves two linked but distinct morphogenetic events: leaflet initiation and boundary establishment for leaflet separation. Although some progress in understanding the regulatory pathways for each event have been made, it is unclear how they are intrinsically coordinated. Here, we identify the PINNATE-LIKE PENTAFOLIATA2 (PINNA2) gene encoding a newly identified GRAS transcription factor in Medicago truncatula. PINNA2 transcripts are preferentially detected at organ boundaries. Its loss-of-function mutations convert trifoliate leaves into a pinnate pentafoliate pattern. PINNA2 directly binds to the promoter region of the LEAFY orthologue SINGLE LEAFLET1 (SGL1), which encodes a key positive regulator of leaflet initiation, and downregulates its expression. Further analysis revealed that PINNA2 synergizes with two other repressors of SGL1 expression, the BEL1-like homeodomain protein PINNA1 and the C2H2 zinc finger protein PALMATE-LIKE PENTAFOLIATA1 (PALM1), to precisely define the spatiotemporal expression of SGL1 in compound leaf primordia, thereby maintaining a proper pattern of leaflet initiation. Moreover, we showed that the enriched expression of PINNA2 at the leaflet-to-leaflet boundaries is positively regulated by the boundary-specific gene MtNAM, which is essential for leaflet boundary formation. Together, these results unveil a pivotal role of the boundary-expressed transcription factor PINNA2 in regulating leaflet initiation, providing molecular insights into the coordination of intricate developmental processes underlying compound leaf pattern formation.

Overview of molecular mechanisms of plant leaf development: a systematic review

Leaf growth initiates in the peripheral region of the meristem at the apex of the stem, eventually forming flat structures. Leaves are pivotal organs in plants, serving as the primary sites for photosynthesis, respiration, and transpiration. Their development is intricately governed by complex regulatory networks. Leaf development encompasses five processes: the leaf primordium initiation, the leaf polarity establishment, leaf size expansion, shaping of leaf, and leaf senescence. The leaf primordia starts from the side of the growth cone at the apex of the stem. Under the precise regulation of a series of genes, the leaf primordia establishes adaxial-abaxial axes, proximal-distal axes and medio-lateral axes polarity, guides the primordia cells to divide and differentiate in a specific direction, and finally develops into leaves of a certain shape and size. Leaf senescence is a kind of programmed cell death that occurs in plants, and as it is the last stage of leaf development. Each of these processes is meticulously coordinated through the intricate interplay among transcriptional regulatory factors, microRNAs, and plant hormones. This review is dedicated to examining the regulatory influences of major regulatory factors and plant hormones on these five developmental aspects of leaves.

Control of compound leaf patterning by MULTI-PINNATE LEAF1 (MPL1) in chickpea

Plant lateral organs are often elaborated through repetitive formation of developmental units, which progress robustly in predetermined patterns along their axes. Leaflets in compound leaves provide an example of such units that are generated sequentially along the longitudinal axis, in species-specific patterns. In this context, we explored the molecular mechanisms underlying an acropetal mode of leaflet initiation in chickpea pinnate compound leaf patterning. By analyzing naturally occurring mutants multi-pinnate leaf1 (mpl1) that develop higher-ordered pinnate leaves with more than forty leaflets, we show that MPL1 encoding a C2H2-zinc finger protein sculpts a morphogenetic gradient along the proximodistal axis of the early leaf primordium, thereby conferring the acropetal leaflet formation. This is achieved by defining the spatiotemporal expression pattern of CaLEAFY, a key regulator of leaflet initiation, and also perhaps by modulating the auxin signaling pathway. Our work provides novel molecular insights into the sequential progression of leaflet formation.

The role of Class Ⅱ KNOX family in controlling compound leaf patterning in Medicago truncatula

Compound leaf development requires the coordination of genetic factors, hormones, and other signals. In this study, we explored the functions of Class Ⅱ KNOTTED-like homeobox (KNOXII) genes in the model leguminous plant Medicago truncatula. Phenotypic and genetic analyses suggest that MtKNOX4, 5 are able to repress leaflet formation, while MtKNOX3, 9, 10 are not involved in this developmental process. Further investigations have shown that MtKNOX4 represses the CK signal transduction, which is downstream of MtKNOXⅠ-mediated CK biosynthesis. Additionally, two boundary genes, FUSED COMPOUND LEAF1 (orthologue of Arabidopsis Class M KNOX) and NO APICAL MERISTEM (orthologue of Arabidopsis CUP-SHAPED COTYLEDON), are necessary for MtKNOX4-mediated compound leaf formation. These findings suggest, that among the members of MtKNOXⅡ, MtKNOX4 plays a crucial role in integrating the CK pathway and boundary regulators, providing new insights into the roles of MtKNOXⅡ in regulating the elaboration of compound leaves in M. truncatula.

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.

SHOOT MERISTEMLESS participates in the heterophylly of Hygrophila difformis (Acanthaceae)

In heterophyllous plants, leaf shape shows remarkable plasticity in response to environmental conditions. However, transgenic studies of heterophylly are lacking and the molecular mechanism remains unclear. Here, we cloned the KNOTTED1-LIKE HOMEOBOX family gene SHOOT MERISTEMLESS (STM) from the heterophyllous plant Hygrophila difformis (Acanthaceae). We used molecular, morphogenetic, and biochemical tools to explore its functions in heterophylly. HdSTM was detected in different organs of H. difformis, and its expression changed with environmental conditions. Heterologous, ectopic expression of HdSTM in Arabidopsis (Arabidopsis thaliana) increased leaf complexity and CUP-SHAPED COTYLEDON (CUC) transcript levels. However, overexpression of HdSTM in H. difformis did not induce the drastic leaf change in the terrestrial condition. Overexpression of HdSTM in H. difformis induced quick leaf variations in submergence, while knockdown of HdSTM led to disturbed leaf development and weakened heterophylly in H. difformis. HdCUC3 had the same spatiotemporal expression pattern as HdSTM. Biochemical analysis revealed a physical interaction between HdSTM and HdCUC3. Our results provide genetic evidence that HdSTM is involved in regulating heterophylly in H. difformis.

The cellular basis for synergy between RCO and KNOX1 homeobox genes in leaf shape diversity

Leaves of seed plants provide an attractive system to study the development and evolution of form. Leaves show varying degrees of margin complexity ranging from simple, as in Arabidopsis thaliana, to fully dissected into leaflets in the closely related species Cardamine hirsuta. Leaflet formation requires actions of Class I KNOTTED1-LIKE HOMEOBOX (KNOX1) and REDUCED COMPLEXITY (RCO) homeobox genes, which are expressed in the leaves of C. hirsuta but not A. thaliana. Evolutionary studies indicate that diversification of KNOX1 and RCO genes was repeatedly associated with increased leaf complexity. However, whether this gene combination represents a developmentally favored avenue for leaflet formation remains unknown, and the cell-level events through which the combined action of these genes drives leaflet formation are also poorly understood. Here we show, through a genetic screen, that when a C. hirsuta RCO transgene is expressed in A. thaliana, then ectopic KNOX1 expression in leaves represents a preferred developmental path for leaflet formation. Using time-lapse growth analysis, we demonstrate that KNOX1 expression in the basal domain of leaves leads to prolonged and anisotropic cell growth. This KNOX1 action, in synergy with local growth repression by RCO, is instrumental in generating rachises and petiolules, the linear geometrical elements, that bear leaflets in complex leaves. Our results show how the combination of cell-level growth analyses and genetics can help us understand how evolutionary modifications in expression of developmentally important genes are translated into diverse leaf shapes.

Developmental Analysis of Compound Leaf Development in Arachis hypogaea

Leaves are the primary photosynthetic structures, while photosynthesis is the direct motivation of crop yield formation. As a legume plant, peanut (Arachis hypogaea) is one of the most economically essential crops as well as an important source of edible oil and protein. The leaves of A. hypogaea are in the tetrafoliate form, which is different from the trifoliate leaf pattern of Medicago truncatula, a model legume species. In A. hypogaea, an even-pinnate leaf with a pair of proximal and distal leaflets was developed; however, only a single terminal leaflet and a pair of lateral leaflets were formed in the odd-pinnate leaf in M. truncatula. In this study, the development of compound leaf in A. hypogaea was investigated. Transcriptomic profiles revealed that the common and unique differentially expressed genes were identified in a proximal leaflet and a distal leaflet, which provided a research route to understand the leaf development in A. hypogaea. Then, a naturally occurring mutant line with leaf developmental defects in A. hypogaea was obtained, which displayed a pentafoliate form with an extra terminal leaflet. The characterization of the mutant indicated that cytokinin and class I KNOTTED-LIKE HOMEOBOX were involved in the control of compound leaf pattern in A. hypogaea. These results expand our knowledge and provide insights into the molecular mechanism underlying the formation of different compound leaf patterns among species.

The Genetic Control of the Compound Leaf Patterning in Medicago truncatula

Simple and compound which are the two basic types of leaves are distinguished by the pattern of the distribution of blades on the petiole. Compared to simple leaves comprising a single blade, compound leaves have multiple blade units and exhibit more complex and diverse patterns of organ organization, and the molecular mechanisms underlying their pattern formation are receiving more and more attention in recent years. Studies in model legume Medicago truncatula have led to an improved understanding of the genetic control of the compound leaf patterning. This review is an attempt to summarize the current knowledge about the compound leaf morphogenesis of M. truncatula, with a focus on the molecular mechanisms involved in pattern formation. It also includes some comparisons of the molecular mechanisms between leaf morphogenesis of different model species and offers useful information for the molecular design of legume crops.

Integration of Core Mechanisms Underlying Plant Aerial Architecture

Over the last decade or so important progress has been made in identifying and understanding a set of patterning mechanisms that have the potential to explain many aspects of plant morphology. These include the feedback loop between mechanical stresses and interphase microtubules, the regulation of plant cell polarity and the role of adaxial and abaxial cell type boundaries. What is perhaps most intriguing is how these mechanisms integrate in a combinatorial manner that provides a means to generate a large variety of commonly seen plant morphologies. Here, I review our current understanding of these mechanisms and discuss the links between them.

From genes to networks: The genetic control of leaf development

Substantial diversity exists for both the size and shape of the leaf, the main photosynthetic organ of flowering plants. The two major forms of leaf are simple leaves, in which the leaf blade is undivided, and compound leaves, which comprise several leaflets. Leaves form at the shoot apical meristem from a group of undifferentiated cells, which first establish polarity, then grow and differentiate. Each of these processes is controlled by a combination of transcriptional regulators, microRNAs and phytohormones. The present review documents recent advances in our understanding of how these various factors modulate the development of both simple leaves (focusing mainly on the model plant Arabidopsis thaliana) and compound leaves (focusing mainly on the model legume species Medicago truncatula).

Leaf Shape Diversity: From Genetic Modules to Computational Models

Plant leaves display considerable variation in shape. Here, we introduce key aspects of leaf development, focusing on the morphogenetic basis of leaf shape diversity. We discuss the importance of the genetic control of the amount, duration, and direction of cellular growth for the emergence of leaf form. We highlight how the combined use of live imaging and computational frameworks can help conceptualize how regulated cellular growth is translated into different leaf shapes. In particular, we focus on the morphogenetic differences between simple and complex leaves and how carnivorous plants form three-dimensional insect traps. We discuss how evolution has shaped leaf diversity in the case of complex leaves, by tinkering with organ-wide growth and local growth repression, and in carnivorous plants, by modifying the relative growth of the lower and upper sides of the leaf primordium to create insect-digesting traps.

GhARF16-1 modulates leaf development by transcriptionally regulating the GhKNOX2-1 gene in cotton

The leaf is a crucial organ evolved with remarkable morphological diversity to maximize plant photosynthesis. The leaf shape is a key trait that affects photosynthesis, flowering rates, disease resistance and yield. Although many genes regulating leaf development have been identified in the past years, the precise regulatory architecture underlying the generation of diverse leaf shapes remains to be elucidated. We used cotton as a reference model to probe the genetic framework underlying divergent leaf forms. Comparative transcriptome analysis revealed that the GhARF16-1 and GhKNOX2-1 genes might be potential regulators of leaf shape. We functionally characterized the auxin-responsive factor ARF16-1 acting upstream of GhKNOX2-1 to determine leaf morphology in cotton. The transcription of GhARF16-1 was significantly higher in lobed-leaved cotton than in smooth-leaved cotton. Furthermore, the overexpression of GhARF16-1 led to the up-regulation of GhKNOX2-1 and resulted in more and deeper serrations in cotton leaves, similar to the leaf shape of cotton plants overexpressing GhKNOX2-1. We found that GhARF16-1 specifically bound to the promoter of GhKNOX2-1 to induce its expression. The heterologous expression of GhARF16-1 and GhKNOX2-1 in Arabidopsis led to lobed and curly leaves, and a genetic analysis revealed that GhKNOX2-1 is epistatic to GhARF16-1 in Arabidopsis, suggesting that the GhARF16-1 and GhKNOX2-1 interaction paradigm also functions to regulate leaf shape in Arabidopsis. To our knowledge, our results uncover a novel mechanism by which auxin, through the key component ARF16-1 and its downstream-activated gene KNOX2-1, determines leaf morphology in eudicots.

Coordination of differentiation rate and local patterning in compound-leaf development

The variability in leaf form in nature is immense. Leaf patterning occurs by differential growth, taking place during a limited window of morphogenetic activity at the leaf marginal meristem. While many regulators have been implicated in the designation of the morphogenetic window and in leaf patterning, how these effectors interact to generate a particular form is still not well understood. We investigated the interaction among different effectors of tomato (Solanum lycopersicum) compound-leaf development, using genetic and molecular analyses. Mutations in the tomato auxin response factor SlARF5/SlMP, which normally promotes leaflet formation, suppressed the increased leaf complexity of mutants with extended morphogenetic window. Impaired activity of the NAC/CUC transcription factor GOBLET (GOB), which specifies leaflet boundaries, also reduced leaf complexity in these backgrounds. Analysis of genetic interactions showed that the patterning factors SlMP, GOB and the MYB transcription factor LYRATE (LYR) coordinately regulate leaf patterning by modulating in parallel different aspects of leaflet formation and shaping. This work places an array of developmental regulators in a morphogenetic context. It reveals how organ-level differentiation rate and local growth are coordinated to sculpture an organ. These concepts are applicable to the coordination of pattering and differentiation in other species and developmental processes.

TCP5 controls leaf margin development by regulating KNOX and BEL-like transcription factors in Arabidopsis

Development of leaf margins is an important process in leaf morphogenesis. CIN-clade TCP (TEOSINTE BRANCHED1/CYCLOIDEA/PCF) transcription factors are known to have redundant roles in specifying leaf margins, but the specific mechanisms through which individual TCP genes function remain elusive. In this study, we report that the CIN-TCP gene TCP5 is involved in repressing the initiation and outgrowth of leaf serrations by activating two key regulators of margin development, the Class II KNOX factor KNAT3 and BEL-like SAW1. Specifically, TCP5 directly promotes the transcription of KNAT3 and indirectly activates the expression of SAW1. We also show that TCP5 regulates KNAT3 and SAW1 in a temporal- and spatial- specific manner that is largely in accordance with the progress of formation of serrations. This regulation might serve as a key mechanism in patterning margin morphogenesis and in sculpting the final form of the leaf.

SAW homeodomain transcription factors regulate initiation of leaf margin serrations

Plant leaves are the main photosynthetic organ of plants and they occur in an array of different shapes. Leaf shape is determined by morphogenesis whereby patterning of the leaf margin can result in interspaced leaf serrations, lobes, or leaflets, depending on the species, developmental stage, and in some instances the environment. In Arabidopsis, mutations in the homeodomain transcription factors SAW1 and SAW2 result in more prominent leaf margin serrations. Here we show that serrations appear precociously in the saw1 saw2 mutant. The pattern of auxin maxima, and of PIN1 and CUC2 expression, which form a feedback loop that drives serration outgrowth, is altered in saw1 saw2 and correlates with precocious serration initiation. SAW1 is not expressed in the outer epidermal cell layer where PIN1 convergence points generate auxin maxima. Instead, SAW1 is expressed on the adaxial side of the leaf and expression in this domain is sufficient for function. We suggest that SAW1 and SAW2 repress serration initiation and outgrowth by promoting the transition to a determinate fate in the leaf margin.

Upregulation of a KN1 homolog by transposon insertion promotes leafy head development in lettuce

Leafy head is a unique type of plant architecture found in some vegetable crops, with leaves bending inward to form a compact head. The genetic and molecular mechanisms underlying leafy head in vegetables remain poorly understood. We genetically fine-mapped and cloned a major quantitative trait locus controlling heading in lettuce. The candidate gene (LsKN1) is a homolog of knotted 1 (KN1) from Zea mays Complementation and CRISPR/Cas9 knockout experiments confirmed the role of LsKN1 in heading. In heading lettuce, there is a CACTA-like transposon inserted into the first exon of LsKN1 (LsKN1▽). The transposon sequences act as a promoter rather than an enhancer and drive high expression of LsKN1▽. The enhanced expression of LsKN1▽ is necessary but not sufficient for heading in lettuce. Data from ChIP-sequencing, electrophoretic mobility shift assays, and dual luciferase assays indicate that the LsKN1▽ protein binds the promoter of LsAS1 and down-regulates its expression to alter leaf dorsoventrality. This study provides insight into plant leaf development and will be useful for studies on heading in other vegetable crops.

Characterization of the KNOTTED1-LIKE HOMEOBOX (KNOX) gene family in Pinus pinaster Ait

KNOTTED1-LIKE HOMEOBOX (KNOX) genes are a family of plant-specific homeobox transcription factors with important roles in plant development that have been classified into two subfamilies with differential expression domains and functions. Studies in angiosperms have shown that class I members are related to the maintenance of meristem homeostasis and leaf development, whereas class II members promote differentiation of tissues and organs. However, little is known about its diversification and function in gymnosperms. By combining PCR-based detection and transcriptome data analysis, we identified four class I and two class II KNOX genes in Pinus pinaster. Expression analyses showed that class I members were mainly expressed in meristematic regions and differentiating tissues, with practically no expression in lateral organs, whereas expression of class II members was restricted to lateral organs. Furthermore, overexpression of P. pinaster KNOX genes in Arabidopsis thaliana caused similar phenotypic effects to those described for their angiosperms counterparts. This is the first time to our knowledge that functional analyses of class II members are reported in a conifer species. These results suggest a high conservation of the KNOX gene family throughout seed plants, as the functional differentiation of both subfamilies observed in angiosperms might be partially conserved in gymnosperms.

Comparative transcriptome and metabolomic profiling reveal the complex mechanisms underlying the developmental dynamics of tobacco leaves

Although the leaf is the most important photosynthetic organ in most plants, many of the molecular mechanisms underlying leaf developmental dynamics remain to be explored. To better understand the transcriptional regulatory mechanisms involved in leaf development, we conducted comparative transcriptomic and metabolomic analysis of leaves from seven positions on tobacco (Nicotiana tabacum) plants. A total of 35,622 unique differentially expressed genes and 79 metabolites were identified. A time-series expression analysis detected two interesting transcriptional profiles, one comprising 10,197 genes that displayed continual up-regulation during leaf development and another comprising 4696 genes that displayed continual down-regulation. Combining these data with co-expression network results identified four important regulatory networks involved in photorespiration and the tricarboxylic acid cycle; these networks may regulate carbon/nitrogen balance during leaf development. We also found that the transcription factor NtGATA5 acts as a hub associated with C and N metabolism and chloroplast development during leaf development through regulation of phytohormones. Furthermore, we investigated the transcriptional dynamics of genes involved in the auxin, cytokinin, and jasmonic acid biosynthesis and signaling pathways during tobacco leaf development. Overall, our study greatly expands the understanding of the regulatory network controlling developmental dynamics in plant leaves.

Patterning a Leaf by Establishing Polarities

Leaves are the major organ for photosynthesis in most land plants, and leaf structure is optimized for the maximum capture of sunlight and gas exchange. Three polarity axes, the adaxial-abaxial axis, the proximal-distal axis, and the medial-lateral axis are established during leaf development to give rise to a flattened lamina with a large area for photosynthesis and blades that are extended on petioles for maximum sunlight. Adaxial cells are elongated, tightly packed cells with many chloroplasts, and their fate is specified by HD-ZIP III and related factors. Abaxial cells are rounder and loosely packed cells and their fate is established and maintained by YABBY family and KANADI family proteins. The activities of adaxial and abaxial regulators are coordinated by ASYMMETRIC LEAVES2 and auxin. Establishment of the proximodistal axis involves the BTB/POZ domain proteins BLADE-ON-PETIOLE1 and 2, whereas homeobox genes PRESSED FLOWER and WUSCHEL-RELATED HOMEOBOX1 mediate leaf development along the mediolateral axis. This review summarizes recent advances in leaf polarity establishment with a focus on the regulatory networks involved.

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.

Evolution: How a Homeobox Gene Cuts the Mustard Leaf

Gene duplication and cis-regulatory divergence created the growth-repressing RCO homeobox gene and facilitated evolution of dissected crucifer leaves. Identification of RCO targets reveals that auto-repression evolved to fine-tune RCO activity and that RCO dissects leaves by increasing cytokinin signalling to inhibit growth locally.

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.

Transcriptional Landscapes of Divergent Sporophyte Development in Two Mosses, Physcomitrium (Physcomitrella) patens and Funaria hygrometrica

Understanding the molecular basis of morphological shifts is a fundamental question of evolutionary biology. New morphologies may arise through the birth/death of genes (gene gain/loss) or by reutilizing existing gene sets. Yet, the relative contribution of these two processes to radical morphological shifts is still poorly understood. Here, we use the model system of two mosses, Funaria hygrometrica and Physcomitrium (Physcomitrella) patens, to investigate the molecular mechanisms underlying contrasting sporophyte architectures. We used comparative analysis of time-series expression data for four stages of sporophyte development in both species to address this question in detail. We found that large-scale differences in sporophytic architecture are mainly governed by orthologous (i.e., shared) genes frequently experiencing temporal gene expression shifts between the two species. While the absolute number of species-specific genes expressed during sporophyte development is somewhat smaller, we observed a significant increase of their proportion in preferentially sporophyte expressed genes, suggesting a fundamental role in the sporophyte phase. However, further functional studies are necessary to determine their contribution to diverging sporophyte morphologies. Our results add to the growing set of studies suggesting that radical changes in morphology may rely on the heterochronic expression of conserved regulators.

Overexpression of a Novel LcKNOX Transcription Factor from Liriodendron chinense Induces Lobed Leaves in Arabidopsis thaliana

Liriodendron chinense is a common ornamental tree that has attractive leaves, which is a valuable trait for use in landscape architecture. In this work, we aimed to identify the potential genes that control and regulate the development of L. chinense leaf lobes. Based on the transcriptome data for the leaf developmental stages we previously generated, two candidate genes were identified in this study. KNOTTED-LIKE HOMEOBOX(KNOX), encoding homeobox family proteins, play a large role in leaf lobe and leaf complexity regulation. Here, two full length KNOX genes from L. chinense were amplified and named LcKNOX1 and LcKNOX6 according to their sequence similarities with the respective Arabidopsis thaliana KNOX family genes. Overexpression vectors were constructed and subsequently transformed into wild type (WT) A. thaliana. Additionally, LcKNOX6 was expressed in tobacco leaves to examine its subcellular localization, and the 35S::LcKNOX6 transgenic A. thaliana leaf cells were imaged with the use of SEM. The expression of several genes that participate in KNOX gene regulation were validated by quantitative real-time PCR. The results show that LcKNOX1 produces almost the same phenotype as that found in WT A. thaliana. Notably, the LcKNOX6-1 lines presented deep leaf lobes that were similar to L. chinense leaf lobes. Two 35S::LcKNOX6 lines induced an abnormal growth phenotype whose seeds were abortive. In short, these results indicate that the LcKNOX6 gene might affect leaf development in A. thaliana and provide insights into the regulation of L. chinense leaf shaping.

Autoregulation of RCO by Low-Affinity Binding Modulates Cytokinin Action and Shapes Leaf Diversity

Mechanisms through which the evolution of gene regulation causes morphological diversity are largely unclear. The tremendous shape variation among plant leaves offers attractive opportunities to address this question. In cruciferous plants, the REDUCED COMPLEXITY (RCO) homeodomain protein evolved via gene duplication and acquired a novel expression domain that contributed to leaf shape diversity. However, the molecular pathways through which RCO regulates leaf growth are unknown. A key question is to identify genome-wide transcriptional targets of RCO and the DNA sequences to which RCO binds. We investigate this question using Cardamine hirsuta, which has complex leaves, and its relative Arabidopsis thaliana, which evolved simple leaves through loss of RCO. We demonstrate that RCO directly regulates genes controlling homeostasis of the hormone cytokinin to repress growth at the leaf base. Elevating cytokinin signaling in the RCO expression domain is sufficient to both transform A. thaliana simple leaves into complex ones and partially bypass the requirement for RCO in C. hirsuta complex leaf development. We also identify RCO as its own target gene. RCO directly represses its own transcription via an array of low-affinity binding sites, which evolved after RCO duplicated from its progenitor sequence. This autorepression is required to limit RCO expression. Thus, evolution of low-affinity binding sites created a negative autoregulatory loop that facilitated leaf shape evolution by defining RCO expression and fine-tuning cytokinin activity. In summary, we identify a transcriptional mechanism through which conflicts between novelty and pleiotropy are resolved during evolution and lead to morphological differences between species.

SOBIR1/EVR prevents precocious initiation of fiber differentiation during wood development through a mechanism involving BP and ERECTA

In plants, secondary growth results in radial expansion of stems and roots, generating large amounts of biomass in the form of wood. Using genome-wide association studies (GWAS)-guided reverse genetics in Arabidopsis thaliana, we discovered SOBIR1/EVR, previously known to control plant immunoresponses and abscission, as a regulator of secondary growth. We present anatomical, genetic, and molecular evidence indicating that SOBIR1/EVR prevents the precocious differentiation of xylem fiber, a key cell type for wood development. SOBIR1/EVR acts through a mechanism that involves BREVIPEDICELLUS (BP) and ERECTA (ER), 2 proteins previously known to regulate xylem fiber development. We demonstrate that BP binds SOBIR1/EVR promoter and that SOBIR1/EVR expression is enhanced in bp mutants, suggesting a direct, negative regulation of BP over SOBIR1/EVR expression. We show that SOBIR1/EVR physically interacts with ER and that defects caused by the sobir1/evr mutation are aggravated by mutating ER, indicating that SOBIR1/EVR and ERECTA act together in the control of the precocious formation of xylem fiber development.

Growth and biomechanics of shoot organs

Plant organs arise through complex interactions between biological and physical factors that control morphogenesis. While there has been tremendous progress in the understanding of the genetics behind development, we know much less about how mechanical forces control growth in plants. In recent years, new multidisciplinary research combining genetics, live-imaging, physics, and computational modeling has begun to fill this gap by revealing the crucial role of biomechanics in the establishment of plant organs. In this review, we provide an overview of our current understanding of growth during initiation, patterning, and expansion of shoot lateral organs. We discuss how growth is controlled by physical forces, and how mechanical stresses generated during growth can control morphogenesis at the level of both cells and tissues. Understanding the mechanical basis of growth and morphogenesis in plants is in its early days, and many puzzling facts are yet to be deciphered.

A Growth-Based Framework for Leaf Shape Development and Diversity

How do genes modify cellular growth to create morphological diversity? We study this problem in two related plants with differently shaped leaves: Arabidopsis thaliana (simple leaf shape) and Cardamine hirsuta (complex shape with leaflets). We use live imaging, modeling, and genetics to deconstruct these organ-level differences into their cell-level constituents: growth amount, direction, and differentiation. We show that leaf shape depends on the interplay of two growth modes: a conserved organ-wide growth mode that reflects differentiation; and a local, directional mode that involves the patterning of growth foci along the leaf edge. Shape diversity results from the distinct effects of two homeobox genes on these growth modes: SHOOTMERISTEMLESS broadens organ-wide growth relative to edge-patterning, enabling leaflet emergence, while REDUCED COMPLEXITY inhibits growth locally around emerging leaflets, accentuating shape differences created by patterning. We demonstrate the predictivity of our findings by reconstructing key features of C. hirsuta leaf morphology in A. thaliana. VIDEO ABSTRACT.

Resolving the backbone of the Brassicaceae phylogeny for investigating trait diversity

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

Overexpression of the KNOX gene Tkn4 affects pollen development and confers sensitivity to gibberellin and auxin in tomato

The knotted1-like homeobox genes not only regulate the formation and differentiation of meristems and vascular system but are also involved in biosynthesis and signal transduction of diverse plant hormones in tomato. Here, we showed that a knotted1-like homeobox gene Tkn4 is required for pollen and pollen tube growth when this gene is overexpressed in tomato. Pollen grains in the Tkn4 overexpressed plants (Tkn4-OX) germinated quicker than those in the wild-type (WT) plant cultured in vitro in germination media. The percentage of fruit set was higher in Tkn4-OX than in WT plants and the transgenic plants showed an ordered inflorescence. Tkn4-OX seedlings also exhibited sensitivity to gibberellins (GA) and auxins. RNA sequencing results showed that the expression of genes related to sugar, cell wall-modification, microtubule-associated vesicular transport for pollen growth, GA and auxin synthesis were significantly changed. Hence, Tkn4 contributes to a function in the development of pollen and pollen tube and the regulation of phytohormones to participate in plant growth. These results provided a potential application value for agricultural improvement to enhance the rate of fruit set in tomato.

LMI1-like and KNOX1 genes coordinately regulate plant leaf development in dicotyledons

This report reveals that the LMI1-like and KNOX1 genes coordinately control the leaf development and different combinations of those genes which produce diverse leaf shapes including broad, lobed and compound leaves. Class I KNOTTED1-like homeobox (KNOX1) genes are involved in compound leaf development and are repressed by the ASYMMETRIC LEAVES1 (AS1)-AS2 complex. Cotton plants have a variety of leaf shapes, including broad leaves and lobed leaves. GhOKRA, a LATE MERISTEM IDENTITY 1 (LMI1)-like gene, controls the development of an okra leaf shape. We cloned the corresponding cotton homologs of Arabidopsis thaliana AS1 and AS2 and seven KNOX1 genes. Through virus-induced gene silencing technology, we found that either GhAS1 or GhAS2-silenced cotton plants showed a great change in leaf shape from okra leaves to trifoliolate dissected leaves. In the shoot tips of these plants, the expression of the cotton ortholog of Knotted in A. thaliana 1 (KNAT1), GhKNOTTED1-LIKE2/3/4 (GhKNL2/3/4), was increased. However, GhKNOX1s-silenced plants maintained the wild-type okra leaves. A novel dissected-like leaf in A. thaliana was further generated by crossing plants constitutively expressing GhOKRA with either as1-101 or as2-101 mutant plants. The dissected-like leaves showed two different leaf vein patterns. This report reveals that the LMI1-like and KNOX1 genes coordinately control leaf development, and different combinations of these genes produce diverse leaf shapes including broad leaves, lobed leaves and compound leaves. This is the first report on the artificial generation of compound leaves from simple leaves in cotton.

Regulatory Diversification of INDEHISCENT in the Capsella Genus Directs Variation in Fruit Morphology

Evolution of gene-regulatory sequences is considered the primary driver of morphological variation [1-3]. In animals, the diversity of body plans between distantly related phyla is due to the differential expression patterns of conserved "toolkit" genes [4]. In plants, variation in expression domains similarly underlie most of the reported diversity of organ shape both in natural evolution and in the domestication of crops [5-9]. The heart-shaped fruit from members of the Capsella genus is a morphological novelty that has evolved after Capsella diverged from Arabidopsis ∼8 mya [10]. Comparative studies of fruit growth in Capsella and Arabidopsis revealed that the difference in shape is caused by local control of anisotropic growth [11]. Here, we show that sequence variation in regulatory domains of the fruit-tissue identity gene, INDEHISCENT (IND), is responsible for expansion of its expression domain in the heart-shaped fruits from Capsella rubella. We demonstrate that expression of this CrIND gene in the apical part of the valves in Capsella contributes to the heart-shaped appearance. While studies on morphological diversity have revealed the importance of cis-regulatory sequence evolution, few examples exist where the downstream effects of such variation have been characterized in detail. We describe here how CrIND exerts its function on Capsella fruit shape by binding sequence elements of auxin biosynthesis genes to activate their expression and ensure auxin accumulation into highly localized maxima in the fruit valves. Thus, our data provide a direct link between changes in expression pattern and altered hormone homeostasis in the evolution of morphological novelty.

Genetic control of compound leaf development in the mungbean (Vigna radiata L.)

Many studies suggest that there are distinct regulatory processes controlling compound leaf development in different clades of legumes. Loss of function of the LEAFY (LFY) orthologs results in a reduction of leaf complexity to different degrees in inverted repeat-lacking clade (IRLC) and non-IRLC species. To further understand the role of LFY orthologs and the molecular mechanism in compound leaf development in non-IRLC plants, we studied leaf development in unifoliate leaf (un) mutant, a classical mutant of mungbean (Vigna radiata L.), which showed a complete conversion of compound leaves into simple leaves. Our analysis revealed that UN encoded the mungbean LFY ortholog (VrLFY) and played a significant role in leaf development. In situ RNA hybridization results showed that STM-like KNOXI genes were expressed in compound leaf primordia in mungbean. Furthermore, increased leaflet number in heptafoliate leaflets1 (hel1) mutants was demonstrated to depend on the function of VrLFY and KNOXI genes in mungbean. Our results suggested that HEL1 is a key factor coordinating distinct processes in the control of compound leaf development in mungbean and its related non-IRLC legumes.

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.

The Analysis of Genes and Phytohormone Metabolic Pathways Associated with Leaf Shape Development in Liriodendron chinense via De Novo Transcriptome Sequencing

The leaf, a photosynthetic organ that plays an indispensable role in plant development and growth, has a certain ability to adapt to the environment and exhibits tremendous diversity among angiosperms. Liriodendron chinense, an ancestral angiosperm species, is very popular in landscaping. The leaf of this species has two lobes and resembles a Qing Dynasty Chinese robe; thus, leaf shape is the most valuable ornamental trait of the tree. In this work, to determine the candidate genes associated with leaf development in L. chinense, scanning electron microscopy (SEM) was employed to distinguish the developmental stages of tender leaves. Four stages were clearly separated, and transcriptome sequencing was performed for two special leaf stages. Altogether, there were 48.23 G clean reads in the libraries of the two leaf developmental stages, and 48,107 assembled unigenes were annotated with five databases. Among four libraries, 3118 differentially expressed genes (DEGs) were enriched in expression profiles. We selected ten DEGs associated with leaf development and validated their expression patterns via quantitative real-time PCR (qRT-PCR) assays. Most validation results were closely correlated with the RNA-sequencing data. Taken together, we examined the dynamic process of leaf development and indicated that several transcription factors and phytohormone metabolism genes may participate in leaf shape development. The transcriptome data analysis presented in this work aims to provide basic insights into the mechanisms mediating leaf development, and the results serve as a reference for the genetic breeding of ornamental traits in L. chinense.

Regulation of xylem fiber differentiation by gibberellins through DELLA-KNAT1 interaction

The thickening of plant organs is supported by secondary growth, a process by which new vascular tissues (xylem and phloem) are produced. Xylem is composed of several cell types, including xylary fibers, parenchyma and vessel elements. In Arabidopsis, it has been shown that fibers are promoted by the class-I KNOX gene KNAT1 and the plant hormones gibberellins, and are repressed by a small set of receptor-like kinases; however, we lack a mechanistic framework to integrate their relative contributions. Here, we show that DELLAs, negative elements of the gibberellin signaling pathway, physically interact with KNAT1 and impair its binding to KNAT1-binding sites. Our analysis also indicates that at least 37% of the transcriptome mobilized by KNAT1 is potentially dependent on this interaction, and includes genes involved in secondary cell wall modifications and phenylpropanoid biosynthesis. Moreover, the promotion by constitutive overexpression of KNAT1 of fiber formation and the expression of genes required for fiber differentiation were still reverted by DELLA accumulation, in agreement with post-translational regulation of KNAT1 by DELLA proteins. These results suggest that gibberellins enhance fiber development by promoting KNAT1 activity.

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.

Gene networks and the evolution of plant morphology

Elaboration of morphology depends on the precise orchestration of gene expression by key regulatory genes. The hierarchy and relationship among the participating genes is commonly known as gene regulatory network (GRN). Therefore, the evolution of morphology ultimately occurs by the rewiring of gene network structures or by the co-option of gene networks to novel domains. The availability of high-resolution expression data combined with powerful statistical tools have opened up new avenues to formulate and test hypotheses on how diverse gene networks influence trait development and diversity. Here we summarize recent studies based on both big-data and genetics approaches to understand the evolution of plant form and physiology. We also discuss recent genome-wide investigations on how studying open-chromatin regions may help study the evolution of gene expression patterns.

Getting leaves into shape: a molecular, cellular, environmental and evolutionary view

Leaves arise from groups of undifferentiated cells as small primordia that go through overlapping phases of morphogenesis, growth and differentiation. These phases are genetically controlled and modulated by environmental cues to generate a stereotyped, yet plastic, mature organ. Over the past couple of decades, studies have revealed that hormonal signals, transcription factors and miRNAs play major roles during leaf development, and more recent findings have highlighted the contribution of mechanical signals to leaf growth. In this Review, we discuss how modulating the activity of some of these regulators can generate diverse leaf shapes during development, in response to a varying environment, or between species during evolution.

Successive duplication-divergence mechanisms at the RCO locus contributed to leaf shape diversity in the Brassicaceae

Gene duplication is a major driver for the increase of biological complexity. The divergence of newly duplicated paralogs may allow novel functions to evolve, while maintaining the ancestral one. Alternatively, partitioning the ancestral function among paralogs may allow parts of that role to follow independent evolutionary trajectories. We studied the REDUCED COMPLEXITY (RCO) locus, which contains three paralogs that have evolved through two independent events of gene duplication, and which underlies repeated events of leaf shape evolution within the Brassicaceae. In particular, we took advantage of the presence of three potentially functional paralogs in Capsella to investigate the extent of functional divergence among them. We demonstrate that the RCO copies control growth in different areas of the leaf. Consequently, the copies that are retained active in the different Brassicaceae lineages contribute to define the leaf dissection pattern. Our results further illustrate how successive gene duplication events and subsequent functional divergence can increase trait evolvability by providing independent evolutionary trajectories to specialized functions that have an additive effect on a given trait.

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.

A common developmental program can produce diverse leaf shapes

Eudicot leaves have astoundingly diverse shapes. The central problem addressed in this paper is the developmental origin of this diversity. To investigate this problem, we propose a computational model of leaf development that generalizes the largely conserved molecular program for the reference plants Arabidopsis thaliana, Cardamine hirsuta and Solanum lycopersicum. The model characterizes leaf development as a product of three interwoven processes: the patterning of serrations, lobes and/or leaflets on the leaf margin; the patterning of the vascular system; and the growth of the leaf blade spanning the main veins. The veins play a significant morphogenetic role as a local determinant of growth directions. We show that small variations of this model can produce diverse leaf shapes, from simple to lobed to compound. It is thus plausible that diverse shapes of eudicot leaves result from small variations of a common developmental program.