QTL analysis of divergent floral morphology traits between Gilia yorkii and G. capitata

Speciation is a complex process typically accompanied by significant genetic and morphological differences between sister populations. In plants, divergent floral morphologies and pollinator differences can result in reproductive isolation between populations. Here, we explore floral trait differences between two recently diverged species, Gilia yorkii and G. capitata. The distributions of floral traits in parental, F1, and F2 populations are compared, and groups of correlated traits are identified. We describe the genetic architecture of floral traits through a quantitative trait locus (QTL) analysis using an F2 population of 187 individuals. While all identified QTLs were of moderate (10-25%) effect, interestingly, most QTL intervals were non-overlapping, suggesting that, in general, traits do not share a common genetic basis. Our results provide a framework for future identification of genes involved in the evolution of floral morphology.

Intra-leaf modeling of Cannabis leaflet shape produces leaf models that predict genetic and developmental identities

The iconic, palmately compound leaves of Cannabis have attracted significant attention in the past. However, investigations into the genetic basis of leaf shape or its connections to phytochemical composition have yielded inconclusive results. This is partly due to prominent changes in leaflet number within a single plant during development, which has so far prevented the proper use of common morphometric techniques. Here, we present a new method that overcomes the challenge of nonhomologous landmarks in palmate, pinnate, and lobed leaves, using Cannabis as an example. We model corresponding pseudo-landmarks for each leaflet as angle-radius coordinates and model them as a function of leaflet to create continuous polynomial models, bypassing the problems associated with variable number of leaflets between leaves. We analyze 341 leaves from 24 individuals from nine Cannabis accessions. Using 3591 pseudo-landmarks in modeled leaves, we accurately predict accession identity, leaflet number, and relative node number. Intra-leaf modeling offers a rapid, cost-effective means of identifying Cannabis accessions, making it a valuable tool for future taxonomic studies, cultivar recognition, and possibly chemical content analysis and sex identification, in addition to permitting the morphometric analysis of leaves in any species with variable numbers of leaflets or lobes.

From buds to shoots: insights into grapevine development from the Witch's Broom bud sport

BACKGROUND

Bud sports occur spontaneously in plants when new growth exhibits a distinct phenotype from the rest of the parent plant. The Witch's Broom bud sport occurs occasionally in various grapevine (Vitis vinifera) varieties and displays a suite of developmental defects, including dwarf features and reduced fertility. While it is highly detrimental for grapevine growers, it also serves as a useful tool for studying grapevine development. We used the Witch's Broom bud sport in grapevine to understand the developmental trajectories of the bud sports, as well as the potential genetic basis. We analyzed the phenotypes of two independent cases of the Witch's Broom bud sport, in the Dakapo and Merlot varieties of grapevine, alongside wild type counterparts. To do so, we quantified various shoot traits, performed 3D X-ray Computed Tomography on dormant buds, and landmarked leaves from the samples. We also performed Illumina and Oxford Nanopore sequencing on the samples and called genetic variants using these sequencing datasets.

RESULTS

The Dakapo and Merlot cases of Witch's Broom displayed severe developmental defects, with no fruit/clusters formed and dwarf vegetative features. However, the Dakapo and Merlot cases of Witch's Broom studied were also phenotypically different from one another, with distinct differences in bud and leaf development. We identified 968-974 unique genetic mutations in our two Witch's Broom cases that are potential causal variants of the bud sports. Examining gene function and validating these genetic candidates through PCR and Sanger-sequencing revealed one strong candidate mutation in Merlot Witch's Broom impacting the gene GSVIVG01008260001.

CONCLUSIONS

The Witch's Broom bud sports in both varieties studied had dwarf phenotypes, but the two instances studied were also vastly different from one another and likely have distinct genetic bases. Future work on Witch's Broom bud sports in grapevine could provide more insight into development and the genetic pathways involved in grapevine.

Herbarium specimens reveal links between Capsella bursa-pastoris leaf shape and climate

Studies into the evolution and development of leaf shape have connected variation in plant form, function, and fitness. For species with consistent leaf margin features, patterns in leaf architecture are related to both biotic and abiotic factors. However, for species with inconsistent leaf margin features, quantifying leaf shape variation and the effects of environmental factors on leaf shape has proven challenging. To investigate leaf shape variation in species with inconsistent shapes, we analyzed approximately 500 digitized Capsella bursa-pastoris specimens collected throughout the continental U.S. over a 100-year period with geometric morphometric modeling and deterministic techniques. We generated a morphospace of C. bursa-pastoris leaf shapes and modeled leaf shape as a function of environment and time. Our results suggest C. bursa-pastoris leaf shape variation is strongly associated with temperature over the C. bursa-pastoris growing season, with lobing decreasing as temperature increases. While we expected to see changes in variation over time, our results show that level of leaf shape variation is consistent over the 100-year period. Our findings showed that species with inconsistent leaf shape variation can be quantified using geometric morphometric modeling techniques and that temperature is the main environmental factor influencing leaf shape variation.

Topological data analysis reveals core heteroblastic and ontogenetic programs embedded in leaves of grapevine (Vitaceae) and maracuyá (Passifloraceae)

Leaves are often described in language that evokes a single shape. However, embedded in that descriptor is a multitude of latent shapes arising from evolutionary, developmental, environmental, and other effects. These confounded effects manifest at distinct developmental time points and evolve at different tempos. Here, revisiting datasets comprised of thousands of leaves of vining grapevine (Vitaceae) and maracuyá (Passifloraceae) species, we apply a technique from the mathematical field of topological data analysis to comparatively visualize the structure of heteroblastic and ontogenetic effects on leaf shape in each group. Consistent with a morphologically closer relationship, members of the grapevine dataset possess strong core heteroblasty and ontogenetic programs with little deviation between species. Remarkably, we found that most members of the maracuyá family also share core heteroblasty and ontogenetic programs despite dramatic species-to-species leaf shape differences. This conservation was not initially detected using traditional analyses such as principal component analysis or linear discriminant analysis. We also identify two morphotypes of maracuyá that deviate from the core structure, suggesting the evolution of new developmental properties in this phylogenetically distinct sub-group. Our findings illustrate how topological data analysis can be used to disentangle previously confounded developmental and evolutionary effects to visualize latent shapes and hidden relationships, even ones embedded in complex, high-dimensional datasets.

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.

A critical analysis of plant science literature reveals ongoing inequities

