Metagenomes provide valuable comparative information on soil microeukaryotes

Ciliate diversity in rodrigo de freitas lagoon (Rio de Janeiro, Brazil) from an integrative standpoint

The Rodrigo de Freitas Lagoon is a highly eutrophic lacustrine system and has one of the longest histories of exploration and anthropic alteration in Brazil. Despite its relevance, limited studies explored the diversity of micro-eukaryotes in the lagoon. Ciliates (Alveolata, Ciliophora) are overlooked in environmental microbiology, especially in tropical and subtropical ecosystems, resulting in limited knowledge about their diversity and functional relevance in South American habitats, particularly in coastal lagoons. To fill this gap, here we investigated the diversity of ciliates in a brackish coastal lagoon in an urban area of Rio de Janeiro, Brazil, applying and comparing the performance of morphological and metabarcoding approaches. The metabarcoding analysis, based on high-throughput sequencing of the hipervariable region V4 of the 18S rRNA genes detected 37 molecular operational taxonomic units (MOTUs) assigned to Ciliophora, representing only about a half (56.9%) of the diversity detected by microscopy, which counted 65 ciliate morphotypes. The most representative classes in both approaches were Spirotrichea and Oligohymenophorea. The metabarcoding analysis revealed that 35.3% of the ciliate MOTUs had less than 97% similarity to available sequences in the NCBI database, indicating that more than one-third of these MOTUs potentially represents still not represented or undescribed ciliate species in current databases. Our findings indicate that metabarcoding techniques can significantly enhance the comprehension of ciliate diversity in tropical environments, but the scarcity of reference sequences of brackish ciliates in molecular databases represents a challenge to the taxonomic assignment of the MOTUs. This study provides new insights into the diversity of ciliates in a threatened coastal lagoon, revealing a vast array of still unknown and rare ciliate taxonomic units in tropical environments.

Microbial dynamics in soils of the Damma glacier forefield show succession in the functional genetic potential

Glacier retreat is a visible consequence of climate change worldwide. Although taxonomic change of the soil microbiomes in glacier forefields have been widely documented, how microbial genetic potential changes along succession is little known. Here, we used shotgun metagenomics to analyse whether the soil microbial genetic potential differed between four stages of soil development (SSD) sampled along three transects in the Damma glacier forefield (Switzerland). The SSDs were characterized by an increasing vegetation cover, from barren soil, to biological soil crust, to sparsely vegetated soil and finally to vegetated soil. Results suggested that SSD significantly influenced microbial genetic potential, with the lowest functional diversity surprisingly occurring in the vegetated soils. Overall, carbohydrate metabolism and secondary metabolite biosynthesis genes overrepresented in vegetated soils, which could be partly attributed to plant-soil feedbacks. For C degradation, glycoside hydrolase genes enriched in vegetated soils, while auxiliary activity and carbohydrate esterases genes overrepresented in barren soils, suggested high labile C degradation potential in vegetated, and high recalcitrant C degradation potential in barren soils. For N-cycling, organic N degradation and synthesis genes dominated along succession, and gene families involved in nitrification were overrepresented in barren soils. Our study provides new insights into how the microbial genetic potential changes during soil formation along the Damma glacier forefield.

Biological Interactions and Environmental Influences Shift Microeukaryotes in Permafrost Active Layer Soil Across the Qinghai-Tibet Plateau

Permafrost active layer soils are harsh environments with thaw/freeze cycles and sub-zero temperatures, harboring diverse microorganisms. However, the distribution patterns, assembly mechanism, and driving forces of soil microeukaryotes in permafrost remain largely unknown. In this study, we investigated microeukaryotes in permafrost active layer across the Qinghai-Tibet Plateau (QTP) using 18S rRNA gene sequencing. The results showed that the microbial eukaryotic communities were dominated by Nematozoa, Ciliophora, Ascomycota, Cercozoa, Arthropoda, and Basidiomycota in terms of relative abundance and operational taxonomic unit (OTU) richness. Nematozoa had the highest relative abundance, while Ciliophora had the highest OTU richness. These phyla had strong interactions between each other. Their alpha diversity and community structure were differently influenced by the factors associated to location, climate, and soil properties, particularly the soil properties. Significant but weak distance-decay relationships with different slopes were established for the communities of these dominant phyla, except for Basidiomycota. According to the null model, community assemblies of Nematozoa and Cercozoa were dominated by heterogeneous selection, Ciliophora and Ascomycota were dominated by dispersal limitation, while Arthropoda and Basidiomycota were highly dominated by non-dominant processes. The assembly mechanisms can be jointly explained by biotic interactions, organism treats, and environmental influences. Modules in the co-occurrence network of the microeukaryotes were composed by members from different taxonomic groups. These modules also had interactions and responded to different environmental factors, within which, soil properties had strong influences on these modules. The results suggested the importance of biological interactions and soil properties in structuring microbial eukaryotic communities in permafrost active layer soil across the QTP.

