Limited intraspecific variation in drought resistance along a pronounced tropical rainfall gradient.

Assessing within-species variation in response to drought is crucial for predicting species' responses to climate change and informing restoration and conservation efforts, yet experimental data are lacking for the vast majority of tropical tree species. We assessed intraspecific variation in response to water availability across a strong rainfall gradient for 16 tropical tree species using reciprocal transplant and common garden field experiments, along with measurements of gene flow and key functional traits linked to drought resistance. Although drought resistance varies widely among species in these forests, we found little evidence for within-species variation in drought resistance. For the majority of functional traits measured, we detected no significant intraspecific variation. The few traits that did vary significantly between drier and wetter origins of the same species all showed relationships opposite to expectations based on drought stress. Furthermore, seedlings of the same species originating from drier and wetter sites performed equally well under drought conditions in the common garden experiment and at the driest transplant site. However, contrary to expectation, wetter-origin seedlings survived better than drier-origin seedlings under wetter conditions in both the reciprocal transplant and common garden experiment, potentially due to lower insect herbivory. Our study provides the most comprehensive picture to date of intraspecific variation in tropical tree species' responses to water availability. Our findings suggest that while drought plays an important role in shaping species composition across moist tropical forests, its influence on within-species variation is limited.

Cryptic species diversity in a widespread neotropical tree genus: The case of Cedrela odorata.

PREMISE: Reconciling the use of taxonomy to partition morphological variation and describe genetic divergence within and among closely related species is a persistent challenge in phylogenetics. We reconstructed phylogenetic relationships among Cedrela odorata (Meliaceae) and five closely allied species to test the genetic basis for the current model of species delimitation in this economically valuable and threatened genus. METHODS: We prepared a nuclear species tree with the program SNPhylo and 16,000 single-nucleotide polymorphisms from 168 Cedrela specimens. Based on clades present and ancestral patterns ADMIXTURE, we designed nine species delimitation models and compared each model to current taxonomy with Bayes factor delimitation. Timing of major lineage divergences was estimated with the program SNAPP. RESULTS: The resulting analysis revealed that modern C. odorata evolved from two genetically distinct ancestral sources. All species delimitation models tested better fit the data than the model representing current taxonomic delimitation. Models with the greatest marginal likelihoods separated Mesoamerican C. odorata and South American C. odorata into two species and lumped C. angustifolia and C. montana as a single species. We estimated that Cedrela diversified in South America within the last 19 million years following one or more dispersal events from Mesoamerican lineages. CONCLUSIONS: Our analyses show that the present taxonomic understanding within the genus obscures divergent lineages in C. odorata due in part to morphological differentiation and taxonomic distinctions that are not predictably associated with genetic divergence. A more accurate application of taxonomy to C. odorata and related species may aid in its conservation, management, and restoration efforts.

Elevated inbreeding in Heliconia tortuosa is determined by tropical forest stand age, isolation and loss of hummingbird functional diversity.

Forest conversion and habitat loss are major threats to biological diversity. Forest regeneration can mitigate the negative effects of old-growth forest loss on species diversity, but less is known about the extent to which forest loss reduces genetic diversity in remnant populations and whether secondary forests play a role in the maintenance of genetic diversity. We quantified genetic diversity in a tropical hummingbird-pollinated understorey herb, Heliconia tortuosa, across a landscape mosaic of primary and secondary forest regrowth. Using microsatellite genotypes from >850 adult and juvenile plants within 33 forest patches and extensive bird surveys, we examined the effect of contemporary and historical landscape features including forest age (primary vs. secondary forest), stand isolation and pollinator assemblages on genetic diversity and levels of inbreeding in H. tortuosa. We found that inbreeding was up to three times higher in secondary forest, and this effect was amplified with reductions in primary forest in the surrounding landscape through reduced observed heterozygosity in isolated fragments. Inbreeding in forest patches was negatively correlated with the local frequency of specialist long-distance foraging traplining hummingbirds. Traplining hummingbirds therefore appear to facilitate mating among unrelated plants-an inference we tested using empirically parameterized simulations. Higher levels of inbreeding in H. tortuosa are therefore associated with reduced functional diversity of hummingbirds in secondary forests and forest patches isolated from primary forests. Our findings suggest a cryptic consequence of primary forest loss and secondary forest regeneration through the disruption of mutualistic interactions resulting in the erosion of genetic diversity in a common understorey plant.