The field of plant science has grown dramatically in the past two decades, but global disparities and systemic inequalities persist. Here, we analyzed ~300,000 papers published over the past two decades to quantify disparities across nations, genders, and taxonomy in the plant science literature. Our analyses reveal striking geographical biases-affluent nations dominate the publishing landscape and vast areas of the globe have virtually no footprint in the literature. Authors in Northern America are cited nearly twice as many times as authors based in Sub-Saharan Africa and Latin America, despite publishing in journals with similar impact factors. Gender imbalances are similarly stark and show remarkably little improvement over time. Some of the most affluent nations have extremely male biased publication records, despite supposed improvements in gender equality. In addition, we find that most studies focus on economically important crop and model species, and a wealth of biodiversity is underrepresented in the literature. Taken together, our analyses reveal a problematic system of publication, with persistent imbalances that poorly capture the global wealth of scientific knowledge and biological diversity. We conclude by highlighting disparities that can be addressed immediately and offer suggestions for long-term solutions to improve equity in the plant sciences.

X-ray imaging of 30 year old wine grape wood reveals cumulative impacts of rootstocks on scion secondary growth and Ravaz index

Annual rings from 30 year old vines in a California rootstock trial were measured to determine the effects of 15 different rootstocks on Chardonnay and Cabernet Sauvignon scions. Viticultural traits measuring vegetative growth, yield, berry quality, and nutrient uptake were collected at the beginning (1995 to 1999) and end (2017 to 2020) of the lifetime of a vineyard initially planted in 1991 and removed in 2021. X-ray Computed Tomography (CT) was used to measure ring widths in 103 vines. Ring width was modeled as a function of ring number using a negative exponential model. Early and late wood ring widths, cambium width, and scion trunk radius were correlated with 27 traits. Modeling of annual ring width shows that scions alter the width of the first rings but that rootstocks alter the decay of later rings, consistently shortening ring width throughout the lifetime of the vine. Ravaz index, juice pH, photosynthetic assimilation and transpiration rates, and instantaneous water use efficiency are correlated with scion trunk radius. Ultimately, our research indicates that rootstocks modulate secondary growth over years, altering physiology and agronomic traits. Rootstocks act in similar but distinct ways from climate to modulate ring width, which borrowing techniques from dendrochronology, can be used to monitor both genetic and environmental effects in woody perennial crop species.

Increases in vein length compensate for leaf area lost to lobing in grapevine

PREMISE

Leaf lobing and leaf size vary considerably across and within species, including among grapevines (Vitis spp.), some of the best-studied leaves. We examined the relationship between leaf lobing and leaf area across grapevine populations that varied in extent of leaf lobing.

METHODS

We used homologous landmarking techniques to measure 2632 leaves across 2 years in 476 unique, genetically distinct grapevines from five biparental crosses that vary primarily in the extent of lobing. We determined to what extent leaf area explained variation in lobing, vein length, and vein to blade ratio.

RESULTS

Although lobing was the primary source of variation in shape across the leaves we measured, leaf area varied only slightly as a function of lobing. Rather, leaf area increases as a function of total major vein length, total branching vein length, and vein to blade ratio. These relationships are stronger for more highly lobed leaves, with the residuals for each model differing as a function of distal lobing.

CONCLUSIONS

For leaves with different extents of lobing but the same area, the more highly lobed leaves have longer veins and higher vein to blade ratios, allowing them to maintain similar leaf areas despite increased lobing. These findings show how more highly lobed leaves may compensate for what would otherwise result in a reduced leaf area, allowing for increased photosynthetic capacity through similar leaf size.

A predicted developmental and evolutionary morphospace for grapevine leaves

Using conventional statistical approaches there exist powerful methods to classify shapes. Embedded in morphospaces is information that allows us to visualise theoretical leaves. These unmeasured leaves are never considered nor how the negative morphospace can inform us about the forces responsible for shaping leaf morphology. Here, we model leaf shape using an allometric indicator of leaf size, the ratio of vein to blade areas. The borders of the observable morphospace are restricted by constraints and define an orthogonal grid of developmental and evolutionary effects which can predict the shapes of possible grapevine leaves. Leaves in the genus Vitis are found to fully occupy morphospace available to them. From this morphospace, we predict the developmental and evolutionary shapes of grapevine leaves that are not only possible, but exist, and argue that rather than explaining leaf shape in terms of discrete nodes or species, that a continuous model is more appropriate.

Multi-dimensional leaf phenotypes reflect root system genotype in grafted grapevine over the growing season

BACKGROUND

Modern biological approaches generate volumes of multi-dimensional data, offering unprecedented opportunities to address biological questions previously beyond reach owing to small or subtle effects. A fundamental question in plant biology is the extent to which below-ground activity in the root system influences above-ground phenotypes expressed in the shoot system. Grafting, an ancient horticultural practice that fuses the root system of one individual (the rootstock) with the shoot system of a second, genetically distinct individual (the scion), is a powerful experimental system to understand below-ground effects on above-ground phenotypes. Previous studies on grafted grapevines have detected rootstock influence on scion phenotypes including physiology and berry chemistry. However, the extent of the rootstock's influence on leaves, the photosynthetic engines of the vine, and how those effects change over the course of a growing season, are still largely unknown.

RESULTS

Here, we investigate associations between rootstock genotype and shoot system phenotypes using 5 multi-dimensional leaf phenotyping modalities measured in a common grafted scion: ionomics, metabolomics, transcriptomics, morphometrics, and physiology. Rootstock influence is ubiquitous but subtle across modalities, with the strongest signature of rootstock observed in the leaf ionome. Moreover, we find that the extent of rootstock influence on scion phenotypes and patterns of phenomic covariation are highly dynamic across the season.

CONCLUSIONS

These findings substantially expand previously identified patterns to demonstrate that rootstock influence on scion phenotypes is complex and dynamic and underscore that broad understanding necessitates volumes of multi-dimensional data previously unmet.

Grapevine rootstocks affect growth-related scion phenotypes

Grape growers use rootstocks to provide protection against pests and pathogens and to modulate viticulture performance such as shoot growth. Our study examined two grapevine scion varieties ('Chardonnay' and 'Cabernet Sauvignon') grafted to 15 different rootstocks and determined the effect of rootstocks on eight traits important to viticulture. We assessed the vines across five years and identified both year and variety as contributing strongly to trait variation. The effect of rootstock was relatively consistent across years and varieties, explaining between 8.99% and 9.78% of the variation in growth-related traits including yield, pruning weight, berry weight and Ravaz index (yield to pruning weight ratio). Increases in yield due to rootstock were generally the result of increases in berry weight, likely due to increased water uptake by vines grafted to a particular rootstock. We demonstrated a greater than 50% increase in yield, pruning weight, or Ravaz index by choosing the optimal rootstock, indicating that rootstock choice is crucial for grape growers looking to improve vine performance.

Vein-to-blade ratio is an allometric indicator of leaf size and plasticity

PREMISE