Soil filtration-sedimentation improves shelled protist recovery in eukaryotic eDNA surveys

A large part of the soil protist diversity is missed in metabarcoding studies based on 0.25 g of soil environmental DNA (eDNA) and universal primers due to ca. 80% co-amplification of non-target plants, animals and fungi. To overcome this problem, enrichment of the substrate used for eDNA extraction is an easily implemented option but its effect has not yet been tested. In this study, we evaluated the effect of a 150 μm mesh size filtration and sedimentation method to improve the recovery of protist eDNA, while reducing the co-extraction of plant, animal and fungal eDNA, using a set of contrasted forest and alpine soils from La Réunion, Japan, Spain and Switzerland. Total eukaryotic diversity was estimated by V4 18S rRNA metabarcoding and classical amplicon sequence variant calling. A 2- to 3-fold enrichment in shelled protists (Euglyphida, Arcellinida and Chrysophyceae) was observed at the sample level with the proposed method, with, at the same time, a 2-fold depletion of Fungi and a 3-fold depletion of Embryophyceae. Protist alpha diversity was slightly lower in filtered samples due to reduced coverage in Variosea and Sarcomonadea, but significant differences were observed in only one region. Beta diversity varied mostly between regions and habitats, which explained the same proportion of variance in bulk soil and filtered samples. The increased resolution in soil protist diversity estimates provided by the filtration-sedimentation method is a strong argument in favour of including it in the standard protocol for soil protist eDNA metabarcoding studies.

Response of soil protozoa to acid mine drainage in a contaminated terrace

Acid mine drainage (AMD) system represents one of the most unfavorable habitats for microorganisms due to its low pH and high concentrations of metals. Compared to bacteria and fungi, our understanding regarding the response of soil protozoa to such extremely acidic environments remains limited. This study characterized the structures of protozoan communities inhabiting a terrace heavily contaminated by AMD. The sharp environmental gradient of this terrace was generated by annual flooding from an AMD lake located below, which provided a natural setting to unravel the environment-protozoa interactions. Previously unrecognized protozoa, such as Apicomplexa and Euglenozoa, dominated the extremely acidic soils, rather than the commonly recognized members (e.g., Ciliophora and Cercozoa). pH was the most important factor regulating the abundance of protozoan taxa. Metagenomic analysis of protozoan metabolic potential showed that many functional genes encoding for the alleviation of acid stress and various metabolic pathways were enriched, which may facilitate the survival and adaptation of protozoa to acidic environments. In addition, numerous co-occurrences between protozoa and bacterial or fungal taxa were observed, suggesting shared environmental preferences or potential bio-interactions among them. Future studies are required to confirm the ecological roles of these previously unrecognized protozoa as being important soil microorganisms.

Effects of Rotations With Legume on Soil Functional Microbial Communities Involved in Phosphorus Transformation

Including legumes in the cereal cropping could improve the crop yield and the uptake of nitrogen (N) and phosphorus (P) of subsequent cereals. The effects of legume-cereal crop rotations on the soil microbial community have been studied in recent years, the impact on soil functional genes especially involved in P cycling is raising great concerns. The metagenomic approach was used to investigate the impacts of crop rotation managements of soybean-wheat (SW) and maize-wheat (MW) lasting 2 and 7years on soil microbial communities and genes involved in P transformation in a field experiment. Results indicated that SW rotation increased the relative abundances of Firmicutes and Bacteroidetes, reduced Actinobacteria, Verrucomicrobia, and Chloroflexi compared to MW rotation. gcd, phoR, phoD, and ppx predominated in genes involved in P transformation in both rotations. Genes of gcd, ppa, and ugpABCE showed higher abundances in SW rotation than in MW rotation, whereas gadAC and pstS showed less abundances. Proteobacteria, Acidobacteria, and Gemmatimonadetes played predominant roles in microbial P cycling. Our study provides a novel insight into crop P, which requires strategy and help to understand the mechanism of improving crop nutrient uptake and productivity in different rotations.

Changes in Soil Prokaryotic Diversity in Response to Land-Use Changes in Sub-Saharan Africa

