LEAFY and WAPO1 jointly regulate spikelet number per spike and floret development in wheat.

In wheat, the transition of the inflorescence meristem to a terminal spikelet (IM→TS) determines the spikelet number per spike (SNS), an important yield component. In this study, we demonstrate that the plant-specific transcription factor LEAFY (LFY) physically and genetically interacts with WHEAT ORTHOLOG OF APO1 (WAPO1) to regulate SNS and floret development. Loss-of-function mutations in either or both genes result in significant and similar reductions in SNS, as a result of a reduction in the rate of spikelet meristem formation per day. SNS is also modulated by significant genetic interactions between LFY and the SQUAMOSA MADS-box genes VRN1 and FUL2, which promote the IM→TS transition. Single-molecule fluorescence in situ hybridization revealed a downregulation of LFY and upregulation of the SQUAMOSA MADS-box genes in the distal part of the developing spike during the IM→TS transition, supporting their opposite roles in the regulation of SNS in wheat. Concurrently, the overlap of LFY and WAPO1 transcription domains in the developing spikelets contributes to normal floret development. Understanding the genetic network regulating SNS is a necessary first step to engineer this important agronomic trait.

Sequencing 4.3 million mutations in wheat promoters to understand and modify gene expression.

Wheat is an important contributor to global food security, and further improvements are required to feed a growing human population. Functional genetics and genomics tools can help us to understand the function of different genes and to engineer beneficial changes. In this study, we used a promoter capture assay to sequence 2-kb regions upstream of all high-confidence annotated genes from 1,513 mutagenized plants from the tetraploid wheat variety Kronos. We identified 4.3 million induced mutations with an accuracy of 99.8%, resulting in a mutation density of 41.9 mutations per kb. We also remapped Kronos exome capture reads to Chinese Spring RefSeq v1.1, identified 4.7 million mutations, and predicted their effects on annotated genes. Using these predictions, we identified 59% more nonsynonymous substitutions and 49% more truncation mutations than in the original study. To show the biological value of the promoter dataset, we selected two mutations within the promoter of the VRN-A1 vernalization gene. Both mutations, located within transcription factor binding sites, significantly altered VRN-A1 expression, and one reduced the number of spikelets per spike. These publicly available sequenced mutant datasets provide rapid and inexpensive access to induced variation in the promoters and coding regions of most wheat genes. These mutations can be used to understand and modulate gene expression and phenotypes for both basic and commercial applications, where limited governmental regulations can facilitate deployment. These mutant collections, together with gene editing, provide valuable tools to accelerate functional genetic studies in this economically important crop.

GIGANTEA accelerates wheat heading time through gene interactions converging on FLOWERING LOCUS T1.

Precise regulation of flowering time is critical for cereal crops to synchronize reproductive development with optimum environmental conditions, thereby maximizing grain yield. The plant-specific gene GIGANTEA (GI) plays an important role in the control of flowering time, with additional functions on the circadian clock and plant stress responses. In this study, we show that GI loss-of-function mutants in a photoperiod-sensitive tetraploid wheat background exhibit significant delays in heading time under both long-day (LD) and short-day photoperiods, with stronger effects under LD. However, this interaction between GI and photoperiod is no longer observed in isogenic lines carrying either a photoperiod-insensitive allele in the PHOTOPERIOD1 (PPD1) gene or a loss-of-function allele in EARLY FLOWERING 3 (ELF3), a known repressor of PPD1. These results suggest that the normal circadian regulation of PPD1 is required for the differential effect of GI on heading time in different photoperiods. Using crosses between mutant or transgenic plants of GI and those of critical genes in the flowering regulation pathway, we show that GI accelerates wheat heading time by promoting FLOWERING LOCUS T1 (FT1) expression via interactions with ELF3, VERNALIZATION 2 (VRN2), CONSTANS (CO), and the age-dependent microRNA172-APETALA2 (AP2) pathway, at both transcriptional and protein levels. Our study reveals conserved GI mechanisms between wheat and Arabidopsis but also identifies specific interactions of GI with the distinctive photoperiod and vernalization pathways of the temperate grasses. These results provide valuable knowledge for modulating wheat heading time and engineering new varieties better adapted to a changing environment.

Use of GRF-GIF chimeras and a ternary vector system to improve maize (Zea mays L.) transformation frequency.

