Mycorrhizal feedbacks influence global forest structure and diversity

One mechanism proposed to explain high species diversity in tropical systems is strong negative conspecific density dependence (CDD), which reduces recruitment of juveniles in proximity to conspecific adult plants. Although evidence shows that plant-specific soil pathogens can drive negative CDD, trees also form key mutualisms with mycorrhizal fungi, which may counteract these effects. Across 43 large-scale forest plots worldwide, we tested whether ectomycorrhizal tree species exhibit weaker negative CDD than arbuscular mycorrhizal tree species. We further tested for conmycorrhizal density dependence (CMDD) to test for benefit from shared mutualists. We found that the strength of CDD varies systematically with mycorrhizal type, with ectomycorrhizal tree species exhibiting higher sapling densities with increasing adult densities than arbuscular mycorrhizal tree species. Moreover, we found evidence of positive CMDD for tree species of both mycorrhizal types. Collectively, these findings indicate that mycorrhizal interactions likely play a foundational role in global forest diversity patterns and structure.

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.

Distribution of biomass dynamics in relation to tree size in forests across the world

Tree size shapes forest carbon dynamics and determines how trees interact with their environment, including a changing climate. Here, we conduct the first global analysis of among-site differences in how aboveground biomass stocks and fluxes are distributed with tree size. We analyzed repeat tree censuses from 25 large-scale (4-52 ha) forest plots spanning a broad climatic range over five continents to characterize how aboveground biomass, woody productivity, and woody mortality vary with tree diameter. We examined how the median, dispersion, and skewness of these size-related distributions vary with mean annual temperature and precipitation. In warmer forests, aboveground biomass, woody productivity, and woody mortality were more broadly distributed with respect to tree size. In warmer and wetter forests, aboveground biomass and woody productivity were more right skewed, with a long tail towards large trees. Small trees (1-10 cm diameter) contributed more to productivity and mortality than to biomass, highlighting the importance of including these trees in analyses of forest dynamics. Our findings provide an improved characterization of climate-driven forest differences in the size structure of aboveground biomass and dynamics of that biomass, as well as refined benchmarks for capturing climate influences in vegetation demographic models.

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.

Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide

Arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) associations are critical for host-tree performance. However, how mycorrhizal associations correlate with the latitudinal tree beta-diversity remains untested. Using a global dataset of 45 forest plots representing 2,804,270 trees across 3840 species, we test how AM and EcM trees contribute to total beta-diversity and its components (turnover and nestedness) of all trees. We find AM rather than EcM trees predominantly contribute to decreasing total beta-diversity and turnover and increasing nestedness with increasing latitude, probably because wide distributions of EcM trees do not generate strong compositional differences among localities. Environmental variables, especially temperature and precipitation, are strongly correlated with beta-diversity patterns for both AM trees and all trees rather than EcM trees. Results support our hypotheses that latitudinal beta-diversity patterns and environmental effects on these patterns are highly dependent on mycorrhizal types. Our findings highlight the importance of AM-dominated forests for conserving global forest biodiversity.

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

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.

Direct and indirect effects of climate on richness drive the latitudinal diversity gradient in forest trees

Climate is widely recognised as an important determinant of the latitudinal diversity gradient. However, most existing studies make no distinction between direct and indirect effects of climate, which substantially hinders our understanding of how climate constrains biodiversity globally. Using data from 35 large forest plots, we test hypothesised relationships amongst climate, topography, forest structural attributes (stem abundance, tree size variation and stand basal area) and tree species richness to better understand drivers of latitudinal tree diversity patterns. Climate influences tree richness both directly, with more species in warm, moist, aseasonal climates and indirectly, with more species at higher stem abundance. These results imply direct limitation of species diversity by climatic stress and more rapid (co-)evolution and narrower niche partitioning in warm climates. They also support the idea that increased numbers of individuals associated with high primary productivity are partitioned to support a greater number of species.

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-30^o) 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.

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.

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.

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.

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.

