The trichome pattern diversity of Cardamine shares genetic mechanisms with Arabidopsis but differs in environmental drivers.

Natural variation in trichome pattern (amount and distribution) is prominent among populations of many angiosperms. However, the degree of parallelism in the genetic mechanisms underlying this diversity and its environmental drivers in different species remain unclear. To address these questions, we analyzed the genomic and environmental bases of leaf trichome pattern diversity in Cardamine hirsuta, a relative of Arabidopsis (Arabidopsis thaliana). We characterized 123 wild accessions for their genomic diversity, leaf trichome patterns at different temperatures, and environmental adjustments. Nucleotide diversities and biogeographical distribution models identified two major genetic lineages with distinct demographic and adaptive histories. Additionally, C. hirsuta showed substantial variation in trichome pattern and plasticity to temperature. Trichome amount in C. hirsuta correlated positively with spring precipitation but negatively with temperature, which is opposite to climatic patterns in A. thaliana. Contrastingly, genetic analysis of C. hirsuta glabrous accessions indicated that, like for A. thaliana, glabrousness is caused by null mutations in ChGLABRA1 (ChGL1). Phenotypic genome-wide association studies (GWAS) further identified a ChGL1 haplogroup associated with low trichome density and ChGL1 expression. Therefore, a ChGL1 series of null and partial loss-of-function alleles accounts for the parallel evolution of leaf trichome pattern in C. hirsuta and A. thaliana. Finally, GWAS also detected other candidate genes (e.g. ChETC3, ChCLE17) that might affect trichome pattern. Accordingly, the evolution of this trait in C. hirsuta and A. thaliana shows partially conserved genetic mechanisms but is likely involved in adaptation to different environments.

Cell-cycle-linked growth reprogramming encodes developmental time into leaf morphogenesis.

How is time encoded into organ growth and morphogenesis? We address this question by investigating heteroblasty, where leaf development and form are modified with progressing plant age. By combining morphometric analyses, fate-mapping through live-imaging, computational analyses, and genetics, we identify age-dependent changes in cell-cycle-associated growth and histogenesis that underpin leaf heteroblasty. We show that in juvenile leaves, cell proliferation competence is rapidly released in a "proliferation burst" coupled with fast growth, whereas in adult leaves, proliferative growth is sustained for longer and at a slower rate. These effects are mediated by the SPL9 transcription factor in response to inputs from both shoot age and individual leaf maturation along the proximodistal axis. SPL9 acts by activating CyclinD3 family genes, which are sufficient to bypass the requirement for SPL9 in the control of leaf shape and in heteroblastic reprogramming of cellular growth. In conclusion, we have identified a mechanism that bridges across cell, tissue, and whole-organism scales by linking cell-cycle-associated growth control to age-dependent changes in organ geometry.

Genealogical Analyses of 3 Cultivated and 1 Wild Specimen of Vitis vinifera from Greece.

Grapevine (Vitis vinifera) has been an important crop with considerable cultural and economic significance for over 2,500 years, and Greece has been an important entry point into Europe for lineages that were domesticated in Western Asia and the Caucasus. However, whole-genome-based investigation of the demographic history of Greek cultivars relative to other European lineages has only started recently. To understand how Greek cultivars relate to Eurasian domesticated and wild populations, we sequenced 3 iconic domesticated strains ('Xinomavro,' 'Agiorgitiko,' 'Mavrotragano') along with 1 wild accession (the vinetree of Pausanias-a historically important wild specimen) and analyzed their genomic diversity together with a large sample of publicly available domesticated and wild strains. We also reconstructed genealogies by leveraging the powerful tsinfer methodology which has not previously been used in this system. We show that cultivated strains from Greece differ genetically from other strains in Europe. Interestingly, all the 3 cultivated Greek strains clustered with cultivated and wild accessions from Transcaucasia, South Asia, and the Levant and are amongst the very few cultivated European strains belonging to this cluster. Furthermore, our results indicate that 'Xinomavro' shares close genealogical proximity with European elite cultivars such as 'Chardonnay,' 'Riesling,' and 'Gamay' but not 'Pinot.' Therefore, the proximity of 'Xinomavro' to Gouais/Heunisch Weiss is confirmed and the utility of ancestral recombination graph reconstruction approaches to study genealogical relationships in crops is highlighted.

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.

From genes to shape in leaf development and evolution.

Plant leaves display tremendous variation in shape. Here, we discuss how information obtained from genetics, live imaging and computational modeling has helped conceptualize the ways in which gene activity is translated into different leaf shapes. In this framework, the action of genes on leaf form can be captured as the sum of their effects on the amount, duration, and direction of cellular growth, which together produce leaf geometry. We use three different examples to illustrate this point. First, the emergence of complex versus simple leaves in eudicots, which arises from differences in organ-wide growth duration as well as local growth repression at the leaf margin. Second, the development of strap-shaped grass leaves with a broad sheathing base versus the typical eudicot leaves with a narrow petiole, where these features of grass leaves emerged through lateral expansion of the zone of leaf progenitor cells, coupled with later remodeling of growth of early domains of the leaf blade. Third, the formation of insect traps on carnivorous plants that arose through constrained directional growth that produced a 3D deformation. In all the above examples, changes in gene expression of different classes of homeobox genes have contributed to the altered growth patterns underlying these different aspects of leaf shape diversity.

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.

