An adaptive teosinte mexicana introgression modulates phosphatidylcholine levels and is associated with maize flowering time.

Native Americans domesticated maize (Zea mays ssp. mays) from lowland teosinte parviglumis (Zea mays ssp. parviglumis) in the warm Mexican southwest and brought it to the highlands of Mexico and South America where it was exposed to lower temperatures that imposed strong selection on flowering time. Phospholipids are important metabolites in plant responses to low-temperature and phosphorus availability and have been suggested to influence flowering time. Here, we combined linkage mapping with genome scans to identify High PhosphatidylCholine 1 (HPC1), a gene that encodes a phospholipase A1 enzyme, as a major driver of phospholipid variation in highland maize. Common garden experiments demonstrated strong genotype-by-environment interactions associated with variation at HPC1, with the highland HPC1 allele leading to higher fitness in highlands, possibly by hastening flowering. The highland maize HPC1 variant resulted in impaired function of the encoded protein due to a polymorphism in a highly conserved sequence. A meta-analysis across HPC1 orthologs indicated a strong association between the identity of the amino acid at this position and optimal growth in prokaryotes. Mutagenesis of HPC1 via genome editing validated its role in regulating phospholipid metabolism. Finally, we showed that the highland HPC1 allele entered cultivated maize by introgression from the wild highland teosinte Zea mays ssp. mexicana and has been maintained in maize breeding lines from the Northern United States, Canada, and Europe. Thus, HPC1 introgressed from teosinte mexicana underlies a large metabolic QTL that modulates phosphatidylcholine levels and has an adaptive effect at least in part via induction of early flowering time.

Allele-specific Expression Reveals Multiple Paths to Highland Adaptation in Maize.

Maize is a staple food of smallholder farmers living in highland regions up to 4,000 m above sea level worldwide. Mexican and South American highlands are two major highland maize growing regions, and population genetic data suggest the maize's adaptation to these regions occurred largely independently, providing a case study for convergent evolution. To better understand the mechanistic basis of highland adaptation, we crossed maize landraces from 108 highland and lowland sites of Mexico and South America with the inbred line B73 to produce F1 hybrids and grew them in both highland and lowland sites in Mexico. We identified thousands of genes with divergent expression between highland and lowland populations. Hundreds of these genes show patterns of convergent evolution between Mexico and South America. To dissect the genetic architecture of the divergent gene expression, we developed a novel allele-specific expression analysis pipeline to detect genes with divergent functional cis-regulatory variation between highland and lowland populations. We identified hundreds of genes with divergent cis-regulation between highland and lowland landrace alleles, with 20 in common between regions, further suggesting convergence in the genes underlying highland adaptation. Further analyses suggest multiple mechanisms contribute to this convergence in gene regulation. Although the vast majority of evolutionary changes associated with highland adaptation were region specific, our findings highlight an important role for convergence at the gene expression and gene regulation levels as well.

Gene regulatory effects of a large chromosomal inversion in highland maize.

Chromosomal inversions play an important role in local adaptation. Inversions can capture multiple locally adaptive functional variants in a linked block by repressing recombination. However, this recombination suppression makes it difficult to identify the genetic mechanisms underlying an inversion's role in adaptation. In this study, we used large-scale transcriptomic data to dissect the functional importance of a 13 Mb inversion locus (Inv4m) found almost exclusively in highland populations of maize (Zea mays ssp. mays). Inv4m was introgressed into highland maize from the wild relative Zea mays ssp. mexicana, also present in the highlands of Mexico, and is thought to be important for the adaptation of these populations to cultivation in highland environments. However, the specific genetic variants and traits that underlie this adaptation are not known. We created two families segregating for the standard and inverted haplotypes of Inv4m in a common genetic background and measured gene expression effects associated with the inversion across 9 tissues in two experimental conditions. With these data, we quantified both the global transcriptomic effects of the highland Inv4m haplotype, and the local cis-regulatory variation present within the locus. We found diverse physiological effects of Inv4m across the 9 tissues, including a strong effect on the expression of genes involved in photosynthesis and chloroplast physiology. Although we could not confidently identify the causal alleles within Inv4m, this research accelerates progress towards understanding this inversion and will guide future research on these important genomic features.

