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.

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.

Species Distribution Modelling: Contrasting presence-only models with plot abundance data

Species distribution models (SDMs) are widely used in ecology and conservation. Presence-only SDMs such as MaxEnt frequently use natural history collections (NHCs) as occurrence data, given their huge numbers and accessibility. NHCs are often spatially biased which may generate inaccuracies in SDMs. Here, we test how the distribution of NHCs and MaxEnt predictions relates to a spatial abundance model, based on a large plot dataset for Amazonian tree species, using inverse distance weighting (IDW). We also propose a new pipeline to deal with inconsistencies in NHCs and to limit the area of occupancy of the species. We found a significant but weak positive relationship between the distribution of NHCs and IDW for 66% of the species. The relationship between SDMs and IDW was also significant but weakly positive for 95% of the species, and sensitivity for both analyses was high. Furthermore, the pipeline removed half of the NHCs records. Presence-only SDM applications should consider this limitation, especially for large biodiversity assessments projects, when they are automatically generated without subsequent checking. Our pipeline provides a conservative estimate of a species’ area of occupancy, within an area slightly larger than its extent of occurrence, compatible to e.g. IUCN red list assessments.

Evolutionary heritage influences Amazon tree ecology

Lineages tend to retain ecological characteristics of their ancestors through time. However, for some traits, selection during evolutionary history may have also played a role in determining trait values. To address the relative importance of these processes requires large-scale quantification of traits and evolutionary relationships among species. The Amazonian tree flora comprises a high diversity of angiosperm lineages and species with widely differing life-history characteristics, providing an excellent system to investigate the combined influences of evolutionary heritage and selection in determining trait variation. We used trait data related to the major axes of life-history variation among tropical trees (e.g. growth and mortality rates) from 577 inventory plots in closed-canopy forest, mapped onto a phylogenetic hypothesis spanning more than 300 genera including all major angiosperm clades to test for evolutionary constraints on traits. We found significant phylogenetic signal (PS) for all traits, consistent with evolutionarily related genera having more similar characteristics than expected by chance. Although there is also evidence for repeated evolution of pioneer and shade tolerant life-history strategies within independent lineages, the existence of significant PS allows clearer predictions of the links between evolutionary diversity, ecosystem function and the response of tropical forests to global change.

RADseq dataset with 90% missing data fully resolves recent radiation of Petalidium (Acanthaceae) in the ultra-arid deserts of Namibia

Deserts, even those at tropical latitudes, often have strikingly low levels of plant diversity, particularly within genera. One remarkable exception to this pattern is the genus Petalidium (Acanthaceae), in which 37 of 40 named species occupy one of the driest environments on Earth, the Namib Desert of Namibia and neighboring Angola. To contribute to understanding this enigmatic diversity, we generated RADseq data for 47 accessions of Petalidium representing 22 species. We explored the impacts of 18 different combinations of assembly parameters in de novo assembly of the data across nine levels of missing data plus a best practice assembly using a reference Acanthaceae genome for a total of 171 sequence datasets assembled. RADseq data assembled at several thresholds of missing data, including 90% missing data, yielded phylogenetic hypotheses of Petalidium that were confidently and nearly fully resolved, which is notable given that divergence time analyses suggest a crown age for African species of 3.6-1.4 Ma. De novo assembly of our data yielded the most strongly supported and well-resolved topologies; in contrast, reference-based assembly performed poorly, perhaps due in part to moderate phylogenetic divergence between the reference genome, Ruellia speciosa, and the ingroup. Overall, we found that Petalidium, despite the harshness of the environment in which species occur, shows a net diversification rate (0.8-2.1 species per my) on par with those of diverse genera in tropical, Mediterranean, and alpine environments.

Aboveground Carbon Storage and Its Links to Stand Structure, Tree Diversity and Floristic Composition in South-Eastern Tanzania

