Spatial transcriptional signatures define margin morphogenesis along the proximal-distal and medio-lateral axes in tomato (Solanum lycopersicum) leaves.

Leaf morphogenesis involves cell division, expansion, and differentiation in the developing leaf, which take place at different rates and at different positions along the medio-lateral and proximal-distal leaf axes. The gene expression changes that control cell fate along these axes remain elusive due to difficulties in precisely isolating tissues. Here, we combined rigorous early leaf characterization, laser capture microdissection, and transcriptomic sequencing to ask how gene expression patterns regulate early leaf morphogenesis in wild-type tomato (Solanum lycopersicum) and the leaf morphogenesis mutant trifoliate. We observed transcriptional regulation of cell differentiation along the proximal-distal axis and identified molecular signatures delineating the classically defined marginal meristem/blastozone region during early leaf development. We describe the role of endoreduplication during leaf development, when and where leaf cells first achieve photosynthetic competency, and the regulation of auxin transport and signaling along the leaf axes. Knockout mutants of BLADE-ON-PETIOLE2 exhibited ectopic shoot apical meristem formation on leaves, highlighting the role of this gene in regulating margin tissue identity. We mapped gene expression signatures in specific leaf domains and evaluated the role of each domain in conferring indeterminacy and permitting blade outgrowth. Finally, we generated a global gene expression atlas of the early developing compound leaf.

The importance of habitat type and historical fire regimes in arthropod community response following large-scale wildfires.

Novel wildfire regimes are rapidly changing global ecosystems and pose significant challenges for biodiversity conservation and ecosystem management. In this study, we used DNA metabarcoding to assess the response of arthropod pollinator communities to large-scale wildfires across diverse habitat types in California. We sampled six reserves within the University of California Natural Reserve System, each of which was partially burned in the 2020 Lightning Complex wildfires in California. Using yellow pan traps to target pollinators, we collected arthropods from burned and unburned sites across multiple habitat types including oak woodland, redwood, scrub, chamise, grassland, forest, and serpentine habitats. We found no significant difference in alpha diversity values between burned and unburned sites; instead, seasonal variations played a significant role in arthropod community dynamics, with the emergence of plant species in Spring promoting increased pollinator richness at all sites. When comparing all sites, we found that burn status was not a significant grouping factor. Instead, compositional differences were largely explained by geographic differences, with distinct communities within each reserve. Within a geographic area, the response of arthropods to fire was dependent on habitat type. While communities in grasslands and oak woodlands exhibited recovery following burn, scrublands experienced substantial changes in community composition. Our study highlights the importance of examining community responses to wildfires across broad spatial scales and diverse habitat types. By understanding the nuanced dynamics of arthropod communities in response to fire disturbances, we can develop effective conservation strategies that promote resilience and maintain biodiversity in the face of increasing wildfire frequency and severity driven by climate change.

Principles for data analysis workflows.

A systematic and reproducible "workflow"-the process that moves a scientific investigation from raw data to coherent research question to insightful contribution-should be a fundamental part of academic data-intensive research practice. In this paper, we elaborate basic principles of a reproducible data analysis workflow by defining 3 phases: the Explore, Refine, and Produce Phases. Each phase is roughly centered around the audience to whom research decisions, methodologies, and results are being immediately communicated. Importantly, each phase can also give rise to a number of research products beyond traditional academic publications. Where relevant, we draw analogies between design principles and established practice in software development. The guidance provided here is not intended to be a strict rulebook; rather, the suggestions for practices and tools to advance reproducible, sound data-intensive analysis may furnish support for both students new to research and current researchers who are new to data-intensive work.

Ten simple rules on writing clean and reliable open-source scientific software.

