The evolution and diversification of oakleaf butterflies.

Oakleaf butterflies in the genus Kallima have a polymorphic wing phenotype, enabling these insects to masquerade as dead leaves. This iconic example of protective resemblance provides an interesting evolutionary paradigm that can be employed to study biodiversity. We integrated multi-omic data analyses and functional validation to infer the evolutionary history of Kallima species and investigate the genetic basis of their variable leaf wing patterns. We find that Kallima butterflies diversified in the eastern Himalayas and dispersed to East and Southeast Asia. Moreover, we find that leaf wing polymorphism is controlled by the wing patterning gene cortex, which has been maintained in Kallima by long-term balancing selection. Our results provide macroevolutionary and microevolutionary insights into a model species originating from a mountain ecosystem.

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.

Adaxial-abaxial bipolar leaf genes encode a putative cytokinin receptor and HD-Zip III, and control the formation of ectopic shoot meristems in rice.

Shoot apical meristems (SAMs) continuously initiate organ formation and maintain pluripotency through dynamic genetic regulations and cell-to-cell communications. The activity of meristems directly affects the plant's structure by determining the number and arrangement of organs and tissues. We have taken a forward genetic approach to dissect the genetic pathway that controls cell differentiation around the SAM. The rice mutants, adaxial-abaxial bipolar leaf 1 and 2 (abl1 and abl2), produce an ectopic leaf that is fused back-to-back with the fourth leaf, the first leaf produced after embryogenesis. The abaxial-abaxial fusion is associated with the formation of an ectopic shoot meristem at the adaxial base of the fourth leaf primordium. We cloned the ABL1 and ABL2 genes of rice by mapping their chromosomal positions. ABL1 encodes OsHK6, a histidine kinase, and ABL2 encodes a transcription factor, OSHB3 (Class III homeodomain leucine zipper). Expression analyses of these mutant genes as well as OSH1, a rice ortholog of the Arabidopsis STM gene, unveiled a regulatory circuit that controls the formation of an ectopic meristem near the SAM at germination.

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.

A phytoplasma effector destabilizes chloroplastic glutamine synthetase inducing chlorotic leaves that attract leafhopper vectors.

Leaf yellowing is a well-known phenotype that attracts phloem-feeding insects. However, it remains unclear how insect-vectored plant pathogens induce host leaf yellowing to facilitate their own transmission by insect vectors. Here, we report that an effector protein secreted by rice orange leaf phytoplasma (ROLP) inhibits chlorophyll biosynthesis and induces leaf yellowing to attract leafhopper vectors, thereby presumably promoting pathogen transmission. This effector, designated secreted ROLP protein 1 (SRP1), first secreted into rice phloem by ROLP, was subsequently translocated to chloroplasts by interacting with the chloroplastic glutamine synthetase (GS2). The direct interaction between SRP1 and GS2 disrupts the decamer formation of the GS2 holoenzyme, attenuating its enzymatic activity, thereby suppressing the synthesis of chlorophyll precursors glutamate and glutamine. Transgenic expression of SRP1 in rice plants decreased GS2 activity and chlorophyll precursor accumulation, finally inducing leaf yellowing. This process is correlated with the previous evidence that the knockout of GS2 expression in rice plants causes a similar yellow chlorosis phenotype. Consistently, these yellowing leaves attracted higher numbers of leafhopper vectors, caused the vectors to probe more frequently, and presumably facilitate more efficient phytoplasma transmission. Together, these results uncover the mechanism used by phytoplasmas to manipulate the leaf color of infected plants for the purpose of enhancing attractiveness to insect vectors.

Dehydration-induced corrugated folding in Rhapis excelsa plant leaves.