Molecular Parallelism Underlies Convergent Highland Adaptation of Maize Landraces.

Convergent phenotypic evolution provides some of the strongest evidence for adaptation. However, the extent to which recurrent phenotypic adaptation has arisen via parallelism at the molecular level remains unresolved, as does the evolutionary origin of alleles underlying such adaptation. Here, we investigate genetic mechanisms of convergent highland adaptation in maize landrace populations and evaluate the genetic sources of recurrently selected alleles. Population branch excess statistics reveal substantial evidence of parallel adaptation at the level of individual single-nucleotide polymorphism (SNPs), genes, and pathways in four independent highland maize populations. The majority of convergently selected SNPs originated via migration from a single population, most likely in the Mesoamerican highlands, while standing variation introduced by ancient gene flow was also a contributor. Polygenic adaptation analyses of quantitative traits reveal that alleles affecting flowering time are significantly associated with elevation, indicating the flowering time pathway was targeted by highland adaptation. In addition, repeatedly selected genes were significantly enriched in the flowering time pathway, indicating their significance in adapting to highland conditions. Overall, our study system represents a promising model to study convergent evolution in plants with potential applications to crop adaptation across environmental gradients.

ROOT PENETRATION INDEX 3, a major quantitative trait locus associated with root system penetrability in Arabidopsis.

Soil mechanical impedance precludes root penetration, confining root system development to shallow soil horizons where mobile nutrients are scarce. Using a two-phase-agar system, we characterized Arabidopsis responses to low and high mechanical impedance at three root penetration stages. We found that seedlings whose roots fail to penetrate agar barriers show a significant reduction in leaf area, root length, and elongation zone and an increment in root diameter, while those capable of penetrating show only minor morphological effects. Analyses using different auxin-responsive reporter lines, exogenous auxins, and inhibitor treatments suggest that auxin responsiveness and PIN-mediated auxin distribution play an important role in regulating root responses to mechanical impedance. The assessment of 21 Arabidopsis accessions revealed that primary root penetrability varies widely among accessions. To search for quantitative trait loci (QTLs) associated to root system penetrability, we evaluated a recombinant inbred population derived from Landsberg erecta (Ler-0, with a high primary root penetrability) and Shahdara (Sha, with a low primary root penetrability) accessions. QTL analysis revealed a major-effect QTL localized in chromosome 3, ROOT PENETRATION INDEX 3 (q-RPI3), which accounted for 29.98% (logarithm of odds=8.82) of the total phenotypic variation. Employing an introgression line (IL-321) with a homozygous q-RPI3 region from Sha in the Ler-0 genetic background, we demonstrated that q-RPI3 plays a crucial role in root penetrability. This multiscale study reveals new insights into root plasticity during the penetration process in hard agar layers, natural variation, and genetic architecture behind primary root penetrability in Arabidopsis.

Malate-dependent Fe accumulation is a critical checkpoint in the root developmental response to low phosphate.

Low phosphate (Pi) availability constrains plant development and seed production in both natural and agricultural ecosystems. When Pi is scarce, modifications of root system architecture (RSA) enhance the soil exploration ability of the plant and lead to an increase in Pi uptake. In Arabidopsis, an iron-dependent mechanism reprograms primary root growth in response to low Pi availability. This program is activated upon contact of the root tip with low-Pi media and induces premature cell differentiation and the arrest of mitotic activity in the root apical meristem, resulting in a short-root phenotype. However, the mechanisms that regulate the primary root response to Pi-limiting conditions remain largely unknown. Here we report on the isolation and characterization of two low-Pi insensitive mutants (lpi5 and lpi6), which have a long-root phenotype when grown in low-Pi media. Cellular, genomic, and transcriptomic analysis of low-Pi insensitive mutants revealed that the genes previously shown to underlie Arabidopsis Al tolerance via root malate exudation, known as SENSITIVE TO PROTON RHIZOTOXICITY (STOP1) and ALUMINUM ACTIVATED MALATE TRANSPORTER 1 (ALMT1), represent a critical checkpoint in the root developmental response to Pi starvation in Arabidopsis thaliana Our results also show that exogenous malate can rescue the long-root phenotype of lpi5 and lpi6 Malate exudation is required for the accumulation of Fe in the apoplast of meristematic cells, triggering the differentiation of meristematic cells in response to Pi deprivation.