Sub-Saharan Africa is one of the most severely affected regions regarding soil degradation, a global issue with the loss of nutrients caused by inappropriate management, leading to low agricultural productivity. Here we asked the question of how soil prokaryotic communities are affected by shifts in land use management and subsequent losses in soil organic carbon. We sampled soils from three sites in Zambia which have neighboring natural and managed sites. After the measurement of soil properties, soil DNA was sequenced, targeting the 16S rRNA gene. As expected, total carbon in soil was decreased in the managed sites, with significant reductions of bacterial biomass. However, the diversity indices in the managed soils were higher than in natural soils. Particularly, the relative abundance of nitrifiers was increased in the managed soils, most likely as a result of fertilization. However also other bacteria, e.g., those which formed tight interactions with the cultivated crops including the genera Balneimonas, and Bacillus, were increased in the managed soils. In contrast bacteria belonging to the family Chloroflexi, which were high in abundance in the natural soil were outcompeted by other prokaryotes in the managed soils most likely as a result of changes in the amount of soil organic carbon. Overall, our results suggest that we need to discuss the trends of prokaryotic diversity separately from those for prokaryotic abundance. Even when bacterial abundances were decreased in the managed soils, nitrifiers’ relative abundance and diversity increased in our experiment, suggesting the possible alteration of the nitrogen cycle in managed soils in sub-Saharan Africa.

All together now: Limitations and recommendations for the simultaneous analysis of all eukaryotic soil sequences

The soil environment contains a large, but historically underexplored, reservoir of biodiversity. Sequencing prokaryotic marker genes has become commonplace for the discovery and characterization of soil bacteria and archaea. Increasingly, this approach is also applied to eukaryotic marker genes to characterize the diversity and distribution of soil eukaryotes. However, understanding the properties and limitations of eukaryotic marker sequences is essential for correctly analysing, interpreting, and synthesizing the resulting data. Here, we illustrate several biases from sequencing data that affect measurements of biodiversity that arise from variation in morphology, taxonomy and phylogeny between organisms, as well as from sampling designs. We recommend analytical approaches to overcome these limitations, and outline how the benchmarking and standardization of sequencing protocols may improve the comparability of the data.

Geochemical and microbial characteristics of seepage water and mineral precipitates in a radwaste disposal facility impacted by seawater intrusion and high alkalinity

The construction of an underground facility can dramatically change the quality, flow direction, and level of groundwater. It may also impact subsurface microbial composition and activity. Groundwater quality was monitored over eight years in two observational wells near an underground disposal facility on the east coast of South Korea. The results showed dramatic increases in dissolved ions such as O₂, Na, Ca, Mg, and SO₄ during facility construction. Seepage water samples downgradient from the silos and tunnels, and precipitates deposited along the seepage water flow path were collected to determine the impact inside the disposal facility. X-ray analysis (powder X-ray diffraction (pXRD) and X-ray absorption fine structure (XAFS)) were used to characterize the mineral precipitates. Microbial community composition was determined by 16S rRNA gene sequencing. The seepage water composition was of two types: Ca-Cl and Ca-Na-HCO₃. The ratio of Cl and δ¹⁸O showed that the Ca-Cl type seepage water was influenced by groundwater mixed with seawater ranging from 2.7% to 15.1%. Various sulfate-reducing bacteria were identified in the Ca-Cl type seepage water, exhibiting relatively high sulfate content from seawater intrusion. Samples from the Ca-Na-HCO₃ type seepage water had an extremely high pH (>10) and abundance of Hydrogenophaga. The precipitates observed along the flow path of the seepage water included calcite, ferrihydrite, green rust, and siderite, depending on seepage water chemistry and microbial activity. This study suggests that the construction of underground structures creates distinct, localized geochemical conditions (e.g., high alkalinity, high salinity, and oxic conditions), which may impact microbial communities. These biogeochemical changes may have undesirable large-scale impacts such as water pump clogging. An understanding of the process and long-term monitoring are essential to assess the safety of underground facilities.

High temperatures enhance the microbial genetic potential to recycle C and N from necromass in high-mountain soils