Functional, usable, and maintainable open-source software is increasingly essential to scientific research, but there is a large variation in formal training for software development and maintainability. Here, we propose 10 "rules" centered on 2 best practice components: clean code and testing. These 2 areas are relatively straightforward and provide substantial utility relative to the learning investment. Adopting clean code practices helps to standardize and organize software code in order to enhance readability and reduce cognitive load for both the initial developer and subsequent contributors; this allows developers to concentrate on core functionality and reduce errors. Clean coding styles make software code more amenable to testing, including unit tests that work best with modular and consistent software code. Unit tests interrogate specific and isolated coding behavior to reduce coding errors and ensure intended functionality, especially as code increases in complexity; unit tests also implicitly provide example usages of code. Other forms of testing are geared to discover erroneous behavior arising from unexpected inputs or emerging from the interaction of complex codebases. Although conforming to coding styles and designing tests can add time to the software development project in the short term, these foundational tools can help to improve the correctness, quality, usability, and maintainability of open-source scientific software code. They also advance the principal point of scientific research: producing accurate results in a reproducible way. In addition to suggesting several tips for getting started with clean code and testing practices, we recommend numerous tools for the popular open-source scientific software languages Python, R, and Julia.

A sister of PIN1 gene in tomato (Solanum lycopersicum) defines leaf and flower organ initiation patterns by maintaining epidermal auxin flux.

The spatiotemporal localization of the plant hormone auxin acts as a positional cue during early leaf and flower organogenesis. One of the main contributors to auxin localization is the auxin efflux carrier PIN-FORMED1 (PIN1). Phylogenetic analysis has revealed that PIN1 genes are split into two sister clades; PIN1 and the relatively uncharacterized Sister-Of-PIN1 (SoPIN1). In this paper we identify entire-2 as a loss-of-function SlSoPIN1a (Solyc10g078370) mutant in Solanum lycopersicum. The entire-2 plants are unable to specify proper leaf initiation leading to a frequent switch from the wild type spiral phyllotactic pattern to distichous and decussate patterns. Leaves in entire-2 are large and less complex and the leaflets display spatial deformities in lamina expansion, vascular development, and margin specification. During sympodial growth in entire-2 the specification of organ position and identity is greatly affected resulting in variable branching patterns on the main sympodial and inflorescence axes. To understand how SlSoPIN1a functions in establishing proper auxin maxima we used the auxin signaling reporter DR5: Venus to visualize differences in auxin localization between entire-2 and wild type. DR5: Venus visualization shows a widening of auxin localization which spreads to subepidermal tissue layers during early leaf and flower organogenesis, showing that SoPIN1 functions to focus auxin signaling to the epidermal layer. The striking spatial deformities observed in entire-2 help provide a mechanistic framework for explaining the function of the SoPIN1 clade in S.lycopersicum.

Resolving distinct genetic regulators of tomato leaf shape within a heteroblastic and ontogenetic context.

Leaf shape is mutable, changing in ways modulated by both development and environment within genotypes. A complete model of leaf phenotype would incorporate the changes in leaf shape during juvenile-to-adult phase transitions and the ontogeny of each leaf. Here, we provide a morphometric description of >33,000 leaflets from a set of tomato (Solanum spp) introgression lines grown under controlled environment conditions. We first compare the shape of these leaves, arising during vegetative development, with >11,000 previously published leaflets from a field setting and >11,000 leaflets from wild tomato relatives. We then quantify the changes in shape, across ontogeny, for successive leaves in the heteroblastic series. Using principal component analysis, we then separate genetic effects modulating (1) the overall shape of all leaves versus (2) the shape of specific leaves in the series, finding the former more heritable than the latter and comparing quantitative trait loci regulating each. Our results demonstrate that phenotype is highly contextual and that unbiased assessments of phenotype, for quantitative genetic or other purposes, would ideally sample the many developmental and environmental factors that modulate it.

Light-Induced Indeterminacy Alters Shade-Avoiding Tomato Leaf Morphology.