Plant leaves, whose remarkable ability for morphogenesis results in a wide range of petal and leaf shapes in response to environmental cues, have inspired scientific studies as well as the development of engineering structures and devices. Although some typical shape changes in plants and the driving force for such shape evolution have been extensively studied, there remain many poorly understood mechanisms, characteristics, and principles associated with the vast array of shape formation of plant leaves in nature. Here, we present a comprehensive study that combines experiment, theory, and numerical simulations of one such topic-the mechanics and mechanisms of corrugated leaf folding induced by differential shrinking in Rhapis excelsa. Through systematic measurements of the dehydration process in sectioned leaves, we identify a linear correlation between change in the leaf-folding angle and water loss. Building on experimental findings, we develop a generalized model that provides a scaling relationship for water loss in sectioned leaves. Furthermore, our study reveals that corrugated folding induced by dehydration in R. excelsa leaves is achieved by the deformation of a structural architecture-the "hinge" cells. Utilizing such connections among structure, morphology, environmental stimuli, and mechanics, we fabricate several biomimetic machines, including a humidity sensor and morphing devices capable of folding in response to dehydration. The mechanisms of corrugated folding in R. excelsa identified in this work provide a general understanding of the interactions between plant leaves and water. The actuation mechanisms identified in this study also provide insights into the rational design of soft machines.

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.

Linking fine root lifespan to root chemical and morphological traits-A global analysis.

Fine root lifespan is a critical trait associated with contrasting root strategies of resource acquisition and protection. Yet, its position within the multidimensional "root economics space" synthesizing global root economics strategies is largely uncertain, and it is rarely represented in frameworks integrating plant trait variations. Here, we compiled the most comprehensive dataset of absorptive median root lifespan (MRL) data including 98 observations from 79 woody species using (mini-)rhizotrons across 40 sites and linked MRL to other plant traits to address questions of the regulators of MRL at large spatial scales. We demonstrate that MRL not only decreases with plant investment in root nitrogen (associated with more metabolically active tissues) but also increases with construction of larger diameter roots which is often associated with greater plant reliance on mycorrhizal symbionts. Although theories linking organ structure and function suggest that root traits should play a role in modulating MRL, we found no correlation between root traits associated with structural defense (root tissue density and specific root length) and MRL. Moreover, fine root and leaf lifespan were globally unrelated, except among evergreen species, suggesting contrasting evolutionary selection between leaves and roots facing contrasting environmental influences above vs. belowground. At large geographic scales, MRL was typically longer at sites with lower mean annual temperature and higher mean annual precipitation. Overall, this synthesis uncovered several key ecophysiological covariates and environmental drivers of MRL, highlighting broad avenues for accurate parametrization of global biogeochemical models and the understanding of ecosystem response to global climate change.

HEXOKINASE1 and glucose-6-phosphate fuel plant growth and development.

As photoautotrophic organisms, plants produce an incredible spectrum of pigments, anti-herbivory compounds, structural materials and energic intermediates. These biosynthetic routes help plants grow, reproduce and mitigate stress. HEXOKINASE1 (HXK1), a metabolic enzyme and glucose sensor, catalyzes the phosphorylation of hexoses, a key introductory step for many of these pathways. However, previous studies have largely focused on the glucose sensing and signaling functions of HXK1, and the importance of the enzyme's catalytic function is only recently being connected to plant development. In this brief Spotlight, we describe the developmental significance of plant HXK1 and its role in plant metabolic pathways, specifically in glucose-6-phosphate production. Furthermore, we describe the emerging connections between metabolism and development and suggest that HXK1 signaling and catalytic activity regulate discrete areas of plant development.

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.

Fossil leaves reveal drivers of herbivore functional diversity during the Cenozoic.

Herbivorous arthropods are the most diverse group of multicellular organisms on Earth. The most discussed drivers of their inordinate taxonomic and functional diversity are high niche availability associated with the diversity of host plants and dense niche packing due to host partitioning among herbivores. However, the relative contributions of these two factors to dynamics in the diversity of herbivores throughout Earth's history remain unresolved. Using fossil data on herbivore-induced leaf damage from across the Cenozoic, we infer quantitative bipartite interaction networks between plants and functional feeding types of herbivores. We fit a general model of diversity to these interaction networks and discover that host partitioning among functional groups of herbivores contributed twice as much to herbivore functional diversity as host diversity. These findings indicate that niche packing primarily shaped the dynamics in the functional diversity of herbivores during the past 66 my. Our study highlights how the fossil record can be used to test fundamental theories of biodiversity and represents a benchmark for assessing the drivers of herbivore functional diversity in modern ecosystems.