As a leaf expands, its shape dynamically changes. Previously, we documented an allometric relationship between vein and blade area in grapevine leaves. Larger leaves have a smaller ratio of primary and secondary vein area relative to blade area compared to smaller leaves. We sought to use allometry as an indicator of leaf size and plasticity.

METHODS

We measured the ratio of vein-to-blade area from the same 208 vines across four growing seasons (2013, 2015, 2016, and 2017). Matching leaves by vine and node, we analyzed the correlation between the size and shape of grapevine leaves as repeated measures with climate variables across years.

RESULTS

The proportion of leaf area occupied by vein and blade exponentially decreased and increased, respectively, during leaf expansion making their ratio a stronger indicator of leaf size than area itself. Total precipitation and leaf wetness hours of the previous year but not the current showed strong negative correlations with vein-to-blade ratio, whereas maximum air temperature from the previous year was positively correlated.

CONCLUSIONS

Our results demonstrate that vein-to-blade ratio is a strong allometric indicator of leaf size and plasticity in grapevines measured across years. Grapevine leaf primordia are initiated in buds the year before they emerge, and we found that total precipitation and maximum air temperature of the previous growing season exerted the largest statistically significant effects on leaf morphology. Vein-to-blade ratio is a promising allometric indicator of relationships between leaf morphology and climate, the robustness of which should be explored further.

Directional phyllotactic bias in calatheas (Goeppertia, Marantaceae): A citizen science approach

Lateral organs arranged in spiral phyllotaxy are separated by the golden angle, ≈137.5°, leading to chirality: either clockwise or counter-clockwise. In some species, leaves are asymmetric such that they are smaller and curved towards the side ascending the phyllotactic spiral. As such, these asymmetries lead to mirroring of leaf shapes in plants of opposite phyllotactic handedness. Previous reports had suggested that the pin-stripe calathea (Goeppertia ornata) may be exclusively of one phyllotactic direction, counter-clockwise, but had limited sampling to a single population. Here, we use a citizen science approach leveraging a social media poll, internet image searches, in-person verification at nurseries in four countries and digitally-curated, research-grade observations to demonstrate that calatheas (Goeppertia spp.) around the world are biased towards counter-clockwise phyllotaxy. The possibility that this bias is genetic and its implications for models of phyllotaxy that assume handedness is stochastically specified in equal proportions is discussed.

Composite modeling of leaf shape along shoots discriminates Vitis species better than individual leaves

PREMISE

Leaf morphology is dynamic, continuously deforming during leaf expansion and among leaves within a shoot. Here, we measured the leaf morphology of more than 200 grapevines (Vitis spp.) over four years and modeled changes in leaf shape along the shoot to determine whether a composite leaf shape comprising all the leaves from a single shoot can better capture the variation and predict species identity compared with individual leaves.

METHODS

Using homologous universal landmarks found in grapevine leaves, we modeled various morphological features as polynomial functions of leaf nodes. The resulting functions were used to reconstruct modeled leaf shapes across the shoots, generating composite leaves that comprehensively capture the spectrum of leaf morphologies present.

RESULTS

We found that composite leaves are better predictors of species identity than individual leaves from the same plant. We were able to use composite leaves to predict the species identity of previously unassigned grapevines, which were verified with genotyping.

DISCUSSION

Observations of individual leaf shape fail to capture the true diversity between species. Composite leaf shape-an assemblage of modeled leaf snapshots across the shoot-is a better representation of the dynamic and essential shapes of leaves, in addition to serving as a better predictor of species identity than individual leaves.

Leaf shape is a predictor of fruit quality and cultivar performance in tomato

Commercial tomato (Solanum lycopersicum) is one of the most widely grown vegetable crops worldwide. Heirloom tomatoes retain extensive genetic diversity and a considerable range of fruit quality and leaf morphological traits. Here the role of leaf morphology was investigated for its impact on fruit quality. Heirloom cultivars were grown in field conditions, and BRIX by yield (BY) and other traits were measured over a 14-wk period. The complex relationships among these morphological and physiological traits were evaluated using partial least-squares path modeling, and a consensus model was developed. Photosynthesis contributed strongly to vegetative biomass and sugar content of fruits but had a negative impact on yield. Conversely leaf shape, specifically rounder leaves, had a strong positive impact on both fruit sugar content and yield. Cultivars such as Stupice and Glacier, with very round leaves, had the highest performance in both fruit sugar and yield. Our model accurately predicted BY for two commercial cultivars using leaf shape data as input. This study revealed the importance of leaf shape to fruit quality in tomato, with rounder leaves having significantly improved fruit quality. This correlation was maintained across a range of diverse genetic backgrounds and shows the importance of leaf morphology in tomato crop improvement.

The shape of things to come: Topological data analysis and biology, from molecules to organisms

Shape is data and data is shape. Biologists are accustomed to thinking about how the shape of biomolecules, cells, tissues, and organisms arise from the effects of genetics, development, and the environment. Less often do we consider that data itself has shape and structure, or that it is possible to measure the shape of data and analyze it. Here, we review applications of topological data analysis (TDA) to biology in a way accessible to biologists and applied mathematicians alike. TDA uses principles from algebraic topology to comprehensively measure shape in data sets. Using a function that relates the similarity of data points to each other, we can monitor the evolution of topological features-connected components, loops, and voids. This evolution, a topological signature, concisely summarizes large, complex data sets. We first provide a TDA primer for biologists before exploring the use of TDA across biological sub-disciplines, spanning structural biology, molecular biology, evolution, and development. We end by comparing and contrasting different TDA approaches and the potential for their use in biology. The vision of TDA, that data are shape and shape is data, will be relevant as biology transitions into a data-driven era where the meaningful interpretation of large data sets is a limiting factor.

The influences of environmental change and development on leaf shape in Vitis

PREMISE

The size and shape (physiognomy) of woody, dicotyledonous angiosperm leaves are correlated with climate. These relationships are the basis for multiple paleoclimate proxies. Here we test whether Vitis exhibits phenotypic plasticity and whether physiognomy varies along the vine.

METHODS

We used Digital Leaf Physiognomy (DiLP) to measure leaf characters of four Vitis species from the USDA Germplasm Repository (Geneva, New York) from the 2012-2013 and 2014-2015 leaf-growing seasons, which had different environmental conditions.

RESULTS

Leaf shape changed allometrically through developmental stages; early stages were more linear than later stages. There were significant differences in physiognomy in the same developmental stage between the growing seasons, and species had significant differences in mean physiognomy between growing seasons. Phenotypic plasticity was defined as changes between growing seasons after controlling for developmental stage or after averaging all developmental stages. Vitis amurensis and V. riparia had the greatest phenotypic plasticity. North American species exhibited significant differences in tooth area:blade area. Intermediate developmental stages were most likely to exhibit phenotypic plasticity, and only V. amurensis exhibited phenotypic plasticity in later developmental stages.