Plants sense the foliar shade of competitors and alter their developmental programs through the shade-avoidance response. Internode and petiole elongation, and changes in overall leaf area and leaf mass per area, are the stereotypical architectural responses to foliar shade in the shoot. However, changes in leaf shape and complexity in response to shade remain incompletely, and qualitatively, described. Using a meta-analysis of more than 18,000 previously published leaflet outlines, we demonstrate that shade avoidance alters leaf shape in domesticated tomato (Solanum lycopersicum) and wild relatives. The effects of shade avoidance on leaf shape are subtle with respect to individual traits but are combinatorially strong. We then seek to describe the developmental origins of shade-induced changes in leaf shape by swapping plants between light treatments. Leaf size is light responsive late into development, but patterning events, such as stomatal index, are irrevocably specified earlier. Observing that shade induces increases in shoot apical meristem size, we then describe gene expression changes in early leaf primordia and the meristem using laser microdissection. We find that in leaf primordia, shade avoidance is not mediated through canonical pathways described in mature organs but rather through the expression of KNOTTED1-LIKE HOMEOBOX and other indeterminacy genes, altering known developmental pathways responsible for patterning leaf shape. We also demonstrate that shade-induced changes in leaf primordium gene expression largely do not overlap with those found in successively initiated leaf primordia, providing evidence against classic hypotheses that shaded leaf morphology results from the prolonged production of juvenile leaf types.

Trifoliate encodes an MYB transcription factor that modulates leaf and shoot architecture in tomato.

Leaf morphology and the pattern of shoot branching determine to a large extent the growth habit of seed plants. Until recently, the developmental processes that led to the establishment of these morphological structures seemed unrelated. Here, we show that the tomato Trifoliate (Tf) gene plays a crucial role in both processes, affecting the formation of leaflets in the compound tomato leaf and the initiation of axillary meristems in the leaf axil. Tf encodes a myeloblastosis oncoprotein (MYB)-like transcription factor related to the Arabidopsis thaliana LATERAL ORGAN FUSION1 (LOF1) and LOF2 proteins. Tf is expressed in the leaf margin, where leaflets are formed, and in the leaf axil, where axillary meristems initiate. During tomato ontogeny, expression of Tf in young leaf primordia increases, correlating with a rise in leaf dissection (heteroblasty). Formation of leaflets and initiation of axillary meristems can be traced back to groups of pluripotent cells. Tf function is required to inhibit differentiation of these cells and thereby to maintain their morphogenetic competence, a fundamental process in plant development. KNOTTED1-LIKE proteins, which are known regulators in tomato leaf dissection, require Tf activity to exert their function in the basal part of the leaf. Similarly, the plant hormone auxin needs Tf activity to initiate the formation of lateral leaflets. Thus, leaf dissection and shoot branching rely on a conserved mechanism that regulates the morphogenetic competence of cells at the leaf margin and in the leaf axil.

A modern ampelography: a genetic basis for leaf shape and venation patterning in grape.

Terroir, the unique interaction between genotype, environment, and culture, is highly refined in domesticated grape (Vitis vinifera). Toward cultivating terroir, the science of ampelography tried to distinguish thousands of grape cultivars without the aid of genetics. This led to sophisticated phenotypic analyses of natural variation in grape leaves, which within a palmate-lobed framework exhibit diverse patterns of blade outgrowth, hirsuteness, and venation patterning. Here, we provide a morphometric analysis of more than 1,200 grape accessions. Elliptical Fourier descriptors provide a global analysis of leaf outlines and lobe positioning, while a Procrustes analysis quantitatively describes venation patterning. Correlation with previous ampelography suggests an important genetic component, which we confirm with estimates of heritability. We further use RNA-Seq of mutant varieties and perform a genome-wide association study to explore the genetic basis of leaf shape. Meta-analysis reveals a relationship between leaf morphology and hirsuteness, traits known to correlate with climate in the fossil record and extant species. Together, our data demonstrate a genetic basis for the intricate diversity present in grape leaves. We discuss the possibility of using grape leaves as a breeding target to preserve terroir in the face of anticipated climate change, a major problem facing viticulture.

Leaf asymmetry as a developmental constraint imposed by auxin-dependent phyllotactic patterning.

