Echoes of ancient introgression punctuate stable genomic lineages in the evolution of figs.

Studies investigating the evolution of flowering plants have long focused on isolating mechanisms such as pollinator specificity. Some recent studies have proposed a role for introgressive hybridization between species, recognizing that isolating processes such as pollinator specialization may not be complete barriers to hybridization. Occasional hybridization may therefore lead to distinct yet reproductively connected lineages. We investigate the balance between introgression and reproductive isolation in a diverse clade using a densely sampled phylogenomic study of fig trees (Ficus, Moraceae). Codiversification with specialized pollinating wasps (Agaonidae) is recognized as a major engine of fig diversity, leading to about 850 species. Nevertheless, some studies have focused on the importance of hybridization in Ficus, highlighting the consequences of pollinator sharing. Here, we employ dense taxon sampling (520 species) throughout Moraceae and 1,751 loci to investigate phylogenetic relationships and the prevalence of introgression among species throughout the history of Ficus. We present a well-resolved phylogenomic backbone for Ficus, providing a solid foundation for an updated classification. Our results paint a picture of phylogenetically stable evolution within lineages punctuated by occasional local introgression events likely mediated by local pollinator sharing, illustrated by clear cases of cytoplasmic introgression that have been nearly drowned out of the nuclear genome through subsequent lineage fidelity. The phylogenetic history of figs thus highlights that while hybridization is an important process in plant evolution, the mere ability of species to hybridize locally does not necessarily translate into ongoing introgression between distant lineages, particularly in the presence of obligate plant-pollinator relationships.

Common spatial patterns of trees in various tropical forests: Small trees are associated with increased diversity at small spatial scales.

Tropical forests are notable for their high species diversity, even on small spatial scales, and right-skewed species and size abundance distributions. The role of individual species as drivers of the spatial organization of diversity in these forests has been explained by several hypotheses and processes, for example, stochastic dilution, negative density dependence, or gap dynamics. These processes leave a signature in spatial distribution of small trees, particularly in the vicinity of large trees, likely having stronger effects on their neighbors. We are exploring species diversity patterns within the framework of various diversity-generating hypotheses using individual species-area relationships. We used the data from three tropical forest plots (Wanang-Papua New Guinea, Barro Colorado Island-Panama, and Sinharaja-Sri Lanka) and included also the saplings (DBH ≥ 1 cm). Resulting cross-size patterns of species richness and evenness reflect the dynamics of saplings affected by the distribution of large trees. When all individuals with DBH ≥1 cm are included, ~50% of all tree species from the 25- or 50-ha plot can be found within 35 m radius of an individual tree. For all trees, 72%-78% of species were identified as species richness accumulators, having more species present in their surroundings than expected by null models. This pattern was driven by small trees as the analysis of DBH >10 cm trees showed much lower proportion of accumulators, 14%-65% of species identified as richness repellers and had low richness of surrounding small trees. Only 11%-26% of species had lower species evenness than was expected by null models. High proportions of species richness accumulators were probably due to gap dynamics and support Janzen-Connell hypothesis driven by competition or top-down control by pathogens and herbivores. Observed species diversity patterns show the importance of including small tree size classes in analyses of the spatial organization of diversity.

Language and ethnobiological skills decline precipitously in Papua New Guinea, the world's most linguistically diverse nation.