Grasses suppress shoot-borne roots to conserve water during drought.

Many important crops are members of the Poaceae family, which develop root systems characterized by a high degree of root initiation from the belowground basal nodes of the shoot, termed the crown. Although this postembryonic shoot-borne root system represents the major conduit for water uptake, little is known about the effect of water availability on its development. Here we demonstrate that in the model C4 grass Setaria viridis, the crown locally senses water availability and suppresses postemergence crown root growth under a water deficit. This response was observed in field and growth room environments and in all grass species tested. Luminescence-based imaging of root systems grown in soil-like media revealed a shift in root growth from crown-derived to primary root-derived branches, suggesting that primary root-dominated architecture can be induced in S. viridis under certain stress conditions. Crown roots of Zea mays and Setaria italica, domesticated relatives of teosinte and S. viridis, respectively, show reduced sensitivity to water deficit, suggesting that this response might have been influenced by human selection. Enhanced water status of maize mutants lacking crown roots suggests that under a water deficit, stronger suppression of crown roots actually may benefit crop productivity.

Co-ordinated Changes in the Accumulation of Metal Ions in Maize (Zea mays ssp. mays L.) in Response to Inoculation with the Arbuscular Mycorrhizal Fungus Funneliformis mosseae.

Arbuscular mycorrhizal symbiosis is an ancient interaction between plants and fungi of the phylum Glomeromycota. In exchange for photosynthetically fixed carbon, the fungus provides the plant host with greater access to soil nutrients via an extensive network of root-external hyphae. Here, to determine the impact of the symbiosis on the host ionome, the concentration of 19 elements was determined in the roots and leaves of a panel of 30 maize varieties, grown under phosphorus-limiting conditions, with or without inoculation with the fungus Funneliformis mosseae. Although the most recognized benefit of the symbiosis to the host plant is greater access to soil phosphorus, the concentration of a number of other elements responded significantly to inoculation across the panel as a whole. In addition, variety-specific effects indicated the importance of plant genotype to the response. Clusters of elements were identified that varied in a co-ordinated manner across genotypes, and that were maintained between non-inoculated and inoculated plants.

Nicotianamine functions in the Phloem-based transport of iron to sink organs, in pollen development and pollen tube growth in Arabidopsis.

The metal chelator nicotianamine promotes the bioavailability of Fe and reduces cellular Fe toxicity. For breeding Fe-efficient crops, we need to explore the fundamental impact of nicotianamine on plant development and physiology. The quadruple nas4x-2 mutant of Arabidopsis thaliana cannot synthesize any nicotianamine, shows strong leaf chlorosis, and is sterile. To date, these phenotypes have not been fully explained. Here, we show that sink organs of this mutant were Fe deficient, while aged leaves were Fe sufficient. Upper organs were also Zn deficient. We demonstrate that transport of Fe to aged leaves relied on citrate, which partially complemented the loss of nicotianamine. In the absence of nicotianamine, Fe accumulated in the phloem. Our results show that rather than enabling the long-distance movement of Fe in the phloem (as is the case for Zn), nicotianamine facilitates the transport of Fe from the phloem to sink organs. We delimit nicotianamine function in plant reproductive biology and demonstrate that nicotianamine acts in pollen development in anthers and pollen tube passage in the carpels. Since Fe and Zn both enhance pollen germination, a lack of either metal may contribute to the reproductive defect. Our study sheds light on the physiological functions of nicotianamine.

Iron-dependent modifications of the flower transcriptome, proteome, metabolome, and hormonal content in an Arabidopsis ferritin mutant.

