Root-associated fungal communities are influenced more by soils than by plant-host root traits in a Chinese tropical forest.

Forest fungal communities are shaped by the interactions between host tree root systems and the associated soil conditions. We investigated how the soil environment, root morphological traits, and root chemistry influence root-inhabiting fungal communities in three tropical forest sites of varying successional status in Xishuangbanna, China. For 150 trees of 66 species, we measured root morphology and tissue chemistry. Tree species identity was confirmed by sequencing rbcL, and root-associated fungal (RAF) communities were determined using high-throughput ITS2 sequencing. Using distance-based redundancy analysis and hierarchical variation partitioning, we quantified the relative importance of two soil variables (site average total phosphorus and available phosphorus), four root traits (dry matter content, tissue density, specific tip abundance, and forks), and three root tissue elemental concentrations (nitrogen, calcium, and manganese) on RAF community dissimilarity. The root and soil environment collectively explained 23% of RAF compositional variation. Soil phosphorus explained 76% of that variation. Twenty fungal taxa differentiated RAF communities among the three sites. Soil phosphorus most strongly affects RAF assemblages in this tropical forest. Variation in root calcium and manganese concentrations and root morphology among tree hosts, principally an architectural trade-off between dense, highly branched vs less-dense, herringbone-type root systems, are important secondary determinants.

Effects of cavity orientation on nesting success inferred from long-term monitoring of the endangered red-cockaded woodpecker.

Animals that create structures often display non-random patterns in the direction of their constructions. This tendency of oriented construction is widely presumed to be an adaptive trait of the constructor's extended phenotype, but there is little empirical support for this hypothesis. Particularly, for cavity nesting-birds there is a lack of studies examining this issue. In this study of a primary cavity excavator, the endangered red-cockaded woodpecker (Dryobates borealis), we show that cavity entrances exhibited a strong westward bias in all 11 of the populations examined throughout the geographic range of the species in the southeastern United States. This species requires cavities in living pine trees for roosting and nesting that often take many years to complete, resulting in many incomplete excavations on the landscape. We used population monitoring data to show that orientation was stronger among completed cavities than incomplete cavities. There was a significant correlation between latitude and average cavity direction among populations, turning northward with increasing latitude, suggesting adaptation to local conditions. Long-term monitoring data showed that cavity orientation and breeding group size are correlated with egg hatching rates, fledging rates, and the total number of fledglings produced per nest. Our results provide empirical evidence from extensive long-term data that directional orientation in animal constructions is an important feature of the extended animal phenotype and have immediate implications for animal ecology and the conservation of endangered species.

Disease-related population declines in bats demonstrate non-exchangeability in generalist predators.

The extent to which persisting species may fill the functional role of extirpated or declining species has profound implications for the structure of biological communities and ecosystem functioning. In North America, arthropodivorous bats are threatened on a continent-wide scale by the spread of white-nose syndrome (WNS), a disease caused by the fungus Pseudogymnoascus destructans. We tested whether bat species that display lower mortality from this disease can partially fill the functional role of other bat species experiencing population declines. Specifically, we performed high-throughput amplicon sequencing of guano from two generalist predators: the little brown bat (Myotis lucifugus) and big brown bat (Eptesicus fuscus). We then compared changes in prey consumption before versus after population declines related to WNS. Dietary niches contracted for both species after large and abrupt declines in little brown bats and smaller declines in big brown bats, but interspecific dietary overlap did not change. Furthermore, the incidence and taxonomic richness of agricultural pest taxa detected in diet samples decreased following bat population declines. Our results suggest that persisting generalist predators do not necessarily expand their dietary niches following population declines in other predators, providing further evidence that the functional roles of different generalist predators are ecologically distinct.

Understudied, underrepresented, and unknown: Methodological biases that limit detection of early diverging fungi from environmental samples.

