Structure and function of the global topsoil microbiome.

Soils harbour some of the most diverse microbiomes on Earth and are essential for both nutrient cycling and carbon storage. To understand soil functioning, it is necessary to model the global distribution patterns and functional gene repertoires of soil microorganisms, as well as the biotic and environmental associations between the diversity and structure of both bacterial and fungal soil communities1-4. Here we show, by leveraging metagenomics and metabarcoding of global topsoil samples (189 sites, 7,560 subsamples), that bacterial, but not fungal, genetic diversity is highest in temperate habitats and that microbial gene composition varies more strongly with environmental variables than with geographic distance. We demonstrate that fungi and bacteria show global niche differentiation that is associated with contrasting diversity responses to precipitation and soil pH. Furthermore, we provide evidence for strong bacterial-fungal antagonism, inferred from antibiotic-resistance genes, in topsoil and ocean habitats, indicating the substantial role of biotic interactions in shaping microbial communities. Our results suggest that both competition and environmental filtering affect the abundance, composition and encoded gene functions of bacterial and fungal communities, indicating that the relative contributions of these microorganisms to global nutrient cycling varies spatially.

Effects of nitrogen deposition on carbon and nutrient cycling along a natural soil acidity gradient as revealed by metagenomics.

Nitrogen (N) deposition and soil acidification are environmental challenges affecting ecosystem functioning, health, and biodiversity, but their effects on functional genes are poorly understood. Here, we utilized metabarcoding and metagenomics to investigate the responses of soil functional genes to N deposition along a natural soil pH gradient. Soil N content was uncorrelated with pH, enabling us to investigate their effects separately. Soil acidity strongly and negatively affected the relative abundances of most cluster of orthologous gene categories of the metabolism supercategory. Similarly, soil acidity negatively affected the diversity of functional genes related to carbon and N but not phosphorus cycling. Multivariate analyses showed that soil pH was the most important factor affecting microbial and functional gene composition, while the effects of N deposition were less important. Relative abundance of KEGG functional modules related to different parts of the studied cycles showed variable responses to soil acidity and N deposition. Furthermore, our results suggested that the diversity-function relationship reported for other organisms also applies to soil microbiomes. Since N deposition and soil pH affected microbial taxonomic and functional composition to a different extent, we conclude that N deposition effects might be primarily mediated through soil acidification in forest ecosystems.

Best practices in metabarcoding of fungi: From experimental design to results.

The development of high-throughput sequencing (HTS) technologies has greatly improved our capacity to identify fungi and unveil their ecological roles across a variety of ecosystems. Here we provide an overview of current best practices in metabarcoding analysis of fungal communities, from experimental design through molecular and computational analyses. By reanalysing published data sets, we demonstrate that operational taxonomic units (OTUs) outperform amplified sequence variants (ASVs) in recovering fungal diversity, a finding that is particularly evident for long markers. Additionally, analysis of the full-length ITS region allows more accurate taxonomic placement of fungi and other eukaryotes compared to the ITS2 subregion. Finally, we show that specific methods for compositional data analyses provide more reliable estimates of shifts in community structure. We conclude that metabarcoding analyses of fungi are especially promising for integrating fungi into the full microbiome and broader ecosystem functioning context, recovery of novel fungal lineages and ancient organisms as well as barcoding of old specimens including type material.

Global patterns in endemicity and vulnerability of soil fungi.

Fungi are highly diverse organisms, which provide multiple ecosystem services. However, compared with charismatic animals and plants, the distribution patterns and conservation needs of fungi have been little explored. Here, we examined endemicity patterns, global change vulnerability and conservation priority areas for functional groups of soil fungi based on six global surveys using a high-resolution, long-read metabarcoding approach. We found that the endemicity of all fungi and most functional groups peaks in tropical habitats, including Amazonia, Yucatan, West-Central Africa, Sri Lanka, and New Caledonia, with a negligible island effect compared with plants and animals. We also found that fungi are predominantly vulnerable to drought, heat and land-cover change, particularly in dry tropical regions with high human population density. Fungal conservation areas of highest priority include herbaceous wetlands, tropical forests, and woodlands. We stress that more attention should be focused on the conservation of fungi, especially root symbiotic arbuscular mycorrhizal and ectomycorrhizal fungi in tropical regions as well as unicellular early-diverging groups and macrofungi in general. Given the low overlap between the endemicity of fungi and macroorganisms, but high conservation needs in both groups, detailed analyses on distribution and conservation requirements are warranted for other microorganisms and soil organisms.

