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.

The GRAS transcription factor CsTL regulates tendril formation in cucumber.

Cucumber (Cucumis sativus, Cs) tendrils are slender vegetative organs that typically require manual removal to ensure orderly growth during greenhouse cultivation. Here, we identified cucumber tendril-less (tl), a Tnt1 retrotransposon-induced insertion mutant lacking tendrils. Map-based cloning identified the mutated gene, CsaV3_3G003590, which we designated as CsTL, which is homologous to Arabidopsis thaliana LATERAL SUPPRESSOR (AtLAS). Knocking out CsTL repressed tendril formation but did not affect branch initiation, whereas overexpression (OE) of CsTL resulted in the formation of two or more tendrils in one leaf axil. Although expression of two cucumber genes regulating tendril formation, Tendril (CsTEN) and Unusual Floral Organs (CsUFO), was significantly decreased in CsTL knockout lines, these two genes were not direct downstream targets of CsTL. Instead, CsTL physically interacted with CsTEN, an interaction that further enhanced CsTEN-mediated expression of CsUFO. In Arabidopsis, the CsTL homolog AtLAS acts upstream of REVOLUTA (REV) to regulate branch initiation. Knocking out cucumber CsREV inhibited branch formation without affecting tendril initiation. Furthermore, genomic regions containing CsTL and AtLAS were not syntenic between the cucumber and Arabidopsis genomes, whereas REV orthologs were found on a shared syntenic block. Our results revealed not only that cucumber CsTL possesses a divergent function in promoting tendril formation but also that CsREV retains its conserved function in shoot branching.

Genetic regulation of self-organizing azimuthal canopy orientations and their impacts on light interception in maize.

The efficiency of solar radiation interception contributes to the photosynthetic efficiency of crop plants. Light interception is a function of canopy architecture, including plant density; leaf number, length, width, and angle; and azimuthal canopy orientation. We report on the ability of some maize (Zea mays) genotypes to alter the orientations of their leaves during development in coordination with adjacent plants. Although the upper canopies of these genotypes retain the typical alternate-distichous phyllotaxy of maize, their leaves grow parallel to those of adjacent plants. A genome-wide association study (GWAS) on this parallel canopy trait identified candidate genes, many of which are associated with shade avoidance syndrome, including phytochromeC2. GWAS conducted on the fraction of photosynthetically active radiation (PAR) intercepted by canopies also identified multiple candidate genes, including liguleless1 (lg1), previously defined by its role in ligule development. Under high plant densities, mutants of shade avoidance syndrome and liguleless genes (lg1, lg2, and Lg3) exhibit altered canopy patterns, viz, the numbers of interrow leaves are greatly reduced as compared to those of nonmutant controls, resulting in dramatically decreased PAR interception. In at least the case of lg2, this phenotype is not a consequence of abnormal ligule development. Instead, liguleless gene functions are required for normal light responses, including azimuth canopy re-orientation.

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.

Nudix hydrolase 23 post-translationally regulates carotenoid biosynthesis in plants.

Carotenoids are essential for photosynthesis and photoprotection. Plants must evolve multifaceted regulatory mechanisms to control carotenoid biosynthesis. However, the regulatory mechanisms and the regulators conserved among plant species remain elusive. Phytoene synthase (PSY) catalyzes the highly regulated step of carotenogenesis and geranylgeranyl diphosphate synthase (GGPPS) acts as a hub to interact with GGPP-utilizing enzymes for the synthesis of specific downstream isoprenoids. Here, we report a function of Nudix hydrolase 23 (NUDX23), a Nudix domain-containing protein, in post-translational regulation of PSY and GGPPS for carotenoid biosynthesis. NUDX23 expresses highly in Arabidopsis (Arabidopsis thaliana) leaves. Overexpression of NUDX23 significantly increases PSY and GGPPS protein levels and carotenoid production, whereas knockout of NUDX23 dramatically reduces their abundances and carotenoid accumulation in Arabidopsis. NUDX23 regulates carotenoid biosynthesis via direct interactions with PSY and GGPPS in chloroplasts, which enhances PSY and GGPPS protein stability in a large PSY-GGPPS enzyme complex. NUDX23 was found to co-migrate with PSY and GGPPS proteins and to be required for the enzyme complex assembly. Our findings uncover a regulatory mechanism underlying carotenoid biosynthesis in plants and offer promising genetic tools for developing carotenoid-enriched food crops.

Transcription factor FveMYB117a inhibits axillary bud outgrowth by regulating cytokinin homeostasis in woodland strawberry.