CONCLUSIONS

Leaves have variable phenotypic plasticity along the vine. Environmental signal was strongest in intermediate developmental stages. This is significant for leaf physiognomic-paleoclimate proxies because these leaves are likely the most common in leaf litter and reflect leaves primarily included in paleoclimate reconstructions. Early season and early developmental stages have the potential to be confounding factors but are unlikely to exert significant influence because of differential preservation potential.

Characterizing 3D inflorescence architecture in grapevine using X-ray imaging and advanced morphometrics: implications for understanding cluster density

Inflorescence architecture provides the scaffold on which flowers and fruits develop, and consequently is a primary trait under investigation in many crop systems. Yet the challenge remains to analyse these complex 3D branching structures with appropriate tools. High information content datasets are required to represent the actual structure and facilitate full analysis of both the geometric and the topological features relevant to phenotypic variation in order to clarify evolutionary and developmental inflorescence patterns. We combined advanced imaging (X-ray tomography) and computational approaches (topological and geometric data analysis and structural simulations) to comprehensively characterize grapevine inflorescence architecture (the rachis and all branches without berries) among 10 wild Vitis species. Clustering and correlation analyses revealed unexpected relationships, for example pedicel branch angles were largely independent of other traits. We identified multivariate traits that typified species, which allowed us to classify species with 78.3% accuracy, versus 10% by chance. Twelve traits had strong signals across phylogenetic clades, providing insight into the evolution of inflorescence architecture. We provide an advanced framework to quantify 3D inflorescence and other branched plant structures that can be used to tease apart subtle, heritable features for a better understanding of genetic and environmental effects on plant phenotypes.

Quantitative Trait Locus Analysis of Leaf Morphology Indicates Conserved Shape Loci in Grapevine

Leaf shape in plants plays important roles in water use, canopy structure, and physiological tolerances to abiotic stresses; all important traits for the future development and sustainability of grapevine cultivation. Historically, researchers have used ampelography, the study of leaf shape in grapevines, to differentiate Vitis species and cultivars based on finite leaf attributes. However, ampelographic measurements have limitations and new methods for quantifying shape are now available. We paired an analysis of finite trait attributes with a 17-point landmark survey and generalized Procrustes analysis (GPA) to reconstruct grapevine leaves digitally from five interspecific hybrid mapping families. Using the reconstructed leaves, we performed three types of quantitative trait loci (QTL) analyses to determine the genetic architecture that defines leaf shape. In the first analysis, we compared several important ampelographic measurements as finite trait QTL. In the second and third analyses, we identified significant shape variation via principal components analysis (PCA) and using a multivariate least squares interval mapping (MLSIM) approach. In total, we identified 271 significant QTL across the three measures of leaf shape and identified specific QTL hotspots in the grape genome which appear to drive major aspects of grapevine leaf shape.

Isolating phyllotactic patterns embedded in the secondary growth of sweet cherry (Prunus avium L.) using magnetic resonance imaging

BACKGROUND

Epicormic branches arise from dormant buds patterned during the growth of previous years. Dormant epicormic buds remain just below the surface of trees, pushed outward from the pith during secondary growth, but maintain vascular connections. Epicormic buds can be activated to elongate into a new shoot, either through natural processes or horticultural intervention, to potentially rejuvenate orchards and restructure tree architecture. Because epicormic structures are embedded within secondary growth, tomographic approaches are a useful method to study them and understand their development.

RESULTS

We apply techniques from image processing to determine the locations of epicormic vascular traces embedded within secondary growth of sweet cherry (Prunus avium L.), revealing the juvenile phyllotactic pattern in the trunk of an adult tree. Techniques include the flood fill algorithm to find the pith of the tree, edge detection to approximate the radius, and a conversion to polar coordinates to threshold and segment phyllotactic features. Intensity values from magnetic resonance imaging (MRI) of the trunk are projected onto the surface of a perfect cylinder to find the locations of traces in the "boundary image". Mathematical phyllotaxy provides a means to capture the patterns in the boundary image by modeling phyllotactic parameters. Our cherry tree specimen has the conspicuous parastichy pair (2,3), phyllotactic fraction 2/5, and divergence angle of approximately 143°.

CONCLUSIONS

The methods described provide a framework not only for studying phyllotaxy, but also for processing of volumetric image data in plants. Our results have practical implications for orchard rejuvenation and directed approaches to influence tree architecture. The study of epicormic structures, which are hidden within secondary growth, using tomographic methods also opens the possibility of studying genetic and environmental influences such structures.

Rootstock effects on scion phenotypes in a 'Chambourcin' experimental vineyard

Understanding how root systems modulate shoot system phenotypes is a fundamental question in plant biology and will be useful in developing resilient agricultural crops. Grafting is a common horticultural practice that joins the roots (rootstock) of one plant to the shoot (scion) of another, providing an excellent method for investigating how these two organ systems affect each other. In this study, we used the French-American hybrid grapevine 'Chambourcin' (Vitis L.) as a model to explore the rootstock-scion relationship. We examined leaf shape, ion concentrations, and gene expression in 'Chambourcin' grown ungrafted as well as grafted to three different rootstocks ('SO4', '1103P' and '3309C') across 2 years and three different irrigation treatments. We found that a significant amount of the variation in leaf shape could be explained by the interaction between rootstock and irrigation. For ion concentrations, the primary source of variation identified was the position of a leaf in a shoot, although rootstock and rootstock by irrigation interaction also explained a significant amount of variation for most ions. Lastly, we found rootstock-specific patterns of gene expression in grafted plants when compared to ungrafted vines. Thus, our work reveals the subtle and complex effect of grafting on 'Chambourcin' leaf morphology, ionomics, and gene expression.

Bundle sheath extensions affect leaf structural and physiological plasticity in response to irradiance

Coordination between structural and physiological traits is key to plants' responses to environmental fluctuations. In heterobaric leaves, bundle sheath extensions (BSEs) increase photosynthetic performance (light-saturated rates of photosynthesis, A_max ) and water transport capacity (leaf hydraulic conductance, K_leaf ). However, it is not clear how BSEs affect these and other leaf developmental and physiological parameters in response to environmental conditions. The obscuravenosa (obv) mutation, found in many commercial tomato varieties, leads to absence of BSEs. We examined structural and physiological traits of tomato heterobaric and homobaric (obv) near-isogenic lines grown at two different irradiance levels. K_leaf , minor vein density, and stomatal pore area index decreased with shading in heterobaric but not in homobaric leaves, which show similarly lower values in both conditions. Homobaric plants, on the other hand, showed increased A_max , leaf intercellular air spaces, and mesophyll surface area exposed to intercellular airspace (S_mes ) in comparison with heterobaric plants when both were grown in the shade. BSEs further affected carbon isotope discrimination, a proxy for long-term water-use efficiency. BSEs confer plasticity in traits related to leaf structure and function in response to irradiance levels and might act as a hub integrating leaf structure, photosynthetic function, and water supply and demand.