African savannas and dry forests represent a large, but poorly quantified store of biomass carbon and biodiversity. Improving this information is hindered by a lack of recent forest inventories, which are necessary for calibrating earth observation data and for evaluating the relationship between carbon stocks and tree diversity in the context of forest conservation (for example, REDD+). Here, we present new inventory data from south-eastern Tanzania, comprising more than 15,000 trees at 25 locations located across a gradient of aboveground woody carbon (AGC) stocks. We find that larger trees disproportionately contribute to AGC, with the largest 3.7% of individuals containing half the carbon. Tree species diversity and carbon stocks were positively related, implying a potential functional relationship between the two, and a ‘win–win’ scenario for conservation; however, lower biomass areas also contain diverse species assemblages meaning that carbon-oriented conservation may miss important areas of biodiversity. Despite these variations, we find that total tree abundance and biomass is skewed towards a few locally dominant species, with eight and nine species (5.7% of the total) accounting for over half the total measured trees and carbon, respectively. This finding implies that carbon production in these areas is channelled through a small number of relatively abundant species. Our results provide key insights into the structure and functioning of these heterogeneous ecosystems and indicate the need for novel strategies for future measurement and monitoring of carbon stocks and biodiversity, including the use for larger plots to capture spatial variations in large tree density and AGC stocks, and to allow the calibration of earth observation data.

Coevolutionary arms race versus host defense chase in a tropical herbivore-plant system

Coevolutionary models suggest that herbivores drive diversification and community composition in plants. For herbivores, many questions remain regarding how plant defenses shape host choice and community structure. We addressed these questions using the tree genus Inga and its lepidopteran herbivores in the Amazon. We constructed phylogenies for both plants and insects and quantified host associations and plant defenses. We found that similarity in herbivore assemblages between Inga species was correlated with similarity in defenses. There was no correlation with phylogeny, a result consistent with our observations that the expression of defenses in Inga is independent of phylogeny. Furthermore, host defensive traits explained 40% of herbivore community similarity. Analyses at finer taxonomic scales showed that different lepidopteran clades select hosts based on different defenses, suggesting taxon-specific histories of herbivore-host plant interactions. Finally, we compared the phylogeny and defenses of Inga to phylogenies for the major lepidopteran clades. We found that closely related herbivores fed on Inga with similar defenses rather than on closely related plants. Together, these results suggest that plant defenses might be more evolutionarily labile than the herbivore traits related to host association. Hence, there is an apparent asymmetry in the evolutionary interactions between Inga and its herbivores. Although plants may evolve under selection by herbivores, we hypothesize that herbivores may not show coevolutionary adaptations, but instead "chase" hosts based on the herbivore's own traits at the time that they encounter a new host, a pattern more consistent with resource tracking than with the arms race model of coevolution.

Biogeographic distributions of neotropical trees reflect their directly measured drought tolerances

High levels of species diversity hamper current understanding of how tropical forests may respond to environmental change. In the tropics, water availability is a leading driver of the diversity and distribution of tree species, suggesting that many tropical taxa may be physiologically incapable of tolerating dry conditions, and that their distributions along moisture gradients can be used to predict their drought tolerance. While this hypothesis has been explored at local and regional scales, large continental-scale tests are lacking. We investigate whether the relationship between drought-induced mortality and distributions holds continentally by relating experimental and observational data of drought-induced mortality across the Neotropics to the large-scale bioclimatic distributions of 115 tree genera. Across the different experiments, genera affiliated to wetter climatic regimes show higher drought-induced mortality than dry-affiliated ones, even after controlling for phylogenetic relationships. This pattern is stronger for adult trees than for saplings or seedlings, suggesting that the environmental filters exerted by drought impact adult tree survival most strongly. Overall, our analysis of experimental, observational, and bioclimatic data across neotropical forests suggests that increasing moisture-stress is indeed likely to drive significant changes in floristic composition.

Dispersal assembly of rain forest tree communities across the Amazon basin

We investigate patterns of historical assembly of tree communities across Amazonia using a newly developed phylogeny for the species-rich neotropical tree genus Inga We compare our results with those for three other ecologically important, diverse, and abundant Amazonian tree lineages, Swartzia, Protieae, and Guatteria Our analyses using phylogenetic diversity metrics demonstrate a clear lack of geographic phylogenetic structure, and show that local communities of Inga and regional communities of all four lineages are assembled by dispersal across Amazonia. The importance of dispersal in the biogeography of Inga and other tree genera in Amazonian and Guianan rain forests suggests that speciation is not driven by vicariance, and that allopatric isolation following dispersal may be involved in the speciation process. A clear implication of these results is that over evolutionary timescales, the metacommunity for any local or regional tree community in the Amazon is the entire Amazon basin.