Shoot branching affects plant architecture. In strawberry (Fragaria L.), short branches (crowns) develop from dormant axillary buds to form inflorescences and flowers. While this developmental transition contributes greatly to perenniality and yield in strawberry, its regulatory mechanism remains unclear and understudied. In the woodland strawberry (Fragaria vesca), we identified and characterized 2 independent mutants showing more crowns. Both mutant alleles reside in FveMYB117a, a R2R3-MYB transcription factor gene highly expressed in shoot apical meristems, axillary buds, and young leaves. Transcriptome analysis revealed that the expression of several cytokinin pathway genes was altered in the fvemyb117a mutant. Consistently, active cytokinins were significantly increased in the axillary buds of the fvemyb117a mutant. Exogenous application of cytokinin enhanced crown outgrowth in the wild type, whereas the cytokinin inhibitors suppressed crown outgrowth in the fvemyb117a mutant. FveMYB117a binds directly to the promoters of the cytokinin homeostasis genes FveIPT2 encoding an isopentenyltransferase and FveCKX1 encoding a cytokinin oxidase to regulate their expression. Conversely, the type-B Arabidopsis response regulators FveARR1 and FveARR2b can directly inhibit the expression of FveMYB117a, indicative of a negative feedback regulation. In conclusion, we identified FveMYB117a as a key repressor of crown outgrowth by inhibiting cytokinin accumulation and provide a mechanistic basis for bud fate transition in an herbaceous perennial plant.

ENHANCER OF SHOOT REGENERATION1 promotes de novo root organogenesis after wounding in Arabidopsis leaf explants.

Plants have an astonishing ability to regenerate new organs after wounding. Here, we report that the wound-inducible transcription factor ENHANCER OF SHOOT REGENERATION1 (ESR1) has a dual mode of action in activating ANTHRANILATE SYNTHASE ALPHA SUBUNIT1 (ASA1) expression to ensure auxin-dependent de novo root organogenesis locally at wound sites of Arabidopsis (Arabidopsis thaliana) leaf explants. In the first mode, ESR1 interacts with HISTONE DEACETYLASE6 (HDA6), and the ESR1-HDA6 complex directly binds to the JASMONATE-ZIM DOMAIN5 (JAZ5) locus, inhibiting JAZ5 expression through histone H3 deacetylation. As JAZ5 interferes with the action of ETHYLENE RESPONSE FACTOR109 (ERF109), the transcriptional repression of JAZ5 at the wound site allows ERF109 to activate ASA1 expression. In the second mode, the ESR1 transcriptional activator directly binds to the ASA1 promoter to enhance its expression. Overall, our findings indicate that the dual biochemical function of ESR1, which specifically occurs near wound sites of leaf explants, maximizes local auxin biosynthesis and de novo root organogenesis in Arabidopsis.

Multiplexed in situ hybridization reveals distinct lineage identities for major and minor vein initiation during maize leaf development.

Leaves of flowering plants are characterized by diverse venation patterns. Patterning begins with the selection of vein-forming procambial initial cells from within the ground meristem of a developing leaf, a process which is considered to be auxin-dependent, and continues until veins are anatomically differentiated with functional xylem and phloem. At present, the mechanisms responsible for leaf venation patterning are primarily characterized in the model eudicot Arabidopsis thaliana which displays a reticulate venation network. However, evidence suggests that vein development may proceed via a different mechanism in monocot leaves where venation patterning is parallel. Here, we employed Molecular Cartography, a multiplexed in situ hybridization technique, to analyze the spatiotemporal localization of a subset of auxin-related genes and candidate regulators of vein patterning in maize leaves. We show how different combinations of auxin influx and efflux transporters are recruited during leaf and vein specification and how major and minor vein ranks develop with distinct identities. The localization of the procambial marker PIN1a and the spatial arrangement of procambial initial cells that give rise to major and minor vein ranks further suggests that vein spacing is prepatterned across the medio-lateral leaf axis prior to accumulation of the PIN1a auxin transporter. In contrast, patterning in the adaxial-abaxial axis occurs progressively, with markers of xylem and phloem gradually becoming polarized as differentiation proceeds. Collectively, our data suggest that both lineage- and position-based mechanisms may underpin vein patterning in maize leaves.

Genetic improvement of phosphate-limited photosynthesis for high yield in rice.

