Molecular characterization of airborne fungal spores in boreal forests of contrasting human disturbance

In this study we present a new approach to characterize fungal diversity with DNA sequencing of mycelium grown from trapped airborne spores. Fungal spores were extracted systematically from air in three boreal forest sites (clear-cut, young and old-growth forests) using an air sampling device. Internal transcribed spacer (ITS) sequences from the nuclear ribosomal DNA (nrDNA) were generated, and the sequences most likely taxon affinities were established through DNA homology searches. Phylogenetic analyses were used to classify similar sequences into operational taxonomic units (OTUs). The analyses indicated that a total of 84 different OTUs had been sampled, 24 basidiomycetes and 60 ascomycetes. OTUs belonging to the ascomycete orders Helotiales and Pleosporales were most frequent (31 and 18 respectively). A total of 54, 29 and 33 OTUs were sampled, respectively, in the old-growth, young and clear-cut forest sites. Although heavy generalization should be avoided due to few replicates, the results could indicate that old-growth boreal forests have significantly higher airborne fungal species richness than recently managed forests. The study shows that the spore-trapping approach has a great potential for targeting and studying anonymous fungi.     

Plant identity strongly structures the root-associated fungal community in a diverse subtropical forest

Revealing the relationship between plants and root-associated fungi is very important in understanding diversity maintenance and community assembly in ecosystems. However, the community assembly of root-associated fungi of focal plant species along a subtropical plant species diversity gradient is less documented. Here, we examined root-associated fungal communities associated with five ectomycorrhizal (EM) plant species (Betula luminifera, Castanea henryi, Castanopsis fargesii, C. sclerophylla, and Quercus serrate) in a Chinese subtropical woody plant species diversity (1, 2, 4, 8, 16 and 24 species) experiment, using paired-end Illumina MiSeq sequencing of the ITS2 region. In total, we detected 1933 root-associated fungal operational taxonomic units (OTUs) at a 97% sequence similarity level. Plant identity had a significant effect on total and saprotrophic fungal OTU richness, but plant species diversity level had a significant effect on saprotrophic and pathogenic fungal OTU richness. The community composition of total, saprotrophic and EM fungi was structured by plant identity and plant species diversity level. However, the community composition of pathogenic fungi was only shaped by plant identity. This study highlights that plant identity has a stronger effect on the root-associated fungal community than plant species diversity level in a diverse subtropical forest ecosystem.     

Smaller and more numerous harvesting gaps emulate natural forest disturbances: a biodiversity test case using rove beetles (Coleoptera,Staphylinidae)

Aim To evaluate changes in the abundance, species richness and community composition of rove beetles (Coleoptera, Staphylinidae) in response to three configurations of experimental gap cuts and to the effects of ground scarification in early succession yellow birch‐dominated boreal forest. In each experimental treatment, total forest removed was held constant (35% removal by partial cutting with a concomitant decrease in gap size) but the total number of gaps was increased (two, four and eight gaps, respectively), resulting in an experimental increase in the total amount of ‘edge’ within each stand. Location Early succession yellow birch‐dominated forests, Quebec, Canada. Methods Pitfall traps, ANOVA, MIXED procedure in sas ®, post hoc Tukey's adjustment, rarefaction estimates, sum‐of‐squares and distance‐based multivariate regression trees (ssMRT, dbMRT). Results Estimates of species richness using rarefaction were highest in clearcut and two‐gap treatments, decreased in smaller and more numerous gaps and were significantly higher in scarified areas than in unscarified areas. ANOVA indicated a significant impact of harvesting on the overall standardized catch. Post hoc Tukey's tests indicated that the total catch of all rove beetles was significantly higher in uncut forests than in the treated areas. Both sum‐of‐squares and distance‐based multivariate regression trees indicated that community structure of rove beetles differed among treatments. Assemblages were grouped into (a) control plots, (b) four‐ and eight‐gap treatments and (c) two‐gap and clearcut treatments. Main conclusions Rove beetle composition responded significantly to increasing gap size. Composition among intermediate and small‐sized gap treatments (four‐ and eight‐gap treatments) was more similar to uncut control forests than were larger gap treatments (two‐gap) and clearcuts. Effects of scarification were nested within the harvested treatments. When the total area of forest removed is held constant, smaller, more numerous gaps are more similar to uncut control stands than to larger gaps and falls more closely within the natural forest heterogeneity.     