Turgor loss point predicts survival responses to experimental and natural drought in tropical tree seedlings.

Identifying key traits that can serve as proxies for species drought resistance is crucial for predicting and mitigating the effects of climate change in diverse plant communities. Turgor loss point (πtlp ) is a recently emerged trait that has been linked to species distributions across gradients of water availability. However, a direct relationship between πtlp and species ability to survive drought has yet to be established for woody species. Using a manipulative field experiment to quantify species drought resistance (i.e., their survival response to drought), combined with measurements of πtlp for 16 tree species, we show a negative relationship between πtlp and seedling drought resistance. Using long-term forest plot data, we also show that πtlp predicts seedling survival responses to a severe El Niño-related drought, although additional factors are clearly also important. Our study demonstrates that species with lower πtlp exhibit higher survival under both experimental and natural drought. These results provide a missing cornerstone in the assessment of the traits underlying drought resistance in woody species and strengthen πtlp as a proxy for evaluating which species will lose or win under projections of exacerbating drought regimes.

Tree species diversity increases with conspecific negative density dependence across an elevation gradient.

Elevational and latitudinal gradients in species diversity may be mediated by biotic interactions that cause density-dependent effects of conspecifics on survival or growth to differ from effects of heterospecifics (i.e. conspecific density dependence), but limited evidence exists to support this. We tested the hypothesis that conspecific density dependence varies with elevation using over 40 years of data on tree survival and growth from 23 old-growth temperate forest stands across a 1,000-m elevation gradient. We found that conspecific-density-dependent effects on survival of small-to-intermediate-sized focal trees were negative in lower elevation, higher diversity forest stands typically characterised by warmer temperatures and greater relative humidity. Conspecific-density-dependent effects on survival were less negative in higher elevation stands and ridges than in lower elevation stands and valley bottoms for small-to-intermediate-sized trees, but were neutral for larger trees across elevations. Conspecific-density-dependent effects on growth were negative across all tree size classes and elevations. These findings reveal fundamental differences in biotic interactions that may contribute to relationships between species diversity, elevation and climate.

Increased mortality of tropical tree seedlings during the extreme 2015-16 El Niño.

As extreme climate events are predicted to become more frequent because of global climate change, understanding their impacts on natural systems is crucial. Tropical forests are vulnerable to droughts associated with extreme El Niño events. However, little is known about how tropical seedling communities respond to El Niño-related droughts, even though patterns of seedling survival shape future forest structure and diversity. Using long-term data from eight tropical moist forests spanning a rainfall gradient in central Panama, we show that community-wide seedling mortality increased by 11% during the extreme 2015-16 El Niño, with mortality increasing most in drought-sensitive species and in wetter forests. These results indicate that severe El Niño-related droughts influence understory dynamics in tropical forests, with effects varying both within and across sites. Our findings suggest that predicted increases in the frequency of extreme El Niño events will alter tropical plant communities through their effects on early life stages.

Co-occurring Fungal Functional Groups Respond Differently to Tree Neighborhoods and Soil Properties Across Three Tropical Rainforests in Panama.