Maize (Zea mays L.) is an important crop that has been widely studied for its agronomic and industrial applications and is one of the main classical model organisms for genetic research. Agrobacterium-mediated transformation of immature maize embryos is a commonly used method to introduce transgenes, but a low transformation frequency remains a bottleneck for many gene-editing applications. Previous approaches to enhance transformation included the improvement of tissue culture media and the use of morphogenic regulators such as BABY BOOM and WUSCHEL2. Here, we show that the frequency can be increased using a pVS1-VIR2 virulence helper plasmid to improve T-DNA delivery, and/or expressing a fusion protein between a GROWTH-REGULATING FACTOR (GRF) and GRF-INTERACTING FACTOR (GIF) protein to improve regeneration. Using hygromycin as a selection agent to avoid escapes, the transformation frequency in the maize inbred line B104 significantly improved from 2.3 to 8.1% when using the pVS1-VIR2 helper vector with no effect on event quality regarding T-DNA copy number. Combined with a novel fusion protein between ZmGRF1 and ZmGIF1, transformation frequencies further improved another 3.5- to 6.5-fold with no obvious impact on plant growth, while simultaneously allowing efficient CRISPR-/Cas9-mediated gene editing. Our results demonstrate how a GRF-GIF chimera in conjunction with a ternary vector system has the potential to further improve the efficiency of gene-editing applications and molecular biology studies in maize.

MiR172-APETALA2-like genes integrate vernalization and plant age to control flowering time in wheat.

Plants possess regulatory mechanisms that allow them to flower under conditions that maximize reproductive success. Selection of natural variants affecting those mechanisms has been critical in agriculture to modulate the flowering response of crops to specific environments and to increase yield. In the temperate cereals, wheat and barley, the photoperiod and vernalization pathways explain most of the natural variation in flowering time. However, other pathways also participate in fine-tuning the flowering response. In this work, we integrate the conserved microRNA miR172 and its targets APETALA2-like (AP2L) genes into the temperate grass flowering network involving VERNALIZATION 1 (VRN1), VRN2 and FLOWERING LOCUS T 1 (FT1 = VRN3) genes. Using mutants, transgenics and different growing conditions, we show that miR172 promotes flowering in wheat, while its target genes AP2L1 (TaTOE1) and AP2L5 (Q) act as flowering repressors. Moreover, we reveal that the miR172-AP2L pathway regulates FT1 expression in the leaves, and that this regulation is independent of VRN2 and VRN1. In addition, we show that the miR172-AP2L module and flowering are both controlled by plant age through miR156 in spring cultivars. However, in winter cultivars, flowering and the regulation of AP2L1 expression are decoupled from miR156 downregulation with age, and induction of VRN1 by vernalization is required to repress AP2L1 in the leaves and promote flowering. Interestingly, the levels of miR172 and both AP2L genes modulate the flowering response to different vernalization treatments in winter cultivars. In summary, our results show that conserved and grass specific gene networks interact to modulate the flowering response, and that natural or induced mutations in AP2L genes are useful tools for fine-tuning wheat flowering time in a changing environment.

Interactions between SQUAMOSA and SHORT VEGETATIVE PHASE MADS-box proteins regulate meristem transitions during wheat spike development.

Inflorescence architecture is an important determinant of crop productivity. The number of spikelets produced by the wheat inflorescence meristem (IM) before its transition to a terminal spikelet (TS) influences the maximum number of grains per spike. Wheat MADS-box genes VERNALIZATION 1 (VRN1) and FRUITFULL 2 (FUL2) (in the SQUAMOSA-clade) are essential to promote the transition from IM to TS and for spikelet development. Here we show that SQUAMOSA genes contribute to spikelet identity by repressing MADS-box genes VEGETATIVE TO REPRODUCTIVE TRANSITION 2 (VRT2), SHORT VEGETATIVE PHASE 1 (SVP1), and SVP3 in the SVP clade. Constitutive expression of VRT2 resulted in leafy glumes and lemmas, reversion of spikelets to spikes, and downregulation of MADS-box genes involved in floret development, whereas the vrt2 mutant reduced vegetative characteristics in spikelets of squamosa mutants. Interestingly, the vrt2 svp1 mutant showed similar phenotypes to squamosa mutants regarding heading time, plant height, and spikelets per spike, but it exhibited unusual axillary inflorescences in the elongating stem. We propose that SQUAMOSA-SVP interactions are important to promote heading, formation of the TS, and stem elongation during the early reproductive phase, and that downregulation of SVP genes is then necessary for normal spikelet and floral development. Manipulating SVP and SQUAMOSA genes can contribute to engineering spike architectures with improved productivity.