Sex-dependent phenological responses to climate vary across species' ranges.

Anthropogenic climate change has significantly altered the flowering times (i.e., phenology) of plants worldwide, affecting their reproduction, survival, and interactions. Recent studies utilizing herbarium specimens have uncovered significant intra- and inter-specific variation in flowering phenology and its response to changes in climate but have mostly been limited to animal-pollinated species. Thus, despite their economic and ecological importance, variation in phenological responses to climate remain largely unexplored among and within wind-pollinated dioecious species and across their sexes. Using both herbarium specimens and volunteer observations of cottonwood (Populus) species, we examined how phenological sensitivity to climate varies across species, their ranges, sexes, and phenophases. The timing of flowering varied significantly across and within species, as did their sensitivity to spring temperature. In particular, male flowering generally happened earlier in the season and was more sensitive to warming than female flowering. Further, the onset of flowering was more sensitive to changes in temperature than leaf out. Increased temporal gaps between male and female flowering time and between the first open flower date and leaf out date were predicted for the future under two climate change scenarios. These shifts will impact the efficacy of sexual reproduction and gene flow among species. Our study demonstrates significant inter- and intra-specific variation in phenology and its responses to environmental cues, across species' ranges, phenophases, and sex, in wind-pollinated species. These variations need to be considered to predict accurately the effects of climate change and assess their ecological and evolutionary consequences.

El Niño-Southern Oscillation forcing on carbon and water cycling in a Bornean tropical rainforest.

Carbon dioxide and water vapor exchanges between tropical forest canopies and the atmosphere through photosynthesis, respiration, and evapotranspiration (ET) influence carbon and water cycling at the regional and global scales. Their inter- and intra-annual variations are sensitive to seasonal rhythms and longer-timescale tropical climatic events. In the present study, we assessed the El Niño-Southern Oscillation (ENSO) influence on ET and on the net ecosystem exchange (NEE), using eddy-covariance flux observations in a Bornean rainforest over a 10-y period (2010-2019) that included several El Niño and La Niña events. From flux model inversions, we inferred ecophysiological properties, notably the canopy stomatal conductance and "big-leaf" maximum carboxylation rate (Vcmax25_BL). Mean ET values were similar between ENSO phases (El Niño, La Niña, and neutral conditions). Conversely, the mean net ecosystem productivity was highest during La Niña events and lowest during El Niño events. Combining Shapley additive explanation calculations for nine controlling factors with a machine-learning algorithm, we determined that the primary factors for ET and NEE in the La Niña and neutral phases were incoming shortwave solar radiation and Vcmax25_BL, respectively, but that canopy stomatal conductance was the most significant factor for both ET and NEE in the El Niño phase. A combined stomatal-photosynthesis model approach further indicated that Vcmax25_BL differences between ENSO phases were the most significant controlling factor for canopy photosynthesis, emphasizing the strong need to account for ENSO-induced ecophysiological factor variations in model projections of the long-term carbon balance in Southeast Asian tropical rainforests.

LMI1 homeodomain protein regulates organ proportions by spatial modulation of endoreduplication.

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.

Multi-omics analyses of sid2 mutant reflect the need of isochorismate synthase ICS1 to cope with sulfur limitation in Arabidopsis thaliana.