Host associations and beta diversity of fungal endophyte communities in New Guinea rainforest trees

Processes shaping the distribution of foliar fungal endophyte species remain poorly understood. Despite increasing evidence that these cryptic fungal symbionts of plants mediate interactions with pathogens and herbivores, there remain basic questions regarding the extent to which dispersal limitation and host specificity might shape fungal endophyte community composition in rainforests. To assess the relative importance of spatial pattern and host specificity, we isolated fungi from a sample of mapped trees in lowland Papua New Guinea. Sequences of the internal transcribed spacer (ITS) region were obtained for 2079 fungal endophytes from three sites and clustered into molecular operational taxonomic units (MOTUs) at 95% similarity. Multivariate analyses suggest that host affinity plays a significant role in structuring endophyte community composition whereas there was no evidence of endophyte spatial pattern at the scale of tens to hundreds of metres. Differences in endophyte communities between sampled trees were weakly correlated with variation in foliar traits but not with tree species relatedness. The dominance of relatively few generalist endophytes and the presence of a large number of rare MOTUs was a consistent observation at three sites separated by hundreds of kilometres and regional turnover was low. Host specificity appears to play a relatively weak but more important role than dispersal limitation in shaping the distribution of fungal endophyte communities in New Guinea forests. Our results suggest that in the absence of strong ecological gradients and host turnover, beta diversity of endophyte communities could be low in large areas of contiguous forest.

Gene duplication and divergence affecting drug content in Cannabis sativa

Cannabis sativa is an economically important source of durable fibers, nutritious seeds, and psychoactive drugs but few economic plants are so poorly understood genetically. Marijuana and hemp were crossed to evaluate competing models of cannabinoid inheritance and to explain the predominance of tetrahydrocannabinolic acid (THCA) in marijuana compared with cannabidiolic acid (CBDA) in hemp. Individuals in the resulting F2 population were assessed for differential expression of cannabinoid synthase genes and were used in linkage mapping. Genetic markers associated with divergent cannabinoid phenotypes were identified. Although phenotypic segregation and a major quantitative trait locus (QTL) for the THCA/CBDA ratio were consistent with a simple model of codominant alleles at a single locus, the diversity of THCA and CBDA synthase sequences observed in the mapping population, the position of enzyme coding loci on the map, and patterns of expression suggest multiple linked loci. Phylogenetic analysis further suggests a history of duplication and divergence affecting drug content. Marijuana is distinguished from hemp by a nonfunctional CBDA synthase that appears to have been positively selected to enhance psychoactivity. An unlinked QTL for cannabinoid quantity may also have played a role in the recent escalation of drug potency.

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.

Forest carbon in lowland Papua New Guinea: Local variation and the importance of small trees

Efforts to incentivize the reduction of carbon emissions from deforestation and forest degradation require accurate carbon accounting. The extensive tropical forest of Papua New Guinea (PNG) is a target for such efforts and yet local carbon estimates are few. Previous estimates, based on models of neotropical vegetation applied to PNG forest plots, did not consider such factors as the unique species composition of New Guinea vegetation, local variation in forest biomass, or the contribution of small trees. We analysed all trees >1 cm in diameter at breast height (DBH) in Melanesia's largest forest plot (Wanang) to assess local spatial variation and the role of small trees in carbon storage. Above-ground living biomass (AGLB) of trees averaged 210.72 Mg  ha⁻¹ at Wanang. Carbon storage at Wanang was somewhat lower than in other lowland tropical forests, whereas local variation among 1-ha subplots and the contribution of small trees to total AGLB were substantially higher. We speculate that these differences may be attributed to the dynamics of Wanang forest where erosion of a recently uplifted and unstable terrain appears to be a major source of natural disturbance. These findings emphasize the need for locally calibrated forest carbon estimates if accurate landscape level valuation and monetization of carbon is to be achieved. Such estimates aim to situate PNG forests in the global carbon context and provide baseline information needed to improve the accuracy of PNG carbon monitoring schemes.