Perspectives and Benefits of High-Throughput Long-Read Sequencing in Microbial Ecology.

Short-read, high-throughput sequencing (HTS) methods have yielded numerous important insights into microbial ecology and function. Yet, in many instances short-read HTS techniques are suboptimal, for example, by providing insufficient phylogenetic resolution or low integrity of assembled genomes. Single-molecule and synthetic long-read (SLR) HTS methods have successfully ameliorated these limitations. In addition, nanopore sequencing has generated a number of unique analysis opportunities, such as rapid molecular diagnostics and direct RNA sequencing, and both Pacific Biosciences (PacBio) and nanopore sequencing support detection of epigenetic modifications. Although initially suffering from relatively low sequence quality, recent advances have greatly improved the accuracy of long-read sequencing technologies. In spite of great technological progress in recent years, the long-read HTS methods (PacBio and nanopore sequencing) are still relatively costly, require large amounts of high-quality starting material, and commonly need specific solutions in various analysis steps. Despite these challenges, long-read sequencing technologies offer high-quality, cutting-edge alternatives for testing hypotheses about microbiome structure and functioning as well as assembly of eukaryote genomes from complex environmental DNA samples.

Temperature and pH define the realised niche space of arbuscular mycorrhizal fungi.

The arbuscular mycorrhizal (AM) fungi are a globally distributed group of soil organisms that play critical roles in ecosystem function. However, the ecological niches of individual AM fungal taxa are poorly understood. We collected > 300 soil samples from natural ecosystems worldwide and modelled the realised niches of AM fungal virtual taxa (VT; approximately species-level phylogroups). We found that environmental and spatial variables jointly explained VT distribution worldwide, with temperature and pH being the most important abiotic drivers, and spatial effects generally occurring at local to regional scales. While dispersal limitation could explain some variation in VT distribution, VT relative abundance was almost exclusively driven by environmental variables. Several environmental and spatial effects on VT distribution and relative abundance were correlated with phylogeny, indicating that closely related VT exhibit similar niche optima and widths. Major clades within the Glomeraceae exhibited distinct niche optima, Acaulosporaceae generally had niche optima in low pH and low temperature conditions, and Gigasporaceae generally had niche optima in high precipitation conditions. Identification of the realised niche space occupied by individual and phylogenetic groups of soil microbial taxa provides a basis for building detailed hypotheses about how soil communities respond to gradients and manipulation in ecosystems worldwide.

Diversity and substrate-specificity of green algae and other micro-eukaryotes colonizing amphibian clutches in Germany, revealed by DNA metabarcoding.

Amphibian clutches are colonized by diverse but poorly studied communities of micro-organisms. One of the most noted ones is the unicellular green alga, Oophila amblystomatis, but the occurrence and role of other micro-organisms in the capsular chamber surrounding amphibian clutches have remained largely unstudied. Here, we undertook a multi-marker DNA metabarcoding study to characterize the community of algae and other micro-eukaryotes associated with agile frog (Rana dalmatina) clutches. Samplings were performed at three small ponds in Germany, from four substrates: water, sediment, tree leaves from the bottom of the pond, and R. dalmatina clutches. Sampling substrate strongly determined the community compositions of algae and other micro-eukaryotes. Therefore, as expected, the frog clutch-associated communities formed clearly distinct clusters. Clutch-associated communities in our study were structured by a plethora of not only green algae, but also diatoms and other ochrophytes. The most abundant operational taxonomic units (OTUs) in clutch samples were taxa from Chlamydomonas, Oophila, but also from Nitzschia and other ochrophytes. Sequences of Oophila "Clade B" were found exclusively in clutches. Based on additional phylogenetic analyses of 18S rDNA and of a matrix of 18 nuclear genes derived from transcriptomes, we confirmed in our samples the existence of two distinct clades of green algae assigned to Oophila in past studies. We hypothesize that "Clade B" algae correspond to the true Oophila, whereas "Clade A" algae are a series of Chlorococcum species that, along with other green algae, ochrophytes and protists, colonize amphibian clutches opportunistically and are often cultured from clutch samples due to their robust growth performance. The clutch-associated communities were subject to filtering by sampling location, suggesting that the taxa colonizing amphibian clutches can drastically differ depending on environmental conditions.