Low phosphate (Pi) availability decreases photosynthesis, with phosphate limitation of photosynthesis occurring particularly during grain filling of cereal crops; however, effective genetic solutions remain to be established. We previously discovered that rice phosphate transporter OsPHO1;2 controls seed (sink) development through Pi reallocation during grain filling. Here, we find that OsPHO1;2 regulates Pi homeostasis and thus photosynthesis in leaves (source). Loss-of-function of OsPHO1;2 decreased Pi levels in leaves, leading to decreased photosynthetic electron transport activity, CO2 assimilation rate, and early occurrence of phosphate-limited photosynthesis. Interestingly, ectopic expression of OsPHO1;2 greatly increased Pi availability, and thereby, increased photosynthetic rate in leaves during grain filling, contributing to increased yield. This was supported by the effect of foliar Pi application. Moreover, analysis of core rice germplasm resources revealed that higher OsPHO1;2 expression was associated with enhanced photosynthesis and yield potential compared to those with lower expression. These findings reveal that phosphate-limitation of photosynthesis can be relieved via a genetic approach, and the OsPHO1;2 gene can be employed to reinforce crop breeding strategies for achieving higher photosynthetic efficiency.

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.

Membralin is required for maize development and defines a branch of the endoplasmic reticulum-associated degradation pathway in plants.

Endoplasmic reticulum (ER)-associated degradation (ERAD) plays key roles in controlling protein levels and quality in eukaryotes. The Ring Finger Protein 185 (RNF185)/membralin ubiquitin ligase complex was recently identified as a branch in mammals and is essential for neuronal function, but its function in plant development is unknown. Here, we report the map-based cloning and characterization of Narrow Leaf and Dwarfism 1 (NLD1), which encodes the ER membrane-localized protein membralin and specifically interacts with maize homologs of RNF185 and related components. The nld1 mutant shows defective leaf and root development due to reduced cell number. The defects of nld1 were largely restored by expressing membralin genes from Arabidopsis thaliana and mice, highlighting the conserved roles of membralin proteins in animals and plants. The excessive accumulation of ß-hydroxy ß-methylglutaryl-CoA reductase in nld1 indicates that the enzyme is a membralin-mediated ERAD target. The activation of bZIP60 mRNA splicing-related unfolded protein response signaling and marker gene expression in nld1, as well as DNA fragment and cell viability assays, indicate that membralin deficiency induces ER stress and cell death in maize, thereby affecting organogenesis. Our findings uncover the conserved, indispensable role of the membralin-mediated branch of the ERAD pathway in plants. In addition, ZmNLD1 contributes to plant architecture in a dose-dependent manner, which can serve as a potential target for genetic engineering to shape ideal plant architecture, thereby enhancing high-density maize yields.

Understanding the Genetic Basis of C4 Kranz Anatomy with a View to Engineering C3 Crops.

One of the most remarkable examples of convergent evolution is the transition from C3 to C4 photosynthesis, an event that occurred on over 60 independent occasions. The evolution of C4 is particularly noteworthy because of the complexity of the developmental and metabolic changes that took place. In most cases, compartmentalized metabolic reactions were facilitated by the development of a distinct leaf anatomy known as Kranz. C4 Kranz anatomy differs from ancestral C3 anatomy with respect to vein spacing patterns across the leaf, cell-type specification around veins, and cell-specific organelle function. Here we review our current understanding of how Kranz anatomy evolved and how it develops, with a focus on studies that are dissecting the underlying genetic mechanisms. This research field has gained prominence in recent years because understanding the genetic regulation of Kranz may enable the C3-to-C4 transition to be engineered, an endeavor that would significantly enhance crop productivity.

The genome formula of a multipartite virus is regulated both at the individual segment and the segment group levels.

Differential accumulation of the distinct genome segments is a common feature of viruses with segmented genomes. The reproducible and specific pattern of genome segment accumulation within the host is referred to as the "genome formula". There is speculation and some experimental support for a functional role of the genome formula by modulating gene expression through copy number variations. However, the mechanisms of genome formula regulation have not yet been identified. In this study, we investigated whether the genome formula of the octopartite nanovirus faba bean necrotic stunt virus (FBNSV) is regulated by processes acting at the individual segment vs. viral population levels. We used a leaf infiltration system to show that the two most accumulated genome segments of the FBNSV possess a greater intrinsic accumulation capacity in Vicia faba tissues than the other segments. Nevertheless, processes acting at the individual segment level are insufficient to generate the genome formula, suggesting the involvement of additional mechanisms acting at the supra-segment level. Indeed, the absence of segments with important functions during systemic infection strongly modifies the relative frequency of the others, indicating that the genome formula is a property of the segment group. Together, these results demonstrate that the FBNSV genome formula is shaped by a complex process acting at both the individual segment and the segment group levels.