Metabarcoding is an important tool for understanding fungal communities. The internal transcribed spacer (ITS) rDNA is the accepted fungal barcode but has known problems. The large subunit (LSU) rDNA has also been used to investigate fungal communities but available LSU metabarcoding primers were mostly designed to target Dikarya (Ascomycota + Basidiomycota) with little attention to early diverging fungi (EDF). However, evidence from multiple studies suggests that EDF comprise a large portion of unknown diversity in community sampling. Here, we investigate how DNA marker choice and methodological biases impact recovery of EDF from environmental samples. We focused on one EDF lineage, Zoopagomycota, as an example. We evaluated three primer sets (ITS1F/ITS2, LROR/LR3, and LR3 paired with new primer LR22F) to amplify and sequence a Zoopagomycota mock community and a set of 146 environmental samples with Illumina MiSeq. We compared two taxonomy assignment methods and created an LSU reference database compatible with AMPtk software. The two taxonomy assignment methods recovered strikingly different communities of fungi and EDF. Target fragment length variation exacerbated PCR amplification biases and influenced downstream taxonomic assignments, but this effect was greater for EDF than Dikarya. To improve identification of LSU amplicons we performed phylogenetic reconstruction and illustrate the advantages of this critical tool for investigating identified and unidentified sequences. Our results suggest much of the EDF community may be missed or misidentified with "standard" metabarcoding approaches and modified techniques are needed to understand the role of these taxa in a broader ecological context.

Preinvasion Assessment of Exotic Bark Beetle-Vectored Fungi to Detect Tree-Killing Pathogens.

Exotic diseases and pests of trees have caused continental-scale disturbances in forest ecosystems and industries, and their invasions are considered largely unpredictable. We tested the concept of preinvasion assessment of not yet invasive organisms, which enables empirical risk assessment of potential invasion and impact. Our example assesses fungi associated with Old World bark and ambrosia beetles and their potential to impact North American trees. We selected 55 Asian and European scolytine beetle species using host use, economic, and regulatory criteria. We isolated 111 of their most consistent fungal associates and tested their effect on four important southeastern American pine and oak species. Our test dataset found no highly virulent pathogens that should be classified as an imminent threat. Twenty-two fungal species were minor pathogens, which may require context-dependent response for their vectors at North American borders, while most of the tested fungi displayed no significant impact. Our results are significant in three ways; they ease the concerns over multiple overseas fungus vectors suspected of heightened potential risk, they provide a basis for the focus on the prevention of introduction and establishment of species that may be of consequence, and they demonstrate that preinvasion assessment, if scaled up, can support practical risk assessment of exotic pathogens.

Discovering the role of Patagonian birds in the dispersal of truffles and other mycorrhizal fungi.

Dispersal is a key process that impacts population dynamics and structures biotic communities. Dispersal limitation influences the assembly of plant and microbial communities, including mycorrhizal fungi and their plant hosts. Mycorrhizal fungi play key ecological roles in forests by feeding nutrients to plants in exchange for sugars, so the dispersal of mycorrhizal fungi spores actively shapes plant communities. Although many fungi rely on wind for spore dispersal, some fungi have lost the ability to shoot their spores into the air and instead produce enclosed belowground fruiting bodies (truffles) that rely on animals for dispersal. The role of mammals in fungal spore dispersal is well documented, but the relevance of birds as dispersal agents of fungi has been understudied, despite the prominence of birds as seed dispersal vectors. Here, we use metagenomics and epifluorescence microscopy to demonstrate that two common, widespread, and endemic Patagonian birds, chucao tapaculos (Scelorchilus rubecula) and black-throated huet-huets (Pteroptochos tarnii), regularly consume mycorrhizal fungi and disperse viable spores via mycophagy. Our metagenomic analysis indicates that these birds routinely consume diverse mycorrhizal fungi, including many truffles, that are symbiotically associated with Nothofagaceae trees that dominate Patagonian forests. Epifluorescence microscopy of fecal samples confirmed that the birds dispersed copious viable spores from truffles and other mycorrhizal fungi. We show that fungi are a common food for both bird species and that this animal-fungi symbiosis is widespread and ecologically important in Patagonia. Our findings indicate that birds may also act as cryptic but critical fungal dispersal agents in other ecosystems.

Effects of Field Fumigation and Inoculation With the Pecan Truffle (Tuber lyonii) on the Fungal Community of Pecan (Carya illinoinensis) Seedlings Over 5 Years.