Surveying for Batrachochytrium dendrobatidis and B. salamandrivorans in wild and captive amphibian populations in Estonia and Latvia.

The pathogenic chytrid fungi Batrachochytrium dendrobatidis (Bd) and B. salamandrivorans (Bsal) cause infections that have become primary drivers of amphibian biodiversity loss. While globally widespread, the distribution margins of Bd and Bsal have not been determined, and the presence of these pathogens has probably gone unnoticed in many areas, especially in northern Eurasia. To better understand the presence and distribution of both pathogens in the northern temperate and boreal forest biomes, 243 individuals were sampled from 8 native amphibian species across Estonia. Additionally, 68 amphibians were sampled from captive collections in Estonia and Latvia. Pathogen infection was assessed using metabarcoding of the ITS2 marker. No positive matches for Bsal infection were found. Bd was detected in 13 specimens, 3 of which were sampled at the Riga Zoo (with a prevalence of 5.2%) and 10 in natural environments in Estonia (3.3%). The infected wild individuals belonged to 6 amphibian species and were detected throughout the mainland of Estonia, but not on islands. Prevalence of infection with Bd ranged between 3.1 and 12.5% among native species. In addition, we found molecular evidence for a potentially new sister species to Bd in nature. Although outbreaks of chytridiomycosis have never been observed in Estonia, it cannot be excluded that the dynamics of local amphibian populations are affected by Bd infections. Therefore, further work, including capture-mark-recapture studies and long-term monitoring, are required to clarify the impact of Bd on amphibians in Northern Europe.

Plant nutrient-acquisition strategies drive topsoil microbiome structure and function.

Plant nutrient-acquisition strategies drive soil processes and vegetation performance, but their effect on the soil microbiome remains poorly understood. This knowledge is important to predict the shifts in microbial diversity and functions due to increasing changes in vegetation traits under global change. Here we documented the topsoil microbiomes of 145 boreal and temperate terrestrial sites in the Baltic region that broadly differed in vegetation type and nutritional traits, such as mycorrhizal types and symbiotic nitrogen-fixation. We found that sites dominated by arbuscular mycorrhizal (AM) vegetation harbor relatively more AM fungi, bacteria, fungal saprotrophs, and pathogens in the topsoil compared with sites dominated by ectomycorrhizal (EM) plants. These differences in microbiome composition reflect the rapid nutrient cycling and negative plant-soil feedback in AM soils. Lower fungal diversity and bacteria : fungi ratios in EM-dominated habitats are driven by monodominance of woody vegetation as well as soil acidification by EM fungi, which are associated with greater diversity and relative abundance of carbohydrate-active enzymes. Our study suggests that shifts in vegetation related to global change and land use may strongly alter the topsoil microbiome structure and function.

Relative Performance of MinION (Oxford Nanopore Technologies) versus Sequel (Pacific Biosciences) Third-Generation Sequencing Instruments in Identification of Agricultural and Forest Fungal Pathogens.

Culture-based molecular identification methods have revolutionized detection of pathogens, yet these methods are slow and may yield inconclusive results from environmental materials. The second-generation sequencing tools have much-improved precision and sensitivity of detection, but these analyses are costly and may take several days to months. Of the third-generation sequencing techniques, the portable MinION device (Oxford Nanopore Technologies) has received much attention because of its small size and possibility of rapid analysis at reasonable cost. Here, we compare the relative performances of two third-generation sequencing instruments, MinION and Sequel (Pacific Biosciences), in identification and diagnostics of fungal and oomycete pathogens from conifer (Pinaceae) needles and potato (Solanum tuberosum) leaves and tubers. We demonstrate that the Sequel instrument is efficient for metabarcoding of complex samples, whereas MinION is not suited for this purpose due to a high error rate and multiple biases. However, we find that MinION can be utilized for rapid and accurate identification of dominant pathogenic organisms and other associated organisms from plant tissues following both amplicon-based and PCR-free metagenomics approaches. Using the metagenomics approach with shortened DNA extraction and incubation times, we performed the entire MinION workflow, from sample preparation through DNA extraction, sequencing, bioinformatics, and interpretation, in 2.5 h. We advocate the use of MinION for rapid diagnostics of pathogens and potentially other organisms, but care needs to be taken to control or account for multiple potential technical biases.IMPORTANCE Microbial pathogens cause enormous losses to agriculture and forestry, but current combined culturing- and molecular identification-based detection methods are too slow for rapid identification and application of countermeasures. Here, we develop new and rapid protocols for Oxford Nanopore MinION-based third-generation diagnostics of plant pathogens that greatly improve the speed of diagnostics. However, due to high error rate and technical biases in MinION, the Pacific BioSciences Sequel platform is more useful for in-depth amplicon-based biodiversity monitoring (metabarcoding) from complex environmental samples.