The SID2 (SA INDUCTION-DEFICIENT2) gene that encodes ICS1 (isochorismate synthase), plays a central role in salicylic acid biosynthesis in Arabidopsis. The sid2 and NahG (encoding a bacterial SA hydroxylase) overexpressing mutants (NahG-OE) have currently been shown to outperform wild type, presenting delayed leaf senescence, higher plant biomass and better seed yield. When grown under sulfate-limited conditions (low-S), sid2 mutants exhibited early leaf yellowing compared to the NahG-OE, the npr1 mutant affected in SA signaling pathway, and WT. This indicated that the hypersensitivity of sid2 to sulfate limitation was independent of the canonical npr1 SA-signaling pathway. Transcriptomic and proteomic analyses revealed that major changes occurred in sid2 when cultivated under low-S, changes that were in good accordance with early senescence phenotype and showed the exacerbation of stress responses. The sid2 mutants displayed a lower sulfate uptake capacity when cultivated under low-S and lower S concentrations in their rosettes. Higher glutathione concentrations in sid2 rosettes under low-S were in good accordance with the higher abundance of proteins involved in glutathione and ascorbate redox metabolism. Amino acid and lipid metabolisms were also strongly modified in sid2 under low-S. Depletion of total fatty acids in sid2 under low-S was consistent with the fact that S-metabolism plays a central role in lipid synthesis. Altogether, our results show that functional ICS1 is important for plants to cope with S limiting conditions.

Group C MAP kinases phosphorylate MBD10 to regulate ABA-induced leaf senescence in Arabidopsis.

Abscisic acid (ABA) is a phytohormone that promotes leaf senescence in response to environmental stress. We previously identified methyl CpG-binding domain 10 (MBD10) as a phosphoprotein that becomes differentially phosphorylated after ABA treatment in Arabidopsis. ABA-induced leaf senescence was delayed in mbd10 knockout plants but accelerated in MBD10-overexpressing plants, suggesting that MBD10 positively regulates ABA-induced leaf senescence. ABA-induced phosphorylation of MBD10 occurs in planta on Thr-89, and our results demonstrated that Thr-89 phosphorylation is essential for MBD10's function in leaf senescence. The in vivo phosphorylation of Thr-89 in MBD10 was significantly downregulated in a quadruple mutant of group C MAPKs (mpk1/2/7/14), and group C MAPKs directly phosphorylated MBD10 in vitro. Furthermore, mpk1/2/7/14 showed a similar phenotype as seen in mbd10 for ABA-induced leaf senescence, suggesting that group C MAPKs are the cognate kinases of MBD10 for Thr-89. Because group C MAPKs have been reported to function downstream of SnRK2s, our results indicate that group C MAPKs and MBD10 constitute a regulatory pathway for ABA-induced leaf senescence.

ORANGE interplays with TCP7 to regulate endoreduplication and leaf size.

Leaf size is a crucial agronomic trait directly affecting crop yield, which is mainly determined by coordinated cell proliferation, growth, and differentiation. Although endoreduplication is known to be correlated with the onset of cell differentiation and leaf size, the underlying molecular mechanisms are largely unclear. The DnaJ-like zinc finger domain-containing protein ORANGE (OR) was initially demonstrated to confer the massive accumulation of carotenoids in cauliflower curds. However, the cauliflower or mutant also possesses other phenotypes such as smaller curds, smaller leaves with elongated petioles, and delayed flowering. Here, we demonstrated that OR physically interacts with the transcription factor TCP7, which promotes endoreduplication by inducing the expression of the cell cycle gene CYCLIN D 1;1 (CYCD1;1). Overexpression of OR resulted in smaller rosette leaves, whereas the OR-silencing plants had larger rosette leaves than wild-type plants. Our microscopic observations and flow cytometry analysis revealed that the variation in leaf size was a result of different endoreduplication levels. Genetic analyses showed that OR functions antagonistically with TCP7 in regulating the endoreduplication levels in leaf cells. While the expression of OR is induced by TCP7, OR represses the transactivation activity of TCP7 by affecting its binding capability to the TCP-binding motif in the promoter region of CYCD1;1. Through this interaction, OR negatively regulates the expression of CYCD1;1 and reduces the nuclear ploidy level in rosette leaf cells. Our findings provide new insights into the regulatory network of leaf size and also reveal a regulatory circuit controlling endoreduplication in leaf cells.

SlPP2C2 interacts with FZY/SAUR and regulates tomato development via signaling crosstalk of ABA and auxin.

