Anthony S. Amend

Quantifying microbial communities with 454 pyrosequencing: does read abundance count?

Pyrosequencing technologies have revolutionized how we describe and compare complex microbial communities. In 454 pyrosequencing data sets, the abundance of reads pertaining to taxa or phylotypes is commonly interpreted as a measure of genic or taxon abundance, useful for quantitative comparisons of community similarity. Potentially systematic biases inherent in sample processing, amplification and sequencing, however, may alter read abundance and reduce the utility of quantitative metrics. Here, we examine the relationship between read abundance and biological abundance in a sample of house dust spiked with known quantities and identities of fungi along a dilution gradient. Our results show one order of magnitude differences in read abundance among species. Precision of quantification within species along the dilution gradient varied from R2 of 0.96–0.54. Read‐quality based processing stringency profoundly affected the abundance of one species containing long homopolymers in a read orientation‐biased manner. Order‐level composition of background environmental fungal communities determined from pyrosequencing data was comparable with that derived from cloning and Sanger sequencing and was not biased by read orientation. We conclude that read abundance is approximately quantitative within species, but between‐species comparisons can be biased by innate sequence structure. Our results showed a trade off between sequence quality stringency and quantification. Careful consideration of sequence processing methods and community analyses are warranted when testing hypotheses using read abundance data.

Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data

Summary 1. The nuclear ribosomal internal transcribed spacer (ITS) region is the primary choice for molecular identification of fungi. Its two highly variable spacers (ITS1 and ITS2) are usually species specific, whereas the intercalary 5.8S gene is highly conserved. For sequence clustering and BLAST searches, it is often advantageous to rely on either one of the variable spacers but not the conserved 5.8S gene. To identify and extract ITS1 and ITS2 from large taxonomic and environmental data sets is, however, often difficult, and many ITS sequences are incorrectly delimited in the public sequence databases. 2. We introduce ITSx, a Perl-based software tool to extract ITS1, 5.8S and ITS2 – as well as full-length ITS sequences – from both Sanger and high-throughput sequencing data sets. ITSx uses hidden Markov models computed from large alignments of a total of 20 groups of eukaryotes, including fungi, metazoans and plants, and the sequence extraction is based on the predicted positions of the ribosomal genes in the sequences. 3. ITSx has a very high proportion of true-positive extractions and a low proportion of false-positive extractions. Additionally, process parallelization permits expedient analyses of very large data sets, such as a one million sequence amplicon pyrosequencing data set. ITSx is rich in features and written to be easily incorporated into automated sequence analysis pipelines. 4. ITSx paves the way for more sensitive BLAST searches and sequence clustering operations for the ITS region in eukaryotes. The software also permits elimination of non-ITS sequences from any data set. This is particularly useful for amplicon-based next-generation sequencing data sets, where insidious non-target sequences are often found among the target sequences. Such non-target sequences are difficult to find by other means and would contribute noise to diversity estimates if left in the data set.

Indoor fungal composition is geographically patterned and more diverse in temperate zones than in the tropics.

Fungi are ubiquitous components of indoor human environments, where most contact between humans and microbes occurs. The majority of these organisms apparently play a neutral role, but some are detrimental to human lifestyles and health. Recent studies that used culture-independent sampling methods demonstrated a high diversity of indoor fungi distinct from that of outdoor environments. Others have shown temporal fluctuations of fungal assemblages in human environments and modest correlations with human activity, but global-scale patterns have not been examined, despite the manifest significance of biogeography in other microbial systems. Here we present a global survey of fungi from indoor environments (n = 72), using both taxonomic and phylogeny-informative molecular markers to determine whether global or local indoor factors determine indoor fungal composition. Contrary to common ecological patterns, we show that fungal diversity is significantly higher in temperate zones than in the tropics, with distance from the equator being the best predictor of phylogenetic community similarity. Fungal composition is significantly auto-correlated at the national and hemispheric spatial scales. Remarkably, building function has no significant effect on indoor fungal composition, despite stark contrasts between architecture and materials of some buildings in close proximity. Distribution of individual taxa is significantly range- and latitude-limited compared with a null model of randomized distribution. Our results suggest that factors driving fungal composition are primarily global rather than mediated by building design or function.

Coral-associated marine fungi form novel lineages and heterogeneous assemblages

