Global diversity and geography of soil fungi
Assessing fungal diversity worldwide
Fungi are hyperdiverse but poorly known, despite their ecological and economic impacts. Tedersoo et al. collected nearly 15,000 topsoil samples from 365 sites worldwide and sequenced their genomes (see the Perspective by Wardle and Lindahl). Overall, they found a striking decline in fungal species richness with distance from the equator. For some specialist groups though, diversity depended more on the abundance of host plants than host diversity or geography. The findings reveal a huge gap between known and described species and the actual numbers of distinct fungi in the world's soils.
Science, this issue 10.1126/science.1256688; see also p. 1052
Structured Abstract
Introduction
The kingdom Fungi is one of the most diverse groups of organisms on Earth, and they are integral ecosystem agents that govern soil carbon cycling, plant nutrition, and pathology. Fungi are widely distributed in all terrestrial ecosystems, but the distribution of species, phyla, and functional groups has been poorly documented. On the basis of 365 global soil samples from natural ecosystems, we determined the main drivers and biogeographic patterns of fungal diversity and community composition.
Rationale
We identified soil-inhabiting fungi using 454 Life Sciences (Branford, CN) pyrosequencing and through comparison against taxonomically and functionally annotated sequence databases. Multiple regression models were used to disentangle the roles of climatic, spatial, edaphic, and floristic parameters on fungal diversity and community composition. Structural equation models were used to determine the direct and indirect effects of climate on fungal diversity, soil chemistry, and vegetation. We also examined whether fungal biogeographic patterns matched paradigms derived from plants and animals—namely, that species’ latitudinal ranges increase toward the poles (Rapoport’s rule) and diversity increases toward the equator. Last, we sought group-specific global biogeographic links among major biogeographic regions and biomes using a network approach and area-based clustering.
Results
Metabarcoding analysis of global soils revealed fungal richness estimates approaching the number of species recorded to date. Distance from equator and mean annual precipitation had the strongest effects on richness of fungi, including most fungal taxonomic and functional groups. Diversity of most fungal groups peaked in tropical ecosystems, but ectomycorrhizal fungi and several fungal classes were most diverse in temperate or boreal ecosystems, and many fungal groups exhibited distinct preferences for specific edaphic conditions (such as pH, calcium, or phosphorus). Consistent with Rapoport’s rule, the geographic range of fungal taxa increased toward the poles. Fungal endemicity was particularly strong in tropical regions, but multiple fungal taxa had cosmopolitan distribution.
Conclusions
Climatic factors, followed by edaphic and spatial patterning, are the best predictors of soil fungal richness and community composition at the global scale. Richness of all fungi and functional groups is causally unrelated to plant diversity, with the exception of ectomycorrhizal root symbionts, suggesting that plant-soil feedbacks do not influence the diversity of soil fungi at the global scale. The plant-to-fungi richness ratio declined exponentially toward the poles, indicating that current predictions—assuming globally constant ratios—overestimate fungal richness by 1.5- to 2.5-fold. Fungi follow similar biogeographic patterns as plants and animals, with the exception of several major taxonomic and functional groups that run counter to overall patterns. Strong biogeographic links among distant continents reflect relatively efficient long-distance dispersal compared with macro-organisms.
Abstract
Fungi play major roles in ecosystem processes, but the determinants of fungal diversity and biogeographic patterns remain poorly understood. Using DNA metabarcoding data from hundreds of globally distributed soil samples, we demonstrate that fungal richness is decoupled from plant diversity. The plant-to-fungus richness ratio declines exponentially toward the poles. Climatic factors, followed by edaphic and spatial variables, constitute the best predictors of fungal richness and community composition at the global scale. Fungi show similar latitudinal diversity gradients to other organisms, with several notable exceptions. These findings advance our understanding of global fungal diversity patterns and permit integration of fungi into a general macroecological framework.
Get full access to this article
View all available purchase options and get full access to this article.
Already a Subscriber?Sign In
Supplementary Material
Summary
Figs. S1 to S17
Tables S1 to S3
Data Files S1 and S2
Resources
References and Notes
1
Blackwell M., The fungi: 1, 2, 3 ... 5.1 million species? Am. J. Bot. 98, 426–438 (2011).
2
Fierer N., Strickland M. S., Liptzin D., Bradford M. A., Cleveland C. C., Global patterns in belowground communities. Ecol. Lett. 12, 1238 (2009).
3
Serna-Chavez H., Fierer N., van Bodegom P. M., Global drivers and patterns of microbial abundance in soil. Glob. Ecol. Biogeogr. 22, 1162–1172 (2013).
4
Xu X., Thornton P., Post W. M., A global analysis of soil microbial biomass carbon, nitrogen and phosphorus in terrestrial ecosystems. Glob. Ecol. Biogeogr. 22, 737–749 (2013).
5
Hillebrand H., On the generality of the latitudinal diversity gradient. Am. Nat. 163, 192–211 (2004).
6
Mittelbach G. G., Schemske D. W., Cornell H. V., Allen A. P., Brown J. M., Bush M. B., Harrison S. P., Hurlbert A. H., Knowlton N., Lessios H. A., McCain C. M., McCune A. R., McDade L. A., McPeek M. A., Near T. J., Price T. D., Ricklefs R. E., Roy K., Sax D. F., Schluter D., Sobel J. M., Turelli M., Evolution and the latitudinal diversity gradient: Speciation, extinction and biogeography. Ecol. Lett. 10, 315–331 (2007).
7
Finlay B. J., Global dispersal of free-living microbial eukaryote species. Science 296, 1061–1063 (2002).
8
Nemergut D. R., Schmidt S. K., Fukami T., O’Neill S. P., Bilinski T. M., Stanish L. F., Knelman J. E., Darcy J. L., Lynch R. C., Wickey P., Ferrenberg S., Patterns and processes of microbial community assembly. Microbiol. Mol. Biol. Rev. 77, 342–356 (2013).
9
Peay K. G., Bidartondo M. I., Arnold A. E., Not every fungus is everywhere: Scaling to the biogeography of fungal-plant interactions across roots, shoots and ecosystems. New Phytol. 185, 878–882 (2010).
10
Talbot J. M., Bruns T. D., Taylor J. W., Smith D. P., Branco S., Glassman S. I., Erlandson S., Vilgalys R., Liao H. L., Smith M. E., Peay K. G., Endemism and functional convergence across the North American soil mycobiome. Proc. Natl. Acad. Sci. U.S.A. 111, 6341–6346 (2014).
11
Lindahl B. D., Nilsson R. H., Tedersoo L., Abarenkov K., Carlsen T., Kjøller R., Kõljalg U., Pennanen T., Rosendahl S., Stenlid J., Kauserud H., Fungal community analysis by high-throughput sequencing of amplified markers—A user’s guide. New Phytol. 199, 288–299 (2013).
12
Kõljalg U., Nilsson R. H., Abarenkov K., Tedersoo L., Taylor A. F., Bahram M., Bates S. T., Bruns T. D., Bengtsson-Palme J., Callaghan T. M., Douglas B., Drenkhan T., Eberhardt U., Dueñas M., Grebenc T., Griffith G. W., Hartmann M., Kirk P. M., Kohout P., Larsson E., Lindahl B. D., Lücking R., Martín M. P., Matheny P. B., Nguyen N. H., Niskanen T., Oja J., Peay K. G., Peintner U., Peterson M., Põldmaa K., Saag L., Saar I., Schüßler A., Scott J. A., Senés C., Smith M. E., Suija A., Taylor D. L., Telleria M. T., Weiss M., Larsson K. H., Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277 (2013).
13
Abarenkov K., et al., Evol. Bioinform. 6, 189 (2010).
14
Adl S. M., Simpson A. G. B., Lane C. E., Lukeš J., Bass D., Bowser S. S., Brown M. W., Burki F., Dunthorn M., Hampl V., Heiss A., Hoppenrath M., Lara E., le Gall L., Lynn D. H., McManus H., Mitchell E. A. D., Mozley-Stanridge S. E., Parfrey L. W., Pawlowski J., Rueckert S., Shadwick L., Schoch C. L., Smirnov A., Spiegel F. W., The revised classification of eukaryotes. J. Eukaryot. Microbiol. 59, 429 (2012).
15
Tedersoo L., Smith M. E., Lineages of ectomycorrhizal fungi revisited: Foraging strategies and novel lineages revealed by sequences from belowground. Fungal Biol. Rev. 27, 83–99 (2013).
16
Tedersoo L., Naadel T., Bahram M., Pritsch K., Buegger F., Leal M., Kõljalg U., Põldmaa K., Enzymatic activities and stable isotope patterns of ectomycorrhizal fungi in relation to phylogeny and exploration types in an afrotropical rain forest. New Phytol. 195, 832–843 (2012).
17
Hiiesalu I., Pärtel M., Davison J., Gerhold P., Metsis M., Moora M., Öpik M., Vasar M., Zobel M., Wilson S. D., Species richness of arbuscular mycorrhizal fungi: Associations with grassland plant richness and biomass. New Phytol. 203, 233–244 (2014).
18
Kreft H., Jetz W., Global patterns and determinants of vascular plant diversity. Proc. Natl. Acad. Sci. U.S.A. 104, 5925–5930 (2007).
19
Taylor D. L., Hollingsworth T. N., McFarland J. W., Lennon N. J., Nusbaum C., Ruess R. W., A first comprehensive census of fungi in soil reveals both hyperdiversity and fine-scale niche partitioning. Ecol. Monogr. 84, 3–20 (2014).
20
O’Brien H. E., Parrent J. L., Jackson J. A., Moncalvo J. M., Vilgalys R., Fungal community analysis by large-scale sequencing of environmental samples. Appl. Environ. Microbiol. 71, 5544–5550 (2005).
21
R Core Team, R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, Vienna 2014).
22
Rangel T. F., Diniz-Filho J. A. F., Bini L. M., SAM: A comprehensive application for Spatial Analysis in Macroecology. Ecography 33, 46–50 (2010).
23
Schoch C. L., Seifert K. A., Huhndorf S., Robert V., Spouge J. L., Levesque C. A., Chen W., Bolchacova E., Voigt K., Crous P. W., Miller A. N., Wingfield M. J., Aime M. C., An K.-D., Bai F.-Y., Barreto R. W., Begerow D., Bergeron M.-J., Blackwell M., Boekhout T., Bogale M., Boonyuen N., Burgaz A. R., Buyck B., Cai L., Cai Q., Cardinali G., Chaverri P., Coppins B. J., Crespo A., Cubas P., Cummings C., Damm U., de Beer Z. W., de Hoog G. S., Del-Prado R., Dentinger B., Dieguez-Uribeondo J., Divakar P. K., Douglas B., Duenas M., Duong T. A., Eberhardt U., Edwards J. E., Elshahed M. S., Fliegerova K., Furtado M., Garcia M. A., Ge Z.-W., Griffith G. W., Griffiths K., Groenewald J. Z., Groenewald M., Grube M., Gryzenhout M., Guo L.-D., Hagen F., Hambleton S., Hamelin R. C., Hansen K., Harrold P., Heller G., Herrera C., Hirayama K., Hirooka Y., Ho H.-M., Hoffmann K., Hofstetter V., Hognabba F., Hollingsworth P. M., Hong S.-B., Hosaka K., Houbraken J., Hughes K., Huhtinen S., Hyde K. D., James T., Johnson E. M., Johnson J. E., Johnston P. R., Jones E. B. G., Kelly L. J., Kirk P. M., Knapp D. G., Koljalg U., Kovacs G. M., Kurtzman C. P., Landvik S., Leavitt S. D., Liggenstoffer A. S., Liimatainen K., Lombard L., Luangsa-ard J. J., Lumbsch H. T., Maganti H., Maharachchikumbura S. S. N., Martin M. P., May T. W., McTaggart A. R., Methven A. S., Meyer W., Moncalvo J.-M., Mongkolsamrit S., Nagy L. G., Nilsson R. H., Niskanen T., Nyilasi I., Okada G., Okane I., Olariaga I., Otte J., Papp T., Park D., Petkovits T., Pino-Bodas R., Quaedvlieg W., Raja H. A., Redecker D., Rintoul T. L., Ruibal C., Sarmiento-Ramirez J. M., Schmitt I., Schussler A., Shearer C., Sotome K., Stefani F. O. P., Stenroos S., Stielow B., Stockinger H., Suetrong S., Suh S.-O., Sung G.-H., Suzuki M., Tanaka K., Tedersoo L., Telleria M. T., Tretter E., Untereiner W. A., Urbina H., Vagvolgyi C., Vialle A., Vu T. D., Walther G., Wang Q.-M., Wang Y., Weir B. S., Weiss M., White M. M., Xu J., Yahr R., Yang Z. L., Yurkov A., Zamora J.-C., Zhang N., Zhuang W.-Y., Schindel D.Fungal Barcoding ConsortiumFungal Barcoding Consortium Author List, Nuclear ribosomal internal transcribed spacer (ITS) region as a universal DNA barcode marker for Fungi. Proc. Natl. Acad. Sci. U.S.A. 109, 6241–6246 (2012).
24
Stevens G. C., The latitudinal gradient in geographical range: How so many species coexist in the tropics. Am. Nat. 133, 240 (1989).
25
Dickie I. A., Insidious effects of sequencing errors on perceived diversity in molecular surveys. New Phytol. 188, 916–918 (2010).
26
Jones M. D. M., Forn I., Gadelha C., Egan M. J., Bass D., Massana R., Richards T. A., Discovery of novel intermediate forms redefines the fungal tree of life. Nature 474, 200–203 (2011).
27
Lauber C. L., Hamady M., Knight R., Fierer N., Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120 (2009).
28
Robeson M. S., King A. J., Freeman K. R., Birky C. W., Martin A. P., Schmidt S. K., Soil rotifer communities are extremely diverse globally but spatially autocorrelated locally. Proc. Natl. Acad. Sci. U.S.A. 108, 4406–4410 (2011).
29
Rosling A., Cox F., Cruz-Martinez K., Ihrmark K., Grelet G. A., Lindahl B. D., Menkis A., James T. Y., Archaeorhizomycetes: Unearthing an ancient class of ubiquitous soil fungi. Science 333, 876–879 (2011).
30
Põlme S., Bahram M., Yamanaka T., Nara K., Dai Y. C., Grebenc T., Kraigher H., Toivonen M., Wang P. H., Matsuda Y., Naadel T., Kennedy P. G., Kõljalg U., Tedersoo L., Biogeography of ectomycorrhizal fungi associated with alders (Alnus spp.) in relation to biotic and abiotic variables at the global scale. New Phytol. 198, 1239–1249 (2013).
31
Reich P. B., Oleksyn J., Modrzynski J., Mrozinski P., Hobbie S. E., Eissenstat D. M., Chorover J., Chadwick O. A., Hale C. M., Tjoelker M. G., Linking litter calcium, earthworms and soil properties: A common garden test with 14 tree species. Ecol. Lett. 8, 811–818 (2005).
32
Arnold A. E., Understanding the diversity of foliar endophytic fungi: Progress, challenges, and frontiers. Fungal Biol. Rev. 21, 51–66 (2007).
33
Berbee M. L., Taylor J. W., Dating the molecular clock in fungi—How close are we? Fungal Biol. Rev. 24, 1–16 (2010).
34
Treseder K. K., Maltz M. R., Hawkins B. A., Fierer N., Stajich J. E., McGuire K. L., Evolutionary histories of soil fungi are reflected in their large-scale biogeography. Ecol. Lett. 17, 1086–1093 (2014).
35
Gilbert G. S., Webb C. O., Phylogenetic signal in plant pathogen-host range. Proc. Natl. Acad. Sci. U.S.A. 104, 4979–4983 (2007).
36
Liu X., Liang M., Etienne R. S., Wang Y., Staehelin C., Yu S., Experimental evidence for a phylogenetic Janzen-Connell effect in a subtropical forest. Ecol. Lett. 15, 111–118 (2012).
37
Bagchi R., Gallery R. E., Gripenberg S., Gurr S. J., Narayan L., Addis C. E., Freckleton R. P., Lewis O. T., Pathogens and insect herbivores drive rainforest plant diversity and composition. Nature 506, 85–88 (2014).
38
Bahram M., Harend H., Tedersoo L., Network perspectives of ectomycorrhizal associations. Fungal Ecol. 7, 70–77 (2013).
39
Tedersoo L., Bahram M., Toots M., Diédhiou A. G., Henkel T. W., Kjøller R., Morris M. H., Nara K., Nouhra E., Peay K. G., Põlme S., Ryberg M., Smith M. E., Kõljalg U., Towards global patterns in the diversity and community structure of ectomycorrhizal fungi. Mol. Ecol. 21, 4160–4170 (2012).
40
Bahram M., Kõljalg U., Courty P.-E., Diédhiou A. G., Kjøller R., Põlme S., Ryberg M., Veldre V., Tedersoo L., The distance decay of similarity in communities of ectomycorrhizal fungi in different ecosystems and scales. J. Ecol. 101, 1335–1344 (2013).
41
Qian H., Chen S., Mao L., Ouyang Z., Drivers of β-diversity along latitudinal gradients revisited. Glob. Ecol. Biogeogr. 22, 659–670 (2013).
42
Qian H., Beta diversity in relation to dispersal ability for vascular plants in North America. Glob. Ecol. Biogeogr. 18, 327–332 (2009).
43
Geml J., Timling I., Robinson C. H., Lennon N., Nusbaum H. C., Brochmann C., Noordeloos M. E., Taylor D. L., An arctic community of symbiotic fungi assembled by long-distance dispersers: Phylogenetic diversity of ectomycorrhizal basidiomycetes in Svalbard based on soil and sporocarp DNA. J. Biogeogr. 39, 74–88 (2012).
44
Timling I., Walker D. A., Nusbaum C., Lennon N. J., Taylor D. L., Rich and cold: Diversity, distribution and drivers of fungal communities in patterned-ground ecosystems of the North American Arctic. Mol. Ecol. 23, 3258–3272 (2014).
45
Rousk J., Bååth E., Brookes P. C., Lauber C. L., Lozupone C., Caporaso J. G., Knight R., Fierer N., Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 4, 1340–1351 (2010).
46
Crisp M. D., Arroyo M. T., Cook L. G., Gandolfo M. A., Jordan G. J., McGlone M. S., Weston P. H., Westoby M., Wilf P., Linder H. P., Phylogenetic biome conservatism on a global scale. Nature 458, 754–756 (2009).
47
Bates S. T., Clemente J. C., Flores G. E., Walters W. A., Parfrey L. W., Knight R., Fierer N., Global biogeography of highly diverse protistan communities in soil. ISME J. 7, 652–659 (2013).
48
Sul W. J., Oliver T. A., Ducklow H. W., Amaral-Zettler L. A., Sogin M. L., Marine bacteria exhibit a bipolar distribution. Proc. Natl. Acad. Sci. U.S.A. 110, 2342–2347 (2013).
49
Sato H., Tsujino R., Kurita K., Yokoyama K., Agata K., Modelling the global distribution of fungal species: New insights into microbial cosmopolitanism. Mol. Ecol. 21, 5599–5612 (2012).
50
Sanmartín I., Ronquist F., Southern Hemisphere biogeography inferred by event-based models: Plant versus animal patterns. Syst. Biol. 53, 216–243 (2004).
51
Altizer S., Ostfeld R. S., Johnson P. T., Kutz S., Harvell C. D., Climate change and infectious diseases: From evidence to a predictive framework. Science 341, 514–519 (2013).
52
Van der Putten W. H., Macel M., Visser M. E., Predicting species distribution and abundance responses to climate change: Why it is essential to include biotic interactions across trophic levels. Philos. Trans. R. Soc. B 365, 2025–2034 (2010).
53
Desprez-Loustau M.-L., Robin C., Buée M., Courtecuisse R., Garbaye J., Suffert F., Sache I., Rizzo D. M., The fungal dimension of biological invasions. Trends Ecol. Evol. 22, 472–480 (2007).
Information & Authors
Information
Published In