Long-Distance and Trans-Generational Stomatal Patterning by CO2 Across Arabidopsis Organs

Stomata control water loss and carbon dioxide uptake by both altering pore aperture and developmental patterning. Stomatal patterning is regulated by environmental factors including atmospheric carbon dioxide (p[CO2]), which is increasing globally at an unprecedented rate. Mature leaves are known to convey developmental cues to immature leaves in response to p[CO2], but the developmental mechanisms are unknown. To characterize changes in stomatal patterning resulting from signals moving from mature to developing leaves, we constructed a dual-chamber growth system in which rosette and cauline leaves of Arabidopsis thaliana were subjected to differing p[CO2]. Young rosette tissue was found to adjust stomatal index (SI, the proportion of stomata to total cell number) in response to both the current environment and the environment experienced by mature rosette tissue, whereas cauline leaves appear to be insensitive to p[CO2] treatment. It is likely that cauline leaves and cotyledons deploy mechanisms for controlling stomatal development that share common but also deploy distinctive mechanisms to that operating in rosette leaves. The effect of p[CO2] on stomatal development is retained in cotyledons of the next generation, however, this effect does not occur in pre-germination stomatal lineage cells but only after germination. Finally, these data suggest that p[CO2] affects regulation of stomatal development specifically through the development of satellite stomata (stomata induced by signals from a neighboring stomate) during spacing divisions and not the basal pathway. To our knowledge, this is the first report identifying developmental steps responsible for altered stomatal patterning to p[CO2] and its trans-generational inheritance.

The Persistent Homology Mathematical Framework Provides Enhanced Genotype-to-Phenotype Associations for Plant Morphology

Efforts to understand the genetic and environmental conditioning of plant morphology are hindered by the lack of flexible and effective tools for quantifying morphology. Here, we demonstrate that persistent-homology-based topological methods can improve measurement of variation in leaf shape, serrations, and root architecture. We apply these methods to 2D images of leaves and root systems in field-grown plants of a domesticated introgression line population of tomato (Solanum pennellii). We find that compared with some commonly used conventional traits, (1) persistent-homology-based methods can more comprehensively capture morphological variation; (2) these techniques discriminate between genotypes with a larger normalized effect size and detect a greater number of unique quantitative trait loci (QTLs); (3) multivariate traits, whether statistically derived from univariate or persistent-homology-based traits, improve our ability to understand the genetic basis of phenotype; and (4) persistent-homology-based techniques detect unique QTLs compared to conventional traits or their multivariate derivatives, indicating that previously unmeasured aspects of morphology are now detectable. The QTL results further imply that genetic contributions to morphology can affect both the shoot and root, revealing a pleiotropic basis to natural variation in tomato. Persistent homology is a versatile framework to quantify plant morphology and developmental processes that complements and extends existing methods.

Erratum to: Morphometric analysis of Passiflora leaves: the relationship between landmarks of the vasculature and elliptical Fourier descriptors of the blade

Background

Leaf shape among Passiflora species is spectacularly diverse. Underlying this diversity in leaf shape are profound changes in the patterning of the primary vasculature and laminar outgrowth. Each of these aspects of leaf morphology—vasculature and blade—provides different insights into leaf patterning.

Results

Here, we morphometrically analyze >3300 leaves from 40 different Passiflora species collected sequentially across the vine. Each leaf is measured in two different ways: using 1) 15 homologous Procrustes-adjusted landmarks of the vasculature, sinuses, and lobes; and 2) Elliptical Fourier Descriptors (EFDs), which quantify the outline of the leaf. The ability of landmarks, EFDs, and both datasets together are compared to determine their relative ability to predict species and node position within the vine. Pairwise correlation of x and y landmark coordinates and EFD harmonic coefficients reveals close associations between traits and insights into the relationship between vasculature and blade patterning.

Conclusions

Landmarks, more reflective of the vasculature, and EFDs, more reflective of the blade contour, describe both similar and distinct features of leaf morphology. Landmarks and EFDs vary in ability to predict species identity and node position in the vine and exhibit a correlational structure (both within landmark or EFD traits and between the two data types) revealing constraints between vascular and blade patterning underlying natural variation in leaf morphology among Passiflora species.

Topological Data Analysis as a Morphometric Method: Using Persistent Homology to Demarcate a Leaf Morphospace

Current morphometric methods that comprehensively measure shape cannot compare the disparate leaf shapes found in seed plants and are sensitive to processing artifacts. We explore the use of persistent homology, a topological method applied as a filtration across simplicial complexes (or more simply, a method to measure topological features of spaces across different spatial resolutions), to overcome these limitations. The described method isolates subsets of shape features and measures the spatial relationship of neighboring pixel densities in a shape. We apply the method to the analysis of 182,707 leaves, both published and unpublished, representing 141 plant families collected from 75 sites throughout the world. By measuring leaves from throughout the seed plants using persistent homology, a defined morphospace comparing all leaves is demarcated. Clear differences in shape between major phylogenetic groups are detected and estimates of leaf shape diversity within plant families are made. The approach predicts plant family above chance. The application of a persistent homology method, using topological features, to measure leaf shape allows for a unified morphometric framework to measure plant form, including shapes, textures, patterns, and branching architectures.

Genetic and Developmental Basis for Increased Leaf Thickness in the Arabidopsis Cvi Ecotype

Leaf thickness is a quantitative trait that is associated with the ability of plants to occupy dry, high irradiance environments. Despite its importance, leaf thickness has been difficult to measure reproducibly, which has impeded progress in understanding its genetic basis, and the associated anatomical mechanisms that pattern it. Here, we used a custom-built dual confocal profilometer device to measure leaf thickness in the Arabidopsis Ler × Cvi recombinant inbred line population and found statistical support for four quantitative trait loci (QTL) associated with this trait. We used publically available data for a suite of traits relating to flowering time and growth responses to light quality and show that three of the four leaf thickness QTL coincide with QTL for at least one of these traits. Using time course photography, we quantified the relative growth rate and the pace of rosette leaf initiation in the Ler and Cvi ecotypes. We found that Cvi rosettes grow slower than Ler, both in terms of the rate of leaf initiation and the overall rate of biomass accumulation. Collectively, these data suggest that leaf thickness is tightly linked with physiological status and may present a tradeoff between the ability to withstand stress and rapid vegetative growth. To understand the anatomical basis of leaf thickness, we compared cross-sections of Cvi and Ler leaves and show that Cvi palisade mesophyll cells elongate anisotropically contributing to leaf thickness. Flow cytometry of whole leaves show that endopolyploidy accompanies thicker leaves in Cvi. Overall, our data suggest that mechanistically, an altered schedule of cellular events affecting endopolyploidy and increasing palisade mesophyll cell length contribute to increase of leaf thickness in Cvi. Ultimately, knowledge of the genetic basis and developmental trajectory leaf thickness will inform the mechanisms by which natural selection acts to produce variation in this adaptive trait.

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'.