Diversity measures in environmental sequences are highly dependent on alignment quality--data from ITS and new LSU primers targeting basidiomycetes

The ribosomal DNA comprised of the ITS1-5.8S-ITS2 regions is widely used as a fungal marker in molecular ecology and systematics but cannot be aligned with confidence across genetically distant taxa. In order to study the diversity of Agaricomycotina in forest soils, we designed primers targeting the more alignable 28S (LSU) gene, which should be more useful for phylogenetic analyses of the detected taxa. This paper compares the performance of the established ITS1F/4B primer pair, which targets basidiomycetes, to that of two new pairs. Key factors in the comparison were the diversity covered, off-target amplification, rarefaction at different Operational Taxonomic Unit (OTU) cutoff levels, sensitivity of the method used to process the alignment to missing data and insecure positional homology, and the congruence of monophyletic clades with OTU assignments and BLAST-derived OTU names. The ITS primer pair yielded no off-target amplification but also exhibited the least fidelity to the expected phylogenetic groups. The LSU primers give complementary pictures of diversity, but were more sensitive to modifications of the alignment such as the removal of difficult-to align stretches. The LSU primers also yielded greater numbers of singletons but also had a greater tendency to produce OTUs containing sequences from a wider variety of species as judged by BLAST similarity. We introduced some new parameters to describe alignment heterogeneity based on Shannon entropy and the extent and contents of the OTUs in a phylogenetic tree space. Our results suggest that ITS should not be used when calculating phylogenetic trees from genetically distant sequences obtained from environmental DNA extractions and that it is inadvisable to define OTUs on the basis of very heterogeneous alignments.     

Metabarcoding of fungal communities associated with bark beetles

Many species of fungi are closely allied with bark beetles, including many tree pathogens, but their species richness and patterns of distribution remain largely unknown. We established a protocol for metabarcoding of fungal communities directly from total genomic DNA extracted from individual beetles, showing that the ITS3/4 primer pair selectively amplifies the fungal ITS. Using three specimens of bark beetle from different species, we assess the fungal diversity associated with these specimens and the repeatability of these estimates in PCRs conducted with different primer tags. The combined replicates produced 727 fungal Operational Taxonomic Units (OTUs) for the specimen of Hylastes ater, 435 OTUs for Tomicus piniperda, and 294 OTUs for Trypodendron lineatum, while individual PCR reactions produced on average only 229, 54, and 31 OTUs for the three specimens, respectively. Yet, communities from PCR replicates were very similar in pairwise comparisons, in particular when considering species abundance, but differed greatly among the three beetle specimens. Different primer tags or the inclusion of amplicons in separate libraries did not impact the species composition. The ITS2 sequences were identified with the Lowest Common Ancestor approach and correspond to diverse lineages of fungi, including Ophiostomaceae and Leotiomycetes widely found to be tree pathogens. We conclude that Illumina MiSeq metabarcoding reliably captures fungal diversity associated with bark beetles, although numerous PCR replicates are recommended for an exhaustive sample. Direct PCR from beetle DNA extractions provides a rapid method for future surveys of fungal species diversity and their associations with bark beetles and environmental variables.     

Host and tissue variations overshadow the response of boreal moss‐associated fungal communities to increased nitrogen load