Iron homeostasis is an important process for flower development and plant fertility. The role of plastids in these processes has been shown to be essential. To document the relationships between plastid iron homeostasis and flower biology further, a global study (transcriptome, proteome, metabolome, and hormone analysis) was performed of Arabidopsis flowers from wild-type and triple atfer1-3-4 ferritin mutant plants grown under iron-sufficient or excess conditions. Some major modifications in specific functional categories were consistently observed at these three omic levels, although no significant overlaps of specific transcripts and proteins were detected. These modifications concerned redox reactions and oxidative stress, as well as amino acid and protein catabolism, this latter point being exemplified by an almost 10-fold increase in urea concentration of atfer1-3-4 flowers from plants grown under iron excess conditions. The mutant background caused alterations in Fe-haem redox proteins located in membranes and in hormone-responsive proteins. Specific effects of excess Fe in the mutant included further changes in these categories, supporting the idea that the mutant is facing a more intense Fe/redox stress than the wild type. The mutation and/or excess Fe had a strong impact at the membrane level, as denoted by the changes in the transporter and lipid metabolism categories. In spite of the large number of genes and proteins responsive to hormones found to be regulated in this study, changes in the hormonal balance were restricted to cytokinins, especially in the mutant plants grown under Fe excess conditions.

Genetic mapping of maize metabolites using high-throughput mass profiling.

Plant metabolites are the basis of human nutrition and have biological relevance in ecology. Farmers selected plants with favorable characteristics since prehistoric times and improved the cultivars, but without knowledge of underlying mechanisms. Understanding the genetic basis of metabolite production can facilitate the successful breeding of plants with augmented nutritional value. To identify genetic factors related to the metabolic composition in maize, we generated mass profiles of 198 recombinant inbred lines (RILs) and their parents (B73 and Mo17) using direct-injection electrospray ionization mass spectrometry (DLI-ESI MS). Mass profiling allowed the correct clustering of samples according to genotype. We quantified 71 mass features from grains and 236 mass features from leaf extracts. For the corresponding ions, we identified tissue-specific metabolic 'Quantitative Trait Loci' (mQTLs) distributed across the maize genome. These genetic regions could regulate multiple metabolite biosynthesis pathways. Our findings demonstrate that DLI-ESI MS has sufficient analytical resolution to map mQTLs. These identified genetic loci will be helpful in metabolite-focused maize breeding. Mass profiling is a powerful tool for detecting mQTLs in maize and enables the high-throughput screening of loci responsible for metabolite biosynthesis.

Oxidative stress is a consequence, not a cause, of aluminum toxicity in the forage legume Lotus corniculatus.

• Aluminum (Al) toxicity is a major limiting factor of crop production on acid soils, but the implication of oxidative stress in this process is controversial. A multidisciplinary approach was used here to address this question in the forage legume Lotus corniculatus. • Plants were treated with low Al concentrations in hydroponic culture, and physiological and biochemical parameters, together with semiquantitative metabolic and proteomic profiles, were determined. • The exposure of plants to 10 µM Al inhibited root and leaf growth, but had no effect on the production of reactive oxygen species or lipid peroxides. By contrast, exposure to 20 µM Al elicited the production of superoxide radicals, peroxide and malondialdehyde. In response to Al, there was a progressive replacement of the superoxide dismutase isoforms in the cytosol, a loss of ascorbate and consistent changes in amino acids, sugars and associated enzymes. • We conclude that oxidative stress is not a causative factor of Al toxicity. The increased contents in roots of two powerful Al chelators, malic and 2-isopropylmalic acids, together with the induction of an Al-activated malate transporter gene, strongly suggest that both organic acids are implicated in Al detoxification. The effects of Al on key proteins involved in cytoskeleton dynamics, protein turnover, transport, methylation reactions, redox control and stress responses underscore a metabolic dysfunction, which affects multiple cellular compartments, particularly in plants exposed to 20 µM Al.

Demonstration of local adaptation in maize landraces by reciprocal transplantation.