Plant diversity patterns in neotropical dry forests and their conservation implications

Seasonally dry tropical forests are distributed across Latin America and the Caribbean and are highly threatened, with less than 10% of their original extent remaining in many countries. Using 835 inventories covering 4660 species of woody plants, we show marked floristic turnover among inventories and regions, which may be higher than in other neotropical biomes, such as savanna. Such high floristic turnover indicates that numerous conservation areas across many countries will be needed to protect the full diversity of tropical dry forests. Our results provide a scientific framework within which national decision-makers can contextualize the floristic significance of their dry forest at a regional and continental scale.

Maximising Synergy among Tropical Plant Systematists, Ecologists, and Evolutionary Biologists

Closer collaboration among ecologists, systematists, and evolutionary biologists working in tropical forests, centred on studies within long-term permanent plots, would be highly beneficial for their respective fields. With a key unifying theme of the importance of vouchered collection and precise identification of species, especially rare ones, we identify four priority areas where improving links between these communities could achieve significant progress in biodiversity and conservation science: (i) increasing the pace of species discovery; (ii) documenting species turnover across space and time; (iii) improving models of ecosystem change; and (iv) understanding the evolutionary assembly of communities and biomes.

Evolutionary patterns of volatile terpene emissions across 202 tropical tree species

Plant responses to natural enemies include formation of secondary metabolites acting as direct or indirect defenses. Volatile terpenes represent one of the most diverse groups of secondary metabolites. We aimed to explore evolutionary patterns of volatile terpene emission. We measured the composition of damage‐induced volatile terpenes from 202 Amazonian tree species, spanning the angiosperm phylogeny. Volatile terpenes were extracted with solid‐phase micro extraction and desorbed in a gas chromatography–mass spectrometry for compound identification. The chemical diversity of the terpene blend showed a strong phylogenetic signal as closely related species emitted a similar number of compounds. Closely related species also tended to have compositionally similar blends, although this relationship was weak. Meanwhile, the ability to emit a given compound showed no significant phylogenetic signal for 200 of 286 compounds, indicating a high rate of diversification in terpene synthesis and/or great variability in their expression. Three lineages (Magnoliales, Laurales, and Sapindales) showed exceptionally high rates of terpene diversification. Of the 70 compounds found in >10% of their species, 69 displayed significant correlated evolution with at least one other compound. These results provide insights into the complex evolutionary history of volatile terpenes in angiosperms, while highlighting the need for further research into this important class of compounds.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae)

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae)

Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.

Ficus insipida subsp. insipida (Moraceae) reveals the role of ecology in the phylogeography of widespread Neotropical rain forest tree species

AIM

To examine the phylogeography of Ficus insipida subsp. insipida in order to investigate patterns of spatial genetic structure across the Neotropics and within Amazonia.

LOCATION

Neotropics.

METHODS

Plastid DNA (trnH-psbA; 410 individuals from 54 populations) and nuclear ribosomal internal transcribed spacer (ITS; 85 individuals from 27 populations) sequences were sampled from Mexico to Bolivia, representing the full extent of the taxon's distribution. Divergence of plastid lineages was dated using a Bayesian coalescent approach. Genetic diversity was assessed with indices of haplotype and nucleotide diversities, and genetic structure was examined using spatial analysis of molecular variance (SAMOVA) and haplotype networks. Population expansion within Amazonia was tested using neutrality and mismatch distribution tests.

RESULTS

trnH-psbA sequences yielded 19 haplotypes restricted to either Mesoamerica or Amazonia; six haplotypes were found among ITS sequences. Diversification of the plastid DNA haplotypes began c. 14.6 Ma. Haplotype diversity for trnH-psbA was higher in Amazonia. Seven genetically differentiated SAMOVA groups were described for trnH-psbA, of which two were also supported by the presence of unique ITS sequences. Population expansion was suggested for both markers for the SAMOVA group that contains most Amazonian populations.