Abiotic and biotic drivers of co-occurring fungal functional guilds across regional-scale environmental gradients remain poorly understood. We characterized fungal communities using Illumina sequencing from soil cores collected across three Neotropical rainforests in Panama that vary in soil properties and plant community composition. We classified each fungal OTU into different functional guilds, namely plant pathogens, saprotrophs, arbuscular mycorrhizal (AM), or ectomycorrhizal (ECM). We measured soil properties and nutrients within each core and determined the tree community composition and richness around each sampling core. Canonical correspondence analyses showed that soil pH and moisture were shared potential drivers of fungal communities for all guilds. However, partial the Mantel tests showed different strength of responses of fungal guilds to composition of trees and soils. Plant pathogens and saprotrophs were more strongly correlated with soil properties than with tree composition; ECM fungi showed a stronger correlation with tree composition than with soil properties; and AM fungi were correlated with soil properties, but not with trees. In conclusion, we show that co-occurring fungal guilds respond differently to abiotic and biotic environmental factors, depending on their ecological function. This highlights the joint role that abiotic and biotic factors play in determining composition of fungal communities, including those associated with plant hosts.

Genomic resources for the Neotropical tree genus Cedrela (Meliaceae) and its relatives.

BACKGROUND: Tree species in the genus Cedrela P. Browne are threatened by timber overexploitation across the Neotropics. Genetic identification of processed timber can be used to supplement wood anatomy to assist in the taxonomic and source validation of protected species and populations of Cedrela. However, few genetic resources exist that enable both species and source identification of Cedrela timber products. We developed several 'omic resources including a leaf transcriptome, organelle genome (cpDNA), and diagnostic single nucleotide polymorphisms (SNPs) that may assist the classification of Cedrela specimens to species and geographic origin and enable future research on this widespread Neotropical tree genus. RESULTS: We designed hybridization capture probes to enrich for thousands of genes from both freshly preserved leaf tissue and from herbarium specimens across eight Meliaceae species. We first assembled a draft de novo transcriptome for C. odorata, and then identified putatively low-copy genes. Hybridization probes for 10,001 transcript models successfully enriched 9795 (98%) of these targets, and analysis of target capture efficiency showed that probes worked effectively for five Cedrela species, with each species showing similar mean on-target sequence yield and depth. The probes showed greater enrichment efficiency for Cedrela species relative to the other three distantly related Meliaceae species. We provide a set of candidate SNPs for species identification of four of the Cedrela species included in this analysis, and present draft chloroplast genomes for multiple individuals of eight species from four genera in the Meliaceae. CONCLUSIONS: Deforestation and illegal logging threaten forest biodiversity globally, and wood screening tools offer enforcement agencies new approaches to identify illegally harvested timber. The genomic resources described here provide the foundation required to develop genetic screening methods for Cedrela species identification and source validation. Due to their transferability across the genus and family as well as demonstrated applicability for both fresh leaves and herbarium specimens, the genomic resources described here provide additional tools for studies examining the ecology and evolutionary history of Cedrela and related species in the Meliaceae.

Source identification of western Oregon Douglas-fir wood cores using mass spectrometry and random forest classification.

PREMISE OF THE STUDY: We investigated whether wood metabolite profiles from direct analysis in real time (time-of-flight) mass spectrometry (DART-TOFMS) could be used to determine the geographic origin of Douglas-fir wood cores originating from two regions in western Oregon, USA. METHODS: Three annual ring mass spectra were obtained from 188 adult Douglas-fir trees, and these were analyzed using random forest models to determine whether samples could be classified to geographic origin, growth year, or growth year and geographic origin. Specific wood molecules that contributed to geographic discrimination were identified. RESULTS: Douglas-fir mass spectra could be differentiated into two geographic classes with an accuracy between 70% and 76%. Classification models could not accurately classify sample mass spectra based on growth year. Thirty-two molecules were identified as key for classifying western Oregon Douglas-fir wood cores to geographic origin. DISCUSSION: DART-TOFMS is capable of detecting minute but regionally informative differences in wood molecules over a small geographic scale, and these differences made it possible to predict the geographic origin of Douglas-fir wood with moderate accuracy. Studies involving DART-TOFMS, alone and in combination with other technologies, will be relevant for identifying the geographic origin of illegally harvested wood.