Human activity has more than doubled the amount of nitrogen entering the global nitrogen cycle, and the boreal forest biome is a nitrogen‐limited ecosystem sensitive to nitrogen load perturbation. Although bryophyte‐associated microbes contribute significantly to boreal forest ecosystem function, particularly in carbon and nitrogen cycling, little is known about their responses to anthropogenic global change. Amplicon pyrosequencing of the ITS2 region of rDNA was used to investigate how fungal communities associated with three bryophyte species responded to increased nitrogen loads in a long‐term fertilization experiment in a boreal Picea abies forest in southern Norway. Overall, OTU richness, community composition and the relative abundance of specific ecological guilds were primarily influenced by host species identity and tissue type. Although not the primary factor affecting fungal communities, nitrogen addition did impact the abundance of specific guilds of fungi and the resulting overall community composition. Increased nitrogen loads decreased ectomycorrhizal abundance, with Amphinema, Cortinarius, Russula and Tylospora OTUs responding negatively to fertilization. Pathogen abundance increased with fertilization, particularly in the moss pathogen Eocronartium. Saprophytic fungi were both positively and negatively impacted by the nitrogen addition, indicating a complex community level response. The overshadowing of the effects of increased nitrogen loads by variation related to host and tissue type highlights the complexity of bryophyte‐associated microbial communities and the intricate nature of their responses to anthropogenic global change.     

Rich and cold: diversity,distribution and drivers of fungal communities in patterned‐ground ecosystems of the North American Arctic

Fungi are abundant and functionally important in the Arctic, yet comprehensive studies of their diversity in relation to geography and environment are not available. We sampled soils in paired plots along the North American Arctic Transect (NAAT), which spans all five bioclimatic subzones of the Arctic. Each pair of plots contrasted relatively bare, cryoturbated patterned‐ground features (PGFs) and adjacent vegetated between patterned‐ground features (bPGFs). Fungal communities were analysed via sequencing of 7834 ITS‐LSU clones. We recorded 1834 OTUs – nearly half the fungal richness previously reported for the entire Arctic. These OTUs spanned eight phyla, 24 classes, 75 orders and 120 families, but were dominated by Ascomycota, with one‐fifth belonging to lichens. Species richness did not decline with increasing latitude, although there was a decline in mycorrhizal taxa that was offset by an increase in lichen taxa. The dominant OTUs were widespread even beyond the Arctic, demonstrating no dispersal limitation. Yet fungal communities were distinct in each subzone and were correlated with soil pH, climate and vegetation. Communities in subzone E were distinct from the other subzones, but similar to those of the boreal forest. Fungal communities on disturbed PGFs differed significantly from those of paired stable areas in bPGFs. Indicator species for PGFs included lichens and saprotrophic fungi, while bPGFs were characterized by ectomycorrhizal and pathogenic fungi. Our results suggest that the Arctic does not host a unique mycoflora, while Arctic fungi are highly sensitive to climate and vegetation, with potential to migrate rapidly as global change unfolds.     

Molecular Identification of Ectomycorrhizal Mycelium in Soil Horizons

Molecular identification techniques based on total DNA extraction provide a unique tool for identification of mycelium in soil. Using molecular identification techniques, the ectomycorrhizal (EM) fungal community under coniferous vegetation was analyzed. Soil samples were taken at different depths from four horizons of a podzol profile. A basidiomycete-specific primer pair (ITS1F-ITS4B) was used to amplify fungal internal transcribed spacer (ITS) sequences from total DNA extracts of the soil horizons. Amplified basidiomycete DNA was cloned and sequenced, and a selection of the obtained clones was analyzed phylogenetically. Based on sequence similarity, the fungal clone sequences were sorted into 25 different fungal groups, or operational taxonomic units (OTUs). Out of 25 basidiomycete OTUs, 7 OTUs showed high nucleotide homology (≥99%) with known EM fungal sequences and 16 were found exclusively in the mineral soil. The taxonomic positions of six OTUs remained unclear. OTU sequences were compared to sequences from morphotyped EM root tips collected from the same sites. Of the 25 OTUs, 10 OTUs had ≥98% sequence similarity with these EM root tip sequences. The present study demonstrates the use of molecular techniques to identify EM hyphae in various soil types. This approach differs from the conventional method of EM root tip identification and provides a novel approach to examine EM fungal communities in soil.     

Molecular identification of ectomycorrhizal mycelium in soil horizons