Truffle fungi are esteemed for their aromatic qualities and are among the most widely cultivated edible ectomycorrhizal fungi. Here we document a successful method for establishing Tuber lyonii, the pecan truffle, on pecan (Carya illinoinensis) seedlings in a field setting. We assessed the impacts of soil fumigation and varying concentrations of truffle spore inoculum on the ectomycorrhizal fungal and the complete fungal communities as well as the colonization of T. lyonii on pecan roots at three nurseries in Georgia, United States. To identify fungal communities on pecan seedlings, we performed high-throughput amplicon sequencing of the fungal ITS1 rDNA region. Our 5-year long field experiment demonstrates that fumigation and inoculation together resulted in the highest persistence of T. lyonii on pecan roots. While fungal OTU numbers fluctuated over the years of our experiments, there was no statistical support to demonstrate diversification of communities when Shannon diversity metrics were used. However, we did find that older seedlings were less likely to be dominated by T. lyonii compared to younger ones, suggesting successional changes in the fungal community over time. This suggests that transplanting inoculated seedlings after 2 or 3 years post-inoculation is optimal for future truffle propagation efforts. Our results demonstrate that T. lyonii can be established in situ with methods that are compatible with current pecan nursery industry practices and that fungal communities on pecan seedlings vary depending on the experimental treatments used during planting. While the pecan truffle is not yet widely cultivated, our results provide insights for future large-scale cultivation of this and perhaps other Tuber species.

Predator preferences shape the diets of arthropodivorous bats more than quantitative local prey abundance.

Although most predators are generalists, the majority of studies on the association between prey availability and prey consumption have focused on specialist predators. To investigate the role of highly generalist predators in a complex food web, we measured the relationships between prey consumption and prey availability in two common arthropodivorous bats. Specifically, we used high-throughput amplicon sequencing coupled with a known mock community to characterize seasonal changes in little brown and big brown bat diets. We then linked spatiotemporal variation in prey consumption with quantitative prey availability estimated from intensive prey community sampling. We found that although quantitative prey availability fluctuated substantially over space and time, the most commonly consumed prey items were consistently detected in bat diets independently of their respective abundance. Positive relationships between prey abundance and probability of consumption were found only among prey groups that were less frequently detected in bat diets. While the probability of prey consumption was largely unrelated to abundance, the community structure of prey detected in bat diets was influenced by the local or regional abundance of prey. Observed patterns suggest that while little brown and big brown bats maintain preferences for particular prey independently of quantitative prey availability, total dietary composition may reflect some degree of opportunistic foraging. Overall, our findings suggest that generalist predators can display strong prey preferences that persist despite quantitative changes in prey availability.

A total crapshoot? Evaluating bioinformatic decisions in animal diet metabarcoding analyses.