Optimization of ATAC-seq in wheat seedling roots using INTACT-isolated nuclei.

BACKGROUND: The genetic information contained in the genome of an organism is organized in genes and regulatory elements that control gene expression. The genomes of multiple plants species have already been sequenced and the gene repertory have been annotated, however, cis-regulatory elements remain less characterized, limiting our understanding of genome functionality. These elements act as open platforms for recruiting both positive- and negative-acting transcription factors, and as such, chromatin accessibility is an important signature for their identification. RESULTS: In this work we developed a transgenic INTACT [isolation of nuclei tagged in specific cell types] system in tetraploid wheat for nuclei purifications. Then, we combined the INTACT system together with the assay for transposase-accessible chromatin with sequencing [ATAC-seq] to identify open chromatin regions in wheat root tip samples. Our ATAC-seq results showed a large enrichment of open chromatin regions in intergenic and promoter regions, which is expected for regulatory elements and that is similar to ATAC-seq results obtained in other plant species. In addition, root ATAC-seq peaks showed a significant overlap with a previously published ATAC-seq data from wheat leaf protoplast, indicating a high reproducibility between the two experiments and a large overlap between open chromatin regions in root and leaf tissues. Importantly, we observed overlap between ATAC-seq peaks and cis-regulatory elements that have been functionally validated in wheat, and a good correlation between normalized accessibility and gene expression levels. CONCLUSIONS: We have developed and validated an INTACT system in tetraploid wheat that allows rapid and high-quality nuclei purification from root tips. Those nuclei were successfully used to performed ATAC-seq experiments that revealed open chromatin regions in the wheat genome that will be useful to identify cis-regulatory elements. The INTACT system presented here will facilitate the development of ATAC-seq datasets in other tissues, growth stages, and under different growing conditions to generate a more complete landscape of the accessible DNA regions in the wheat genome.

Evolutionary Footprints Reveal Insights into Plant MicroRNA Biogenesis.

MicroRNAs (miRNAs) are endogenous small RNAs that recognize target sequences by base complementarity and play a role in the regulation of target gene expression. They are processed from longer precursor molecules that harbor a fold-back structure. Plant miRNA precursors are quite variable in size and shape, and are recognized by the processing machinery in different ways. However, ancient miRNAs and their binding sites in target genes are conserved during evolution. Here, we designed a strategy to systematically analyze MIRNAs from different species generating a graphical representation of the conservation of the primary sequence and secondary structure. We found that plant MIRNAs have evolutionary footprints that go beyond the small RNA sequence itself, yet their location along the precursor depends on the specific MIRNA We show that these conserved regions correspond to structural determinants recognized during the biogenesis of plant miRNAs. Furthermore, we found that the members of the miR166 family have unusual conservation patterns and demonstrated that the recognition of these precursors in vivo differs from other known miRNAs. Our results describe a link between the evolutionary conservation of plant MIRNAs and the mechanisms underlying the biogenesis of these small RNAs and show that the MIRNA pattern of conservation can be used to infer the mode of miRNA biogenesis.

Identification of the VERNALIZATION 4 gene reveals the origin of spring growth habit in ancient wheats from South Asia.

Wheat varieties with a winter growth habit require long exposures to low temperatures (vernalization) to accelerate flowering. Natural variation in four vernalization genes regulating this requirement has favored wheat adaptation to different environments. The first three genes (VRN1-VRN3) have been cloned and characterized before. Here we show that the fourth gene, VRN-D4, originated by the insertion of a ∼290-kb region from chromosome arm 5AL into the proximal region of chromosome arm 5DS. The inserted 5AL region includes a copy of VRN-A1 that carries distinctive mutations in its coding and regulatory regions. Three lines of evidence confirmed that this gene is VRN-D4: it cosegregated with VRN-D4 in a high-density mapping population; it was expressed earlier than other VRN1 genes in the absence of vernalization; and induced mutations in this gene resulted in delayed flowering. VRN-D4 was found in most accessions of the ancient subspecies Triticum aestivum ssp. sphaerococcum from South Asia. This subspecies showed a significant reduction of genetic diversity and increased genetic differentiation in the centromeric region of chromosome 5D, suggesting that VRN-D4 likely contributed to local adaptation and was favored by positive selection. Three adjacent SNPs in a regulatory region of the VRN-D4 first intron disrupt the binding of GLYCINE-RICH RNA-BINDING PROTEIN 2 (TaGRP2), a known repressor of VRN1 expression. The same SNPs were identified in VRN-A1 alleles previously associated with reduced vernalization requirement. These alleles can be used to modulate vernalization requirements and to develop wheat varieties better adapted to different or changing environments.