Abscisic acid (ABA) signaling interacts frequently with auxin signaling when it regulates plant development, affecting multiple physiological processes; however, to the best of our knowledge, their interaction during tomato development has not yet been reported. Here, we found that type 2C protein phosphatase (SlPP2C2) interacts with both flavin monooxygenase FZY, an indole-3-acetic acid (IAA) biosynthetic enzyme, and small auxin upregulated RNA (SAUR) of an IAA signaling protein and regulates their activity, thereby affecting the expression of IAA-responsive genes. The expression level of SlPP2C2 was increased by exogenous ABA, IAA, NaCl, or dehydration treatment of fruits, leaves, and seeds, and it decreased in imbibed seeds. Manipulating SlPP2C2 with overexpression, RNA interference, and CRISPR/Cas9-mediated genome editing resulted in pleiotropic changes, such as morphological changes in leaves, stem trichomes, floral organs and fruits, accompanied by alterations in IAA and ABA levels. Furthermore, the RNA-seq analysis indicated that SlPP2C2 regulates the expression of auxin-/IAA-responsive genes in different tissues of tomato. The results demonstrate that SlPP2C2-mediated ABA signaling regulates the development of both vegetative and reproductive organs via interaction with FZY/SAUR, which integrates the cross-talk of ABA and auxin signals during development and affects the expressions of development-related genes in tomato.

Trait association and prediction through integrative k-mer analysis.

Genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) has been widely used to explore genetic controls of phenotypic traits. Alternatively, GWAS can use counts of substrings of length k from longer sequencing reads, k-mers, as genotyping data. Using maize cob and kernel color traits, we demonstrated that k-mer GWAS can effectively identify associated k-mers. Co-expression analysis of kernel color k-mers and genes directly found k-mers from known causal genes. Analyzing complex traits of kernel oil and leaf angle resulted in k-mers from both known and candidate genes. A gene encoding a MADS transcription factor was functionally validated by showing that ectopic expression of the gene led to less upright leaves. Evolution analysis revealed most k-mers positively correlated with kernel oil were strongly selected against in maize populations, while most k-mers for upright leaf angle were positively selected. In addition, genomic prediction of kernel oil, leaf angle, and flowering time using k-mer data resulted in a similarly high prediction accuracy to the standard SNP-based method. Collectively, we showed k-mer GWAS is a powerful approach for identifying trait-associated genetic elements. Further, our results demonstrated the bridging role of k-mers for data integration and functional gene discovery.

The genetic basis of leaf hair development in upland cotton (Gossypium hirsutum).

Trichomes, which originate from the epidermal cell of aerial organs, provide plants with defense and secretion functions. Although numerous genes have been implicated in trichome development, the molecular mechanisms underlying trichome cell formation in plants remain incompletely understood. Here, we using genome-wide association study (GWAS) across 1037 diverse accessions in upland cotton (Gossypium hirsutum) to identify three loci associated with leaf pubescence (hair) amount, located on chromosome A06 (LPA1), A08 (LPA2) and A11 (LPA3), respectively. GhHD1, a previously characterized candidate gene, was identified on LPA1 and encodes an HD-Zip transcription factor. For LPA2 and LPA3, we identified two candidate genes, GhGIR1 and GhGIR2, both encoding proteins with WD40 and RING domains that act as inhibitors of leaf hair formation. Expression analysis revealed that GhHD1 was predominantly expressed in hairy accessions, whereas GhGIR1 and GhGIR2 were expressed in hairless accessions. Silencing GhHD1 or overexpressing GhGIR1 in hairy accessions induced in a hairless phenotype, whereas silencing GhGIR2 in hairless accessions resulted in a hairy phenotype. We also demonstrated that GhHD1 interact with both GhGIR1 and GhGIR2, and GhGIR1 can interact with GhGIR2. Further investigation indicated that GhHD1 functions as a transcriptional activator, binding to the promoters of the GhGIR1 and GhGIR2 to active their expression, whereas GhGIR1 and GhGIR2 can suppress the transcriptional activation of GhHD1. Our findings shed light on the intricate regulatory network involving GhHD1, GhGIR1 and GhGIR2 in the initiation and development of plant epidermal hairs in cotton.