Metabarcoding studies provide a powerful approach to estimate the diversity and abundance of organisms in mixed communities in nature. While strategies exist for optimizing sample and sequence library preparation, best practices for bioinformatic processing of amplicon sequence data are lacking in animal diet studies. Here we evaluate how decisions made in core bioinformatic processes, including sequence filtering, database design, and classification, can influence animal metabarcoding results. We show that denoising methods have lower error rates compared to traditional clustering methods, although these differences are largely mitigated by removing low-abundance sequence variants. We also found that available reference datasets from GenBank and BOLD for the animal marker gene cytochrome oxidase I (COI) can be complementary, and we discuss methods to improve existing databases to include versioned releases. Taxonomic classification methods can dramatically affect results. For example, the commonly used Barcode of Life Database (BOLD) Classification API assigned fewer names to samples from order through species levels using both a mock community and bat guano samples compared to all other classifiers (vsearch-SINTAX and q2-feature-classifier's BLAST + LCA, VSEARCH + LCA, and Naive Bayes classifiers). The lack of consensus on bioinformatics best practices limits comparisons among studies and may introduce biases. Our work suggests that biological mock communities offer a useful standard to evaluate the myriad computational decisions impacting animal metabarcoding accuracy. Further, these comparisons highlight the need for continual evaluations as new tools are adopted to ensure that the inferences drawn reflect meaningful biology instead of digital artifacts.

Major histocompatibility complex variation is similar in little brown bats before and after white-nose syndrome outbreak.

White-nose syndrome (WNS), caused by the fungal pathogen Pseudogymnoascus destructans (Pd), has driven alarming declines in North American hibernating bats, such as little brown bat (Myotis lucifugus). During hibernation, infected little brown bats are able to initiate anti-Pd immune responses, indicating pathogen-mediated selection on the major histocompatibility complex (MHC) genes. However, such immune responses may not be protective as they interrupt torpor, elevate energy costs, and potentially lead to higher mortality rates. To assess whether WNS drives selection on MHC genes, we compared the MHC DRB gene in little brown bats pre- (Wisconsin) and post- (Michigan, New York, Vermont, and Pennsylvania) WNS (detection spanning 2014-2015). We genotyped 131 individuals and found 45 nucleotide alleles (27 amino acid alleles) indicating a maximum of 3 loci (1-5 alleles per individual). We observed high allelic admixture and a lack of genetic differentiation both among sampling sites and between pre- and post-WNS populations, indicating no signal of selection on MHC genes. However, post-WNS populations exhibited decreased allelic richness, reflecting effects from bottleneck and drift following rapid population declines. We propose that mechanisms other than adaptive immunity are more likely driving current persistence of little brown bats in affected regions.

Relationships among wood-boring beetles, fungi, and the decomposition of forest biomass.

A prevailing paradigm in forest ecology is that wood-boring beetles facilitate wood decay and carbon cycling, but empirical tests have yielded mixed results. We experimentally determined the effects of wood borers on fungal community assembly and wood decay within pine trunks in the southeastern United States. Pine trunks were made either beetle-accessible or inaccessible. Fungal communities were compared using culturing and high-throughput amplicon sequencing (HTAS) of DNA and RNA. Prior to beetle infestation, living pines had diverse fungal endophyte communities. Endophytes were displaced by beetle-associated fungi in beetle-accessible trees, whereas some endophytes persisted as saprotrophs in beetle-excluded trees. Beetles increased fungal diversity several fold. Over forty taxa of Ascomycota were significantly associated with beetles, but beetles were not consistently associated with any known wood-decaying fungi. Instead, increasing ambrosia beetle infestations caused reduced decay, consistent with previous in vitro experiments that showed beetle-associated fungi reduce decay rates by competing with decay fungi. No effect of bark-inhabiting beetles on decay was detected. Platypodines carried significantly more fungal taxa than scolytines. Molecular results were validated by synthetic and biological mock communities and were consistent across methodologies. RNA sequencing confirmed that beetle-associated fungi were biologically active in the wood. Metabarcode sequencing of the LSU/28S marker recovered important fungal symbionts that were missed by ITS2, though community-level effects were similar between markers. In contrast to the current paradigm, our results indicate ambrosia beetles introduce diverse fungal communities that do not extensively decay wood, but instead reduce decay rates by competing with wood decay fungi.

A selective fungal transport organ (mycangium) maintains coarse phylogenetic congruence between fungus-farming ambrosia beetles and their symbionts.

Thousands of species of ambrosia beetles excavate tunnels in wood to farm fungi. They maintain associations with particular lineages of fungi, but the phylogenetic extent and mechanisms of fidelity are unknown. We test the hypothesis that selectivity of their mycangium enforces fidelity at coarse phylogenetic scales, while permitting promiscuity among closely related fungal mutualists. We confirm a single evolutionary origin of the Xylosandrus complex-a group of several xyleborine genera that farm fungi in the genus Ambrosiella. Multi-level co-phylogenetic analysis revealed frequent symbiont switching within major Ambrosiella clades, but not between clades. The loss of the mycangium in Diuncus, a genus of evolutionary cheaters, was commensurate with the loss of fidelity to fungal clades, supporting the hypothesis that the mycangium reinforces fidelity. Finally, in vivo experiments tracked symbiotic compatibility throughout the symbiotic life cycle of Xylosandrus compactus and demonstrated that closely related Ambrosiella symbionts are interchangeable, but the probability of fungal uptake in the mycangium was significantly lower in more phylogenetically distant species of symbionts. Symbiont loads in experimental subjects were similar to wild-caught beetles. We conclude that partner choice in ambrosia beetles is achieved in the mycangium, and co-phylogenetic inferences can be used to predict the likelihood of specific symbiont switches.

An improved method for utilizing high-throughput amplicon sequencing to determine the diets of insectivorous animals.

DNA analysis of predator faeces using high-throughput amplicon sequencing (HTS) enhances our understanding of predator-prey interactions. However, conclusions drawn from this technique are constrained by biases that occur in multiple steps of the HTS workflow. To better characterize insectivorous animal diets, we used DNA from a diverse set of arthropods to assess PCR biases of commonly used and novel primer pairs for the mitochondrial gene, cytochrome oxidase C subunit 1 (COI). We compared diversity recovered from HTS of bat guano samples using a commonly used primer pair "ZBJ" to results using the novel primer pair "ANML." To parameterize our bioinformatics pipeline, we created an arthropod mock community consisting of single-copy (cloned) COI sequences. To examine biases associated with both PCR and HTS, mock community members were combined in equimolar amounts both pre- and post-PCR. We validated our system using guano from bats fed known diets and using composite samples of morphologically identified insects collected in pitfall traps. In PCR tests, the ANML primer pair amplified 58 of 59 arthropod taxa (98%), whereas ZBJ amplified 24-40 of 59 taxa (41%-68%). Furthermore, in an HTS comparison of field-collected samples, the ANML primers detected nearly fourfold more arthropod taxa than the ZBJ primers. The additional arthropods detected include medically and economically relevant insect groups such as mosquitoes. Results revealed biases at both the PCR and sequencing levels, demonstrating the pitfalls associated with using HTS read numbers as proxies for abundance. The use of an arthropod mock community allowed for improved bioinformatics pipeline parameterization.

Identifying the "Mushroom of Immortality": Assessing the Ganoderma Species Composition in Commercial Reishi Products.

Species of Ganoderma, commonly called reishi (in Japan) or lingzhi (in China), have been used in traditional medicine for thousands of years, and their use has gained interest from pharmaceutical industries in recent years. Globally, the taxonomy of Ganoderma species is chaotic, and the taxon name Ganoderma lucidum has been used for most laccate (shiny) Ganoderma species. However, it is now known that G. lucidum sensu stricto has a limited native distribution in Europe and some parts of China. It is likely that differences in the quality and quantity of medicinally relevant chemicals occur among Ganoderma species. To determine what species are being sold in commercially available products, twenty manufactured products (e.g., pills, tablets, teas, etc.) and seventeen grow your own (GYO) kits labeled as containing G. lucidum were analyzed. DNA was extracted, and the internal transcribed spacer (ITS) region and translation elongation factor 1-alpha (tef1α) were sequenced with specific fungal primers. The majority (93%) of the manufactured reishi products and almost half of the GYO kits were identified as Ganoderma lingzhi. G. lingzhi is native to Asia and is the most widely cultivated and studied taxon for medicinal use. Illumina MiSeq sequencing of the ITS1 region was performed to determine if multiple Ganoderma species were present. None of the manufactured products tested contained G. lucidum sensu stricto, and it was detected in only one GYO kit. G. lingzhi was detected in most products, but other Ganoderma species were also present, including G. applanatum, G. australe, G. gibbosum, G. sessile, and G. sinense. Our results indicate that the content of these products vary and that better labeling is needed to inform consumers before these products are ingested or marketed as medicine. Of the 17 GYO kits tested, 11 kits contained Ganoderma taxa that are not native to the United States. If fruiting bodies of exotic Ganoderma taxa are cultivated, these GYO kits will likely end up in the environment. The effects of these exotic species to natural ecosystems needs investigation.

Non-biological synthetic spike-in controls and the AMPtk software pipeline improve mycobiome data.

High-throughput amplicon sequencing (HTAS) of conserved DNA regions is a powerful technique to characterize microbial communities. Recently, spike-in mock communities have been used to measure accuracy of sequencing platforms and data analysis pipelines. To assess the ability of sequencing platforms and data processing pipelines using fungal internal transcribed spacer (ITS) amplicons, we created two ITS spike-in control mock communities composed of cloned DNA in plasmids: a biological mock community, consisting of ITS sequences from fungal taxa, and a synthetic mock community (SynMock), consisting of non-biological ITS-like sequences. Using these spike-in controls we show that: (1) a non-biological synthetic control (e.g., SynMock) is the best solution for parameterizing bioinformatics pipelines, (2) pre-clustering steps for variable length amplicons are critically important, (3) a major source of bias is attributed to the initial polymerase chain reaction (PCR) and thus HTAS read abundances are typically not representative of starting values. We developed AMPtk, a versatile software solution equipped to deal with variable length amplicons and quality filter HTAS data based on spike-in controls. While we describe herein a non-biological SynMock community for ITS sequences, the concept and AMPtk software can be widely applied to any HTAS dataset to improve data quality.

Detecting Symbioses in Complex Communities: the Fungal Symbionts of Bark and Ambrosia Beetles Within Asian Pines.

Separating symbioses from incidental associations is a major obstacle in symbiosis research. In this survey of fungi associated with Asian bark and ambrosia beetles, we used quantitative culture and DNA barcode identification to characterize fungal communities associated with co-infesting beetle species in pines (Pinus) of China and Vietnam. To quantitatively discern likely symbioses from coincidental associations, we used multivariate analysis and multilevel pattern analysis (a type of indicator species analysis). Nearly half of the variation in fungal community composition in beetle galleries and on beetle bodies was explained by beetle species. We inferred a spectrum of ecological strategies among beetle-associated fungi: from generalist multispecies associates to highly specialized single-host symbionts that were consistently dominant within the mycangia of their hosts. Statistically significant fungal associates of ambrosia beetles were typically only found with one beetle species. In contrast, bark beetle-associated fungi were often associated with multiple beetle species. Ambrosia beetles and their galleries were frequently colonized by low-prevalence ambrosia fungi, suggesting that facultative ambrosial associations are commonplace, and ecological mechanisms such as specialization and competition may be important in these dynamic associations. The approach used here could effectively delimit symbiotic interactions in any system where symbioses are obscured by frequent incidental associations. It has multiple advantages including (1) powerful statistical tests for non-random associations among potential symbionts, (2) simultaneous evaluation of multiple co-occurring host and symbiont associations, and (3) identifying symbionts that are significantly associated with multiple host species.

Amplicon-Based Sequencing of Soil Fungi from Wood Preservative Test Sites.

Soil samples were collected from field sites in two AWPA (American Wood Protection Association) wood decay hazard zones in North America. Two field plots at each site were exposed to differing preservative chemistries via in-ground installations of treated wood stakes for approximately 50 years. The purpose of this study is to characterize soil fungal species and to determine if long term exposure to various wood preservatives impacts soil fungal community composition. Soil fungal communities were compared using amplicon-based DNA sequencing of the internal transcribed spacer 1 (ITS1) region of the rDNA array. Data show that soil fungal community composition differs significantly between the two sites and that long-term exposure to different preservative chemistries is correlated with different species composition of soil fungi. However, chemical analyses using ICP-OES found levels of select residual preservative actives (copper, chromium and arsenic) to be similar to naturally occurring levels in unexposed areas. A list of indicator species was compiled for each treatment-site combination; functional guild analyses indicate that long-term exposure to wood preservatives may have both detrimental and stimulatory effects on soil fungal species composition. Fungi with demonstrated capacity to degrade industrial pollutants were found to be highly correlated with areas that experienced long-term exposure to preservative testing.

Experimental evidence of a symbiosis between red-cockaded woodpeckers and fungi.

Primary cavity excavators, such as woodpeckers, are ecosystem engineers in many systems. Associations between cavity excavators and fungi have long been hypothesized to facilitate cavity excavation, but these relationships have not been experimentally verified. Fungi may help excavators by softening wood, while excavators may facilitate fungal dispersal. Here we demonstrate that excavators facilitate fungal dispersal and thus we report the first experimental evidence of a symbiosis between fungi and a cavity excavator, the red-cockaded woodpecker (RCW,Picoides borealis). Swab samples of birds showed that RCWs carry fungal communities similar to those found in their completed excavations. A 26-month field experiment using human-made aseptically drilled excavations in live trees, half of which were inaccessible to RCWs, demonstrated that RCWs directly alter fungal colonization and community composition. Experimental excavations that were accessible to RCWs contained fungal communities similar to natural RCW excavations, whereas inaccessible experimental excavations contained significantly different fungal communities. Our work demonstrates a complex symbiosis between cavity excavators and communities of fungi, with implications for forest ecology, wildlife management, and conservation.