Populations are locally adapted when they exhibit higher fitness than foreign populations in their native habitat. Maize landrace adaptations to highland and lowland conditions are of interest to researchers and breeders. To determine the prevalence and strength of local adaptation in maize landraces, we performed a reciprocal transplant experiment across an elevational gradient in Mexico. We grew 120 landraces, grouped into four populations (Mexican Highland, Mexican Lowland, South American Highland, South American Lowland), in Mexican highland and lowland common gardens and collected phenotypes relevant to fitness and known highland-adaptive traits such as anthocyanin pigmentation and macrohair density. 67k DArTseq markers were generated from field specimens to allow comparisons between phenotypic patterns and population genetic structure. We found phenotypic patterns consistent with local adaptation, though these patterns differ between the Mexican and South American populations. Quantitative trait differentiation (Q ST) was greater than neutral allele frequency differentiation (F ST) for many traits, signaling directional selection between pairs of populations. All populations exhibited higher fitness metric values when grown at their native elevation, and Mexican landraces had higher fitness than South American landraces when grown in these Mexican sites. As environmental distance between landraces' native collection sites and common garden sites increased, fitness values dropped, suggesting landraces are adapted to environmental conditions at their natal sites. Correlations between fitness and anthocyanin pigmentation and macrohair traits were stronger in the highland site than the lowland site, supporting their status as highland-adaptive. These results give substance to the long-held presumption of local adaptation of New World maize landraces to elevation and other environmental variables across North and South America.

A B73×Palomero Toluqueño mapping population reveals local adaptation in Mexican highland maize.

Generations of farmer selection in the central Mexican highlands have produced unique maize varieties adapted to the challenges of the local environment. In addition to possessing great agronomic and cultural value, Mexican highland maize represents a good system for the study of local adaptation and acquisition of adaptive phenotypes under cultivation. In this study, we characterize a recombinant inbred line population derived from the B73 reference line and the Mexican highland maize variety Palomero Toluqueño. B73 and Palomero Toluqueño showed classic rank-changing differences in performance between lowland and highland field sites, indicative of local adaptation. Quantitative trait mapping identified genomic regions linked to effects on yield components that were conditionally expressed depending on the environment. For the principal genomic regions associated with ear weight and total kernel number, the Palomero Toluqueño allele conferred an advantage specifically in the highland site, consistent with local adaptation. We identified Palomero Toluqueño alleles associated with expression of characteristic highland traits, including reduced tassel branching, increased sheath pigmentation and the presence of sheath macrohairs. The oligogenic architecture of these three morphological traits supports their role in adaptation, suggesting they have arisen from consistent directional selection acting at distinct points across the genome. We discuss these results in the context of the origin of phenotypic novelty during selection, commenting on the role of de novo mutation and the acquisition of adaptive variation by gene flow from endemic wild relatives.

Genome sequencing reveals evidence of adaptive variation in the genus Zea.

Maize is a globally valuable commodity and one of the most extensively studied genetic model organisms. However, we know surprisingly little about the extent and potential utility of the genetic variation found in wild relatives of maize. Here, we characterize a high-density genomic variation map from 744 genomes encompassing maize and all wild taxa of the genus Zea, identifying over 70 million single-nucleotide polymorphisms. The variation map reveals evidence of selection within taxa displaying novel adaptations. We focus on adaptive alleles in highland teosinte and temperate maize, highlighting the key role of flowering-time-related pathways in their adaptation. To show the utility of variants in these data, we generate mutant alleles for two flowering-time candidate genes. This work provides an extensive sampling of the genetic diversity of Zea, resolving questions on evolution and identifying adaptive variants for direct use in modern breeding.

The arches and spandrels of maize domestication, adaptation, and improvement.

People living in the Balsas River basin in southwest México domesticated maize from the bushy grass teosinte. Nine thousand years later, in 2021, Ms. Deb Haaland - a member of the Pueblo of Laguna tribe of New Mexico - wore a dress adorned with a cornstalk when she was sworn in as the Secretary of Interior of the United States of America. This choice of garment highlights the importance of the coevolution of maize and the farmers who, through careful selection over thousands of years, domesticated maize and adapted the physiology and shoot architecture of maize to fit local environments and growth habits. Some traits such as tillering were directly selected on (arches), and others such as tassel size are the by-products (spandrels) of maize evolution. Here, we review current knowledge of the underlying cellular, developmental, physiological, and metabolic processes that were selected by farmers and breeders, which have positioned maize as a top global staple crop.

Identification of a tri-iron(III), tri-citrate complex in the xylem sap of iron-deficient tomato resupplied with iron: new insights into plant iron long-distance transport.