Coral stress tolerance is intricately tied to the animals association with microbial symbionts. The most well-known of these symbioses is that between corals and their dinoflagellate photobionts (Symbiodinium spp.), whose genotype indirectly affects whether a coral can survive cyclical and anthropogenic warming events. Fungi comprise a lesser-known coral symbiotic community whose taxonomy, stability and function is largely un-examined. To assess how fungal communities inside a coral host correlate with water temperature and the genotype of co-occurring Symbiodinium, we sampled Acropora hyacinthus coral colonies from adjacent natural pools with different water temperatures and Symbiodinium identities. Phylogenetic analysis of coral-associated fungal ribosomal DNA amplicons showed a high diversity of Basidiomycetes and Ascomycetes, including several clades separated from known fungal taxa by long and well-supported branches. Community similarity did not correlate with any measured variables, and total fungal community composition was highly variable among A. hyacinthus coral colonies. Colonies in the warmer pool contained more phylogenetically diverse fungal communities than the colder pool and contained statistically significant ‘indicator’ species. Four taxa were present in all coral colonies sampled, and may represent obligate associates. Messenger RNA sequenced from a subset of these same colonies contained an abundance of transcripts involved in metabolism of complex biological molecules. Coincidence between the taxonomic diversity found in the DNA and RNA analysis indicates a metabolically active and diverse resident marine fungal community.

From Dandruff to Deep-Sea Vents: Malassezia-like Fungi Are Ecologically Hyper-diverse

As the dominant component of the mycobiota on human skin [1] —both healthy and diseased [2] —the genus Malassezia has received a fair amount of attention. Since the middle of the 19th century, researchers have linked these fungi with skin maladies such as dandruff and eczema [3], but their difficulty to culture axenically long hampered studies of their systematics and diversity [4]. Malassezia is the sole genus within the fungal order Malasseziales, contained within the proposed subphylum Malasseziomycetes (anonymous reviewer; personal communication). Although Malassezia is sister to the so-called “smut” plant pathogens, they are markedly divergent in ecological terms. A hallmark of Malassezia species is their incomplete fatty acids synthesis metabolic pathway, and reliance, instead, on a suite of extracellular lipases, phospholipases, and acid sphingomyelinases [5]. In fact, only a single species, M. pachydermatis, is able to survive in axenic culture lacking lipid amendment [6]. Until recently, it was assumed that Malassezia evolved into a specialized and narrow niche associated with the skin of mammalian hosts. However, culture-independent studies of fungi from environmental samples show that Malassezia are exceedingly widespread and ecologically diverse [3]. Recent studies in little-characterized marine environments point to extensive diversification of Malassezia-like organisms, providing exciting opportunities to explore the ecology, evolution and diversity of this enigmatic group.

Microbial response to simulated global change is phylogenetically conserved and linked with functional potential.

The high diversity of microbial communities hampers predictions about their responses to global change. Here we investigate the potential for using a phylogenetic, trait-based framework to capture the response of bacteria and fungi to global change manipulations. Replicated grassland plots were subjected to 3+ years of drought and nitrogen fertilization. The responses of leaf litter bacteria and fungi to these simulated changes were significantly phylogenetically conserved. Proportional changes in abundance were highly correlated among related organisms, such that relatives with approximately 5% ribosomal DNA genetic distance showed similar responses to the treatments. A microbe’s change in relative abundance was significantly correlated between the treatments, suggesting a compromise between numerical abundance in undisturbed environments and resistance to change in general, independent of disturbance type. Lineages in which at least 90% of the microbes shared the same response were circumscribed at a modest phylogenetic depth (τD 0.014–0.021), but significantly larger than randomized simulations predict. In several clades, phylogenetic depth of trait consensus was higher. Fungal response to drought was more conserved than was response to nitrogen fertilization, whereas bacteria responded equally to both treatments. Finally, we show that a bacterium’s response to the manipulations is correlated with its potential functional traits (measured here as the number of glycoside hydrolase genes encoding the capacity to degrade different types of carbohydrates). Together, these results suggest that a phylogenetic, trait-based framework may be useful for predicting shifts in microbial composition and functioning in the face of global change.

Temporal variation overshadows the response of leaf litter microbial communities to simulated global change

Bacteria and fungi drive the decomposition of dead plant biomass (litter), an important step in the terrestrial carbon cycle. Here we investigate the sensitivity of litter microbial communities to simulated global change (drought and nitrogen addition) in a California annual grassland. Using 16S and 28S rDNA amplicon pyrosequencing, we quantify the response of the bacterial and fungal communities to the treatments and compare these results to background, temporal (seasonal and interannual) variability of the communities. We found that the drought and nitrogen treatments both had significant effects on microbial community composition, explaining 2–6% of total compositional variation. However, microbial composition was even more strongly influenced by seasonal and annual variation (explaining 14–39%). The response of microbial composition to drought varied by season, while the effect of the nitrogen addition treatment was constant through time. These compositional responses were similar in magnitude to those seen in microbial enzyme activities and the surrounding plant community, but did not correspond to a consistent effect on leaf litter decomposition rate. Overall, these patterns indicate that, in this ecosystem, temporal variability in the composition of leaf litter microorganisms largely surpasses that expected in a short-term global change experiment. Thus, as for plant communities, future microbial communities will likely be determined by the interplay between rapid, local background variability and slower, global changes.