SHORT ROOT and INDETERMINATE DOMAIN family members govern PIN-FORMED expression to regulate minor vein differentiation in rice.

C3 and C4 grasses directly and indirectly provide the vast majority of calories to the human diet, yet our understanding of the molecular mechanisms driving photosynthetic productivity in grasses is largely unexplored. Ground meristem cells divide to form mesophyll or vascular initial cells early in leaf development in C3 and C4 grasses. Here we define a genetic circuit composed of SHORT ROOT (SHR), INDETERMINATE DOMAIN (IDD), and PIN-FORMED (PIN) family members that specifies vascular identify and ground cell proliferation in leaves of both C3 and C4 grasses. Ectopic expression and loss-of-function mutant studies of SHR paralogs in the C3 plant Oryza sativa (rice) and the C4 plant Setaria viridis (green millet) revealed the roles of these genes in both minor vein formation and ground cell differentiation. Genetic and in vitro studies further suggested that SHR regulates this process through its interactions with IDD12 and 13. We also revealed direct interactions of these IDD proteins with a putative regulatory element within the auxin transporter gene PIN5c. Collectively, these findings indicate that a SHR-IDD regulatory circuit mediates auxin transport by negatively regulating PIN expression to modulate minor vein patterning in the grasses.

Topological data analysis reveals a core gene expression backbone that defines form and function across flowering plants.

Since they emerged approximately 125 million years ago, flowering plants have evolved to dominate the terrestrial landscape and survive in the most inhospitable environments on earth. At their core, these adaptations have been shaped by changes in numerous, interconnected pathways and genes that collectively give rise to emergent biological phenomena. Linking gene expression to morphological outcomes remains a grand challenge in biology, and new approaches are needed to begin to address this gap. Here, we implemented topological data analysis (TDA) to summarize the high dimensionality and noisiness of gene expression data using lens functions that delineate plant tissue and stress responses. Using this framework, we created a topological representation of the shape of gene expression across plant evolution, development, and environment for the phylogenetically diverse flowering plants. The TDA-based Mapper graphs form a well-defined gradient of tissues from leaves to seeds, or from healthy to stressed samples, depending on the lens function. This suggests that there are distinct and conserved expression patterns across angiosperms that delineate different tissue types or responses to biotic and abiotic stresses. Genes that correlate with the tissue lens function are enriched in central processes such as photosynthetic, growth and development, housekeeping, or stress responses. Together, our results highlight the power of TDA for analyzing complex biological data and reveal a core expression backbone that defines plant form and function.

The strength and pattern of natural selection on gene expression in rice.

Levels of gene expression underpin organismal phenotypes1,2, but the nature of selection that acts on gene expression and its role in adaptive evolution remain unknown1,2. Here we assayed gene expression in rice (Oryza sativa)3, and used phenotypic selection analysis to estimate the type and strength of selection on the levels of more than 15,000 transcripts4,5. Variation in most transcripts appears (nearly) neutral or under very weak stabilizing selection in wet paddy conditions (with median standardized selection differentials near zero), but selection is stronger under drought conditions. Overall, more transcripts are conditionally neutral (2.83%) than are antagonistically pleiotropic6 (0.04%), and transcripts that display lower levels of expression and stochastic noise7-9 and higher levels of plasticity9 are under stronger selection. Selection strength was further weakly negatively associated with levels of cis-regulation and network connectivity9. Our multivariate analysis suggests that selection acts on the expression of photosynthesis genes4,5, but that the efficacy of selection is genetically constrained under drought conditions10. Drought selected for earlier flowering11,12 and a higher expression of OsMADS18 (Os07g0605200), which encodes a MADS-box transcription factor and is a known regulator of early flowering13-marking this gene as a drought-escape gene11,12. The ability to estimate selection strengths provides insights into how selection can shape molecular traits at the core of gene action.

A plant genetic network for preventing dysbiosis in the phyllosphere.