The identification of Fe transport forms in plant xylem sap is crucial to the understanding of long-distance Fe transport processes in plants. Previous studies have proposed that Fe may be transported as an Fe-citrate complex in plant xylem sap, but such a complex has never been detected. In this study we report the first direct and unequivocal identification of a natural Fe complex in plant xylem sap. A tri-Fe(III), tri-citrate complex (Fe(3)Cit(3)) was found in the xylem sap of Fe-deficient tomato (Solanum lycopersicum Mill. cv. 'Tres Cantos') resupplied with Fe, by using an integrated mass spectrometry approach based on exact molecular mass, isotopic signature and Fe determination and retention time. This complex has been modeled as having an oxo-bridged tri-Fe core. A second complex, a di-Fe(III), di-citrate complex was also detected in Fe-citrate standards along with Fe(3)Cit(3), with the allocation of Fe between the two complexes depending on the Fe to citrate ratio. These results provide evidence for Fe-citrate complex xylem transport in plants. The consequences for the role of Fe to citrate ratio in long-distance transport of Fe in xylem are also discussed.

Changes in the proteomic and metabolic profiles of Beta vulgaris root tips in response to iron deficiency and resupply.

BACKGROUND: Plants grown under iron deficiency show different morphological, biochemical and physiological changes. These changes include, among others, the elicitation of different strategies to improve the acquisition of Fe from the rhizosphere, the adjustment of Fe homeostasis processes and a reorganization of carbohydrate metabolism. The application of modern techniques that allow the simultaneous and untargeted analysis of multiple proteins and metabolites can provide insight into multiple processes taking place in plants under Fe deficiency. The objective of this study was to characterize the changes induced in the root tip proteome and metabolome of sugar beet plants in response to Fe deficiency and resupply. RESULTS: Root tip extract proteome maps were obtained by 2-D isoelectric focusing polyacrylamide gel electrophoresis, and approximately 140 spots were detected. Iron deficiency resulted in changes in the relative amounts of 61 polypeptides, and 22 of them were identified by mass spectrometry (MS). Metabolites in root tip extracts were analyzed by gas chromatography-MS, and more than 300 metabolites were resolved. Out of 77 identified metabolites, 26 changed significantly with Fe deficiency. Iron deficiency induced increases in the relative amounts of proteins and metabolites associated to glycolysis, tri-carboxylic acid cycle and anaerobic respiration, confirming previous studies. Furthermore, a protein not present in Fe-sufficient roots, dimethyl-8-ribityllumazine (DMRL) synthase, was present in high amounts in root tips from Fe-deficient sugar beet plants and gene transcript levels were higher in Fe-deficient root tips. Also, a marked increase in the relative amounts of the raffinose family of oligosaccharides (RFOs) was observed in Fe-deficient plants, and a further increase in these compounds occurred upon short term Fe resupply. CONCLUSIONS: The increases in DMRL synthase and in RFO sugars were the major changes induced by Fe deficiency and resupply in root tips of sugar beet plants. Flavin synthesis could be involved in Fe uptake, whereas RFO sugars could be involved in the alleviation of oxidative stress, C trafficking or cell signalling. Our data also confirm the increase in proteins and metabolites related to carbohydrate metabolism and TCA cycle pathways.

Inoculation with the mycorrhizal fungus Rhizophagus irregularis modulates the relationship between root growth and nutrient content in maize (Zea mays ssp. mays L.).

Plant root systems play a fundamental role in nutrient and water acquisition. In resource-limited soils, modification of root system architecture is an important strategy to optimize plant performance. Most terrestrial plants also form symbiotic associations with arbuscular mycorrhizal fungi to maximize nutrient uptake. In addition to direct delivery of nutrients, arbuscular mycorrhizal fungi benefit the plant host by promoting root growth. Here, we aimed to quantify the impact of arbuscular mycorrhizal symbiosis on root growth and nutrient uptake in maize. Inoculated plants showed an increase in both biomass and the total content of twenty quantified elements. In addition, image analysis showed mycorrhizal plants to have denser, more branched root systems. For most of the quantified elements, the increase in content in mycorrhizal plants was proportional to root and overall plant growth. However, the increase in boron, calcium, magnesium, phosphorus, sulfur, and strontium was greater than predicted by root system size alone, indicating fungal delivery to be supplementing root uptake.

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.