Using positional information to provide context for biological image analysis with MorphoGraphX 2.0.

Positional information is a central concept in developmental biology. In developing organs, positional information can be idealized as a local coordinate system that arises from morphogen gradients controlled by organizers at key locations. This offers a plausible mechanism for the integration of the molecular networks operating in individual cells into the spatially coordinated multicellular responses necessary for the organization of emergent forms. Understanding how positional cues guide morphogenesis requires the quantification of gene expression and growth dynamics in the context of their underlying coordinate systems. Here, we present recent advances in the MorphoGraphX software (Barbier de Reuille et al., 2015⁠) that implement a generalized framework to annotate developing organs with local coordinate systems. These coordinate systems introduce an organ-centric spatial context to microscopy data, allowing gene expression and growth to be quantified and compared in the context of the positional information thought to control them.

Cytokinin promotes growth cessation in the Arabidopsis root.

The Arabidopsis root offers good opportunities to investigate how regulated cellular growth shapes different tissues and organs, a key question in developmental biology. Along the root's longitudinal axis, cells sequentially occupy different developmental states. Proliferative meristematic cells give rise to differentiating cells, which rapidly elongate in the elongation zone, then mature and stop growing in the differentiation zone. The phytohormone cytokinin contributes to this zonation by positioning the boundary between the meristem and the elongation zone, called the transition zone. However, the cellular growth profile underlying root zonation is not well understood, and the cellular mechanisms that mediate growth cessation remain unclear. By using time-lapse imaging, genetics, and computational analysis, we analyze the effect of cytokinin on root zonation and cellular growth. We found that cytokinin promotes growth cessation in the distal (shootward) elongation zone in conjunction with accelerating the transition from elongation to differentiation. We estimated cell-wall stiffness by using osmotic treatment experiments and found that cytokinin-mediated growth cessation is associated with cell-wall stiffening and requires the action of an auxin influx carrier, AUX1. Our measurement of growth and cell-wall mechanical properties at a cellular resolution reveal mechanisms via which cytokinin influences cell behavior to shape tissue patterns.

The annotation and analysis of complex 3D plant organs using 3DCoordX.

A fundamental question in biology concerns how molecular and cellular processes become integrated during morphogenesis. In plants, characterization of 3D digital representations of organs at single-cell resolution represents a promising approach to addressing this problem. A major challenge is to provide organ-centric spatial context to cells of an organ. We developed several general rules for the annotation of cell position and embodied them in 3DCoordX, a user-interactive computer toolbox implemented in the open-source software MorphoGraphX. 3DCoordX enables rapid spatial annotation of cells even in highly curved biological shapes. Using 3DCoordX, we analyzed cellular growth patterns in organs of several species. For example, the data indicated the presence of a basal cell proliferation zone in the ovule primordium of Arabidopsis (Arabidopsis thaliana). Proof-of-concept analyses suggested a preferential increase in cell length associated with neck elongation in the archegonium of Marchantia (Marchantia polymorpha) and variations in cell volume linked to central morphogenetic features of a trap of the carnivorous plant Utricularia (Utricularia gibba). Our work demonstrates the broad applicability of the developed strategies as they provide organ-centric spatial context to cellular features in plant organs of diverse shape complexity.

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.

Gene expression evolution in pattern-triggered immunity within Arabidopsis thaliana and across Brassicaceae species.

Plants recognize surrounding microbes by sensing microbe-associated molecular patterns (MAMPs) to activate pattern-triggered immunity (PTI). Despite their significance for microbial control, the evolution of PTI responses remains largely uncharacterized. Here, by employing comparative transcriptomics of six Arabidopsis thaliana accessions and three additional Brassicaceae species to investigate PTI responses, we identified a set of genes that commonly respond to the MAMP flg22 and genes that exhibit species-specific expression signatures. Variation in flg22-triggered transcriptome responses across Brassicaceae species was incongruent with their phylogeny, while expression changes were strongly conserved within A. thaliana. We found the enrichment of WRKY transcription factor binding sites in the 5'-regulatory regions of conserved and species-specific responsive genes, linking the emergence of WRKY-binding sites with the evolution of gene expression patterns during PTI. Our findings advance our understanding of the evolution of the transcriptome during biotic stress.

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.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species.

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

A WOX/Auxin Biosynthesis Module Controls Growth to Shape Leaf Form.

