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.

Pupylation-Based Proximity-Tagging of FERONIA-Interacting Proteins in Arabidopsis.

The plasma membrane-localized receptor kinase FERONIA (FER) plays critical roles in a remarkable variety of biological processes throughout the life cycle of Arabidopsis thaliana. Revealing the molecular connections of FER that underlie these processes starts with identifying the proteins that interact with FER. We applied pupylation-based interaction tagging (PUP-IT) to survey cellular proteins in proximity to FER, encompassing weak and transient interactions that can be difficult to capture for membrane proteins. We reproducibly identified 581, 115, and 736 specific FER-interacting protein candidates in protoplasts, seedlings, and flowers, respectively. We also confirmed 14 previously characterized FER-interacting proteins. Protoplast transient gene expression expedited the testing of new gene constructs for PUP-IT analyses and the validation of candidate proteins. We verified the proximity labeling of five selected candidates that were not previously characterized as FER-interacting proteins. The PUP-IT method could be a valuable tool to survey and validate protein-protein interactions for targets of interest in diverse subcellular compartments in plants.

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.

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.

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.

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.

Functional characterisation of BnaA02.TOP1α and BnaC02.TOP1α involved in true leaf biomass accumulation in Brassica napus L.

Leaves, as primary photosynthetic organs essential for high crop yield and quality, have attracted significant attention. The functions of DNA topoisomerase 1α (TOP1α) in various biological processes, including leaf development, in Brassica napus remain unknown. Here, four paralogs of BnaTOP1α, namely BnaA01.TOP1α, BnaA02.TOP1α, BnaC01.TOP1α and BnaC02.TOP1α, were identified and cloned in the B. napus inbred line 'K407'. Expression pattern analysis revealed that BnaA02.TOP1α and BnaC02.TOP1α, but not BnaA01.TOP1α and BnaC01.TOP1α, were persistently and highly expressed in B. napus true leaves. Preliminary analysis in Arabidopsis thaliana revealed that BnaA02.TOP1α and BnaC02.TOP1α paralogs, but not BnaA01.TOP1α and BnaC01.TOP1α, performed biological functions. Targeted mutations of four BnaTOP1α paralogs in B. napus using the CRISPR-Cas9 system revealed that BnaA02.TOP1α and BnaC02.TOP1α served as functional paralogs and redundantly promoted true leaf number and size, thereby promoting true leaf biomass accumulation. Moreover, BnaA02.TOP1α modulated the levels of endogenous gibberellins, cytokinins and auxins by indirectly regulating several genes related to their metabolism processes. BnaA02.TOP1α directly activated BnaA03.CCS52A2 and BnaC09.AN3 by facilitating the recruitment of RNA polymerase II and modulating H3K27me3, H3K36me2 and H3K36me3 levels at these loci and indirectly activated the BnaA08.PARL1 expression, thereby positively controlling the true leaf size in B. napus. Additionally, BnaA02.TOP1α indirectly activated the BnaA07.PIN1 expression to positively regulate the true leaf number. These results reveal the important functions of BnaTOP1α and provide insights into the regulatory network controlling true leaf biomass accumulation in B. napus.

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.

Wind speed affects the rate and kinetics of stomatal conductance.

Understanding the relationship between wind speed and gas exchange in plants is a longstanding challenge. Our aim was to investigate the impact of wind speed on maximum rates of gas exchange and the kinetics of stomatal responses. We conducted experiments in different angiosperm and fern species using an infrared gas analyzer equipped with a controlled leaf fan, enabling precise control of the boundary layer conductance. We first showed that the chamber was adequately mixed even at extremely low wind speed (<0.005 m s-1) and evaluated the link between fan speed, wind speed, and boundary layer conductance. We observed that higher wind speeds led to increased gas exchange of both water vapor and CO2, primarily due to the increase in boundary layer conductance. This increase in transpiration subsequently reduced epidermal pressure, leading to stomatal opening. We documented that stomatal opening in response to light was 2.5 times faster at a wind speed of 2 m s-1 compared to minimal wind speed in Vicia faba, while epidermal peels in a buffer with no transpiration exhibited a similar opening rate. The increase in stomatal conductance under high wind was also observed in four angiosperm species under field conditions, but it was not observed in Boston fern (Nephrolepis exaltata), which lacks epidermal mechanical advantage. Our findings highlight the significant impact of boundary layer conductance on determining gas exchange rates and the kinetics of gas exchange responses to environmental changes.

Novel resources to investigate leaf plasmodesmata formation in C3 and C4 monocots.

Plasmodesmata (PD) are nanochannels that facilitate cell-to-cell transport in plants. More productive and photosynthetically efficient C4 plants form more PD at the mesophyll (M)-bundle sheath (BS) interface in their leaves than their less efficient C3 relatives. In C4 leaves, PD play an essential role in facilitating the rapid metabolite exchange between the M and BS cells to operate a biochemical CO2 concentrating mechanism, which increases the CO2 partial pressure at the site of Rubisco in the BS cells and hence photosynthetic efficiency. The genetic mechanism controlling PD formation in C3 and C4 leaves is largely unknown, especially in monocot crops, due to the technical challenge of quantifying these nanostructures with electron microscopy. To address this issue, we have generated stably transformed lines of Oryza sativa (rice, C3) and Setaria viridis (setaria, C4) with fluorescent protein-tagged PD to build the first spatiotemporal atlas of leaf pit field (cluster of PD) density in monocots without the need for electron microscopy. Across leaf development, setaria had consistently more PD connections at the M-BS wall interface than rice while the difference in M-M pit field density varied. While light was a critical trigger of PD formation, cell type and function determined leaf pit field density. Complementary temporal mRNA sequencing and gene co-expression network analysis revealed that the pattern of pit field density correlated with differentially expressed PD-associated genes and photosynthesis-related genes. PD-associated genes identified from our co-expression network analysis are related to cell wall expansion, translation and chloroplast signalling.

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.