Molecular identification techniques based on total DNA extraction provide a unique tool for identification of mycelium in soil. Using molecular identification techniques, the ectomycorrhizal (EM) fungal community under coniferous vegetation was analyzed. Soil samples were taken at different depths from four horizons of a podzol profile. A basidiomycete-specific primer pair (ITS1F-ITS4B) was used to amplify fungal internal transcribed spacer (ITS) sequences from total DNA extracts of the soil horizons. Amplified basidiomycete DNA was cloned and sequenced, and a selection of the obtained clones was analyzed phylogenetically. Based on sequence similarity, the fungal clone sequences were sorted into 25 different fungal groups, or operational taxonomic units (OTUs). Out of 25 basidiomycete OTUs, 7 OTUs showed high nucleotide homology (> or = 99%) with known EM fungal sequences and 16 were found exclusively in the mineral soil. The taxonomic positions of six OTUs remained unclear. OTU sequences were compared to sequences from morphotyped EM root tips collected from the same sites. Of the 25 OTUs, 10 OTUs had > or = 98% sequence similarity with these EM root tip sequences. The present study demonstrates the use of molecular techniques to identify EM hyphae in various soil types. This approach differs from the conventional method of EM root tip identification and provides a novel approach to examine EM fungal communities in soil.     

Fungal endophyte communities differ between chestnut galls and surrounding foliar tissues

Foliar endophytic fungi are present in almost all vascular plants. The composition of endophyte communities varies among plant individuals. Likely, but understudied, sources of this variation are the species composition of the plant community and initial attacks by insect herbivores. We addressed these issues by characterizing fungal endophyte communities on leaves of chestnut (Castanea sativa) grown in pure vs. mixed stands. We used ITS metabarcoding methods to identify endophytic fungi associated with galls caused by the invasive gall wasp, Dryocosmus kuriphilus, and with surrounding chestnut leaf tissues. We found 1378 different OTUs. The richness, diversity and composition of endophyte communities differed between galls and surrounding leaf tissues but were independent of forest stand composition. Fungal endophyte richness was lower in galls than in surrounding leaf tissues. Most differences in the composition of fungal endophyte communities between galls and foliar tissues were due to OTU turnover. These results suggest that insect-induced galls provide a particular habitat condition for endophytic microorganisms, regardless of forest species composition. A better understanding of endophyte biology is important to improve their use as biocontrol agents of galling insects.     

Species Composition of Saproxylic Fungal Communities on Decaying Logs in the Boreal Forest

Coarse woody debris supports large numbers of saproxylic fungal species. However, most of the current knowledge comes from Scandinavia and studies relating the effect of stand or log characteristics on the diversity and composition of decomposer fungi have not been conducted in Northeastern Canada. Logs from five tree species were sampled along a decomposition gradient in nine stands representing three successional stages of the boreal mixed forest of Northwestern Quebec, Canada. Using a molecular fingerprinting technique, we assessed fungal community Shannon–Weaver diversity index, richness, and composition. We used linear mixed models and multivariate analyses to link changes in fungal communities to log and stand characteristics. We found a total of 33 operational taxonomic units (OTUs) including an indicator species for balsam fir (similar to Athelia sp.) and one found only in aspen stands (similar to Calocera cornea). Spruce logs supported the highest fungal Shannon–Weaver diversity index and OTU number. Our results support the hypothesis that log species influences fungal richness and diversity. However, log decay class does not. Stand composition, volume of coarse woody debris, and log chemical composition were all involved in structuring fungal communities. Maintaining the diversity of wood-decomposing communities therefore requires the presence of dead wood from diverse log species.     

Characterization of fungal communities associated with the bark beetle Ips typographus varies depending on detection method, location, and beetle population levels