Phylogenetic signal variation in the genomes of Medicago (Fabaceae)

Genome-scale data offer the opportunity to clarify phylogenetic relationships that are difficult to resolve with few loci, but they can also identify genomic regions with evolutionary history distinct from that of the species history. We collected whole-genome sequence data from 29 taxa in the legume genus Medicago, then aligned these sequences to the Medicago truncatula reference genome to confidently identify 87 596 variable homologous sites. We used this data set to estimate phylogenetic relationships among Medicago species, to investigate the number of sites needed to provide robust phylogenetic estimates and to identify specific genomic regions supporting topologies in conflict with the genome-wide phylogeny. Our full genomic data set resolves relationships within the genus that were previously intractable. Subsampling the data reveals considerable variation in phylogenetic signal and power in smaller subsets of the data. Even when sampling 5000 sites, no random sample of the data supports a topology identical to that of the genome-wide phylogeny. Phylogenetic relationships estimated from 500-site sliding windows revealed genome regions supporting several alternative species relationships among recently diverged taxa, consistent with the expected effects of deep coalescence or introgression in the recent history of Medicago.

An extreme case of plant-insect codiversification: figs and fig-pollinating wasps

It is thought that speciation in phytophagous insects is often due to colonization of novel host plants, because radiations of plant and insect lineages are typically asynchronous. Recent phylogenetic comparisons have supported this model of diversification for both insect herbivores and specialized pollinators. An exceptional case where contemporaneous plant-insect diversification might be expected is the obligate mutualism between fig trees (Ficus species, Moraceae) and their pollinating wasps (Agaonidae, Hymenoptera). The ubiquity and ecological significance of this mutualism in tropical and subtropical ecosystems has long intrigued biologists, but the systematic challenge posed by >750 interacting species pairs has hindered progress toward understanding its evolutionary history. In particular, taxon sampling and analytical tools have been insufficient for large-scale cophylogenetic analyses. Here, we sampled nearly 200 interacting pairs of fig and wasp species from across the globe. Two supermatrices were assembled: on an average, wasps had sequences from 77% of 6 genes (5.6 kb), figs had sequences from 60% of 5 genes (5.5 kb), and overall 850 new DNA sequences were generated for this study. We also developed a new analytical tool, Jane 2, for event-based phylogenetic reconciliation analysis of very large data sets. Separate Bayesian phylogenetic analyses for figs and fig wasps under relaxed molecular clock assumptions indicate Cretaceous diversification of crown groups and contemporaneous divergence for nearly half of all fig and pollinator lineages. Event-based cophylogenetic analyses further support the codiversification hypothesis. Biogeographic analyses indicate that the present-day distribution of fig and pollinator lineages is consistent with a Eurasian origin and subsequent dispersal, rather than with Gondwanan vicariance. Overall, our findings indicate that the fig-pollinator mutualism represents an extreme case among plant-insect interactions of coordinated dispersal and long-term codiversification. [Biogeography; coevolution; cospeciation; host switching; long-branch attraction; phylogeny.].

Estimating global arthropod species richness: refining probabilistic models using probability bounds analysis

A key challenge in the estimation of tropical arthropod species richness is the appropriate management of the large uncertainties associated with any model. Such uncertainties had largely been ignored until recently, when we attempted to account for uncertainty associated with model variables, using Monte Carlo analysis. This model is restricted by various assumptions. Here, we use a technique known as probability bounds analysis to assess the influence of assumptions about (1) distributional form and (2) dependencies between variables, and to construct probability bounds around the original model prediction distribution. The original Monte Carlo model yielded a median estimate of 6.1 million species, with a 90 % confidence interval of [3.6, 11.4]. Here we found that the probability bounds (p-bounds) surrounding this cumulative distribution were very broad, owing to uncertainties in distributional form and dependencies between variables. Replacing the implicit assumption of pure statistical independence between variables in the model with no dependency assumptions resulted in lower and upper p-bounds at 0.5 cumulative probability (i.e., at the median estimate) of 2.9-12.7 million. From here, replacing probability distributions with probability boxes, which represent classes of distributions, led to even wider bounds (2.4-20.0 million at 0.5 cumulative probability). Even the 100th percentile of the uppermost bound produced (i.e., the absolutely most conservative scenario) did not encompass the well-known hyper-estimate of 30 million species of tropical arthropods. This supports the lower estimates made by several authors over the last two decades.