Divergent leaf shapes among Passiflora species arise from a shared juvenile morphology

Not only does leaf shape vary between Passiflora species, but between sequential nodes of the vine. The profound changes in leaf shape within Passiflora vines reflect the temporal development of the shoot apical meristem from which leaves are derived and patterned, a phenomenon known as heteroblasty. We perform a morphometric analysis of more than 3,300 leaves from 40 different Passiflora species using two different methods: homologous landmarks and Elliptical Fourier Descriptors (EFDs). Changes in leaf shape across the vine are first quantified in allometric terms; that is, changes in the relative area of leaf subregions expressed in terms of overall leaf area. Shape is constrained to strict linear relationships as a function of size that vary between species. Statistical analysis of leaf shape, using landmarks and EFDs, reveals that species effects (regardless of node) are the strongest, followed by interaction effects between species and heteroblasty (i.e., species-specific patterns in leaf shape across nodes) and that heteroblasty effects across nodes (regardless of species) are negligible. The ability of different nodes to predictively discriminate species and the variability of landmark and EFD traits at each node is then analyzed. Heteroblastic trajectories, the changes in leaf shape between the first and last measured leaves in a vine, are then compared between species in a multivariate space. Leaf shape diversity among Passiflora species is expressed in a heteroblastic-dependent manner, unique to each species. Leaf shape is constrained by linear, allometric relationships related to leaf size that vary between species. There is a strong species × heteroblasty interaction effect for leaf shape, suggesting that different leaf shapes between species arise through changes in shape across nodes specific to each species. The first leaves in the series are not only more like each other, but are also less variable across species. From this similar, shared leaf shape, subsequent leaves in the heteroblastic series follow divergent morphological trajectories. The disparate leaf shapes characteristic of Passiflora species arise from a shared, juvenile morphology.

Genetic Architecture and Molecular Networks Underlying Leaf Thickness in Desert-Adapted Tomato Solanum pennellii

Thicker leaves allow plants to grow in water-limited conditions. However, our understanding of the genetic underpinnings of this highly functional leaf shape trait is poor. We used a custom-built confocal profilometer to directly measure leaf thickness in a set of introgression lines (ILs) derived from the desert tomato Solanum pennellii and identified quantitative trait loci. We report evidence of a complex genetic architecture of this trait and roles for both genetic and environmental factors. Several ILs with thick leaves have dramatically elongated palisade mesophyll cells and, in some cases, increased leaf ploidy. We characterized the thick IL2-5 and IL4-3 in detail and found increased mesophyll cell size and leaf ploidy levels, suggesting that endoreduplication underpins leaf thickness in tomato. Next, we queried the transcriptomes and inferred dynamic Bayesian networks of gene expression across early leaf ontogeny in these lines to compare the molecular networks that pattern leaf thickness. We show that thick ILs share S. pennellii-like expression profiles for putative regulators of cell shape and meristem determinacy as well as a general signature of cell cycle-related gene expression. However, our network data suggest that leaf thickness in these two lines is patterned at least partially by distinct mechanisms. Consistent with this hypothesis, double homozygote lines combining introgression segments from these two ILs show additive phenotypes, including thick leaves, higher ploidy levels, and larger palisade mesophyll cells. Collectively, these data establish a framework of genetic, anatomical, and molecular mechanisms that pattern leaf thickness in desert-adapted tomato.

Morphometric analysis of Passiflora leaves: the relationship between landmarks of the vasculature and elliptical Fourier descriptors of the blade

Background

Leaf shape among Passiflora species is spectacularly diverse. Underlying this diversity in leaf shape are profound changes in the patterning of the primary vasculature and laminar outgrowth. Each of these aspects of leaf morphology-vasculature and blade-provides different insights into leaf patterning.

Results

Here, we morphometrically analyze >3300 leaves from 40 different Passiflora species collected sequentially across the vine. Each leaf is measured in two different ways: using 1) 15 homologous Procrustes-adjusted landmarks of the vasculature, sinuses, and lobes; and 2) Elliptical Fourier Descriptors (EFDs), which quantify the outline of the leaf. The ability of landmarks, EFDs, and both datasets together are compared to determine their relative ability to predict species and node position within the vine. Pairwise correlation of x and y landmark coordinates and EFD harmonic coefficients reveals close associations between traits and insights into the relationship between vasculature and blade patterning.

Conclusions

Landmarks, more reflective of the vasculature, and EFDs, more reflective of the blade contour, describe both similar and distinct features of leaf morphology. Landmarks and EFDs vary in ability to predict species identity and node position in the vine and exhibit a correlational structure (both within landmark or EFD traits and between the two data types) revealing constraints between vascular and blade patterning underlying natural variation in leaf morphology among Passiflora species.

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.

Gene expression atlas for the food security crop cassava

Cassava (Manihot esculenta) feeds c. 800 million people world-wide. Although this crop displays high productivity under drought and poor soil conditions, it is susceptible to disease, postharvest deterioration and the roots contain low nutritional content. Here, we provide molecular identities for 11 cassava tissue/organ types through RNA-sequencing and develop an open access, web-based interface for further interrogation of the data. Through this dataset, we consider the physiology of cassava. Specifically, we focus on identification of the transcriptional signatures that define the massive, underground storage roots used as a food source and the favored target tissue for transgene integration and genome editing, friable embryogenic callus (FEC). Further, we identify promoters able to drive strong expression in multiple tissue/organs. The information gained from this study is of value for both conventional and biotechnological improvement programs.

Modifications to a LATE MERISTEM IDENTITY1 gene are responsible for the major leaf shapes of Upland cotton (Gossypium hirsutum L.)