Science
Volume 346 | Issue 6213
28 November 2014
28 November 2014
Copyright
Copyright © 2014, American Association for the Advancement of Science.
Submission history
Received: 29 May 2014
Accepted: 16 October 2014
Published in print: 28 November 2014
Acknowledgments
The sequence data and metadata are deposited in the Short Read Archive (accession SRP043706) and UNITE databases. Data used for analyses are available as supplementary materials, data files S1 and S2. We thank A. Corrales, H. Mann, D. Sveshnikov, F. O. P. Stefani, A. Voitk, and Y. Wu for supplying single soil samples; R. Puusepp, M. Haugas, and M. Nõukas for sample preparation; H. Kreft for providing interpolated plant diversity data; S. Jüris for designing the printed figure; M. I. Bidartondo, K. G. Peay, and three anonymous reviewers for constructive comments on the manuscript; and relevant institutions of multiple countries for issuing permissions for sampling and delivery. The bulk of this project was funded from Estonian Science Foundation grants 9286, 171PUT, and IUT20-30; EMP265; Frontiers in Biodiversity Research; European Research Council; and in part by numerous funding sources that facilitated co-author efforts in collecting and preprocessing samples.
Authors
Metrics & Citations
Metrics
Article Usage
Altmetrics
Citations
Export citation
Select the format you want to export the citation of this publication.
Cited by
- Drivers of soil biodiversity vary with organism type along an extensive aridity gradient, Applied Soil Ecology, 170, (104271), (2022).https://doi.org/10.1016/j.apsoil.2021.104271
- Plant and soil elemental C:N:P ratios are linked to soil microbial diversity during grassland restoration on the Loess Plateau, China, Science of The Total Environment, 806, (150557), (2022).https://doi.org/10.1016/j.scitotenv.2021.150557
- N-induced root exudates mediate the rhizosphere fungal assembly and affect species coexistence, Science of The Total Environment, 804, (150148), (2022).https://doi.org/10.1016/j.scitotenv.2021.150148
- Microbial communities are associated with indicators of soil surface condition across a continental gradient, Geoderma, 405, (115439), (2022).https://doi.org/10.1016/j.geoderma.2021.115439
- Enzymic moderations of bacterial and fungal communities on short- and long-term warming impacts on soil organic carbon, Science of The Total Environment, 804, (150197), (2022).https://doi.org/10.1016/j.scitotenv.2021.150197
- Composition and potential functional roles of soil fungal communities on arid farms in Arequipa (Southern Peru) characterized using SMRT sequencing, Applied Soil Ecology, 169, (104228), (2022).https://doi.org/10.1016/j.apsoil.2021.104228
- Perception of fungi by farmers in the Cerrado, Brazilian Journal of Biology, 82, (2022).https://doi.org/10.1590/1519-6984.236219
- Soil ecology and restoration science, Soils and Landscape Restoration, (39-62), (2021).https://doi.org/10.1016/B978-0-12-813193-0.00002-3
- The factors and scales shaping fungal assemblages in fallen spruce trunks: A DNA metabarcoding study, Forest Ecology and Management, 495, (119381), (2021).https://doi.org/10.1016/j.foreco.2021.119381
- Global Diversity and Biogeography of the Zostera marina Mycobiome , Applied and Environmental Microbiology, 87, 12, (2021).https://doi.org/10.1128/AEM.02795-20
- See more
Loading...
View Options
Get Access
Log in to view the full text
AAAS login provides access to Science for AAAS Members, and access to other journals in the Science family to users who have purchased individual subscriptions.
- Become a AAAS Member
- Activate your AAAS ID
- Purchase Access to Other Journals in the Science Family
- Account Help
Log in via OpenAthens.
Log in via Shibboleth.
More options
Register for free to read this article
As a service to the community, this article is available for free. Login or register for free to read this article.
Buy a single issue of Science for just $15 USD.
View options
PDF format
Download this article as a PDF file
Download PDF