Pollinator-mediated reproductive isolation among dioecious fig species (Ficus, Moraceae)

The extent of isolation among closely related sympatric plant species engaged in obligate pollination mutualisms depends on the fitness consequences of interspecies floral visitation. In figs (Ficus), interspecific gene flow may occur when pollinating wasps (Agaonidae) visit species other than their natal fig species. We studied reproductive isolation in a clade of six sympatric dioecious fig species in New Guinea. Microsatellite genotyping and Bayesian clustering analysis of the fig community indicated strong reproductive barriers among sympatric species. A total of 1-2% of fig populations consisted of hybrid individuals. A new experimental method of manipulating fig wasps investigated the reproductive consequences of conspecific and heterospecific pollinator visitation for both mutualists. Fig wasps introduced to Ficus hispidioides pollinated and oviposited in receptive figs. Seed development and seedling growth were largely comparable between conspecific and heterospecific crosses. Heterospecific pollinator fitness, however, was significantly less than that of conspecific pollinators. Heterospecific pollinators induced gall formation but offspring did not develop to maturity in the new host. Selection on pollinators maintaining host specificity appears to be an important mechanism of contemporary reproductive isolation among these taxa that could potentially influence their diversification.

Why are there more arboreal ant species in primary than in secondary tropical forests?

1. Species diversity of arboreal arthropods tends to increase during rainforest succession so that primary forest communities comprise more species than those from secondary vegetation, but it is not well understood why. Primary forests differ from secondary forests in a wide array of factors whose relative impacts on arthropod diversity have not yet been quantified. 2. We assessed the effects of succession-related determinants on a keystone ecological group, arboreal ants, by conducting a complete census of 1332 ant nests from all trees with diameter at breast height ≥ 5 cm occurring within two (unreplicated) 0·32-ha plots, one in primary and one in secondary lowland forest in New Guinea. Specifically, we used a novel rarefaction-based approach to match number, size distribution and taxonomic structure of trees in primary forest communities to those in secondary forest and compared the resulting numbers of ant species. 3. In total, we recorded 80 nesting ant species from 389 trees in primary forest but only 42 species from 295 trees in secondary forest. The two habitats did not differ in the mean number of ant species per tree or in the relationship between ant diversity and tree size. However, the between-tree similarity of ant communities was higher in secondary forest than in primary forest, as was the between-tree nest site similarity, suggesting that secondary trees were more uniform in providing nesting microhabitats. 4. Using our rarefaction method, the difference in ant species richness between two forest types was partitioned according to the effects of higher tree density (22·6%), larger tree size (15·5%) and higher taxonomic diversity of trees (14·3%) in primary than in secondary forest. The remaining difference (47·6%) was because of higher beta diversity of ant communities between primary forest trees. In contrast, difference in nest density was explained solely by difference in tree density. 5. Our study shows that reduction in plant taxonomic diversity in secondary forests is not the main driver of the reduction in canopy ant species richness. We suggest that the majority of arboreal species losses in secondary tropical forests are attributable to simpler vegetation structure, combined with lower turnover of nesting microhabitats between trees.

Insects on plants: explaining the paradox of low diversity within specialist herbivore guilds

