Diversity and divergence: evolution of secondary metabolism in the tropical tree genus Inga.

Plants are widely recognized as chemical factories, with each species producing dozens to hundreds of unique secondary metabolites. These compounds shape the interactions between plants and their natural enemies. We explore the evolutionary patterns and processes by which plants generate chemical diversity, from evolving novel compounds to unique chemical profiles. We characterized the chemical profile of one-third of the species of tropical rainforest trees in the genus Inga (c. 100, Fabaceae) using ultraperformance liquid chromatography-mass spectrometry-based metabolomics and applied phylogenetic comparative methods to understand the mode of chemical evolution. We show: each Inga species contain structurally unrelated compounds and high levels of phytochemical diversity; closely related species have divergent chemical profiles, with individual compounds, compound classes, and chemical profiles showing little-to-no phylogenetic signal; at the evolutionary time scale, a species' chemical profile shows a signature of divergent adaptation. At the ecological time scale, sympatric species were the most divergent, implying it is also advantageous to maintain a unique chemical profile from community members; finally, we integrate these patterns with a model for how chemical diversity evolves. Taken together, these results show that phytochemical diversity and divergence are fundamental to the ecology and evolution of plants.

Genome-wide data reveal cryptic diversity and hybridization in a group of tree ferns.

Discovery of cryptic diversity is essential to understanding both the process of speciation and the conservation of species. Determining species boundaries in fern lineages represents a major challenge due to lack of morphologically diagnostic characters and frequent hybridization. Genomic data has substantially enhanced our understanding of the speciation process, increased the resolution of species delimitation studies, and led to the discovery of cryptic diversity. Here, we employed restriction-site-associated DNA sequencing (RAD-seq) and integrated phylogenomic and population genomic analyses to investigate phylogenetic relationships and evolutionary history of 16 tree ferns with marginate scales (Cyatheaceae) from China and Vietnam. We conducted multiple species delimitation analyses using the multispecies coalescent (MSC) model and novel approaches based on genealogical divergence index (gdi) and isolation by distance (IBD). In addition, we inferred species trees using concatenation and several coalescent-based methods, and assessed hybridization patterns and rate of gene flow across the phylogeny. We obtained highly supported and generally congruent phylogenies inferred from concatenated and summary-coalescent methods, and the monophyly of all currently recognized species were strongly supported. Our results revealed substantial evidence of cryptic diversity in three widely distributed Gymnosphaera species, each of which was composite of two highly structure lineages that may correspond to cryptic species. We found that hybridization was fairly common between not only closely related species, but also distantly related species. Collectively, it appears that scaly tree ferns may contain cryptic diversity and hybridization has played an important role throughout the evolutionary history of this group.

Genomes shed light on the evolution of Begonia, a mega-diverse genus.

Clarifying the evolutionary processes underlying species diversification and adaptation is a key focus of evolutionary biology. Begonia (Begoniaceae) is one of the most species-rich angiosperm genera with c. 2000 species, most of which are shade-adapted. Here, we present chromosome-scale genome assemblies for four species of Begonia (B. loranthoides, B. masoniana, B. darthvaderiana and B. peltatifolia), and whole genome shotgun data for an additional 74 Begonia representatives to investigate lineage evolution and shade adaptation of the genus. The four genome assemblies range in size from 331.75 Mb (B. peltatifolia) to 799.83 Mb (B. masoniana), and harbor 22 059-23 444 protein-coding genes. Synteny analysis revealed a lineage-specific whole-genome duplication (WGD) that occurred just before the diversification of Begonia. Functional enrichment of gene families retained after WGD highlights the significance of modified carbohydrate metabolism and photosynthesis possibly linked to shade adaptation in the genus, which is further supported by expansions of gene families involved in light perception and harvesting. Phylogenomic reconstructions and genomics studies indicate that genomic introgression has also played a role in the evolution of Begonia. Overall, this study provides valuable genomic resources for Begonia and suggests potential drivers underlying the diversity and adaptive evolution of this mega-diverse clade.

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event.