Leaf shape varies spectacularly among plants. Leaves are the primary source of photoassimilate in crop plants, and understanding the genetic basis of variation in leaf morphology is critical to improving agricultural productivity. Leaf shape played a unique role in cotton improvement, as breeders have selected for entire and lobed leaf morphs resulting from a single locus, okra (l-D1), which is responsible for the major leaf shapes in cotton. The l-D1 locus is not only of agricultural importance in cotton, but through pioneering chimeric and morphometric studies, it has contributed to fundamental knowledge about leaf development. Here we show that an HD-Zip transcription factor homologous to the LATE MERISTEM IDENTITY1 (LMI1) gene of Arabidopsis is the causal gene underlying the l-D1 locus. The classical okra leaf shape allele has a 133-bp tandem duplication in the promoter, correlated with elevated expression, whereas an 8-bp deletion in the third exon of the presumed wild-type normal allele causes a frame-shifted and truncated coding sequence. Our results indicate that subokra is the ancestral leaf shape of tetraploid cotton that gave rise to the okra allele and that normal is a derived mutant allele that came to predominate and define the leaf shape of cultivated cotton. Virus-induced gene silencing (VIGS) of the LMI1-like gene in an okra variety was sufficient to induce normal leaf formation. The developmental changes in leaves conferred by this gene are associated with a photosynthetic transcriptomic signature, substantiating its use by breeders to produce a superior cotton ideotype.

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.

Morphometrics Reveals Complex and Heritable Apple Leaf Shapes

Apple (Malus spp.) is a widely grown and valuable fruit crop. Leaf shape is important for flowering in apple and may also be an early indicator for other agriculturally valuable traits. We examined 9,000 leaves from 869 unique apple accessions using linear measurements and comprehensive morphometric techniques. We identified allometric variation as the result of differing length-to-width aspect ratios between accessions and species of apple. The allometric variation was due to variation in the width of the leaf blade, not the length. Aspect ratio was highly correlated with the first principal component (PC1) of morphometric variation quantified using elliptical Fourier descriptors (EFDs) and persistent homology (PH). While the primary source of variation was aspect ratio, subsequent PCs corresponded to complex shape variation not captured by linear measurements. After linking the morphometric information with over 122,000 genome-wide single nucleotide polymorphisms (SNPs), we found high SNP heritability values even at later PCs, indicating that comprehensive morphometrics can capture complex, heritable phenotypes. Thus, techniques such as EFDs and PH are capturing heritable biological variation that would be missed using linear measurements alone.

Morphometric analysis of Passiflora leaves: the relationship between landmarks of the vasculature and elliptical Fourier descriptors of the blade

BACKGROUND

Leaf shape among Passiflora species is spectacularly diverse. Underlying this diversity in leaf shape are profound changes in the patterning of the primary vasculature and laminar outgrowth. Each of these aspects of leaf morphology-vasculature and blade-provides different insights into leaf patterning.

RESULTS

Here, we morphometrically analyze >3300 leaves from 40 different Passiflora species collected sequentially across the vine. Each leaf is measured in two different ways: using 1) 15 homologous Procrustes-adjusted landmarks of the vasculature, sinuses, and lobes; and 2) Elliptical Fourier Descriptors (EFDs), which quantify the outline of the leaf. The ability of landmarks, EFDs, and both datasets together are compared to determine their relative ability to predict species and node position within the vine. Pairwise correlation of x and y landmark coordinates and EFD harmonic coefficients reveals close associations between traits and insights into the relationship between vasculature and blade patterning.

CONCLUSIONS

Landmarks, more reflective of the vasculature, and EFDs, more reflective of the blade contour, describe both similar and distinct features of leaf morphology. Landmarks and EFDs vary in ability to predict species identity and node position in the vine and exhibit a correlational structure (both within landmark or EFD traits and between the two data types) revealing constraints between vascular and blade patterning underlying natural variation in leaf morphology among Passiflora species.

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'.

Evolutionary and Environmental Forces Sculpting Leaf Development

Leaf shape is spectacularly diverse. As a major component of plant architecture and an interface for light capture, gas exchange, and thermoregulation, the potential contributions of leaves to plant fitness are innumerable. Particularly because of their intimate association and interaction with the surrounding environment, both the plasticity of leaf shape during the lifetime of a plant and the evolution of leaf shape over geologic time are revealing with respect to leaf function. Leaf shapes arise within a developmental context that constrains both their evolution and environmental plasticity. Quantitative models capturing genetic diversity, developmental context, and environmental plasticity will be required to fully understand the evolution and development of leaf shape and its response to environmental pressures. In this review, we discuss recent literature demonstrating that distinct molecular pathways are modulated by specific environmental inputs, the output of which regulates leaf dissection. We propose a synthesis explaining both historical patterns in the paleorecord and conserved plastic responses in extant plants. Understanding the potential adaptive value of leaf shape, and how to molecularly manipulate it, will prove to be invaluable in designing crops optimized for future climates.

A New Advanced Backcross Tomato Population Enables High Resolution Leaf QTL Mapping and Gene Identification

Quantitative Trait Loci (QTL) mapping is a powerful technique for dissecting the genetic basis of traits and species differences. Established tomato mapping populations between domesticated tomato (Solanum lycopersicum) and its more distant interfertile relatives typically follow a near isogenic line (NIL) design, such as the S. pennellii Introgression Line (IL) population, with a single wild introgression per line in an otherwise domesticated genetic background. Here, we report on a new advanced backcross QTL mapping resource for tomato, derived from a cross between the M82 tomato cultivar and S. pennellii This so-called Backcrossed Inbred Line (BIL) population is comprised of a mix of BC2 and BC3 lines, with domesticated tomato as the recurrent parent. The BIL population is complementary to the existing S. pennellii IL population, with which it shares parents. Using the BILs, we mapped traits for leaf complexity, leaflet shape, and flowering time. We demonstrate the utility of the BILs for fine-mapping QTL, particularly QTL initially mapped in the ILs, by fine-mapping several QTL to single or few candidate genes. Moreover, we confirm the value of a backcrossed population with multiple introgressions per line, such as the BILs, for epistatic QTL mapping. Our work was further enabled by the development of our own statistical inference and visualization tools, namely a heterogeneous hidden Markov model for genotyping the lines, and by using state-of-the-art sparse regression techniques for QTL mapping.

eQTL Regulating Transcript Levels Associated with Diverse Biological Processes in Tomato