Classical niche theory explains the coexistence of species through their exploitation of different resources. Assemblages of herbivores coexisting on a particular plant species are thus expected to be dominated by species from host-specific guilds with narrow, coexistence-facilitating niches rather than by species from generalist guilds. Exactly the opposite pattern is observed for folivores feeding on trees in New Guinea. The least specialized mobile chewers were the most species rich, followed by the moderately specialized semiconcealed and exposed chewers. The highly specialized miners and mesophyll suckers were the least species-rich guilds. The Poisson distribution of herbivore species richness among plant species in specialized guilds and the absence of a negative correlation between species richness in different guilds on the same plant species suggest that these guilds are not saturated with species. We show that herbivore assemblages are enriched with generalists because these are more completely sampled from regional species pools. Herbivore diversity increases as a power function of plant diversity, and the rate of increase is inversely related to host specificity. The relative species diversity among guilds is thus scale dependent, as the importance of specialized guilds increases with plant diversity. Specialized insect guilds may therefore comprise a larger component of overall diversity in the tropics (where they are also poorly known taxonomically) than in the temperate zone, which has lower plant diversity.

Guild-specific patterns of species richness and host specialization in plant-herbivore food webs from a tropical forest

1. The extent to which plant-herbivore feeding interactions are specialized is key to understand the processes maintaining the diversity of both tropical forest plants and their insect herbivores. However, studies documenting the full complexity of tropical plant-herbivore food webs are lacking. 2. We describe a complex, species-rich plant-herbivore food web for lowland rain forest in Papua New Guinea, resolving 6818 feeding links between 224 plant species and 1490 herbivore species drawn from 11 distinct feeding guilds. By standardizing sampling intensity and the phylogenetic diversity of focal plants, we are able to make the first rigorous and unbiased comparisons of specificity patterns across feeding guilds. 3. Specificity was highly variable among guilds, spanning almost the full range of theoretically possible values from extreme trophic generalization to monophagy. 4. We identify guilds of herbivores that are most likely to influence the composition of tropical forest vegetation through density-dependent herbivory or apparent competition. 5. We calculate that 251 herbivore species (48 of them unique) are associated with each rain forest tree species in our study site so that the ∼200 tree species coexisting in the lowland rain forest community are involved in ∼50,000 trophic interactions with ∼9600 herbivore species of insects. This is the first estimate of total herbivore and interaction number in a rain forest plant-herbivore food web. 6. A comprehensive classification of insect herbivores into 24 guilds is proposed, providing a framework for comparative analyses across ecosystems and geographical regions.

Development and characterization of microsatellite loci in dioecious figs (Ficus, Moraceae)

PREMISE OF THE STUDY

Microsatellite markers for dioecious fig species, Ficus bernaysii King (section Sycocarpus), Ficus pachyrrhachis K. Schum. & Lauterb. (section Sycocarpus), and Ficus copiosa Steud. (section Sycidium) were developed to estimate gene flow among closely related fig species.

METHODS AND RESULTS

Loci were developed from tri- and tetranucleotide-enriched genomic libraries. The six most repeatable and scorable loci for F. bernaysii and F. pachyrrhachis were characterized in 50 individuals each of five species from section Sycocarpus. Four primer pairs from F. copiosa were screened in 50 individuals each of two species from section Sycidium. Observed heterozygosity ranged from 0.53 to 0.93 in loci developed from F. bernaysii, 0.25 to 0.86 in loci developed from F. pachyrrhachis, and 0.68 to 0.87 in loci developed from F. copiosa.

CONCLUSIONS

Amplification across the two Ficus sections was observed in seven of the ten total primer pairs.

Quantifying uncertainty in estimation of tropical arthropod species richness

There is a bewildering range of estimates for the number of arthropods on Earth. Several measures are based on extrapolation from species specialized to tropical rain forest, each using specific assumptions and justifications. These approaches have not provided any sound measure of uncertainty associated with richness estimates. We present two models that account for parameter uncertainty by replacing point estimates with probability distributions. The models predict medians of 3.7 million and 2.5 million tropical arthropod species globally, with 90% confidence intervals of [2.0, 7.4] million and [1.1, 5.4] million, respectively. Estimates of 30 million or greater are predicted to have <0.00001 probability. Sensitivity analyses identified uncertainty in the proportion of canopy arthropod species that are beetles as the most influential parameter, although uncertainties associated with three other parameters were also important. Using the median estimates suggests that in spite of 250 years of taxonomy and around 855,000 species of arthropods already described, approximately 70% await description.