PacBio metabarcoding of Fungi and other eukaryotes: errors, biases and perspectives.

Second-generation, high-throughput sequencing methods have greatly improved our understanding of the ecology of soil microorganisms, yet the short barcodes (< 500 bp) provide limited taxonomic and phylogenetic information for species discrimination and taxonomic assignment. Here, we utilized the third-generation Pacific Biosciences (PacBio) RSII and Sequel instruments to evaluate the suitability of full-length internal transcribed spacer (ITS) barcodes and longer rRNA gene amplicons for metabarcoding Fungi, Oomycetes and other eukaryotes in soil samples. Metabarcoding revealed multiple errors and biases: Taq polymerase substitution errors and mis-incorporating indels in sequencing homopolymers constitute major errors; sequence length biases occur during PCR, library preparation, loading to the sequencing instrument and quality filtering; primer-template mismatches bias the taxonomic profile when using regular and highly degenerate primers. The RSII and Sequel platforms enable the sequencing of amplicons up to 3000 bp, but the sequence quality remains slightly inferior to Illumina sequencing especially in longer amplicons. The full ITS barcode and flanking rRNA small subunit gene greatly improve taxonomic identification at the species and phylum levels, respectively. We conclude that PacBio sequencing provides a viable alternative for metabarcoding of organisms that are of relatively low diversity, require > 500-bp barcode for reliable identification or when phylogenetic approaches are intended.

The endangered northern bettong, Bettongia tropica, performs a unique and potentially irreplaceable dispersal function for ectomycorrhizal truffle fungi.

Organisms that are highly connected in food webs often perform unique and vital functions within ecosystems. Understanding the unique ecological roles played by highly connected organisms and the consequences of their loss requires a comprehensive understanding of the functional redundancy amongst organisms. One important, yet poorly understood, food web is that between truffle-forming ectomycorrhizal fungi and their mammalian consumers and dispersers. Mammalian fungal specialists rely on fungi as a food source, and they consume and disperse a higher diversity and abundance of fungi than do mycophagous mammals with generalist diets. Therefore, we hypothesize that mammalian fungal specialists are functionally distinct because they disperse a set of fungal taxa not fully nested within the set consumed by the combined generalist mammalian community (i.e., functional redundancy of fungal dispersal is limited). Using high-throughput sequencing, we compared the fungal composition of 93 scats from the endangered fungal specialist northern bettong (Bettongia tropica) and 120 scats from nine co-occurring generalist mammal species across three sites and three seasons. Compared with other generalist mammals, B. tropica consumed a more diverse fungal diet with more unique taxa. This aligns with our hypothesis that B. tropica performs a unique dispersal function for ectomycorrhizal truffle fungi. Additionally, modelling of mammalian extinctions predicted rapid loss of food web connections which could result in loss of gene flow for truffle taxa. Our results suggest that this system is sensitive to the extinction of highly connected specialist species like B. tropica and their loss could have consequences for ectomycorrhizal truffle fungal diversity. This suggests that the conservation of fungal specialists is imperative to maintaining ectomycorrhizal fungal diversity and healthy plant-mycorrhizal relationships.

Tree species identity and diversity drive fungal richness and community composition along an elevational gradient in a Mediterranean ecosystem.

Ecological and taxonomic knowledge is important for conservation and utilization of biodiversity. Biodiversity and ecology of fungi in Mediterranean ecosystems is poorly understood. Here, we examined the diversity and spatial distribution of fungi along an elevational gradient in a Mediterranean ecosystem, using DNA metabarcoding. This study provides novel information about diversity of all ecological and taxonomic groups of fungi along an elevational gradient in a Mediterranean ecosystem. Our analyses revealed that among all biotic and abiotic variables tested, host species identity is the main driver of the fungal richness and fungal community composition. Fungal richness was strongly associated with tree richness and peaked in Quercus-dominated habitats and Cistus-dominated habitats. The highest taxonomic richness of ectomycorrhizal fungi was observed under Quercus ilex, whereas the highest taxonomic richness of saprotrophs was found under Pinus. Our results suggest that the effect of plant diversity on fungal richness and community composition may override that of abiotic variables across environmental gradients.