The consequences of the Cretaceous-Paleogene (K-Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous-Paleogene (K-Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events].

Large-scale genomic sequence data resolve the deepest divergences in the legume phylogeny and support a near-simultaneous evolutionary origin of all six subfamilies.

Phylogenomics is increasingly used to infer deep-branching relationships while revealing the complexity of evolutionary processes such as incomplete lineage sorting, hybridization/introgression and polyploidization. We investigate the deep-branching relationships among subfamilies of the Leguminosae (or Fabaceae), the third largest angiosperm family. Despite their ecological and economic importance, a robust phylogenetic framework for legumes based on genome-scale sequence data is lacking. We generated alignments of 72 chloroplast genes and 7621 homologous nuclear-encoded proteins, for 157 and 76 taxa, respectively. We analysed these with maximum likelihood, Bayesian inference, and a multispecies coalescent summary method, and evaluated support for alternative topologies across gene trees. We resolve the deepest divergences in the legume phylogeny despite lack of phylogenetic signal across all chloroplast genes and the majority of nuclear genes. Strongly supported conflict in the remainder of nuclear genes is suggestive of incomplete lineage sorting. All six subfamilies originated nearly simultaneously, suggesting that the prevailing view of some subfamilies as 'basal' or 'early-diverging' with respect to others should be abandoned, which has important implications for understanding the evolution of legume diversity and traits. Our study highlights the limits of phylogenetic resolution in relation to rapid successive speciation.

Hybrid capture of 964 nuclear genes resolves evolutionary relationships in the mimosoid legumes and reveals the polytomous origins of a large pantropical radiation.

PREMISE: Targeted enrichment methods facilitate sequencing of hundreds of nuclear loci to enhance phylogenetic resolution and elucidate why some parts of the "tree of life" are difficult (if not impossible) to resolve. The mimosoid legumes are a prominent pantropical clade of ~3300 species of woody angiosperms for which previous phylogenies have shown extensive lack of resolution, especially among the species-rich and taxonomically challenging ingoids. METHODS: We generated transcriptomes to select low-copy nuclear genes, enrich these via hybrid capture for representative species of most mimosoid genera, and analyze the resulting data using de novo assembly and various phylogenomic tools for species tree inference. We also evaluate gene tree support and conflict for key internodes and use phylogenetic network analysis to investigate phylogenetic signal across the ingoids. RESULTS: Our selection of 964 nuclear genes greatly improves phylogenetic resolution across the mimosoid phylogeny and shows that the ingoid clade can be resolved into several well-supported clades. However, nearly all loci show lack of phylogenetic signal for some of the deeper internodes within the ingoids. CONCLUSIONS: Lack of resolution in the ingoid clade is most likely the result of hyperfast diversification, potentially causing a hard polytomy of six or seven lineages. The gene set for targeted sequencing presented here offers great potential to further enhance the phylogeny of mimosoids and the wider Caesalpinioideae with denser taxon sampling, to provide a framework for taxonomic reclassification, and to study the ingoid radiation.

Chemocoding as an identification tool where morphological- and DNA-based methods fall short: Inga as a case study.

The need for species identification and taxonomic discovery has led to the development of innovative technologies for large-scale plant identification. DNA barcoding has been useful, but fails to distinguish among many species in species-rich plant genera, particularly in tropical regions. Here, we show that chemical fingerprinting, or 'chemocoding', has great potential for plant identification in challenging tropical biomes. Using untargeted metabolomics in combination with multivariate analysis, we constructed species-level fingerprints, which we define as chemocoding. We evaluated the utility of chemocoding with species that were defined morphologically and subject to next-generation DNA sequencing in the diverse and recently radiated neotropical genus Inga (Leguminosae), both at single study sites and across broad geographic scales. Our results show that chemocoding is a robust method for distinguishing morphologically similar species at a single site and for identifying widespread species across continental-scale ranges. Given that species are the fundamental unit of analysis for conservation and biodiversity research, the development of accurate identification methods is essential. We suggest that chemocoding will be a valuable additional source of data for a quick identification of plants, especially for groups where other methods fall short.

Transcriptome mining for phylogenetic markers in a recently radiated genus of tropical plants (Renealmia L.f., Zingiberaceae).

The reconstruction of relationships within species-rich groups that have recently evolved in biodiversity hotspots is hampered by a lack of phylogenetically informative markers. It is also made difficult by the lack of sampling necessary to reconstruct a species-level phylogeny. We use transcriptome mining to search for markers to reconstruct a phylogeny of the amphi-Atlantic genus Renealmia L. f. (Zingiberaceae). We recover seven introns from single copy genes and use them to reconstruct the phylogeny of the genus together with a commonly used phylogenetic marker, internal transcribed spacers of ribosomal DNA (ITS) that has previously been used to reconstruct the phylogeny of the genus. We targeted genes with low numbers of base pairs that improves sequencing success using highly degraded DNA from herbarium specimens. The use of herbarium specimens greatly increased the number of species in the study as these were readily available in historical collections. Data were obtained for 14 of the 17 African species and 54 of the 65 Neotropical species. The phylogeny was well-supported for a number of Renealmia subgroups although relationships among those clades remained poorly supported.

Ultrastructure and development of non-contiguous stomatal clusters and helicocytic patterning in Begonia.

Background and Aims: Helicocytic stomata are characterized by an inward spiral of mesogenous cells surrounding a central stomatal pore. They represent a relatively rare feature that occurs in some drought-tolerant angiosperm species. In some Begonia species with thick leaves, the stomata are not only helicocytic but also clustered into groups that are spaced apart by at least one cell. This paper presents a detailed ontogenetic study of this characteristic non-contiguous stomatal patterning in a developmental and phylogenetic context. Methods: Light microscopy and both scanning and transmission electron microscopy were used to examine stomatal development in several species of Begonia. Published reports of stomatal development in Begonia and other angiosperms were reviewed to provide a comprehensive discussion of the evolution of stomatal patterning. Key Results: Helicocytic stomata develop from meristemoids that undergo a series of oriented asymmetric divisions to produce a spiral of mesogene stomatal lineage ground cells (SLGCs) surrounding a stoma. A clear developmental similarity between anisocytic and helicocytic stomata is positively correlated with the number of iterations of amplifying divisions that result in SLGCs. Stomatal clusters develop from asymmetric divisions in neighbouring SLGCs. Within each cluster, non-contiguous spacing of meristemoids is maintained by asymmetric divisions oriented away from each developing meristemoid. Conclusions: Formation of non-contiguous stomatal clusters in Begonia relies on two primary developmental factors in the epidermis: an inwardly spiralling series of amplifying divisions that result in helicocytic stomata, and the development of a variable number of meristemoids from neighbouring SLGCs within each cluster. Optimization of these features on an angiosperm phylogeny indicates that the occurrence of amplifying divisions could be pre-adaptive for these factors. Both factors have been thoroughly studied in terms of developmental genetics in Arabidopsis, suggesting gene orthologues that could be implicated in Begonia stomatal patterning.

Maintenance of species boundaries in a Neotropical radiation of Begonia.

A major goal of evolutionary biology is to determine the mechanisms generating biodiversity. In Begonia, one of the largest plant genera (1900+ species), it has been postulated that the high number of endemic species is a by-product of low gene flow among populations, which predisposes the group to speciation. However, this model of divergence requires that reproductive barriers accumulate rapidly among diverging species that overlap in their geographic ranges, otherwise speciation will be opposed by homogenizing gene flow in zones of secondary contact. Here, we test the outcomes of secondary contact in Begonia by genotyping multiple sympatric sites with 12 nuclear and seven plastid loci. We show that three sites of secondary contact between B. heracleifolia and B. nelumbiifolia are highly structured, mostly containing parental genotypes, with few F1 hybrids. A sympatric site between B. heracleifolia and B. sericoneura contains a higher proportion of F1s, but little evidence of introgression. The lack of later-generation hybrids contrasts with that documented in many other plant taxa, where introgression is extensive. Our results, in conjunction with previous genetic work, show that Begonia demonstrate properties making them exceptionally prone to speciation, at multiple stages along the divergence continuum. Not only are populations weakly connected by gene flow, promoting allopatric speciation, but species often show strong reproductive barriers in secondary contact. Whether similar mechanisms contribute to diversification in other large genera remains to be tested.

The evolution of sex ratio differences and inflorescence architectures in Begonia (Begoniaceae).

PREMISE OF THE STUDY: A major benefit conferred by monoecy is the ability to alter floral sex ratio in response to selection. In monoecious species that produce flowers of a given sex at set positions on the inflorescence, floral sex ratio may be related to inflorescence architecture. We studied the loci underlying differences in inflorescence architecture between two monoecious Begonia species and related this to floral sex ratios. METHODS: We performed trait comparisons and quantitative trait locus (QTL) mapping in a segregating backcross population between Central American Begonia plebeja and B. conchifolia. We focused on traits related to inflorescence architecture, sex ratios, and other reproductive traits. KEY RESULTS: The inflorescence branching pattern of B. conchifolia was more asymmetric than B. plebeja, which in turn affects the floral sex ratio. Colocalizing QTLs of moderate effect influenced both the number of male flowers and the fate decisions of axillary meristems, demonstrating the close link between inflorescence architecture and sex ratio. Additional QTLs were found for stamen number (30% variance explained, VE) and pollen sterility (12.3% VE). CONCLUSIONS: One way in which Begonia species develop different floral sex ratios is through modifications of their inflorescence architecture. The potential pleiotropic action of QTL on inflorescence branching and floral sex ratios may have major implications for trait evolution and responses to selection. The presence of a single QTL of large effect on stamen number may allow rapid divergence for this key floral trait. We propose candidate loci for stamen number and inflorescence branching for future characterization.

Target capture and genome skimming for plant diversity studies.

Recent technological advances in long-read high-throughput sequencing and assembly methods have facilitated the generation of annotated chromosome-scale whole-genome sequence data for evolutionary studies; however, generating such data can still be difficult for many plant species. For example, obtaining high-molecular-weight DNA is typically impossible for samples in historical herbarium collections, which often have degraded DNA. The need to fast-freeze newly collected living samples to conserve high-quality DNA can be complicated when plants are only found in remote areas. Therefore, short-read reduced-genome representations, such as target capture and genome skimming, remain important for evolutionary studies. Here, we review the pros and cons of each technique for non-model plant taxa. We provide guidance related to logistics, budget, the genomic resources previously available for the target clade, and the nature of the study. Furthermore, we assess the available bioinformatic analyses, detailing best practices and pitfalls, and suggest pathways to combine newly generated data with legacy data. Finally, we explore the possible downstream analyses allowed by the type of data generated using each technique. We provide a practical guide to help researchers make the best-informed choice regarding reduced genome representation for evolutionary studies of non-model plants in cases where whole-genome sequencing remains impractical.

Developing the Protocol Infrastructure for DNA Sequencing Natural History Collections.

Intentionally preserved biological material in natural history collections represents a vast repository of biodiversity. Advances in laboratory and sequencing technologies have made these specimens increasingly accessible for genomic analyses, offering a window into the genetic past of species and often permitting access to information that can no longer be sampled in the wild. Due to their age, preparation and storage conditions, DNA retrieved from museum and herbarium specimens is often poor in yield, heavily fragmented and biochemically modified. This not only poses methodological challenges in recovering nucleotide sequences, but also makes such investigations susceptible to environmental and laboratory contamination. In this paper, we review the practical challenges associated with making the recovery of DNA sequence data from museum collections more routine. We first review key operational principles and issues to address, to guide the decision-making process and dialogue between researchers and curators about when and how to sample museum specimens for genomic analyses. We then outline the range of steps that can be taken to reduce the likelihood of contamination including laboratory set-ups, workflows and working practices. We finish by presenting a series of case studies, each focusing on protocol practicalities for the application of different mainstream methodologies to museum specimens including: (i) shotgun sequencing of insect mitogenomes, (ii) whole genome sequencing of insects, (iii) genome skimming to recover plant plastid genomes from herbarium specimens, (iv) target capture of multi-locus nuclear sequences from herbarium specimens, (v) RAD-sequencing of bird specimens and (vi) shotgun sequencing of ancient bovid bone samples.

A highly contiguous, scaffold-level nuclear genome assembly for the fever tree (Cinchona pubescens Vahl) as a novel resource for Rubiaceae research.

The Andean fever tree (Cinchona L.; Rubiaceae) is a source of bioactive quinine alkaloids used to treat malaria. C. pubescens Vahl is a valuable cash crop within its native range in northwestern South America, however, genomic resources are lacking. Here we provide the first highly contiguous and annotated nuclear and plastid genome assemblies using Oxford Nanopore PromethION-derived long-read and Illumina short-read data. Our nuclear genome assembly comprises 603 scaffolds with a total length of 904 Mbp (∼82% of the full genome based on a genome size of 1.1 Gbp/1C). Using a combination of de novo and reference-based transcriptome assemblies we annotated 72,305 coding sequences comprising 83% of the BUSCO gene set and 4.6% fragmented sequences. Using additional plastid and nuclear datasets we place C. pubescens in the Gentianales order. This first genomic resource for C. pubescens opens new research avenues, including the analysis of alkaloid biosynthesis in the fever tree.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family.

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

A complex case of simple leaves: indeterminate leaves co-express ARP and KNOX1 genes.

The mutually exclusive relationship between ARP and KNOX1 genes in the shoot apical meristem and leaf primordia in simple leaved plants such as Arabidopsis has been well characterized. Overlapping expression domains of these genes in leaf primordia have been described for many compound leaved plants such as Solanum lycopersicum and Cardamine hirsuta and are regarded as a characteristic of compound leaved plants. Here, we present several datasets illustrating the co-expression of ARP and KNOX1 genes in the shoot apical meristem, leaf primordia, and developing leaves in plants with simple leaves and simple primordia. Streptocarpus plants produce unequal cotyledons due to the continued activity of a basal meristem and produce foliar leaves termed "phyllomorphs" from the groove meristem in the acaulescent species Streptocarpus rexii and leaves from a shoot apical meristem in the caulescent Streptocarpus glandulosissimus. We demonstrate that the simple leaves in both species possess a greatly extended basal meristematic activity that persists over most of the leaf's growth. The area of basal meristem activity coincides with the co-expression domain of ARP and KNOX1 genes. We suggest that the co-expression of ARP and KNOX1 genes is not exclusive to compound leaved plants but is associated with foci of meristematic activity in leaves.

Spatially restricted microRNA directs leaf polarity through ARGONAUTE1.

Gene regulation by RNA interference requires the functions of the PAZ domain protein Argonaute. In plants, mutations in ARGONAUTE1 (AGO1) are associated with distinctive developmental defects that suggest a role for microRNA (miRNA) in organ polarity. Potential targets of miRNA regulation are the homeodomain/leucine zipper genes PHABULOSA (PHB) and PHAVOLUTA (PHV). These genes are expressed in a polar fashion in leaf primordia and are required for adaxial cell fate. Here we show that a 21-nucleotide miRNA that directs cleavage of PHB/PHV messenger RNA accumulates first in the embryonic meristem, and then in the abaxial domain of the developing leaf. miRNA distribution is disrupted by mutations in AGO1, indicating that AGO1 affects the regulation of miRNA. In addition, interactions between homeodomain/leucine zipper genes and an allelic series of ago1 indicate that miRNA acts as a signal to specify leaf polarity.

Mixing and matching pathways in leaf polarity.

Leaves extend a large, porous surface to the environment to catch light and exchange gasses. The extension of the lamina is produced by the interaction of an upper (adaxial) and a lower (abaxial) domain in the developing leaf primordium. Recent studies have revealed that conserved genetic pathways, involving small regulatory RNAs and several distinct transcription factor families, have key roles in adaxial-abaxial patterning, suggesting candidate signals that convey positional information within the shoot to the newly initiated leaf. The interactions of the polarity pathways are distinguished by mutual antagonism and by redundancies. Analysis of these pathways in different model organisms reveals a surprising diversity in the genetic control of such a fundamental developmental process.

Plant stem cells: divergent pathways and common themes in shoots and roots.

Stem cells in plant shoot and root meristems are maintained throughout the life of the plant and produce somatic daughter cells that make up the body of the plant. Plant stem cells can also be derived from somatic cells in vivo and in vitro. Recent findings are refining our knowledge of signaling pathways that define stem cell fate and specify either shoot or root stem cell function. New evidence also highlights a role for epigenetic mechanisms in controlling stem cell fate.

Tracking of Host Defenses and Phylogeny During the Radiation of Neotropical Inga-Feeding Sawflies (Hymenoptera; Argidae).

Coevolutionary theory has long predicted that the arms race between plants and herbivores is a major driver of host selection and diversification. At a local scale, plant defenses contribute significantly to the structure of herbivore assemblages and the high alpha diversity of plants in tropical rain forests. However, the general importance of plant defenses in host associations and divergence at regional scales remains unclear. Here, we examine the role of plant defensive traits and phylogeny in the evolution of host range and species divergence in leaf-feeding sawflies of the family Argidae associated with Neotropical trees in the genus Inga throughout the Amazon, the Guiana Shield and Panama. Our analyses show that the phylogenies of both the sawfly herbivores and their Inga hosts are congruent, and that sawflies radiated at approximately the same time, or more recently than their Inga hosts. Analyses controlling for phylogenetic effects show that the evolution of host use in the sawflies associated with Inga is better correlated with Inga chemistry than with Inga phylogeny, suggesting a pattern of delayed host tracking closely tied to host chemistry. Finally, phylogenetic analyses show that sister species of Inga-sawflies are dispersed across the Neotropics, suggesting a role for allopatric divergence and vicariance in Inga diversification. These results are consistent with the idea that host defensive traits play a key role not only in structuring the herbivore assemblages at a single site, but also in the processes shaping host association and species divergence at a regional scale.