In a majority of species, leaf development is thought to proceed in a bilaterally symmetric fashion without systematic asymmetries. This is despite the left and right sides of an initiating primordium occupying niches that differ in their distance from sinks and sources of auxin. Here, we revisit an existing model of auxin transport sufficient to recreate spiral phyllotactic patterns and find previously overlooked asymmetries between auxin distribution and the centers of leaf primordia. We show that it is the direction of the phyllotactic spiral that determines the side of the leaf these asymmetries fall on. We empirically confirm the presence of an asymmetric auxin response using a DR5 reporter and observe morphological asymmetries in young leaf primordia. Notably, these morphological asymmetries persist in mature leaves, and we observe left-right asymmetries in the superficially bilaterally symmetric leaves of tomato (Solanum lycopersicum) and Arabidopsis thaliana that are consistent with modeled predictions. We further demonstrate that auxin application to a single side of a leaf primordium is sufficient to recapitulate the asymmetries we observe. Our results provide a framework to study a previously overlooked developmental axis and provide insights into the developmental constraints imposed upon leaf morphology by auxin-dependent phyllotactic patterning.

Gradual disintegration of the floral symmetry gene network is implicated in the evolution of a wind-pollination syndrome.

Angiosperms exhibit staggering diversity in floral form, and evolution of floral morphology is often correlated with changes in pollination syndrome. The showy, bilaterally symmetrical flowers of the model species Antirrhinum majus (Plantaginaceae) are highly specialized for bee pollination. In A. majus, Cycloidea (CYC), Dichotoma (DICH), Radialis (RAD), and Divaricata (DIV) specify the development of floral bilateral symmetry. However, it is unclear to what extent evolution of these genes has resulted in flower morphological divergence among closely related members of Plantaginaceae differing in pollination syndrome. We compared floral symmetry genes from insect-pollinated Digitalis purpurea, which has bilaterally symmetrical flowers, with those from closely related Aragoa abietina and wind-pollinated Plantago major, both of which have radially symmetrical flowers. We demonstrate that Plantago, but not Aragoa, species have lost a dorsally expressed CYC-like gene and downstream targets RAD and DIV. Furthermore, the single P. major CYC-like gene is expressed across all regions of the flower, similar to expression of its ortholog in closely related Veronica serpyllifolia. We propose that changes in the expression of duplicated CYC-like genes led to the evolution of radial flower symmetry in Aragoa/Plantago, and that further disintegration of the symmetry gene pathway resulted in the wind-pollination syndrome of Plantago. This model underscores the potential importance of gene loss in the evolution of ecologically important traits.

Whole Genome Sequences of 23 Species from the Drosophila montium Species Group (Diptera: Drosophilidae): A Resource for Testing Evolutionary Hypotheses.

Large groups of species with well-defined phylogenies are excellent systems for testing evolutionary hypotheses. In this paper, we describe the creation of a comparative genomic resource consisting of 23 genomes from the species-rich Drosophila montium species group, 22 of which are presented here for the first time. The montium group is well-positioned for clade genomics. Within the montium clade, evolutionary distances are such that large numbers of sequences can be accurately aligned while also recovering strong signals of divergence; and the distance between the montium group and D. melanogaster is short enough so that orthologous sequence can be readily identified. All genomes were assembled from a single, small-insert library using MaSuRCA, before going through an extensive post-assembly pipeline. Estimated genome sizes within the montium group range from 155 Mb to 223 Mb (mean = 196 Mb). The absence of long-distance information during the assembly process resulted in fragmented assemblies, with the scaffold NG50s varying widely based on repeat content and sample heterozygosity (min = 18 kb, max = 390 kb, mean = 74 kb). The total scaffold length for most assemblies is also shorter than the estimated genome size, typically by 5-15%. However, subsequent analysis showed that our assemblies are highly complete. Despite large differences in contiguity, all assemblies contain at least 96% of known single-copy Dipteran genes (BUSCOs, n = 2,799). Similarly, by aligning our assemblies to the D. melanogaster genome and remapping coordinates for a large set of transcriptional enhancers (n = 3,457), we showed that each montium assembly contains orthologs for at least 91% of D. melanogaster enhancers. Importantly, the genic and enhancer contents of our assemblies are comparable to that of far more contiguous Drosophila assemblies. The alignment of our own D. serrata assembly to a previously published PacBio D. serrata assembly also showed that our longest scaffolds (up to 1 Mb) are free of large-scale misassemblies. Our genome assemblies are a valuable resource that can be used to further resolve the montium group phylogeny; study the evolution of protein-coding genes and cis-regulatory sequences; and determine the genetic basis of ecological and behavioral adaptations.