DNA-metabarcoding supports trophic flexibility and reveals new prey species for the Galapagos sea lion.

Tropical ecosystems are challenging for pinnipeds due to fluctuating food availability. According to previous research, the Galapagos sea lion (GSL, Zalophus wollebaeki) adopts trophic flexibility to face such conditions. However, this hypothesis comes from studies using traditional methods (hard-parts analysis of scat and isotopic analysis from tissue). We studied the diet of five rookeries in the southeastern Galapagos bioregion (which harbors the highest GSL density), via DNA-metabarcoding of scat samples. The DNA-metabarcoding approach may identify consumed prey with a higher taxonomic resolution than isotopic analysis, while not depending on hard-parts remaining through digestion. Our study included five different rookeries to look for evidence of trophic flexibility at the bioregional level. We detected 98 prey OTUs (124 scats), mostly assigned to bony-fish taxa; we identified novel prey items, including a shark, rays, and several deep-sea fish. Our data supported the trophic flexibility of GSL throughout the studied bioregion since different individuals from the same rookery consumed prey coming from different habitats and trophic levels. Significant diet differentiations were found among rookeries, particularly between Punta Pitt and Santa Fe. Punta Pitt rookery, with a more pronounced bathymetry and lower productivity, was distinguished by a high trophic level and consumption of a high proportion of deep-sea prey; meanwhile, Santa Fe, located in more productive, shallow waters over the shelf, consumed a high proportion of epipelagic planktivorous fish. Geographic location and heterogeneous bathymetry of El Malecon, Española, and Floreana rookeries would allow the animals therein to access both, epipelagic prey over the shelf, and deep-sea prey out of the shelf; this would lead to a higher prey richness and diet variability there. These findings provide evidence of GSL adopting a trophic flexibility to tune their diets to different ecological contexts. This strategy would be crucial for this endangered species to overcome the challenges faced in a habitat with fluctuating foraging conditions.

Dataset from RNAseq analysis of differential gene expression among developmental stages of two non-marine ostracodes.

We contribute transcriptomic data for two species of Ostracoda, an early-diverged group of small-sized pancrustaceans. Data include new reference transcriptomes for two asexual non-marine species (Dolerocypris sinensis and Heterocypris aff. salina), as well as single-specimen transcriptomic data that served to analyse gene expression across four developmental stages in D. sinensis. Data are evaluated by computing gene expression profiles of the different developmental stages which consistently placed eggs and small larvae (at the stage of instar A-8) similar to each other, and apart from adults which were distinct from all other developmental stages but closest to large larvae (instar A-4). We further evaluated the transcriptomic data with two newly sequenced low-coverage genomes of the target species. The new data thus document the feasibility of obtaining reliable transcriptomic data from single specimens - even eggs - of these small metazoans.

The influence of tree genus, phylogeny, and richness on the specificity, rarity, and diversity of ectomycorrhizal fungi.

Partner specificity is a well-documented phenomenon in biotic interactions, yet the factors that determine specificity in plant-fungal associations remain largely unknown. By utilizing composite soil samples, we identified the predictors that drive partner specificity in both plants and fungi, with a particular focus on ectomycorrhizal associations. Fungal guilds exhibited significant differences in overall partner preference and avoidance, richness, and specificity to specific tree genera. The highest level of specificity was observed in root endophytic and ectomycorrhizal associations, while the lowest was found in arbuscular mycorrhizal associations. The majority of ectomycorrhizal fungal species showed a preference for one of their partner trees, primarily at the plant genus level. Specialist ectomycorrhizal fungi were dominant in belowground communities in terms of species richness and relative abundance. Moreover, all tree genera (and occasionally species) demonstrated a preference for certain fungal groups. Partner specificity was not related to the rarity of fungi or plants or environmental conditions, except for soil pH. Depending on the partner tree genus, specific fungi became more prevalent and relatively more abundant with increasing stand age, tree dominance, and soil pH conditions optimal for the partner tree genus. The richness of partner tree species and increased evenness of ectomycorrhizal fungi in multi-host communities enhanced the species richness of ectomycorrhizal fungi. However, it was primarily the partner-generalist fungi that contributed to the high diversity of ectomycorrhizal fungi in mixed forests.

EUKARYOME: the rRNA gene reference database for identification of all eukaryotes.