The Eurasian spruce bark beetle Ips typographus and their fungal associates can cause severe damage to Norway spruce forests. In this paper, by using both molecular and cultural methods, we compared fungal assemblages on bark beetles from different locations, characterized by different beetle population levels. Ips typographus was trapped in the western Alps in two outbreak and in two control areas. Sequencing of clone libraries of Internal Transcribed Spacer (ITS) identified 31 fungal Operational Taxonomic Units (OTUs), while fungal isolations yielded 55 OTUs. Only three OTUs were detected by both molecular and cultural methods indicating that both methods are necessary to adequately describe fungal richness. Fungal assemblages on insects from these four and from an additional 12 study sites differed among stands in response to varying ecological conditions and to the limited spreading ability of I. typographus. Ophiostomatoid fungi showed higher diversity in outbreak areas; the pathogenic Ophiostoma polonicum was relatively uncommon, while O. bicolor was the most abundant species. This result was not unexpected, as insects were trapped not at the peak but at the end of the outbreaks and supports the hypothesis of a temporal succession among Ophiostoma species. Ubiquitous endophytes of trees or common airborne fungi were present both in outbreak and in control areas. Wood decaying basidiomycetes were almost never detected on beetles. Yeasts were detected only by molecular analysis, and the OTUs detected matched those reported elsewhere in Europe and in the world, suggesting a very long association between some yeasts and bark beetles.     

Deep Sequencing and Ecological Characterization of Gut Microbial Communities of Diverse Bumble Bee Species

Gut bacterial communities of bumble bees are correlated with defense against pathogens. Further understanding this host-microbe association is vitally important as bumble bees are currently experiencing global population declines, potentially due in part to emergent diseases. In this study, we used pyrosequencing and community fingerprinting (ARISA) to characterize the gut microbial communities of nine bumble species from across the Bombus phylogeny. Overall, we delimited 74 bacterial taxa (operational taxonomic units or OTUs) belonging to Betaproteobacteria, Gammaproteobacteria, Bacilli, Actinobacteria, Flavobacteria and Alphaproteobacteria. Each bacterial community was taxonomically simple, containing an average of 1.9 common (relative abundance per sample > 5%) bacterial OTUs. The most abundant and prevalent (occurring in 92% of the samples) bacterial OTU, based on 16S rRNA sequences, closely matched that of the previously described Betaproteobacteria species Snodgrassella alvi. Bacteria that were first described in bee-related external environments dominated a number of gut bacterial communities, suggesting that they are not strictly dependent on the internal gut environment. The ARISA data showed a correlation between bacterial community structures and the geographic locations where the bees were sampled, suggesting that at least a subset of the bacterial species may be transmitted environmentally. Using light and fluorescent microscopy, we demonstrated that the gut bacteria form a biofilm on the internal epithelial surface of the ileum, corroborating results obtained from Apis mellifera.     

Comparison of the diversity,composition, and host recurrence of xylariaceous endophytes in subtropical,cool temperate,and subboreal regions in Japan

The diversity, composition, and host recurrence of endophytic fungi in the Xylariaceae were compared in subtropical (ST), cool temperate (CT), and subboreal forests (SB) in Japan based on the 28S ribosomal DNA sequences from fungal isolates. A total of 610 isolates were obtained from the leaves of 167 tree species in three sites, which were classified into 42 operational taxonomic units (OTUs) at the 99 % similarity level of the 28S rDNA sequence. ST, CT, and SB yielded 31, 13, and three OTUs, respectively. The OTU richness, diversity, and evenness of fungal communities were in the order: ST > CT > SB. The 42 OTUs were assigned to nine genera in the Xylariaceae: Xylaria, Annulohypoxylon, Anthostomella, Biscogniauxia, Nemania, Hypoxylon, Muscodor, Daldinia, and Rosellinia. Xylarioid isolates in the subfamily Xylarioideae outnumbered Hypoxyloid isolates in the subfamily Hypoxyloideae in ST and CT, whereas the opposite was found in SB. Sørensen’s quotient of similarity was generally low between the three sites. Host recurrence of fungal OTUs was evaluated with the degree of specialization of interaction network between xylariaceous endophytes and plant species and compared between the three sites. We found that the networks in the three sites showed a significantly higher degree of specialization than simulated networks, where partners were associated randomly. Permutational multivariate analyses of variance indicated that plant family and leaf trait significantly affected the OTU composition in ST, which can account for the specialization of interaction network and host recurrence of xylariaceous endophytes.     