Seedling growth responses to phosphorus reflect adult distribution patterns of tropical trees.

Soils influence tropical forest composition at regional scales. In Panama, data on tree communities and underlying soils indicate that species frequently show distributional associations to soil phosphorus. To understand how these associations arise, we combined a pot experiment to measure seedling responses of 15 pioneer species to phosphorus addition with an analysis of the phylogenetic structure of phosphorus associations of the entire tree community. Growth responses of pioneers to phosphorus addition revealed a clear tradeoff: species from high-phosphorus sites grew fastest in the phosphorus-addition treatment, while species from low-phosphorus sites grew fastest in the low-phosphorus treatment. Traits associated with growth performance remain unclear: biomass allocation, phosphatase activity and phosphorus-use efficiency did not correlate with phosphorus associations; however, phosphatase activity was most strongly down-regulated in response to phosphorus addition in species from high-phosphorus sites. Phylogenetic analysis indicated that pioneers occur more frequently in clades where phosphorus associations are overdispersed as compared with the overall tree community, suggesting that selection on phosphorus acquisition and use may be strongest for pioneer species with high phosphorus demand. Our results show that phosphorus-dependent growth rates provide an additional explanation for the regional distribution of tree species in Panama, and possibly elsewhere.

Relating belowground microbial composition to the taxonomic, phylogenetic, and functional trait distributions of trees in a tropical forest.

The complexities of the relationships between plant and soil microbial communities remain unresolved. We determined the associations between plant aboveground and belowground (root) distributions and the communities of soil fungi and bacteria found across a diverse tropical forest plot. Soil microbial community composition was correlated with the taxonomic and phylogenetic structure of the aboveground plant assemblages even after controlling for differences in soil characteristics, but these relationships were stronger for fungi than for bacteria. In contrast to expectations, the species composition of roots in our soil core samples was a poor predictor of microbial community composition perhaps due to the patchy, ephemeral, and highly overlapping nature of fine root distributions. Our ability to predict soil microbial composition was not improved by incorporating information on plant functional traits suggesting that the most commonly measured plant traits are not particularly useful for predicting the plot-level variability in belowground microbial communities.

Do global diversity patterns of vertebrates reflect those of monocots?

Few studies of global diversity gradients in plants exist, largely because the data are not available for all species involved. Instead, most global studies have focussed on vertebrates, as these taxa have historically been associated with the most complete data. Here, we address this shortfall by first investigating global diversity gradients in monocots, a morphologically and functionally diverse clade representing a quarter of flowering plant diversity, and then assessing congruence between monocot and vertebrate diversity patterns. To do this, we create a new dataset that merges biome-level associations for all monocot genera with country-level associations for almost all ∼70,000 species. We then assess the evidence for direct versus indirect effects of this plant diversity on vertebrate diversity using a combination of linear regression and structural equation modelling (SEM). Finally, we also calculate overlap of diversity hotspots for monocots and each vertebrate taxon. Monocots follow a latitudinal gradient although with pockets of extra-tropical diversity, mirroring patterns in vertebrates. Monocot diversity is positively associated with vertebrate diversity, but the strength of correlation varies depending on the clades being compared. Monocot diversity explains marginal amounts of variance (<10%) after environmental factors have been accounted for. However, correlations remain among model residuals, and SEMs apparently reveal some direct effects of monocot richness. Our results suggest that collinear responses to environmental gradients are behind much of the congruence observed, but that there is some evidence for direct effects of producer diversity on consumer diversity. Much remains to be done before broad-scale diversity gradients among taxa are fully explained. Our dataset of monocot distributions will aid in this endeavour.

Barro Colorado Island's phylogenetic assemblage structure across fine spatial scales and among clades of different ages.

