Morphometrics and Phylogenomics of Coca (Erythroxylum spp.) Illuminate Its Reticulate Evolution, With Implications for Taxonomy.

South American coca (Erythroxylum coca and E. novogranatense) has been a keystone crop for many Andean and Amazonian communities for at least 8,000 years. However, over the last half-century, global demand for its alkaloid cocaine has driven intensive agriculture of this plant and placed it in the center of armed conflict and deforestation. To monitor the changing landscape of coca plantations, the United Nations Office on Drugs and Crime collects annual data on their areas of cultivation. However, attempts to delineate areas in which different varieties are grown have failed due to limitations around identification. In the absence of flowers, identification relies on leaf morphology, yet the extent to which this is reflected in taxonomy is uncertain. Here, we analyze the consistency of the current naming system of coca and its four closest wild relatives (the "coca clade"), using morphometrics, phylogenomics, molecular clocks, and population genomics. We include name-bearing type specimens of coca's closest wild relatives E. gracilipes and E. cataractarum. Morphometrics of 342 digitized herbarium specimens show that leaf shape and size fail to reliably discriminate between species and varieties. However, the statistical analyses illuminate that rounder and more obovate leaves of certain varieties could be associated with the subtle domestication syndrome of coca. Our phylogenomic data indicate extensive gene flow involving E. gracilipes which, combined with morphometrics, supports E. gracilipes being retained as a single species. Establishing a robust evolutionary-taxonomic framework for the coca clade will facilitate the development of cost-effective genotyping methods to support reliable identification.

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.

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.

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.

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.

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

The genetics and physiology of seed dormancy, a crucial trait in common bean domestication.

BACKGROUND: Physical seed dormancy is an important trait in legume domestication. Although seed dormancy is beneficial in wild ecosystems, it is generally considered to be an undesirable trait in crops due to reduction in yield and / or quality. The physiological mechanism and underlying genetic factor(s) of seed dormancy is largely unknown in several legume species. Here we employed an integrative approach to understand the mechanisms controlling physical seed dormancy in common bean (Phaseolus vulgaris L.). RESULTS: Using an innovative CT scan imaging system, we were able to track water movements inside the seed coat. We found that water uptake initiates from the bean seed lens. Using a scanning electron microscopy (SEM) we further identified several micro-cracks on the lens surface of non-dormant bean genotypes. Bulked segregant analysis (BSA) was conducted on a bi-parental RIL (recombinant inbred line) population, segregating for seed dormancy. This analysis revealed that the seed water uptake is associated with a single major QTL on Pv03. The QTL region was fine-mapped to a 118 Kb interval possessing 11 genes. Coding sequence analysis of candidate genes revealed a 5-bp insertion in an ortholog of pectin acetylesterase 8 that causes a frame shift, loss-of-function mutation in non-dormant genotype. Gene expression analysis of the candidate genes in the seed coat of contrasting genotypes indicated 21-fold lower expression of pectin acetylesterase 8 in non-dormant genotype. An analysis of mutational polymorphism was conducted among wild and domesticated beans. Although all the wild beans possessed the functional allele of pectin acetylesterase 8, the majority (77%) of domesticated beans had the non-functional allele suggesting that this variant was under strong selection pressure through domestication. CONCLUSIONS: In this study, we identified the physiological mechanism of physical seed dormancy and have identified a candidate allele causing variation in this trait. Our findings suggest that a 5-bp insertion in an ortholog of pectin acetylesterase 8 is likely a major causative mutation underlying the loss of seed dormancy during domestication. Although the results of current study provide strong evidences for the role of pectin acetylesterase 8 in seed dormancy, further confirmations seem necessary by employing transgenic approaches.

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.

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.

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.

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.