60 million years of co-divergence in the fig-wasp symbiosis

Figs (Ficus; ca 750 species) and fig wasps (Agaoninae) are obligate mutualists: all figs are pollinated by agaonines that feed exclusively on figs. This extraordinary symbiosis is the most extreme example of specialization in a plant-pollinator interaction and has fuelled much speculation about co-divergence. The hypothesis that pollinator specialization led to the parallel diversification of fig and pollinator lineages (co-divergence) has so far not been tested due to the lack of robust and comprehensive phylogenetic hypotheses for both partners. We produced and combined the most comprehensive molecular phylogenetic trees to date with fossil data to generate independent age estimates for fig and pollinator lineages, using both non-parametric rate smoothing and penalized likelihood dating methods. Molecular dating of ten pairs of interacting lineages provides an unparalleled example of plant-insect co-divergence over a geological time frame spanning at least 60 million years.

Molecular dating and biogeography of fig-pollinating wasps

Figs and fig-pollinating wasps are obligate mutualists that have coevolved for over 60 million years. But when and where did pollinating fig wasps (Agaonidae) originate? Some studies suggest that agaonids arose in the Late Cretaceous and the current distribution of fig-wasp faunas can be explained by the break-up of the Gondwanan landmass. However, recent molecular-dating studies suggest divergence time estimates that are inconsistent with the Gondwanan vicariance hypothesis and imply that long distance oceanic dispersal could have been an important process for explaining the current distribution of both figs and fig wasps. Here, we use a combination of phylogenetic and biogeographical data to infer the age, the major period of diversification, and the geographic origin of pollinating fig wasps. Age estimates ranged widely depending on the molecular-dating method used and even when using the same method but with slightly different constraints, making it difficult to assess with certainty a Gondwanan origin of agaonids. The reconstruction of ancestral areas suggests that the most recent common ancestor of all extant fig-pollinating wasps was most likely Asian, although a southern Gondwana origin cannot be rejected. Our analysis also suggests that dispersal has played a more important role in the development of the fig-wasp biota than previously assumed.

Identification of candidate genes affecting Delta9-tetrahydrocannabinol biosynthesis in Cannabis sativa

RNA isolated from the glands of a Delta(9)-tetrahydrocannabinolic acid (THCA)-producing strain of Cannabis sativa was used to generate a cDNA library containing over 100 000 expressed sequence tags (ESTs). Sequencing of over 2000 clones from the library resulted in the identification of over 1000 unigenes. Candidate genes for almost every step in the biochemical pathways leading from primary metabolites to THCA were identified. Quantitative PCR analysis suggested that many of the pathway genes are preferentially expressed in the glands. Hexanoyl-CoA, one of the metabolites required for THCA synthesis, could be made via either de novo fatty acids synthesis or via the breakdown of existing lipids. qPCR analysis supported the de novo pathway. Many of the ESTs encode transcription factors and two putative MYB genes were identified that were preferentially expressed in glands. Given the similarity of the Cannabis MYB genes to those in other species with known functions, these Cannabis MYBs may play roles in regulating gland development and THCA synthesis. Three candidates for the polyketide synthase (PKS) gene responsible for the first committed step in the pathway to THCA were characterized in more detail. One of these was identical to a previously reported chalcone synthase (CHS) and was found to have CHS activity. All three could use malonyl-CoA and hexanoyl-CoA as substrates, including the CHS, but reaction conditions were not identified that allowed for the production of olivetolic acid (the proposed product of the PKS activity needed for THCA synthesis). One of the PKS candidates was highly and specifically expressed in glands (relative to whole leaves) and, on the basis of these expression data, it is proposed to be the most likely PKS responsible for olivetolic acid synthesis in Cannabis glands.