Effects of operational taxonomic unit inference methods on soil microeukaryote community analysis using long‐read metabarcoding

Long amplicon metabarcoding has opened the door for phylogenetic analysis of the largely unknown communities of microeukaryotes in soil. Here, we amplified and sequenced the ITS and LSU regions of the rDNA operon (around 1500 bp) from grassland soils using PacBio SMRT sequencing. We tested how three different methods for generation of operational taxonomic units (OTUs) effected estimated richness and identified taxa, and how well large‐scale ecological patterns associated with shifting environmental conditions were recovered in data from the three methods. The field site at Kungsängen Nature Reserve has drawn frequent visitors since Linnaeus''s time, and its species rich vegetation includes the largest population of Fritillaria meleagris in Sweden. To test the effect of different OTU generation methods, we sampled soils across an abrupt moisture transition that divides the meadow community into a Carex acuta dominated plant community with low species richness in the wetter part, which is visually distinct from the mesic‐dry part that has a species rich grass‐dominated plant community including a high frequency of F. meleagris. We used the moisture and plant community transition as a framework to investigate how detected belowground microeukaryotic community composition was influenced by OTU generation methods. Soil communities in both moisture regimes were dominated by protists, a large fraction of which were taxonomically assigned to Ciliophora (Alveolata) while 30%–40% of all reads were assigned to kingdom Fungi. Ecological patterns were consistently recovered irrespective of OTU generation method used. However, different methods strongly affect richness estimates and the taxonomic and phylogenetic resolution of the characterized community with implications for how well members of the microeukaryotic communities can be recognized in the data.     

Global patterns of diversity and community structure in marine bacterioplankton

Because of their small size, great abundance and easy dispersal, it is often assumed that marine planktonic microorganisms have a ubiquitous distribution that prevents any structured assembly into local communities. To challenge this view, marine bacterioplankton communities from coastal waters at nine locations distributed world-wide were examined through the use of comprehensive clone libraries of 16S ribosomal RNA genes, used as operational taxonomic units (OTU). Our survey and analyses show that there were marked differences in the composition and richness of OTUs between locations. Remarkably, the global marine bacterioplankton community showed a high degree of endemism, and conversely included few cosmopolitan OTUs. Our data were consistent with a latitudinal gradient of OTU richness. We observed a positive relationship between the relative OTU abundances and their range of occupation, i.e. cosmopolitans had the largest population sizes. Although OTU richness differed among locations, the distributions of the major taxonomic groups represented in the communities were analogous, and all local communities were similarly structured and dominated by a few OTUs showing variable taxonomic affiliations. The observed patterns of OTU richness indicate that similar evolutionary and ecological processes structured the communities. We conclude that marine bacterioplankton share many of the biogeographical and macroecological features of macroscopic organisms. The general processes behind those patterns are likely to be comparable across taxa and major global biomes.     

The rpb2 gene represents a viable alternative molecular marker for the analysis of environmental fungal communities

Although the commonly used internal transcribed spacer region of rDNA (ITS) is well suited for taxonomic identification of fungi, the information on the relative abundance of taxa and diversity is negatively affected by the multicopy nature of rDNA and the existence of ITS paralogues. Moreover, due to high variability, ITS sequences cannot be used for phylogenetic analyses of unrelated taxa. The part of single‐copy gene encoding the second largest subunit of RNA polymerase II (rpb2) was thus compared with first spacer of ITS as an alternative marker for the analysis of fungal communities in spruce forest topsoil, and their applicability was tested on a comprehensive mock community. In soil, rpb2 exhibited broad taxonomic coverage of the entire fungal tree of life including basal fungal lineages. The gene exhibited sufficient variation for the use in phylogenetic analyses and taxonomic assignments, although it amplifies also paralogues. The fungal taxon spectra obtained with rbp2 region and ITS1 corresponded, but sequence abundance differed widely, especially in the basal lineages. The proportions of OTU counts and read counts of major fungal groups were close to the reality when rpb2 was used as a molecular marker while they were strongly biased towards the Basidiomycota when using the ITS primers ITS1/ITS4. Although the taxonomic placement of rbp2 sequences is currently more difficult than that of the ITS sequences, its discriminative power, quantitative representation of community composition and suitability for phylogenetic analyses represent significant advantages.     