The aboveground parts of terrestrial plants, collectively called the phyllosphere, have a key role in the global balance of atmospheric carbon dioxide and oxygen. The phyllosphere represents one of the most abundant habitats for microbiota colonization. Whether and how plants control phyllosphere microbiota to ensure plant health is not well understood. Here we show that the Arabidopsis quadruple mutant (min7 fls2 efr cerk1; hereafter, mfec)1, simultaneously defective in pattern-triggered immunity and the MIN7 vesicle-trafficking pathway, or a constitutively activated cell death1 (cad1) mutant, carrying a S205F mutation in a membrane-attack-complex/perforin (MACPF)-domain protein, harbour altered endophytic phyllosphere microbiota and display leaf-tissue damage associated with dysbiosis. The Shannon diversity index and the relative abundance of Firmicutes were markedly reduced, whereas Proteobacteria were enriched in the mfec and cad1S205F mutants, bearing cross-kingdom resemblance to some aspects of the dysbiosis that occurs in human inflammatory bowel disease. Bacterial community transplantation experiments demonstrated a causal role of a properly assembled leaf bacterial community in phyllosphere health. Pattern-triggered immune signalling, MIN7 and CAD1 are found in major land plant lineages and are probably key components of a genetic network through which terrestrial plants control the level and nurture the diversity of endophytic phyllosphere microbiota for survival and health in a microorganism-rich environment.

Transcriptional regulation of strigolactone signalling in Arabidopsis.

Plant hormones known as strigolactones control plant development and interactions between host plants and symbiotic fungi or parasitic weeds1-4. In Arabidopsis thaliana and rice, the proteins DWARF14 (D14), MORE AXILLARY GROWTH 2 (MAX2), SUPPRESSOR OF MAX2-LIKE 6, 7 and 8 (SMXL6, SMXL7 and SMXL8) and their orthologues form a complex upon strigolactone perception and play a central part in strigolactone signalling5-10. However, whether and how strigolactones activate downstream transcription remains largely unknown. Here we use a synthetic strigolactone to identify 401 strigolactone-responsive genes in Arabidopsis, and show that these plant hormones regulate shoot branching, leaf shape and anthocyanin accumulation mainly through transcriptional activation of the BRANCHED 1, TCP DOMAIN PROTEIN 1 and PRODUCTION OF ANTHOCYANIN PIGMENT 1 genes. We find that SMXL6 targets 729 genes in the Arabidopsis genome and represses the transcription of SMXL6, SMXL7 and SMXL8 by binding directly to their promoters, showing that SMXL6 serves as an autoregulated transcription factor to maintain the homeostasis of strigolactone signalling. These findings reveal an unanticipated mechanism through which a transcriptional repressor of hormone signalling can directly recognize DNA and regulate transcription in higher plants.

Injury Activates Ca2+/Calmodulin-Dependent Phosphorylation of JAV1-JAZ8-WRKY51 Complex for Jasmonate Biosynthesis.

Insect herbivory causes severe damage to plants and threatens the world's food production. During evolutionary adaptation, plants have evolved sophisticated mechanisms to rapidly accumulate a key defense hormone, jasmonate (JA), that triggers plant defense against herbivory. However, little is known about how plants initially activate JA biosynthesis at encounter with herbivory. Here, we uncover that a novel JAV1-JAZ8-WRKY51 (JJW) complex controls JA biosynthesis to defend against insect attack. In healthy plants, the JJW complex represses JA biosynthesis to restrain JA at a low basal level to ensure proper plant growth. When plants are injured by insect attack, injury rapidly triggers calcium influxes to activate calmodulin-dependent phosphorylation of JAV1, which disintegrates JJW complex and activates JA biosynthesis, giving rise to the rapid burst of JA for plant defense. Our findings offer new insights into the highly sophisticated defense systems evolved by plants to defend against herbivory.

Age-dependent seasonal growth cessation in Populus.

In temperate and boreal regions, perennial plants adapt their annual growth cycle to the change of seasons. In natural forests, juvenile seedlings usually display longer growth seasons compared to adult trees to ensure their establishment and survival under canopy shade. However, how trees adjust their annual growth according to their age is not known. In this study, we show that age-dependent seasonal growth cessation is genetically controlled and found that the miR156-SPL3/5 module, a key regulon of vegetative phase change (VPC), also triggers age-dependent growth cessation in Populus trees. We show that miR156 promotes shoot elongation during vegetative growth, and its targets SPL3/5s function in the same pathway but as repressors. We find that the miR156-SPL3/5s regulon controls growth cessation in both leaves and shoot apices and through multiple pathways, but with a different mechanism compared to how the miR156-SPL regulon controls VPC in annual plants. Taken together, our results reveal an age-dependent genetic network in mediating seasonal growth cessation, a key phenological process in the climate adaptation of perennial trees.