Low beta diversity of herbivorous insects in tropical forests

Recent advances in understanding insect communities in tropical forests have contributed little to our knowledge of large-scale patterns of insect diversity, because incomplete taxonomic knowledge of many tropical species hinders the mapping of their distribution records. This impedes an understanding of global biodiversity patterns and explains why tropical insects are under-represented in conservation biology. Our study of approximately 500 species from three herbivorous guilds feeding on foliage (caterpillars, Lepidoptera), wood (ambrosia beetles, Coleoptera) and fruit (fruitflies, Diptera) found a low rate of change in species composition (beta diversity) across 75,000 square kilometres of contiguous lowland rainforest in Papua New Guinea, as most species were widely distributed. For caterpillars feeding on large plant genera, most species fed on multiple host species, so that even locally restricted plant species did not support endemic herbivores. Large plant genera represented a continuously distributed resource easily colonized by moths and butterflies over hundreds of kilometres. Low beta diversity was also documented in groups with differing host specificity (fruitflies and ambrosia beetles), suggesting that dispersal limitation does not have a substantial role in shaping the distribution of insect species in New Guinea lowland rainforests. Similar patterns of low beta diversity can be expected in other tropical lowland rainforests, as they are typically situated in the extensive low basins of major tropical rivers similar to the Sepik-Ramu region of New Guinea studied here.

Phylogenetic dispersion of host use in a tropical insect herbivore community

Theory has long predicted that insect community structure should be related to host plant phylogeny. We examined the distribution of insect herbivore associations with respect to host plant phylogeny for caterpillars (Lepidoptera), beetles (Coleoptera), and grasshoppers and relatives (orthopteroids) in a New Guinea rain forest. We collected herbivores from three lineages of closely related woody plants and from more distantly related plant lineages in the same locality to examine the phylogenetic scale at which host specificity can be detected in a community sample. By grafting molecular phylogenies inferred from three different genes into a supertree, we developed a phylogenetic hypothesis for the host community. Feeding experiments were performed on more than 100 000 live insects collected from the 62 host species. We examined patterns of host use with respect to the host plant phylogeny. As predicted, we found a negative relationship between faunal similarity, defined as the proportion of all herbivores feeding on two hosts that are shared between the hosts, and the phylogenetic distance between hosts based on DNA sequence divergence. Host phylogenetic distance explained a significant fraction of the variance (25%) in herbivore community similarity, in spite of the many ecological factors that probably influence feeding patterns. Herbivore community similarity among congeneric hosts was high (50% on average) compared to overlap among host families (20-30% on average). We confirmed this pattern using the nearest taxon index (NTI) and net relatedness index (NRI) to quantify the extent of phylogenetic clustering in particular herbivore associations and to test whether patterns are significantly different from chance expectations. We found that 40% of caterpillar species showed significant phylogenetic clustering with respect to host plant associations, somewhat more so than for beetles or orthopteroids. We interpret this as evidence that a substantial fraction of tropical forest insect herbivores are clade specialists.

Why Are There So Many Species of Herbivorous Insects in Tropical Rainforests?

Despite recent progress in understanding mechanisms of tree species coexistence in tropical forests, a simple explanation for the even more extensive diversity of insects feeding on these plants has been missing. We compared folivorous insects from temperate and tropical trees to test the hypothesis that herbivore species coexistence in more diverse communities could reflect narrow host specificity relative to less diverse communities. Temperate and tropical tree species of comparable phylogenetic distribution supported similar numbers of folivorous insect species, 29.0 +/- 2.2 and 23.5 +/- 1.8 per 100 square meters of foliage, respectively. Host specificity did not differ significantly between community samples, indicating that food resources are not more finely partitioned among folivorous insects in tropical than in temperate forests. These findings suggest that the latitudinal gradient in insect species richness could be a direct function of plant diversity, which increased sevenfold from our temperate to tropical study sites.