Repression of growth regulating factors by the microRNA396 inhibits cell proliferation by UV-B radiation in Arabidopsis leaves.

Because of their sessile lifestyle, plants are continuously exposed to solar UV-B radiation. Inhibition of leaf growth is one of the most consistent responses of plants upon exposure to UV-B radiation. In this work, we investigated the role of Growth-Regulating Factors (GRFs) and of microRNA miR396 in UV-B-mediated inhibition of leaf growth in Arabidopsis thaliana plants. We demonstrate that miRNA396 is upregulated by UV-B radiation in proliferating tissues and that this induction is correlated with a decrease in GRF1, GRF2, and GRF3 transcripts. Induction of miR396 results in inhibition of cell proliferation, and this outcome is independent of the UV-B photoreceptor UV resistance locus 8, as well as ATM AND RAD3-related and the mitogen-activated protein kinase MPK6, but is dependent on MPK3. Transgenic plants expressing an artificial target mimic directed against miR396 (MIM396) with a decrease in the endogenous microRNA activity or plants expressing miR396-resistant copies of several GRFs are less sensitive to this inhibition. Consequently, at intensities that can induce DNA damage in Arabidopsis plants, UV-B radiation limits leaf growth by inhibiting cell division in proliferating tissues, a process mediated by miR396 and GRFs.

Post-transcriptional control of GRF transcription factors by microRNA miR396 and GIF co-activator affects leaf size and longevity.

The growth-regulating factors (GRFs) are plant-specific transcription factors. They form complexes with GRF-interacting factors (GIFs), a small family of transcriptional co-activators. In Arabidopsis thaliana, seven out of the nine GRFs are controlled by microRNA miR396. Analysis of Arabidopsis plants carrying a GRF3 allele insensitive to miR396 revealed a strong boost in the number of cells in leaves, which was further enhanced synergistically by an additional increase of GIF1 levels. Genetic experiments revealed that GRF3 can still increase cell number in gif1 mutants, albeit to a much lesser extent. Genome-wide transcript profiling indicated that the simultaneous increase of GRF3 and GIF1 levels causes additional effects in gene expression compared to either of the transgenes alone. We observed that GIF1 interacts in vivo with GRF3, as well as with chromatin-remodeling complexes, providing a mechanistic explanation for the synergistic activities of a GRF3-GIF1 complex. Interestingly, we found that, in addition to the leaf size, the GRF system also affects the organ longevity. Genetic and molecular analysis revealed that the functions of GRFs in leaf growth and senescence can be uncoupled, demonstrating that the miR396-GRF-GIF network impinges on different stages of leaf development. Our results integrate the post-transcriptional control of the GRF transcription factors with the progression of leaf development.

Functional specialization of the plant miR396 regulatory network through distinct microRNA-target interactions.

MicroRNAs (miRNAs) are ∼21 nt small RNAs that regulate gene expression in animals and plants. They can be grouped into families comprising different genes encoding similar or identical mature miRNAs. Several miRNA families are deeply conserved in plant lineages and regulate key aspects of plant development, hormone signaling, and stress response. The ancient miRNA miR396 regulates conserved targets belonging to the GROWTH-REGULATING FACTOR (GRF) family of transcription factors, which are known to control cell proliferation in Arabidopsis leaves. In this work, we characterized the regulation of an additional target for miR396, the transcription factor bHLH74, that is necessary for Arabidopsis normal development. bHLH74 homologs with a miR396 target site could only be detected in the sister families Brassicaceae and Cleomaceae. Still, bHLH74 repression by miR396 is required for margin and vein pattern formation of Arabidopsis leaves. MiR396 contributes to the spatio-temporal regulation of GRF and bHLH74 expression during leaf development. Furthermore, a survey of miR396 sequences in different species showed variations in the 5' portion of the miRNA, a region known to be important for miRNA activity. Analysis of different miR396 variants in Arabidopsis thaliana revealed that they have an enhanced activity toward GRF transcription factors. The interaction between the GRF target site and miR396 has a bulge between positions 7 and 8 of the miRNA. Our data indicate that such bulge modulates the strength of the miR396-mediated repression and that this modulation is essential to shape the precise spatio-temporal pattern of GRF2 expression. The results show that ancient miRNAs can regulate conserved targets with varied efficiency in different species, and we further propose that they could acquire new targets whose control might also be biologically relevant.