MAIN CONCLUSIONS

Our results show marked population genetic structure in F. insipida between Mesoamerica and Amazonia, implying that the Andes and seasonally dry areas of northern South America are eco-climatic barriers to its migration. This pattern is shared with other widespread pioneer species affiliated to wet habitats, indicating that the ecological characteristics of species may impact upon large-scale phylogeography. Ficus insipida also shows genetic structure in north-western Amazonia potentially related to pre-Pleistocene historical events. In contrast, evident population expansion elsewhere in Amazonia, in particular the presence of genetically uniform populations across the south-west, indicate recent colonization. Our findings are consistent with palaeoecological data that suggest recent post-glacial expansion of Amazonian forests in the south.

Fast demographic traits promote high diversification rates of Amazonian trees

The Amazon rain forest sustains the world's highest tree diversity, but it remains unclear why some clades of trees are hyperdiverse, whereas others are not. Using dated phylogenies, estimates of current species richness and trait and demographic data from a large network of forest plots, we show that fast demographic traits – short turnover times – are associated with high diversification rates across 51 clades of canopy trees. This relationship is robust to assuming that diversification rates are either constant or decline over time, and occurs in a wide range of Neotropical tree lineages. This finding reveals the crucial role of intrinsic, ecological variation among clades for understanding the origin of the remarkable diversity of Amazonian trees and forests.

Hyperdominance in the Amazonian tree flora

The vast extent of the Amazon Basin has historically restricted the study of its tree communities to the local and regional scales. Here, we provide empirical data on the commonness, rarity, and richness of lowland tree species across the entire Amazon Basin and Guiana Shield (Amazonia), collected in 1170 tree plots in all major forest types. Extrapolations suggest that Amazonia harbors roughly 16,000 tree species, of which just 227 (1.4%) account for half of all trees. Most of these are habitat specialists and only dominant in one or two regions of the basin. We discuss some implications of the finding that a small group of species--less diverse than the North American tree flora--accounts for half of the world's most diverse tree community.

Historical effects on beta diversity and community assembly in Amazonian trees

We present a unique perspective on the role of historical processes in community assembly by synthesizing analyses of species turnover among communities with environmental data and independent, population genetic-derived estimates of among-community dispersal. We sampled floodplain and terra firme communities of the diverse tree genus Inga (Fabaceae) across a 250-km transect in Amazonian Peru and found patterns of distance-decay in compositional similarity in both habitat types. However, conventional analyses of distance-decay masked a zone of increased species turnover present in the middle of the transect. We estimated past seed dispersal among the same communities by examining geographic plastid DNA variation for eight widespread Inga species and uncovered a population genetic break in the majority of species that is geographically coincident with the zone of increased species turnover. Analyses of these and 12 additional Inga species shared between two communities located on opposite sides of the zone showed that the populations experienced divergence 42,000-612,000 y ago. Our results suggest that the observed distance decay is the result not of environmental gradients or dispersal limitation coupled with ecological drift--as conventionally interpreted under neutral ecological theory--but rather of secondary contact between historically separated communities. Thus, even at this small spatial scale, historical processes seem to significantly impact species' distributions and community assembly. Other documented zones of increased species turnover found in the western Amazon basin or elsewhere may be related to similar historical processes.

Phylogenetic density dependence and environmental filtering predict seedling mortality in a tropical forest

Negative density dependence (NDD) and environmental filtering (EF) shape community assembly, but their relative importance is poorly understood. Recent studies have shown that seedling's mortality risk is positively related to the phylogenetic relatedness of neighbours. However, natural enemies, whose depredations often cause NDD, respond to functional traits of hosts rather than phylogenetic relatedness per se. To understand the roles of NDD and EF in community assembly, we assessed the effects on seedling mortality of functional similarity, phylogenetic relatedness and stem density of neighbouring seedlings and adults in a species-rich tropical forest. Mortality risks increased for common species when their functional traits departed substantially from the neighbourhood mean, and for all species when surrounded by close relatives. This indicates that NDD affects community assembly more broadly than does EF, and leads to the tentative conclusion that natural enemies respond to phylogenetically correlated traits. Our results affirm the prominence of NDD in structuring species-rich communities.