Climate change is strongly affecting high-mountain soils and warming in particular is associated with pronounced changes in microbe-mediated C and N cycling, affecting plant-soil interactions and greenhouse gas balances and therefore feedbacks to global warming. We used shotgun metagenomics to assess changes in microbial community structures, as well as changes in microbial C- and N-cycling potential and stress response genes and we linked these data with changes in soil C and N pools and temperature-dependent measurements of bacterial growth rates. We did so by incubating high-elevation soil from the Swiss Alps at 4°C, 15°C, 25°C, or 35°C for 1 month. We found no shift with increasing temperature in the C-substrate-degrader community towards taxa more capable of degrading recalcitrant organic matter. Conversely, at 35°C, we found an increase in genes associated with the degradation and modification of microbial cell walls, together with high bacterial growth rates. Together, these findings suggest that the rapidly growing high-temperature community is fueled by necromass from heat-sensitive taxa. This interpretation was further supported by a shift in the microbial N-cycling potential towards N mineralization and assimilation under higher temperatures, along with reduced potential for conversions among inorganic N forms. Microbial stress-response genes reacted inconsistently to increasing temperature, suggesting that the high-temperature community was not severely stressed by these conditions. Rather, soil microbes were able to acclimate by changing the thermal properties of membranes and cell walls as indicated by an increase in genes involved in membrane and cell wall modifications as well as a shift in the optimum temperature for bacterial growth towards the treatment temperature. Overall, our results suggest that high temperatures, as they may occur with heat waves under global warming, promote a highly active microbial community capable of rapid mineralization of microbial necromass, which may transiently amplify warming effects.

Different Degrees of Niche Differentiation for Bacteria, Fungi, and Myxomycetes Within an Elevational Transect in the German Alps

We used direct DNA amplification from soil extracts to analyze microbial communities from an elevational transect in the German Alps by parallel metabarcoding of bacteria (16S rRNA), fungi (ITS2), and myxomycetes (18S rRNA). For the three microbial groups, 5710, 6133, and 261 operational taxonomic units (OTU) were found. For the latter group, we can relate OTUs to barcodes from fruit bodies sampled over a 4-year period. The alpha diversity of myxomycetes was positively correlated with that of bacteria. Vegetation type was found to be the main explanatory parameter for the community composition of all three groups and a substantial species turnover with elevation was observed. Bacteria and fungi display similar community responses, driven by symbiont species and plant substrate quality. Myxamoebae show a more patchy distribution, though still clearly stratified between taxa, which seems to be a response to both structural properties of the habitat and interaction with specific bacterial and fungal taxa. Finally, we report a high number of myxomycete OTUs not represented in a reference database from fructifications, which might represent novel species.

Bacterial and protozoan dynamics upon thawing and freezing of an active layer permafrost soil

The active layer of soil overlaying permafrost in the Arctic is subjected to annual changes in temperature and soil chemistry, which we hypothesize to affect the overall soil microbial community. We investigated changes in soil microorganisms at different temperatures during warming and freezing of the active layer soil from Svalbard, Norway. Soil community data were obtained by direct shotgun sequencing of total extracted RNA. No changes in soil microbial communities were detected when warming from -10 to -2 °C or when freezing from -2 to -10 °C. In contrast, within a few days we observed changes when warming from -2 to +2 °C with a decrease in fungal rRNA and an increase in several OTUs belonging to Gemmatimonadetes, Bacteroidetes and Betaproteobacteria. Even more substantial changes occurred when incubating at 2 °C for 16 days, with declines in total fungal potential activity and decreases in oligotrophic members from Actinobacteria and Acidobacteria. Additionally, we detected an increase in transcriptome sequences of bacterial phyla Bacteriodetes, Firmicutes, Betaproteobacteria and Gammaproteobacteria-collectively presumed to be copiotrophic. Furthermore, we detected an increase in putative bacterivorous heterotrophic flagellates, likely due to predation upon the bacterial community via grazing. Although this grazing activity may explain relatively large changes in the bacterial community composition, no changes in total 16S rRNA gene copy number were observed and the total RNA level remained stable during the incubation. Together, these results are showing the first comprehensive ecological evaluation across prokaryotic and eukaryotic microbial communities on thawing and freezing of soil by application of the TotalRNA technique.

Long-term soil metal exposure impaired temporal variation in microbial metatranscriptomes and enriched active phages

BACKGROUND

It remains unclear whether adaptation and changes in diversity associated to a long-term perturbation are sufficient to ensure functional resilience of soil microbial communities. We used RNA-based approaches (16S rRNA gene transcript amplicon coupled to shotgun mRNA sequencing) to study the legacy effects of a century-long soil copper (Cu) pollution on microbial activity and composition, as well as its effect on the capacity of the microbial community to react to temporal fluctuations.

RESULTS