Variation in gene expression, in addition to sequence polymorphisms, is known to influence developmental, physiological, and metabolic traits in plants. Genetic mapping populations have facilitated identification of expression quantitative trait loci (eQTL), the genetic determinants of variation in gene expression patterns. We used an introgression population developed from the wild desert-adapted Solanum pennellii and domesticated tomato (Solanum lycopersicum) to identify the genetic basis of transcript level variation. We established the effect of each introgression on the transcriptome and identified approximately 7,200 eQTL regulating the steady-state transcript levels of 5,300 genes. Barnes-Hut t-distributed stochastic neighbor embedding clustering identified 42 modules revealing novel associations between transcript level patterns and biological processes. The results showed a complex genetic architecture of global transcript abundance pattern in tomato. Several genetic hot spots regulating a large number of transcript level patterns relating to diverse biological processes such as plant defense and photosynthesis were identified. Important eQTL regulating transcript level patterns were related to leaf number and complexity as well as hypocotyl length. Genes associated with leaf development showed an inverse correlation with photosynthetic gene expression, but eQTL regulating genes associated with leaf development and photosynthesis were dispersed across the genome. This comprehensive eQTL analysis details the influence of these loci on plant phenotypes and will be a valuable community resource for investigations on the genetic effects of eQTL on phenotypic traits in tomato.

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.

Plant chimeras: The good, the bad, and the 'Bizzaria'

Chimeras - organisms that are composed of cells of more than one genotype - captured the human imagination long before they were formally described and used in the laboratory. These organisms owe their namesake to a fire-breathing monster from Greek mythology that has the head of a lion, the body of a goat, and the tail of a serpent. The first description of a non-fictional chimera dates back to the middle of the seventeenth century when the Florentine gardener Pietro Nati discovered an adventitious shoot growing from the graft junction between sour orange (Citrus aurantium) and citron (Citrus medica). This perplexing chimera that grows with sectors phenotypically resembling each of the citrus progenitors inspired discussion and wonder from the scientific community and was fittingly named the 'Bizzaria'. Initially, the 'Bizzaria' was believed to be an asexual hybrid that formed from a cellular fusion between the grafted parents; however, in-depth cellular analyses carried out centuries later demonstrated that the 'Bizzaria', along with other chimeras, owe their unique sectored appearance to a conglomeration of cells from the two donors. Since this pivotal discovery at the turn of the twentieth century, chimeras have served both as tools and as unique biological phenomena that have contributed to our understanding of plant development at the cellular, tissue, and organismal level. Rapid advancements in genome sequencing technologies have enabled the establishment of new model species with novel morphological and developmental features that enable the generation of chimeric organisms. In this review, we show that genetic mosaic and chimera studies provide a technologically simple way to delve into the organismal, genetic, and genomic inner workings underlying the development of diverse model organisms. Moreover, we discuss the unique opportunity that chimeras present to explore universal principles governing intercellular communication and the coordination of organismal biology in a heterogenomic landscape.

Rootstocks: Diversity, Domestication, and Impacts on Shoot Phenotypes

Grafting is an ancient agricultural practice that joins the root system (rootstock) of one plant to the shoot (scion) of another. It is most commonly employed in woody perennial crops to indirectly manipulate scion phenotype. While recent research has focused on scions, here we investigate rootstocks, the lesser-known half of the perennial crop equation. We review natural grafting, grafting in agriculture, rootstock diversity and domestication, and developing areas of rootstock research, including molecular interactions and rootstock microbiomes. With growing interest in perennial crops as valuable components of sustainable agriculture, rootstocks provide one mechanism by which to improve and expand woody perennial cultivation in a range of environmental conditions.

Latent developmental and evolutionary shapes embedded within the grapevine leaf

Across plants, leaves exhibit profound diversity in shape. As a single leaf expands, its shape is in constant flux. Plants may also produce leaves with different shapes at successive nodes. In addition, leaf shape varies among individuals, populations and species as a result of evolutionary processes and environmental influences. Because leaf shape can vary in many different ways, theoretically, the effects of distinct developmental and evolutionary processes are separable, even within the shape of a single leaf. Here, we measured the shapes of > 3200 leaves representing > 270 vines from wild relatives of domesticated grape (Vitis spp.) to determine whether leaf shapes attributable to genetics and development are separable from each other. We isolated latent shapes (multivariate signatures that vary independently from each other) embedded within the overall shape of leaves. These latent shapes can predict developmental stages independent from species identity and vice versa. Shapes predictive of development were then used to stage leaves from 1200 varieties of domesticated grape (Vitis vinifera), revealing that changes in timing underlie leaf shape diversity. Our results indicate that distinct latent shapes combine to produce a composite morphology in leaves, and that developmental and evolutionary contributions to shape vary independently from each other.

Latent developmental and evolutionary shapes embedded within the grapevine leaf

Across plants, leaves exhibit profound diversity in shape. As a single leaf expands, its shape is in constant flux. Plants may also produce leaves with different shapes at successive nodes. In addition, leaf shape varies among individuals, populations and species as a result of evolutionary processes and environmental influences. Because leaf shape can vary in many different ways, theoretically, the effects of distinct developmental and evolutionary processes are separable, even within the shape of a single leaf. Here, we measured the shapes of > 3200 leaves representing > 270 vines from wild relatives of domesticated grape (Vitis spp.) to determine whether leaf shapes attributable to genetics and development are separable from each other. We isolated latent shapes (multivariate signatures that vary independently from each other) embedded within the overall shape of leaves. These latent shapes can predict developmental stages independent from species identity and vice versa. Shapes predictive of development were then used to stage leaves from 1200 varieties of domesticated grape (Vitis vinifera), revealing that changes in timing underlie leaf shape diversity. Our results indicate that distinct latent shapes combine to produce a composite morphology in leaves, and that developmental and evolutionary contributions to shape vary independently from each other.

Transcriptomic analysis suggests a key role for SQUAMOSA PROMOTER BINDING PROTEIN LIKE, NAC and YUCCA genes in the heteroblastic development of the temperate rainforest tree Gevuina avellana (Proteaceae)

Heteroblasty, the temporal development of the meristem, can produce diverse leaf shapes within a plant. Gevuina avellana, a tree from the South American temperate rainforest shows strong heteroblasty affecting leaf shape, transitioning from juvenile simple leaves to highly pinnate adult leaves. Light availability within the forest canopy also modulates its leaf size and complexity. Here we studied how the interaction between the light environment and the heteroblastic progression of leaves is coordinated in this species. We used RNA-seq on the Illumina platform to compare the range of transcriptional responses in leaf primordia of G. avellana at different heteroblastic stages and growing under different light environments. We found a steady up-regulation of SQUAMOSA PROMOTER BINDING PROTEIN LIKE (SPL), NAC, YUCCA and AGAMOUS-LIKE genes associated with increases in age, leaf complexity, and light availability. In contrast, expression of TCP, TPR and KNOTTED1 homeobox genes showed a sustained down-regulation. Additionally, genes involved in auxin synthesis/transport and jasmonate activity were differentially expressed, indicating an active regulation of processes controlled by these hormones. Our large-scale transcriptional analysis of the leaf primordia of G. avellana sheds light on the integration of internal and external cues during heteroblastic development in this species.