Nutrient enrichment increased species richness of leaf litter fungal assemblages in a tropical forest

Microbial communities play a major role in terrestrial ecosystem functioning, but the determinates of their diversity and functional interactions are not well known. In this study, we explored leaf litter fungal diversity in a diverse Panama lowland tropical forest in which a replicated factorial N, P, K and micronutrient fertilization experiment of 40 × 40 m plots had been ongoing for nine years. We extracted DNA from leaf litter samples and used fungal‐specific amplification and a 454 pyrosequencing approach to sequence two loci, the nuclear ribosomal internal transcribed spacer (ITS) region and the nuclear ribosomal large subunit (LSU) D1 region. Using a 95% sequence similarity threshold for ITS1 spacer recovered a total of 2523 OTUs, and the number of unique ITS1 OTUs per 0.5–1.0 g leaf litter sample ranged from 55 to 177. Ascomycota were the dominant phylum among the leaf litter fungi (71% of the OTUs), followed by Basidiomycota (26% of the OTUs). In contrast to our expectations based on temperate ecosystems, long‐term addition of nutrients increased, rather than decreased, species richness relative to controls. Effect of individual nutrients was more subtle and seen primarily as changes in community compositions especially at lower taxonomic levels, rather than as significant changes in species richness. For example, plots receiving P tended to show a greater similarity in community composition compared to the other nutrient treatments, the +PK, +NK and +NPK plots appeared to be more dominated by the Nectriaceae than other treatments, and indicator species for particular nutrient combinations were identified.

A Unique Signal Distorts the Perception of Species Richness and Composition in High-Throughput Sequencing Surveys of Microbial Communities: a Case Study of Fungi in Indoor Dust

Sequence-based surveys of microorganisms in varied environments have found extremely diverse assemblages. A standard practice in current high-throughput sequence (HTS) approaches in microbial ecology is to sequence the composition of many environmental samples at once by pooling amplicon libraries at a common concentration before processing on one run of a sequencing platform. Biomass of the target taxa, however, is not typically determined prior to HTS, and here, we show that when abundances of the samples differ to a large degree, this standard practice can lead to a perceived bias in community richness and composition. Fungal signal in settled dust of five university teaching laboratory classrooms, one of which was used for a mycology course, was surveyed. The fungal richness and composition in the dust of the nonmycology classrooms were remarkably similar to each other, while the mycology classroom was dominated by abundantly sporulating specimen fungi, particularly puffballs, and appeared to have a lower overall richness based on rarefaction curves and richness estimators. The fungal biomass was three to five times higher in the mycology classroom than the other classrooms, indicating that fungi added to the mycology classroom swamped the background fungi present in indoor air. Thus, the high abundance of a few taxa can skew the perception of richness and composition when samples are sequenced to an even depth. Next, we used in silico manipulations of the observed data to confirm that a unique signature can be identified with HTS approaches when the source is abundant, whether or not the taxon identity is distinct. Lastly, aerobiology of indoor fungi is discussed.

Measuring carbon gains from fungal networks in understory plants from the tribe Pyroleae (Ericaceae): a field manipulation and stable isotope approach.

Partial mycoheterotrophy, a newly discovered form of mixotrophy in plants, has been described in at least two major lineages of angiosperms, the orchids and ericaceous plants in the tribe Pyroleae. Partial mycoheterotrophy entails carbon gains both directly from photosynthesis and via symbiotic mycorrhizal fungi, but determining the degree of plant dependence on fungal carbon is challenging. The purpose of this study was to determine if two chlorophyllous species of Pyroleae, Chimaphila umbellata and Pyrola picta, were receiving carbon via mycorrhizal networks and, if so, if their proportional dependency on fungal carbon gains increased under reduced light conditions. This was accomplished by a field experiment that manipulated light and plants’ access to mycorrhizal networks, and by using the stable carbon isotope composition (δ13C) of leaf soluble sugars as a marker for the level of mycoheterotrophy. Based on leaf soluble sugars δ13C values, we calculated a site-independent isotope enrichment factor as a measure of fungal contributions to plant C. We found that, under each treatment and over time, the two test species demonstrated different isotopic responses caused by their different intrinsic physiologies. Our data, along with previously published studies, suggest that Chimaphila umbellata is primarily an autotrophic understory plant, while Pyrola picta may be capable of partial mycoheterotrophy. However, in this study, a 50% decrease in light availability did not significantly change the relative dependency of P. picta on carbon gains via mycoheterotrophy.