Despite evidence of microbial adaptation (e.g., iron homeostasis and avoidance/resistance strategies), increased heterogeneity and richness loss in transcribed gene pools were observed with increasing soil Cu, together with an unexpected predominance of phage mRNA signatures. Apparently, phage activation was either triggered directly by Cu, or indirectly via enhanced expression of DNA repair/SOS response systems in Cu-exposed bacteria. Even though total soil carbon and nitrogen had accumulated with increasing Cu, a reduction in temporally induced mRNA functions was observed. Microbial temporal response groups (TRGs, groups of microbes with a specific temporal response) were heavily affected by Cu, both in abundance and phylogenetic composition.

CONCLUSION

Altogether, results point toward a Cu-mediated "decoupling" between environmental fluctuations and microbial activity, where Cu-exposed microbes stopped fulfilling their expected contributions to soil functioning relative to the control. Nevertheless, some functions remained active in February despite Cu, concomitant with an increase in phage mRNA signatures, highlighting that somehow, microbial activity is still happening under these adverse conditions.

A four year survey reveals a coherent pattern between occurrence of fruit bodies and soil amoebae populations for nivicolous myxomycetes

Among soil-inhabiting protists, myxomycetes stand out by their macroscopic fructifications which have allowed studies on their ecology and distribution for more than two hundred years. One of the most distinct ecological guilds in myxomycetes are the nivicolous or "snowbank" myxomycete species, which produce fruit bodies at the edge of melting snowbanks in spring. Relationship between the occurrence of fructifications and myxamoebae remain unknown. In this study we used modern molecular techniques, by direct DNA amplification from soil extracts (NGS metabarcoding) to compare the distribution of soil-inhabiting myxamoebae found in 2016 with fructifications from the same sites collected over the course of four years (2013, 2015-17) along an elevational transect in the northern German Alps. A coherent community composition between fructification and soil myxamoebae, though with species-specific differences in relative abundance, was revealed. Although patterns varied among species, myxamoebae were found at both low and high elevations, whereas fruit bodies were mainly found at higher elevations, likely explained by the presence of a stable and long-lasting snow cover. In addition, a year to year comparison of fructification records support the hypothesis that the abundance of fructifications strongly depends on the onset of snowfall in the previous autumn and the soil temperature regime throughout the winter.

Precise species detection of traditional Chinese patent medicine by shotgun metagenomic sequencing

BACKGROUND

Current quality control methods for traditional Chinese patent medicines (TCPMs), e.g., microscopy, thin-layer chromatography (TLC), and high-performance liquid chromatography (HPLC), cannot detect herbal species composition with adequate precision. To address this issue, more effective detection methods should be explored.

HYPOTHESIS/PURPOSE

We hypothesized that shotgun metagenomic sequencing can fulfill the requirements for the species detection of multi-ingredient TCPMs.

METHODS

Longdan Xiegan Wan (LDXGW), once thought to be the chief culprit in aristolochic acid nephropathy (AAN), was selected to establish the method. It was used for both reference and commercial LDXGW samples. The precision authentication of herbal species contained in multi-ingredient TCPM is based on the shotgun metagenomic sequencing of genomic DNA without PCR amplification. Chemical analyses were also conducted as a contrast test.

RESULTS

Over 100 G of raw data was obtained, and this value represented more than 0.75 billion reads. After assembling and filtering all the reads, a total of 261 contigs were obtained, which belonged to the ITS2, psbA-trnH, and matK regions of the reference and commercial samples. Because the homology of the rbcL region was high, it was not analyzed in the HTS data. Bioinformatics analysis indicated that the ITS2 region, as a DNA barcode, showed the highest identification efficiency. It could successfully detect all prescribed species, including four processed herbal ingredients, in the lab-made reference samples. The commercial samples all met the requirements of the Chinese Pharmacopoeia according to the TLC and HPLC tests. However, the shotgun metagenomic sequencing detected the substitution of Akebiae Caulis (Mutong) in the commercial samples, while the chemical analyses could not distinguish.

CONCLUSION

The results highlight that shotgun metagenomic sequencing is a complementary method for the precise species detection of TCPMs.

Drought Legacy Effects on the Composition of Soil Fungal and Prokaryote Communities