An efficient protoplast-based genome editing protocol for Vitis species.

CRISPR-Cas technologies allow for precise modifications in plant genomes and promise to revolutionize agriculture. These technologies depend on the delivery of editing components into plant cells and the regeneration of fully edited plants. In vegetatively propagated plants, such as grape, protoplast culture provides one of the best avenues for producing non-chimeric and transgene-free genome-edited plants. However, poor regeneration of plants from protoplasts has hindered their implementation for genome editing. Here, we report an efficient protocol for regenerating plants from protoplasts from multiple grape varieties. By encapsulating the protoplasts in calcium alginate beads and co-culturing them with feeder cultures, the protoplasts divide to form callus colonies that regenerate into embryos and ultimately plants. This protocol worked successfully in wine and table grape (Vitis vinifera) varieties, as well as grape rootstocks and the grapevine wild relative Vitis arizonica. Moreover, by transfecting protoplasts with CRISPR-plasmid or ribonucleoprotein (RNP) complexes, we regenerated albino plants with edits in VvPHYTOENE DESATURASE gene in three varieties and in V. arizonica. The results reveal the potential of this platform to facilitate genome editing in Vitis species.

Control of cell proliferation in Arabidopsis thaliana by microRNA miR396.

Cell proliferation is an important determinant of plant form, but little is known about how developmental programs control cell division. Here, we describe the role of microRNA miR396 in the coordination of cell proliferation in Arabidopsis leaves. In leaf primordia, miR396 is expressed at low levels that steadily increase during organ development. We found that miR396 antagonizes the expression pattern of its targets, the GROWTH-REGULATING FACTOR (GRF) transcription factors. miR396 accumulates preferentially in the distal part of young developing leaves, restricting the expression of GRF2 to the proximal part of the organ. This, in turn, coincides with the activity of the cell proliferation marker CYCLINB1;1. We show that miR396 attenuates cell proliferation in developing leaves, through the repression of GRF activity and a decrease in the expression of cell cycle genes. We observed that the balance between miR396 and the GRFs controls the final number of cells in leaves. Furthermore, overexpression of miR396 in a mutant lacking GRF-INTERACTING FACTOR 1 severely compromises the shoot meristem. We found that miR396 is expressed at low levels throughout the meristem, overlapping with the expression of its target, GRF2. In addition, we show that miR396 can regulate cell proliferation and the size of the meristem. Arabidopsis plants with an increased activity of the transcription factor TCP4, which reduces cell proliferation in leaves, have higher miR396 and lower GRF levels. These results implicate miR396 as a significant module in the regulation of cell proliferation in plants.

A second-generation capture panel for cost-effective sequencing of genome regulatory regions in wheat and relatives.

As genome resources for wheat (Triticum L.) expand at a rapid pace, it is important to update targeted sequencing tools to incorporate improved sequence assemblies and regions of previously unknown significance. Here, we developed an updated regulatory region enrichment capture for wheat and other Triticeae species. The core target space includes sequences from 2-Kbp upstream of each gene predicted in the Chinese Spring wheat genome (IWGSC RefSeq Annotation v1.0) and regions of open chromatin identified with an assay for transposase-accessible chromatin using sequencing from wheat leaf and root samples. To improve specificity, we aggressively filtered candidate repetitive sequences using a combination of nucleotide basic local alignment search tool (BLASTN) searches to the Triticeae Repetitive Sequence Database (TREP), identification of regions with read over-coverage from previous target enrichment experiments, and k-mer frequency analyses. The final design comprises 216.5 Mbp of predicted hybridization space in hexaploid wheat and showed increased specificity and coverage of targeted sequences relative to previous protocols. Test captures on hexaploid and tetraploid wheat and other diploid cereals show that the assay has broad potential utility for cost-effective promoter and open chromatin resequencing and general-purpose genotyping of various Triticeae species.