A key challenge in biology is to understand how the regional control of cell growth gives rise to final organ forms. Plant leaves must coordinate growth along both the proximodistal and mediolateral axes to produce their final shape. However, the cell-level mechanisms controlling this coordination remain largely unclear. Here, we show that, in A. thaliana, WOX5, one of the WUSCHEL-RELATED HOMEOBOX (WOX) family of homeobox genes, acts redundantly with WOX1 and WOX3 (PRESSED FLOWER [PRS]) to control leaf shape. Through genetics and hormone measurements, we find that these WOXs act in part through the regional control of YUCCA (YUC) auxin biosynthetic gene expression along the leaf margin. The requirement for WOX-mediated YUC expression in patterning of leaf shape cannot be bypassed by the epidermal expression of YUC, indicating that the precise domain of auxin biosynthesis is important for leaf form. Using time-lapse growth analysis, we demonstrate that WOX-mediated auxin biosynthesis organizes a proximodistal growth gradient that promotes lateral growth and consequently the characteristic ellipsoid A. thaliana leaf shape. We also provide evidence that WOX proteins shape the proximodistal gradient of differentiation by inhibiting differentiation proximally in the leaf blade and promoting it distally. This regulation allows sustained growth of the blade and enables a leaf to attain its final form. In conclusion, we show that the WOX/auxin regulatory module shapes leaf form by coordinating growth along the proximodistal and mediolateral leaf axes.

Accurate and versatile 3D segmentation of plant tissues at cellular resolution.

Quantitative analysis of plant and animal morphogenesis requires accurate segmentation of individual cells in volumetric images of growing organs. In the last years, deep learning has provided robust automated algorithms that approach human performance, with applications to bio-image analysis now starting to emerge. Here, we present PlantSeg, a pipeline for volumetric segmentation of plant tissues into cells. PlantSeg employs a convolutional neural network to predict cell boundaries and graph partitioning to segment cells based on the neural network predictions. PlantSeg was trained on fixed and live plant organs imaged with confocal and light sheet microscopes. PlantSeg delivers accurate results and generalizes well across different tissues, scales, acquisition settings even on non plant samples. We present results of PlantSeg applications in diverse developmental contexts. PlantSeg is free and open-source, with both a command line and a user-friendly graphical interface.

Oil Body Formation in Marchantia polymorpha Is Controlled by MpC1HDZ and Serves as a Defense against Arthropod Herbivores.

The origin of a terrestrial flora in the Ordovician required adaptation to novel biotic and abiotic stressors. Oil bodies, a synapomorphy of liverworts, accumulate secondary metabolites, but their function and development are poorly understood. Oil bodies of Marchantia polymorpha develop within specialized cells as one single large organelle. Here, we show that a class I homeodomain leucine-zipper (C1HDZ) transcription factor controls the differentiation of oil body cells in two different ecotypes of the liverwort M. polymorpha, a model genetic system for early divergent land plants. In flowering plants, these transcription factors primarily modulate responses to abiotic stress, including drought. However, loss-of-function alleles of the single ortholog gene, MpC1HDZ, in M. polymorpha did not exhibit phenotypes associated with abiotic stress. Rather, Mpc1hdz mutant plants were more susceptible to herbivory, and total plant extracts of the mutant exhibited reduced antibacterial activity. Transcriptomic analysis of the mutant revealed a reduction in expression of genes related to secondary metabolism that was accompanied by a specific depletion of oil body terpenoid compounds. Through time-lapse imaging, we observed that MpC1HDZ expression maxima precede oil body formation, indicating that MpC1HDZ mediates differentiation of oil body cells. Our results indicate that M. polymorpha oil bodies, and MpC1HDZ, are critical for defense against herbivory, but not for abiotic stress tolerance. Thus, C1HDZ genes were co-opted to regulate separate responses to biotic and abiotic stressors in two distinct land plant lineages.

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.

Phyllotaxis: is the golden angle optimal for light capture?

Phyllotactic patterns are some of the most conspicuous in nature. To create these patterns plants must control the divergence angle between the appearance of successive organs, sometimes to within a fraction of a degree. The most common angle is the Fibonacci or golden angle, and its prevalence has led to the hypothesis that it has been selected by evolution as optimal for plants with respect to some fitness benefits, such as light capture. We explore arguments for and against this idea with computer models. We have used both idealized and scanned leaves from Arabidopsis thaliana and Cardamine hirsuta to measure the overlapping leaf area of simulated plants after varying parameters such as leaf shape, incident light angles, and other leaf traits. We find that other angles generated by Fibonacci-like series found in nature are equally optimal for light capture, and therefore should be under similar evolutionary pressure. Our findings suggest that the iterative mechanism for organ positioning itself is a more likely target for evolutionary pressure, rather than a specific divergence angle, and our model demonstrates that the heteroblastic progression of leaf shape in A. thaliana can provide a potential fitness benefit via light capture.

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.