Robust manipulation of the behavior of Drosophila melanogaster by a fungal pathogen in the laboratory.

Many microbes induce striking behavioral changes in their animal hosts, but how they achieve this is poorly understood, especially at the molecular level. Mechanistic understanding has been largely constrained by the lack of an experimental system amenable to molecular manipulation. We recently discovered a strain of the behavior-manipulating fungal pathogen Entomophthora muscae infecting wild Drosophila, and established methods to infect D. melanogaster in the lab. Lab-infected flies manifest the moribund behaviors characteristic of E. muscae infection: hours before death, they climb upward, extend their proboscides, affixing in place, then raise their wings, clearing a path for infectious spores to launch from their abdomens. We found that E. muscae invades the nervous system, suggesting a direct means by which the fungus could induce behavioral changes. Given the vast molecular toolkit available for D. melanogaster, we believe this new system will enable rapid progress in understanding how E. muscae manipulates host behavior.

Morphological Plant Modeling: Unleashing Geometric and Topological Potential within the Plant Sciences.

The geometries and topologies of leaves, flowers, roots, shoots, and their arrangements have fascinated plant biologists and mathematicians alike. As such, plant morphology is inherently mathematical in that it describes plant form and architecture with geometrical and topological techniques. Gaining an understanding of how to modify plant morphology, through molecular biology and breeding, aided by a mathematical perspective, is critical to improving agriculture, and the monitoring of ecosystems is vital to modeling a future with fewer natural resources. In this white paper, we begin with an overview in quantifying the form of plants and mathematical models of patterning in plants. We then explore the fundamental challenges that remain unanswered concerning plant morphology, from the barriers preventing the prediction of phenotype from genotype to modeling the movement of leaves in air streams. We end with a discussion concerning the education of plant morphology synthesizing biological and mathematical approaches and ways to facilitate research advances through outreach, cross-disciplinary training, and open science. Unleashing the potential of geometric and topological approaches in the plant sciences promises to transform our understanding of both plants and mathematics.

Left-right leaf asymmetry in decussate and distichous phyllotactic systems.

Leaves in plants with spiral phyllotaxy exhibit directional asymmetries, such that all the leaves originating from a meristem of a particular chirality are similarly asymmetric relative to each other. Models of auxin flux capable of recapitulating spiral phyllotaxis predict handed auxin asymmetries in initiating leaf primordia with empirically verifiable effects on superficially bilaterally symmetric leaves. Here, we extend a similar analysis of leaf asymmetry to decussate and distichous phyllotaxy. We found that our simulation models of these two patterns predicted mirrored asymmetries in auxin distribution in leaf primordia pairs. To empirically verify the morphological consequences of asymmetric auxin distribution, we analysed the morphology of a tomato sister-of-pin-formed1a (sopin1a) mutant, entire-2, in which spiral phyllotaxy consistently transitions to a decussate state. Shifts in the displacement of leaflets on the left and right sides of entire-2 leaf pairs mirror each other, corroborating predicted model results. We then analyse the shape of more than 800 common ivy (Hedera helix) and more than 3000 grapevine (Vitis and Ampelopsis spp.) leaf pairs and find statistical enrichment of predicted mirrored asymmetries. Our results demonstrate that left-right auxin asymmetries in models of decussate and distichous phyllotaxy successfully predict mirrored asymmetric leaf morphologies in superficially symmetric leaves.This article is part of the themed issue 'Provocative questions in left-right asymmetry'.