Polymerase matters: non-proofreading enzymes inflate fungal community richness estimates by up to 15 %

Rare taxa overwhelm metabarcoding data generated using next-generation sequencing (NGS). Low frequency Operational Taxonomic Units (OTUs) may be artifacts generated by PCR-amplification errors resulting from polymerase mispairing. We analyzed two Internal Transcribed Spacer 2 (ITS2) MiSeq libraries generated with proofreading (ThermoScientific Phusion^®) and non-proofreading (ThermoScientific Phire^®) polymerases from the same MiSeq reaction, the same samples, using the same DNA tags, and with two different clustering methods to evaluate the effect of polymerase and clustering tool choices on the estimates of richness, diversity and community composition. Our data show that, while the overall communities are comparable, OTU richness is exaggerated by the use of the non-proofreading polymerase–up to 15 % depending on the clustering method, and on the threshold of low frequency OTU removal. The overestimation of richness also consistently led to underestimation of community evenness, a result of increase in the low frequency OTUs. Stringent thresholds of eliminating the rare reads remedy this issue; exclusion of reads that occurred ≤10 times reduced overestimated OTU numbers to <0.3 %. As a result of these findings, we strongly recommend the use of proofreading polymerases to improve the data integrity as well as the use of stringent culling thresholds for rare sequences to minimize overestimation of community richness.     

Molecular and microscopic analysis of the gut contents of abundant rove beetle species (Coleoptera,?Staphylinidae) in the boreal balsam fir forest of Quebec,?Canada

Experimental research on beetle responses to removal of logging residues following clearcut harvesting in the boreal balsam fir forest of Quebec revealed several abundant rove beetle (Staphylinidae) species potentially important for long-term monitoring. To understand the trophic affiliations of these species in forest ecosystems, it was necessary to analyze their gut contents. We used microscopic and molecular (DNA) methods to identify the gut contents of the following rove beetles: Atheta capsularis Klimaszewski, Atheta klagesi Bernhauer, Oxypoda grandipennis (Casey), Bryophacis smetanai Campbell, Ischnosoma longicorne (Mäklin), Mycetoporus montanus Luze, Tachinus frigidus Erichson, Tachinus fumipennis (Say), Tachinus quebecensis Robert, and Pseudopsis subulata Herman. We found no apparent arthropod fragments within the guts; however, a number of fungi were identified by DNA sequences, including filamentous fungi and budding yeasts [Ascomycota: Candida derodonti Suh & Blackwell (accession number FJ623605), Candida mesenterica (Geiger) Diddens & Lodder (accession number FM178362), Candida railenensis Ramirez and Gonzáles (accession number JX455763), Candida sophie-reginae Ramirez & González (accession number HQ652073), Candida sp. (accession number AY498864), Pichia delftensis Beech (accession number AY923246), Pichia membranifaciens Hansen (accession number JQ26345), Pichia misumaiensis Y. Sasaki and Tak. Yoshida ex Kurtzman 2000 (accession number U73581), Pichia sp. (accession number AM261630), Cladosporium sp. (accession number KF367501), Acremoniumpsammosporum W. Gams (accession number GU566287), Alternaria sp. (accession number GU584946), Aspergillus versicolor Bubak (accession number AJ937750), and Aspergillusamstelodami (L. Mangin) Thom and Church (accession number HQ728257)]. In addition, two species of bacteria [Bradyrhizobium japonicum (Kirchner) Jordan (accession number BA000040) and Serratia marcescens Bizio accession number CP003942] were found in the guts. These results not only provide evidence of the consumer-resource relations of these beetles but also clarify the relationship between rove beetles, woody debris and fungi. Predominance of yeast-feeding by abundant rove beetles suggests that it may play an important role in their dietary requirements.