Molecular identification of micro- and macroorganisms based on nuclear markers has revolutionized our understanding of their taxonomy, phylogeny and ecology. Today, research on the diversity of eukaryotes in global ecosystems heavily relies on nuclear ribosomal RNA (rRNA) markers. Here, we present the research community-curated reference database EUKARYOME for nuclear ribosomal 18S rRNA, internal transcribed spacer (ITS) and 28S rRNA markers for all eukaryotes, including metazoans (animals), protists, fungi and plants. It is particularly useful for the identification of arbuscular mycorrhizal fungi as it bridges the four commonly used molecular markers-ITS1, ITS2, 18S V4-V5 and 28S D1-D2 subregions. The key benefits of this database over other annotated reference sequence databases are that it is not restricted to certain taxonomic groups and it includes all rRNA markers. EUKARYOME also offers a number of reference long-read sequences that are derived from (meta)genomic and (meta)barcoding-a unique feature that can be used for taxonomic identification and chimera control of third-generation, long-read, high-throughput sequencing data. Taxonomic assignments of rRNA genes in the database are verified based on phylogenetic approaches. The reference datasets are available in multiple formats from the project homepage, http://www.eukaryome.org.

Structure of plant-associated microeukaryotes in roots and leaves of aquatic and terrestrial plants revealed by blocking peptide-nucleic acid (PNA) amplification.

Studies of plant-microbe interactions, including mutualistic, antagonistic, parasitic, or commensal microbes, have greatly benefited our understanding of ecosystem functioning. New molecular identification tools have increasingly revealed the association patterns between microorganisms and plants. Here, we integrated long-read PacBio single-molecule sequencing technology with a blocking protein-nucleic acid (PNA) approach to minimise plant amplicons in a survey of plant-eukaryotic microbe relationships in roots and leaves of different aquatic and terrestrial plants to determine patterns of organ, host, and habitat preferences. The PNA approach reduced the samples' relative amounts of plant reads and did not distort the fungal and other microeukaryotic composition. Our analyses revealed that the eukaryotic microbiomes associated with leaves and roots of aquatic plants exhibit a much larger proportion of non-fungal microorganisms than terrestrial plants, and leaf and root microbiomes are similar. Terrestrial plants had much stronger differentiation of leaf and root microbiomes and stronger partner specificity than aquatic plants.

How, not if, is the question mycologists should be asking about DNA-based typification.

Fungal metabarcoding of substrates such as soil, wood, and water is uncovering an unprecedented number of fungal species that do not seem to produce tangible morphological structures and that defy our best attempts at cultivation, thus falling outside the scope of the International Code of Nomenclature for algae, fungi, and plants. The present study uses the new, ninth release of the species hypotheses of the UNITE database to show that species discovery through environmental sequencing vastly outpaces traditional, Sanger sequencing-based efforts in a strongly increasing trend over the last five years. Our findings challenge the present stance of some in the mycological community - that the current situation is satisfactory and that no change is needed to "the code" - and suggest that we should be discussing not whether to allow DNA-based descriptions (typifications) of species and by extension higher ranks of fungi, but what the precise requirements for such DNA-based typifications should be. We submit a tentative list of such criteria for further discussion. The present authors hope for a revitalized and deepened discussion on DNA-based typification, because to us it seems harmful and counter-productive to intentionally deny the overwhelming majority of extant fungi a formal standing under the International Code of Nomenclature for algae, fungi, and plants.

Global diversity and distribution of nitrogen-fixing bacteria in the soil.

Our knowledge of microbial biogeography has advanced in recent years, yet we lack knowledge of the global diversity of some important functional groups. Here, we used environmental DNA from 327 globally collected soil samples to investigate the biodiversity patterns of nitrogen-fixing bacteria by focusing on the nifH gene but also amplifying the general prokaryotic 16S SSU region. Globally, N-fixing prokaryotic communities are driven mainly by climatic conditions, with most groups being positively correlated with stable hot or seasonally humid climates. Among soil parameters, pH, but also soil N content were most often shown to correlate with the diversity of N-fixer groups. However, specific groups of N-fixing prokaryotes show contrasting responses to the same variables, notably in Cyanobacteria that were negatively correlated with stable hot climates, and showed a U-shaped correlation with soil pH, contrary to other N-fixers. Also, the non-N-fixing prokaryotic community composition was differentially correlated with the diversity and abundance of N-fixer groups, showing the often-neglected impact of biotic interactions among bacteria.