It is increasingly acknowledged that climate change is influencing terrestrial ecosystems by increased drought and rainfall intensities. Soil microbes are key drivers of many processes in terrestrial systems and rely on water in soil pores to fulfill their life cycles and functions. However, little is known on how drought and rainfall fluctuations, which affect the composition and structure of microbial communities, persist once original moisture conditions have been restored. Here, we study how simulated short-term drying and re-wetting events shape the community composition of soil fungi and prokaryotes. In a mesocosm experiment, soil was exposed to an extreme drought, then re-wetted to optimal moisture (50% WHC, water holding capacity) or to saturation level (100% WHC). Composition, community structure and diversity of microbes were measured by sequencing ITS and 16S rRNA gene amplicons 3 weeks after original moisture content had been restored. Drying and extreme re-wetting decreased richness of microbial communities, but not evenness. Abundance changes were observed in only 8% of prokaryote OTUs, and 25% of fungal OTUs, whereas all other OTUs did not differ between drying and re-wetting treatments. Two specific legacy response groups (LRGs) were observed for both prokaryotes and fungi. OTUs belonging to the first LRG decreased in relative abundance in soil with a history of drought, whereas OTUs that increased in soil with a history of drought formed a second LRG. These microbial responses were spread among different phyla. Drought appeared to be more important for the microbial community composition than the following extreme re-wetting. 16S profiles were correlated with both inorganic N concentration and basal respiration and ITS profiles correlated with fungal biomass. We conclude that a drying and/or an extreme re-wetting history can persist in soil microbial communities via specific response groups composed of members with broad phylogenetic origins, with possible functional consequences on soil processes and plant species. As a large fraction of OTUs responding to drying and re-wetting belonged to the rare biosphere, our results suggest that low abundant microbial species are potentially important for ecosystem responses to extreme weather events.

Long-term industrial metal contamination unexpectedly shaped diversity and activity response of sediment microbiome

Metal contamination poses serious biotoxicity and bioaccumulation issues, affecting both abiotic conditions and biological activity in ecosystem trophic levels, especially sediments. The MetalEurop foundry released metals directly into the French river "la Deûle" during a century, contaminating sediments with a 30-fold increase compared to upstream unpolluted areas (Férin, Sensée canal). Previous metaproteogenomic work revealed phylogenetically analogous, but functionally different microbial communities between the two locations. However, their potential activity status in situ remains unknown. The present study respectively compares the structures of both total and active fractions of sediment prokaryotic microbiomes by coupling DNA and RNA-based sequencing approaches at the polluted MetalEurop site and its upstream control. We applied the innovative ecological concept of Functional Response Groups (FRGs) to decipher the adaptive tolerance range of the communities through characterization of microbial lifestyles and strategists. The complementing use of DNA and RNA sequencing revealed indications that metals selected for mechanisms such as microbial facilitation via "public-good" providing bacteria, Horizontal Gene Transfer (HGT) and community coalescence, overall resulting in an unexpected higher microbial diversity at the polluted site.

Microbiomes associated with infective stages of root-knot and lesion nematodes in soil

Endoparasitic root-knot (Meloidogyne spp.) and lesion (Pratylenchus spp.) nematodes cause considerable damage in agriculture. Before they invade roots to complete their life cycle, soil microbes can attach to their cuticle or surface coat and antagonize the nematode directly or by induction of host plant defenses. We investigated whether the nematode-associated microbiome in soil differs between infective stages of Meloidogyne incognita and Pratylenchus penetrans, and whether it is affected by variation in the composition of microbial communities among soils. Nematodes were incubated in suspensions of five organically and two integrated horticultural production soils, recovered by sieving and analyzed for attached bacteria and fungi after washing off loosely adhering microbes. Significant effects of the soil type and nematode species on nematode-associated fungi and bacteria were revealed as analyzed by community profiling using denaturing gradient gel electrophoresis. Attached microbes represented a small specific subset of the soil microbiome. Two organic soils had very similar bacterial and fungal community profiles, but one of them was strongly suppressive towards root-knot nematodes. They were selected for deep amplicon sequencing of bacterial 16S rRNA genes and fungal ITS. Significant differences among the microbiomes associated with the two species in both soils suggested specific surface epitopes. Among the 28 detected bacterial classes, Betaproteobacteria, Bacilli and Actinobacteria were the most abundant. The most frequently detected fungal genera were Malassezia, Aspergillus and Cladosporium. Attached microbiomes did not statistically differ between these two soils. However, Malassezia globosa and four fungal species of the family Plectosphaerellaceae, and the bacterium Neorhizobium galegae were strongly enriched on M. incognita in the suppressive soil. In conclusion, the highly specific attachment of microbes to infective stages of phytonematodes in soil suggested an ecological role of this association and might be involved in soil suppressiveness towards them.