Phylogenetic analyses of assemblage membership provide insight into how ecological communities are structured. However, despite the scale-dependency of many ecological processes, little is known about how assemblage and source pool size definitions can be altered, either alone or together, to provide insight into how ecological diversity is maintained. Moreover, although studies have acknowledged that different clades within an assemblage may be structured by different forces, there has been no attempt to relate the age of a clade to its community phylogenetic structure. Using assemblage phylogenies and spatially explicit data for trees from Barro Colorado Island (BCI), we show that larger assemblages, and assemblages with larger source pools, are more phylogenetically clustered. We argue that this reflects competition, the influence of pathogens, and chance assembly at smaller spatial scales, all operating within the context of wider-scale habitat filtering. A community phylogenetic measure that is based on a null model derived explicitly from trait evolution theory, D, is better able to detect these differences than commonly used measures such as SES(MPD) and SES(MNTD). We also detect a moderate tendency for stronger phylogenetic clustering in younger clades, which suggests that coarse analyses of diverse assemblages may be missing important variation among clades. Our results emphasize the importance of spatial and phylogenetic scale in community phylogenetics and show how varying these scales can help to untangle complex assembly processes.

The roots of diversity: below ground species richness and rooting distributions in a tropical forest revealed by DNA barcodes and inverse modeling.

BACKGROUND: Plants interact with each other, nutrients, and microbial communities in soils through extensive root networks. Understanding these below ground interactions has been difficult in natural systems, particularly those with high plant species diversity where morphological identification of fine roots is difficult. We combine DNA-based root identification with a DNA barcode database and above ground stem locations in a floristically diverse lowland tropical wet forest on Barro Colorado Island, Panama, where all trees and lianas >1 cm diameter have been mapped to investigate richness patterns below ground and model rooting distributions. METHODOLOGY/PRINCIPAL FINDINGS: DNA barcode loci, particularly the cpDNA locus trnH-psba, can be used to identify fine and small coarse roots to species. We recovered 33 species of roots from 117 fragments sequenced from 12 soil cores. Despite limited sampling, we recovered a high proportion of the known species in the focal hectare, representing approximately 14% of the measured woody plant richness. This high value is emphasized by the fact that we would need to sample on average 13 m(2) at the seedling layer and 45 m(2) for woody plants >1 cm diameter to obtain the same number of species above ground. Results from inverse models parameterized with the locations and sizes of adults and the species identifications of roots and sampling locations indicates a high potential for distal underground interactions among plants. CONCLUSIONS: DNA barcoding techniques coupled with modeling approaches should be broadly applicable to studying root distributions in any mapped vegetation plot. We discuss the implications of our results and outline how second-generation sequencing technology and environmental sampling can be combined to increase our understanding of how root distributions influence the potential for plant interactions in natural ecosystems.

Plant DNA barcodes and a community phylogeny of a tropical forest dynamics plot in Panama.

The assembly of DNA barcode libraries is particularly relevant within species-rich natural communities for which accurate species identifications will enable detailed ecological forensic studies. In addition, well-resolved molecular phylogenies derived from these DNA barcode sequences have the potential to improve investigations of the mechanisms underlying community assembly and functional trait evolution. To date, no studies have effectively applied DNA barcodes sensu strictu in this manner. In this report, we demonstrate that a three-locus DNA barcode when applied to 296 species of woody trees, shrubs, and palms found within the 50-ha Forest Dynamics Plot on Barro Colorado Island (BCI), Panama, resulted in >98% correct identifications. These DNA barcode sequences are also used to reconstruct a robust community phylogeny employing a supermatrix method for 281 of the 296 plant species in the plot. The three-locus barcode data were sufficient to reliably reconstruct evolutionary relationships among the plant taxa in the plot that are congruent with the broadly accepted phylogeny of flowering plants (APG II). Earlier work on the phylogenetic structure of the BCI forest dynamics plot employing less resolved phylogenies reveals significant differences in evolutionary and ecological inferences compared with our data and suggests that unresolved community phylogenies may have increased type I and type II errors. These results illustrate how highly resolved phylogenies based on DNA barcode sequence data will enhance research focused on the interface between community ecology and evolution.