Papua New Guinea is home to >10% of the world's languages and rich and varied biocultural knowledge, but the future of this diversity remains unclear. We measured language skills of 6,190 students speaking 392 languages (5.5% of the global total) and modeled their future trends using individual-level variables characterizing family language use, socioeconomic conditions, students' skills, and language traits. This approach showed that only 58% of the students, compared to 91% of their parents, were fluent in indigenous languages, while the trends in key drivers of language skills (language use at home, proportion of mixed-language families, urbanization, students' traditional skills) predicted accelerating decline of fluency to an estimated 26% in the next generation of students. Ethnobiological knowledge declined in close parallel with language skills. Varied medicinal plant uses known to the students speaking indigenous languages are replaced by a few, mostly nonnative species for the students speaking English or Tok Pisin, the national lingua franca. Most (88%) students want to teach indigenous language to their children. While crucial for keeping languages alive, this intention faces powerful external pressures as key factors (education, cash economy, road networks, and urbanization) associated with language attrition are valued in contemporary society.

Spatial scaling of plant and bird diversity from 50 to 10,000 ha in a lowland tropical rainforest.

While there are numerous studies of diversity patterns both within local communities and at regional scales, the intermediate scale of tens to thousands of km2 is often neglected. Here we present detailed local data on plant communities (using 20 × 20 m plots) and bird communities (using point counts) for a 50 ha ForestGEO plot in lowland rainforest at Wanang, Papua New Guinea. We compare these local diversity patterns with those documented in the surrounding 10,000 ha of lowland rainforest. Woody plant species richness was lower within 50 ha (88% of 10,000 ha richness), even when both were surveyed with identical sampling effort. In contrast, bird communities exhibited identical species accumulation patterns at both spatial scales. Similarity in species composition (Chao-Jaccard) remained constant while similarity in dominance structure (Bray-Curtis) decreased with increased distance between samples across the range from < 1 to 13.8 km for both plant and bird communities. The similarity decay was more rapid in plants, but in both cases was slow. The results indicate low to zero beta-diversity at the spatial scale represented here, particularly for birds but also for woody plants. A 50 ha plot provided a highly accurate representation of broader-scale diversity and community composition within 10,000 ha for birds, and a relatively good representation for woody plants. This suggests potential for wider generalization of data from ForestGEO plots which are almost always locally unreplicated, at least for those in lowland tropical forest.

A new Cannabis genome assembly associates elevated cannabidiol (CBD) with hemp introgressed into marijuana.

Demand for cannabidiol (CBD), the predominant cannabinoid in hemp (Cannabis sativa), has favored cultivars producing unprecedented quantities of CBD. We investigated the ancestry of a new cultivar and cannabinoid synthase genes in relation to cannabinoid inheritance. A nanopore-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for CBDRx, a potent CBD-type cultivar. We measured cannabinoid synthase expression by cDNA sequencing and conducted a population genetic analysis of diverse Cannabis accessions. Quantitative trait locus mapping of cannabinoids in a hemp × marijuana segregating population was also performed. Cannabinoid synthase paralogs are arranged in tandem arrays embedded in long terminal repeat retrotransposons on chromosome 7. Although CBDRx is predominantly of marijuana ancestry, the genome has cannabidiolic acid synthase (CBDAS) introgressed from hemp and lacks a complete sequence for tetrahydrocannabinolic acid synthase (THCAS). Three additional genomes, including one with complete THCAS, confirmed this genomic structure. Only cannabidiolic acid synthase (CBDAS) was expressed in CBD-type Cannabis, while both CBDAS and THCAS were expressed in a cultivar with an intermediate tetrahydrocannabinol (THC) : CBD ratio. Although variation among cannabinoid synthase loci might affect the THC : CBD ratio, variability among cultivars in overall cannabinoid content (potency) was also associated with other chromosomes.

Validating a predictive model of cannabinoid inheritance with feral, clinical, and industrial Cannabis sativa.

PREMISE: How genetic variation within a species affects phytochemical composition is a fundamental question in botany. The ratio of two specialized metabolites in Cannabis sativa, tetrahydrocannabinol (THC) and cannabidiol (CBD), can be grouped into three main classes (THC-type, CBD-type, and intermediate type). We tested a genetic model associating these three groups with functional and nonfunctional alleles of the cannabidiolic acid synthase gene (CBDAS). METHODS: We characterized cannabinoid content and assayed CBDAS genotypes of >300 feral C. sativa plants in Minnesota, United States. We performed a test cross to assess CBDAS inheritance. Twenty clinical cultivars obtained blindly from the National Institute on Drug Abuse and 12 Canadian-certified grain cultivars were also examined. RESULTS: Frequencies of CBD-type, intermediate-type, and THC-type feral plants were 0.88, 0.11, and 0.01, respectively. Although total cannabinoid content varied substantially, the three groupings were perfectly correlated with CBDAS genotypes. Genotype frequencies observed in the test cross were consistent with codominant Mendelian inheritance of the THC:CBD ratio. Despite significant mean differences in total cannabinoid content, CBDAS genotypes blindly predicted the THC:CBD ratio among clinical cultivars, and the same was true for industrial grain cultivars when plants exhibited >0.5% total cannabinoid content. CONCLUSIONS: Our results extend the generality of the inheritance model for THC:CBD to diverse C. sativa accessions and demonstrate that CBDAS genotyping can predict the ratio in a variety of practical applications. Cannabinoid profiles and associated CBDAS segregation patterns suggest that feral C. sativa populations are potentially valuable experimental systems and sources of germplasm.

Spatial covariance of herbivorous and predatory guilds of forest canopy arthropods along a latitudinal gradient.

In arthropod community ecology, species richness studies tend to be prioritised over those investigating patterns of abundance. Consequently, the biotic and abiotic drivers of arboreal arthropod abundance are still relatively poorly known. In this cross-continental study, we employ a theoretical framework in order to examine patterns of covariance among herbivorous and predatory arthropod guilds. Leaf-chewing and leaf-mining herbivores, and predatory ants and spiders, were censused on > 1000 trees in nine 0.1 ha forest plots. After controlling for tree size and season, we found no negative pairwise correlations between guild abundances per plot, suggestive of weak signals of both inter-guild competition and top-down regulation of herbivores by predators. Inter-guild interaction strengths did not vary with mean annual temperature, thus opposing the hypothesis that biotic interactions intensify towards the equator. We find evidence for the bottom-up limitation of arthropod abundances via resources and abiotic factors, rather than for competition and predation.

Soil nitrogen concentration mediates the relationship between leguminous trees and neighbor diversity in tropical forests.

Legumes provide an essential service to ecosystems by capturing nitrogen from the atmosphere and delivering it to the soil, where it may then be available to other plants. However, this facilitation by legumes has not been widely studied in global tropical forests. Demographic data from 11 large forest plots (16-60 ha) ranging from 5.25° S to 29.25° N latitude show that within forests, leguminous trees have a larger effect on neighbor diversity than non-legumes. Where soil nitrogen is high, most legume species have higher neighbor diversity than non-legumes. Where soil nitrogen is low, most legumes have lower neighbor diversity than non-legumes. No facilitation effect on neighbor basal area was observed in either high or low soil N conditions. The legume-soil nitrogen positive feedback that promotes tree diversity has both theoretical implications for understanding species coexistence in diverse forests, and practical implications for the utilization of legumes in forest restoration.

Compound Specific Trends of Chemical Defences in Ficus Along an Elevational Gradient Reflect a Complex Selective Landscape.

Elevational gradients affect the production of plant secondary metabolites through changes in both biotic and abiotic conditions. Previous studies have suggested both elevational increases and decreases in host-plant chemical defences. We analysed the correlation of alkaloids and polyphenols with elevation in a community of nine Ficus species along a continuously forested elevational gradient in Papua New Guinea. We sampled 204 insect species feeding on the leaves of these hosts and correlated their community structure to the focal compounds. Additionally, we explored species richness of folivorous mammals along the gradient. When we accounted for Ficus species identity, we found a general elevational increase in flavonoids and alkaloids. Elevational trends in non-flavonol polyphenols were less pronounced or showed non-linear correlations with elevation. Polyphenols responded more strongly to changes in temperature and humidity than alkaloids. The abundance of insect herbivores decreased with elevation, while the species richness of folivorous mammals showed an elevational increase. Insect community structure was affected mainly by alkaloid concentration and diversity. Although our results show an elevational increase in several groups of metabolites, the drivers behind these trends likely differ. Flavonoids may provide figs with protection against abiotic stressors. In contrast, alkaloids affect insect herbivores and may provide protection against mammalian herbivores and pathogens. Concurrent analysis of multiple compound groups alongside ecological data is an important approach for understanding the selective landscape that shapes plant defences.

Quantitative assessment of plant-arthropod interactions in forest canopies: A plot-based approach.

Research on canopy arthropods has progressed from species inventories to the study of their interactions and networks, enhancing our understanding of how hyper-diverse communities are maintained. Previous studies often focused on sampling individual tree species, individual trees or their parts. We argue that such selective sampling is not ideal when analyzing interaction network structure, and may lead to erroneous conclusions. We developed practical and reproducible sampling guidelines for the plot-based analysis of arthropod interaction networks in forest canopies. Our sampling protocol focused on insect herbivores (leaf-chewing insect larvae, miners and gallers) and non-flying invertebrate predators (spiders and ants). We quantitatively sampled the focal arthropods from felled trees, or from trees accessed by canopy cranes or cherry pickers in 53 0.1 ha forest plots in five biogeographic regions, comprising 6,280 trees in total. All three methods required a similar sampling effort and provided good foliage accessibility. Furthermore, we compared interaction networks derived from plot-based data to interaction networks derived from simulated non-plot-based data focusing either on common tree species or a representative selection of tree families. All types of non-plot-based data showed highly biased network structure towards higher connectance, higher web asymmetry, and higher nestedness temperature when compared with plot-based data. Furthermore, some types of non-plot-based data showed biased diversity of the associated herbivore species and specificity of their interactions. Plot-based sampling thus appears to be the most rigorous approach for reconstructing realistic, quantitative plant-arthropod interaction networks that are comparable across sites and regions. Studies of plant interactions have greatly benefited from a plot-based approach and we argue that studies of arthropod interactions would benefit in the same way. We conclude that plot-based studies on canopy arthropods would yield important insights into the processes of interaction network assembly and dynamics, which could be maximised via a coordinated network of plot-based study sites.

Faster speciation of fig-wasps than their host figs leads to decoupled speciation dynamics: Snapshots across the speciation continuum.

Even though speciation involving multiple interacting partners, such as plants and their pollinators, has attracted much research, most studies focus on isolated phases of the process. This currently precludes an integrated understanding of the mechanisms leading to cospeciation. Here, we examine population genetic structure across six species-pairs of figs and their pollinating wasps along an elevational gradient in New Guinea. Specifically, we test three hypotheses on the genetic structure within the examined species-pairs and find that the hypothesized genetic structures represent different phases of a single continuum, from incipient cospeciation to the full formation of new species. Our results also illuminate the mechanisms governing cospeciation, namely that fig wasps tend to accumulate population genetic differences faster than their figs, which initially decouples the speciation dynamics between the two interacting partners and breaks down their one-to-one matching. This intermediate phase is followed by genetic divergence of both partners, which may eventually restore the one-to-one matching among the fully formed species. Together, these findings integrate current knowledge on the mechanisms operating during different phases of the cospeciation process. They also reveal that the increasingly reported breakdowns in one-to-one matching may be an inherent part of the cospeciation process. Mechanistic understanding of this process is needed to explain how the extraordinary diversity of species, especially in the tropics, has emerged. Knowing which breakdowns in species interactions are a natural phase of cospeciation and which may endanger further generation of diversity seems critical in a constantly changing world.

Response to Comment on "Plant diversity increases with the strength of negative density dependence at the global scale".

Hülsmann and Hartig suggest that ecological mechanisms other than specialized natural enemies or intraspecific competition contribute to our estimates of conspecific negative density dependence (CNDD). To address their concern, we show that our results are not the result of a methodological artifact and present a null-model analysis that demonstrates that our original findings-(i) stronger CNDD at tropical relative to temperate latitudes and (ii) a latitudinal shift in the relationship between CNDD and species abundance-persist even after controlling for other processes that might influence spatial relationships between adults and recruits.

Response to Comment on "Plant diversity increases with the strength of negative density dependence at the global scale".

Chisholm and Fung claim that our method of estimating conspecific negative density dependence (CNDD) in recruitment is systematically biased, and present an alternative method that shows no latitudinal pattern in CNDD. We demonstrate that their approach produces strongly biased estimates of CNDD, explaining why they do not detect a latitudinal pattern. We also address their methodological concerns using an alternative distance-weighted approach, which supports our original findings of a latitudinal gradient in CNDD and a latitudinal shift in the relationship between CNDD and species abundance.

Community structure of insect herbivores is driven by conservatism, escalation and divergence of defensive traits in Ficus.

Escalation (macroevolutionary increase) or divergence (disparity between relatives) in trait values are two frequent outcomes of the plant-herbivore arms race. We studied the defences and caterpillars associated with 21 sympatric New Guinean figs. Herbivore generalists were concentrated on hosts with low protease and oxidative activity. The distribution of specialists correlated with phylogeny, protease and trichomes. Additionally, highly specialised Asota moths used alkaloid rich plants. The evolution of proteases was conserved, alkaloid diversity has escalated across the studied species, oxidative activity has escalated within one clade, and trichomes have diverged across the phylogeny. Herbivore specificity correlated with their response to host defences: escalating traits largely affected generalists and divergent traits specialists; but the effect of escalating traits on extreme specialists was positive. In turn, the evolution of defences in Ficus can be driven towards both escalation and divergence in individual traits, in combination providing protection against a broad spectrum of herbivores.

Spatial scale changes the relationship between beta diversity, species richness and latitude.

The relationship between ß-diversity and latitude still remains to be a core question in ecology because of the lack of consensus between studies. One hypothesis for the lack of consensus between studies is that spatial scale changes the relationship between latitude and ß-diversity. Here, we test this hypothesis using tree data from 15 large-scale forest plots (greater than or equal to 15 ha, diameter at breast height ≥ 1 cm) across a latitudinal gradient (3-30o) in the Asia-Pacific region. We found that the observed ß-diversity decreased with increasing latitude when sampling local tree communities at small spatial scale (grain size ≤0.1 ha), but the observed ß-diversity did not change with latitude when sampling at large spatial scales (greater than or equal to 0.25 ha). Differences in latitudinal ß-diversity gradients across spatial scales were caused by pooled species richness (γ-diversity), which influenced observed ß-diversity values at small spatial scales, but not at large spatial scales. Therefore, spatial scale changes the relationship between ß-diversity, γ-diversity and latitude, and improving sample representativeness avoids the γ-dependence of ß-diversity.

Variably hungry caterpillars: predictive models and foliar chemistry suggest how to eat a rainforest.

A long-term goal in evolutionary ecology is to explain the incredible diversity of insect herbivores and patterns of host plant use in speciose groups like tropical Lepidoptera. Here, we used standardized food-web data, multigene phylogenies of both trophic levels and plant chemistry data to model interactions between Lepidoptera larvae (caterpillars) from two lineages (Geometridae and Pyraloidea) and plants in a species-rich lowland rainforest in New Guinea. Model parameters were used to make and test blind predictions for two hectares of an exhaustively sampled forest. For pyraloids, we relied on phylogeny alone and predicted 54% of species-level interactions, translating to 79% of all trophic links for individual insects, by sampling insects from only 15% of local woody plant diversity. The phylogenetic distribution of host-plant associations in polyphagous geometrids was less conserved, reducing accuracy. In a truly quantitative food web, only 40% of pair-wise interactions were described correctly in geometrids. Polyphenol oxidative activity (but not protein precipitation capacity) was important for understanding the occurrence of geometrids (but not pyraloids) across their hosts. When both foliar chemistry and plant phylogeny were included, we predicted geometrid-plant occurrence with 89% concordance. Such models help to test macroevolutionary hypotheses at the community level.

Phylogenetic trophic specialization: a robust comparison of herbivorous guilds.

Resource specialization is a key concept in ecology, but it is unexpectedly difficult to parameterize. Differences in resource availability, sampling effort and abundances preclude comparisons of incompletely sampled biotic interaction webs. Here, we extend the distance-based specialization index (DSI) that measures trophic specialization by taking resource phylogenetic relatedness and availability into account into a rescaled version, DSI*. It is a versatile metric of specialization that expands considerably the scope and applicability, hence the usefulness, of DSI. The new metric also accounts for differences in abundance and sampling effort of consumers, which enables robust comparisons among distinct guilds of consumers. It also provides an abundance threshold for the reliability of the metric for rare species, a very desirable property given the difficulty of assessing any aspect of rare species accurately. We apply DSI* to an extensive dataset on interactions between insect herbivores from four folivorous guilds and their host plants in Papua New Guinean rainforests. We demonstrate that DSI*, contrary to the original DSI, is largely independent of sample size and weakly and non-linearly related with several host specificity measures that do not adjust for plant phylogeny. Thus, DSI* provides further insights into host specificity patterns; moreover, it is robust to the number and phylogenetic diversity of plant species selected to be sampled for herbivores. DSI* can be used for a broad range of comparisons of distinct feeding guilds, geographical locations and ecological conditions. This is a key advance in elucidating the interaction structure and evolution of highly diversified systems.

Plant diversity increases with the strength of negative density dependence at the global scale.

Theory predicts that higher biodiversity in the tropics is maintained by specialized interactions among plants and their natural enemies that result in conspecific negative density dependence (CNDD). By using more than 3000 species and nearly 2.4 million trees across 24 forest plots worldwide, we show that global patterns in tree species diversity reflect not only stronger CNDD at tropical versus temperate latitudes but also a latitudinal shift in the relationship between CNDD and species abundance. CNDD was stronger for rare species at tropical versus temperate latitudes, potentially causing the persistence of greater numbers of rare species in the tropics. Our study reveals fundamental differences in the nature of local-scale biotic interactions that contribute to the maintenance of species diversity across temperate and tropical communities.

Trophic phylogenetics: evolutionary influences on body size, feeding, and species associations in grassland arthropods.

Contemporary animal-plant interactions such as herbivory are widely understood to be shaped by evolutionary history. Yet questions remain about the role of plant phylogenetic diversity in generating and maintaining herbivore diversity, and whether evolutionary relatedness of producers might predict the composition of consumer communities. We tested for evidence of evolutionary associations among arthropods and the plants on which they were found, using phylogenetic analysis of naturally occurring arthropod assemblages sampled from a plant-diversity manipulation experiment. Considering phylogenetic relationships among more than 900 arthropod consumer taxa and 29 plant species in the experiment, we addressed several interrelated questions. First, our results support the hypothesis that arthropod functional traits such as body size and trophic role are phylogenetically conserved in community ecological samples. Second, herbivores tended to cooccur with closer phylogenetic relatives than would be expected at random, whereas predators and parasitoids did not show phylogenetic association patterns. Consumer specialization, as measured by association through time with monocultures of particular host plant species, showed significant phylogenetic signal, although the. strength of this association varied among plant species. Polycultures of phylogenetically dissimilar plant species supported more phylogenetically dissimilar consumer communities than did phylogenetically similar polycultures. Finally, we separated the effects of plant species richness and relatedness in predicting the phylogenetic distribution of the arthropod assemblages in this experiment. The phylogenetic diversity of plant communities predicted the phylogenetic diversity of herbivore communities even after accounting for plant species richness. The phylogenetic diversity of secondary consumers differed by guild, with predator phylogenetic diversity responding to herbivore relatedness, while parasitoid phylogenetic diversity was driven by plant relatedness. Evolutionary associations between plants and their consumers are apparent in plots only meters apart in a single field, indicating a strong role for host-plant phylogenetic diversity in sustaining landscape consumer biodiversity.