An efficient transformation method for genome editing of elite bread wheat cultivars.

An efficient genetic transformation protocol is necessary to edit genes for trait improvement directly in elite bread wheat cultivars. We used a protein fusion between a wheat growth-regulating factor 4 (GRF4) and its interacting factor (GIF1) to develop a reproducible genetic transformation and regeneration protocol, which we then used to successfully transform elite bread wheat cultivars Baj, Kachu, Morocco, Reedling, RL6077, and Sujata in addition to the experimental cultivar Fielder. Immature embryos were transformed with the vector using particle bombardment method. Transformation frequency increased nearly 60-fold with the GRF4-GIF1-containing vectors as compared to the control vector and ranged from ~5% in the cultivar Kachu to 13% in the cultivar RL6077. We then edited two genes that confer resistance against leaf rust and powdery mildew directly in the aforementioned elite cultivars. A wheat promoter, TaU3 or TaU6, to drive the expression of guide RNA was effective in gene editing whereas the OsU3 promoter failed to generate any edits. Editing efficiency was nearly perfect with the wheat promoters. Our protocol has made it possible to edit genes directly in elite wheat cultivars and would be useful for gene editing in other wheat varieties, which have been recalcitrant to transformation thus far.

Recent advances in crop transformation technologies.

Agriculture is experiencing a technological inflection point in its history, while also facing unprecedented challenges posed by human population growth and global climate changes. Key advancements in precise genome editing and new methods for rapid generation of bioengineered crops promise to both revolutionize the speed and breadth of breeding programmes and increase our ability to feed and sustain human population growth. Although genome editing enables targeted and specific modifications of DNA sequences, several existing barriers prevent the widespread adoption of editing technologies for basic and applied research in established and emerging crop species. Inefficient methods for the transformation and regeneration of recalcitrant species and the genotype dependency of the transformation process remain major hurdles. These limitations are frequent in monocotyledonous crops, which alone provide most of the calories consumed by human populations. Somatic embryogenesis and de novo induction of meristems - pluripotent groups of stem cells responsible for plant developmental plasticity - are essential strategies to quickly generate transformed plants. Here we review recent discoveries that are rapidly advancing nuclear transformation technologies and promise to overcome the obstacles that have so far impeded the widespread adoption of genome editing in crop species.

A GRF-GIF chimeric protein improves the regeneration efficiency of transgenic plants.

The potential of genome editing to improve the agronomic performance of crops is often limited by low plant regeneration efficiencies and few transformable genotypes. Here, we show that expression of a fusion protein combining wheat GROWTH-REGULATING FACTOR 4 (GRF4) and its cofactor GRF-INTERACTING FACTOR 1 (GIF1) substantially increases the efficiency and speed of regeneration in wheat, triticale and rice and increases the number of transformable wheat genotypes. GRF4-GIF1 transgenic plants were fertile and without obvious developmental defects. Moreover, GRF4-GIF1 induced efficient wheat regeneration in the absence of exogenous cytokinins, which facilitates selection of transgenic plants without selectable markers. We also combined GRF4-GIF1 with CRISPR-Cas9 genome editing and generated 30 edited wheat plants with disruptions in the gene Q (AP2L-A5). Finally, we show that a dicot GRF-GIF chimera improves regeneration efficiency in citrus, suggesting that this strategy can be applied to dicot crops.

Morphogenesis of simple leaves: regulation of leaf size and shape.

Plants produce new organs throughout their life span. Leaves first initiate as rod-like structures protruding from the shoot apical meristem, while they need to pass through different developmental stages to become the flat organ specialized in photosynthesis. Leaf morphogenesis is an active process regulated by many genes and pathways that can generate organs with a wide variety of sizes and shapes. Important differences in leaf architecture can be seen among different species, but also in single individuals. A key aspect of leaf morphogenesis is the precise control of cell proliferation. Modification or manipulation of this process may lead to leaves with different sizes and shapes, and changes in the organ margins and curvature. Many genes required for leaf development have been identified in Arabidopsis thaliana, and the mechanisms underlying leaf morphogenesis are starting to be unraveled at the molecular level.