Philippa Kaur

1.Loreau M, Naeem S, Inchausti P, Bengtsson J, Grime JP, Hector A, et al. Ecology: biodiversity and ecosystem functioning: current knowledge and future challenges. Science (80-). 2001;294:804–8.CAS 
Article 

Google Scholar 
2.Cardinale BJ, Duffy JE, Gonzalez A, Hooper DU, Perrings C, Venail P, et al. Biodiversity loss and its impact on humanity. Nature. 2012;486:59–67.CAS 
PubMed 
Article 

Google Scholar 
3.Jochum M, Fischer M, Isbell F, Roscher C, van der Plas F, Boch S, et al. The results of biodiversity-ecosystem functioning experiments are realistic. Nat Ecol Evol. 2020;4:1485–94.PubMed 
Article 

Google Scholar 
4.Balvanera P, Pfisterer AB, Buchmann N, He JS, Nakashizuka T, Raffaelli D, et al. Quantifying the evidence for biodiversity effects on ecosystem functioning and services. Ecol Lett. 2006;9:1146–56.PubMed 
Article 

Google Scholar 
5.Cardinale BJ, Matulich KL, Hooper DU, Byrnes JE, Duffy E, Gamfeldt L, et al. The functional role of producer diversity in ecosystems. Am J Bot. 2011;98:572–92.PubMed 
Article 

Google Scholar 
6.Weißbecker C, Heintz-Buschart A, Bruelheide H, Buscot F, Wubet T. Linking soil fungal generality to tree richness in young subtropical Chinese forests. Microorganisms. 2019; https://doi.org/10.3390/microorganisms7110547.7.Prada-Salcedo LD, Wambsganss J, Bauhus J, Buscot F, Goldmann K. Low root functional dispersion enhances functionality of plant growth by influencing bacterial activities in European forest soils. Env Microbiol. 2020; https://doi.org/10.1111/1462-2920.15244.8.Prada-Salcedo LD, Goldmann K, Heintz-Buschart A, Reitz T, Wambsganss J, Bauhus J, et al. Fungal guilds and soil functionality respond to tree community traits rather than to tree diversity in European forests. Mol Ecol. 2021;30:572–91.CAS 
PubMed 
Article 

Google Scholar 
9.Baldrian P. The known and the unknown in soil microbial ecology. FEMS Microbiol Ecol. 2019; https://doi.org/10.1093/femsec/fiz005.10.van Dijk EL, Auger H, Jaszczyszyn Y, Thermes C. Ten years of next-generation sequencing technology. Trends Genet. 2014;30:418–26.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
11.Smith SE, Read DJ. Mycorrhizal symbiosis. 2010. Academic press.12.Chen W, Koide RT, Eissenstat DM, Field K. Nutrient foraging by mycorrhizas: from species functional traits to ecosystem processes. Funct Ecol. 2018;32:858–69.Article 

Google Scholar 
13.Bahadur A, Batool A, Nasir F, Jiang S, Mingsen Q, Zhang Q, et al. Mechanistic insights into arbuscular mycorrhizal fungi-mediated drought stress tolerance in plants. Int J Mol Sci. 2019; https://doi.org/10.3390/ijms20174199.14.Pena R, Polle A. Attributing functions to ectomycorrhizal fungal identities in assemblages for nitrogen acquisition under stress. ISME J. 2014;8:321–30.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
15.He X-H, Critchley C, Bledsoe C. Nitrogen transfer within and between plants through common mycorrhizal networks (CMNs). CRC Crit Rev Plant Sci. 2003;22:531–67.Article 

Google Scholar 
16.Aerts R. The role of various types of mycorrhizal fungi in nutrient cycling and plant competition. In: Mycorrhizal Ecol. Springer; 2003. p. 117–33.17.Phillips RP, Brzostek E, Midgley MG. The mycorrhizal-associated nutrient economy: a new framework for predicting carbon-nutrient couplings in temperate forests. New Phytol. 2013;199:41–51.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
18.Hodge A, Campbell CD, Fitter AH. An arbuscular mycorrhizal fungus accelerates decomposition and acquires nitrogen directly from organic material. Nature. 2001;413:297–9.CAS 
PubMed 
Article 

Google Scholar 
19.Koide RT, Kabir Z. Extraradical hyphae of the mycorrhizal fungus Glomus intraradices can hydrolyse organic phosphate. New Phytol. 2000;148:511–7.CAS 
PubMed 
Article 

Google Scholar 
20.Martin F, Kohler A, Murat C, Veneault-Fourrey C, Hibbett DS. Unearthing the roots of ectomycorrhizal symbioses. Nat Rev Microbiol. 2016;14:760–73.CAS 
PubMed 
Article 

Google Scholar 
21.Teste FP, Jones MD, Dickie IA. Dual-mycorrhizal plants: their ecology and relevance. New Phytol. 2020;225:1835–51.PubMed 
Article 

Google Scholar 
22.Heklau H, Schindler N, Buscot F, Eisenhauer N, Ferlian O, Prada Salcedo LD, et al. Mixing tree species associated with arbuscular or ectotrophic mycorrhizae reveals dual mycorrhization and interactive effects on the fungal partners. Ecol Evol. 2021.23.Regvar M, Likar M, Piltaver A, Kugonič N, Smith JE. Fungal community structure under goat willows (Salix caprea L.) growing at metal polluted site: the potential of screening in a model phytostabilisation study. Plant Soil. 2010;330:345–56.CAS 
Article 

Google Scholar 
24.Waldrop MP, Zak DR, Blackwood CB, Curtis CD, Tilman D. Resource availability controls fungal diversity across a plant diversity gradient. Ecol Lett. 2006;9:1127–35.PubMed 
Article 

Google Scholar 
25.Hooper DU, Bignell DE, Brown VK, Brussard L, Mark Dangerfield J, Wall DH, et al. Interactions between aboveground and belowground biodiversity in terrestrial ecosystems: patterns, mechanisms, and feedbacks: we assess the evidence for correlation between aboveground and belowground diversity and conclude that a variety of mechanisms co. Bioscience. 2000;50:1049–61.Article 

Google Scholar 
26.Montesinos‐Navarro A, Segarra‐Moragues JG, Valiente‐Banuet A, Verdú M. The network structure of plant–arbuscular mycorrhizal fungi. New Phytol. 2012;194:536–47.PubMed 
Article 

Google Scholar 
27.Bahram M, Harend H, Tedersoo L. Network perspectives of ectomycorrhizal associations. Fungal Ecol. 2014;7:70–7.Article 

Google Scholar 
28.Querejeta J, Egerton-Warburton LM, Allen MF. Topographic position modulates the mycorrhizal response of oak trees to interannual rainfall variability. Ecology. 2009;90:649–62.PubMed 
Article 

Google Scholar 
29.Bergmann J, Weigelt A, van der Plas F, Laughlin DC, Kuyper TW, Guerrero-Ramirez N, et al. The fungal collaboration gradient dominates the root economics space in plants. Sci Adv. 2020;6:eaba3756.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
30.Weißbecker C, Wubet T, Lentendu G, Kühn P, Scholten T, Bruelheide H, et al. Experimental evidence of functional group-dependent effects of tree diversity on soil fungi in subtropical forests. Front Microbiol. 2018;9:2312.PubMed 
PubMed Central 
Article 

Google Scholar 
31.Ferlian O, Cesarz S, Craven D, Hines J, Barry KE, Bruelheide H, et al. Mycorrhiza in tree diversity–ecosystem function relationships: conceptual framework and experimental implementation. Ecosphere. 2018;9:e02226.PubMed 
PubMed Central 
Article 

Google Scholar 
32.Tedersoo L, May TW, Smith ME. Ectomycorrhizal lifestyle in fungi: global diversity, distribution, and evolution of phylogenetic lineages. Mycorrhiza. 2010;20:217–63.PubMed 
Article 

Google Scholar 
33.Kolaříková Z, Kohout P, Krüger C, Janoušková M, Mrnka L, Rydlová J. Root-associated fungal communities along a primary succession on a mine spoil: distinct ecological guilds assemble differently. Soil Biol Biochem. 2017;113:143–52.Article 
CAS 

Google Scholar 
34.Dang P, Vu NH, Shen Z, Liu J, Zhao F, Zhu H, et al. Changes in soil fungal communities and vegetation following afforestation with Pinus tabulaeformis on the Loess Plateau. Ecosphere. 2018;9:e02401.Article 

Google Scholar 
35.Kalucka IL, Jagodzinski AM. Successional traits of ectomycorrhizal fungi in forest reclamation after surface mining and agricultural disturbances: A review. Dendrobiology. 2016;76:91–104.Article 

Google Scholar 
36.Jones MD, Durall DM, Cairney JWG. Ectomycorrhizal fungal communities in young forest stands regenerating after clearcut logging. New Phytol. 2003;157:399–422.PubMed 
Article 

Google Scholar 
37.Rog I, Rosenstock NP, Korner C, Klein T. Share the wealth: trees with greater ectomycorrhizal species overlap share more carbon. Mol Ecol. 2020;29:2321–33.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
38.Chagnon PL, Bradley RL, Maherali H, Klironomos JN. A trait-based framework to understand life history of mycorrhizal fungi. Trends Plant Sci. 2013;18:484–91.CAS 
PubMed 
Article 

Google Scholar 
39.Ohsowski BM, Zaitsoff PD, Öpik M, Hart MM. Where the wild things are: looking for uncultured Glomeromycota. New Phytol. 2014;204:171–9.PubMed 
Article 
PubMed Central 

Google Scholar 
40.Öpik M, Metsis M, Daniell TJ, Zobel M, Moora M. Large-scale parallel 454 sequencing reveals host ecological group specificity of arbuscular mycorrhizal fungi in a boreonemoral forest. New Phytol. 2009;184:424–37.PubMed 
Article 
CAS 

Google Scholar 
41.Buscot F. Implication of evolution and diversity in arbuscular and ectomycorrhizal symbioses. J Plant Physiol. 2015;172:55–61.CAS 
PubMed 
Article 

Google Scholar 
42.Hiiesalu I, Pärtel M, Davison J, Gerhold P, Metsis M, Moora M, et al. Species richness of arbuscular mycorrhizal fungi: associations with grassland plant richness and biomass. New Phytol. 2014;203:233–44.CAS 
PubMed 
Article 

Google Scholar 
43.Nguyen NH, Williams LJ, Vincent JB, Stefanski A, Cavender-Bares J, Messier C, et al. Ectomycorrhizal fungal diversity and saprotrophic fungal diversity are linked to different tree community attributes in a field‐based tree experiment. Mol Ecol. 2016;25:4032–46.PubMed 
Article 

Google Scholar 
44.Burrows RL, Pfleger FL. Arbuscular mycorrhizal fungi respond to increasing plant diversity. Can J Bot. 2002;80:120–30.Article 

Google Scholar 
45.Eisenhauer N, Lanoue A, Strecker T, Scheu S, Steinauer K, Thakur MP, et al. Root biomass and exudates link plant diversity with soil bacterial and fungal biomass. Sci Rep. 2017;7:44641.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
46.Lange M, Eisenhauer N, Sierra CA, Bessler H, Engels C, Griffiths RI, et al. Plant diversity increases soil microbial activity and soil carbon storage. Nat Commun. 2015;6:1–8.
Google Scholar 
47.Klironomos JN, McCune J, Hart M, Neville J. The influence of arbuscular mycorrhizae on the relationship between plant diversity and productivity. Ecol Lett. 2000;3:137–41.Article 

Google Scholar 
48.Saks Ü, Davison J, Öpik M, Vasar M, Moora M, Zobel M. Root-colonizing and soil-borne communities of arbuscular mycorrhizal fungi in a temperate forest understorey. Botany. 2013;92:277–85.Article 

Google Scholar 
49.Molina R, Horton TR. Mycorrhiza specificity: its role in the development and function of common mycelial networks BT – Mycorrhizal Networks. In: Horton TR, editor. Springer Netherlands, Dordrecht; 2015. p. 1–39.50.van der Linde S, Suz LM, Orme C, Cox F, Andreae H, Asi E, et al. Environment and host as large-scale controls of ectomycorrhizal fungi. Nature. 2018;558:243–8.PubMed 
Article 
CAS 

Google Scholar 
51.Rasmussen AL, Busby RR, Hoeksema JD. Host preference of ectomycorrhizal fungi in mixed pine–oak woodlands. Can J For Res. 2017;48:153–9.Article 
CAS 

Google Scholar 
52.Soudzilovskaia NA, Vaessen S, van’t Zelfde M, Raes N. Global patterns of mycorrhizal distribution and their environmental drivers. In: Biogeogr. mycorrhizal symbiosis. Springer; 2017. p. 223–35.53.Simard SW, Jones MD, Durall DM. Carbon and nutrient fluxes within and between mycorrhizal plants BT – Mycorrhizal Ecology. In: van der Heijden MGA, Sanders IR, editors. Springer Berlin Heidelberg, Berlin, Heidelberg; 2003. p. 33–74.54.Allen EB, Allen MF, Helm DJ, Trappe JM, Molina R, Rincon E. Patterns and regulation of mycorrhizal plant and fungal diversity. Plant Soil. 1995;170:47–62.CAS 
Article 

Google Scholar 
55.Dickie IA, Koide RT, Fayish AC. Vesicular–arbuscular mycorrhizal infection of Quercus rubra seedlings. New Phytol. 2001;151:257–64.PubMed 
Article 
PubMed Central 

Google Scholar 
56.Davison J, Öpik M, Daniell TJ, Moora M, Zobel M. Arbuscular mycorrhizal fungal communities in plant roots are not random assemblages. FEMS Microbiol Ecol. 2011;78:103–15.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
57.Singavarapu B, Beugnon R, Bruelheide H, Cesarz S, Du J, Eisenhauer N, et al. Tree mycorrhizal type and tree diversity shape the forest soil microbiota. Environ Microbiol. 2021.58.Altermann M, Rinklebe J, Merbach I, Körschens M, Langer U, Hofmann B. Chernozem—soil of the year 2005. J Plant Nutr Soil Sci. 2005;168:725–40.CAS 
Article 

Google Scholar 
59.Wang B, Qiu Y-L. Phylogenetic distribution and evolution of mycorrhizas in land plants. Mycorrhiza. 2006;16:299–363.CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
60.Vierheilig H, Schweiger P, Brundrett M. An overview of methods for the detection and observation of arbuscular mycorrhizal fungi in roots. Physiol Plant. 2005;125:393–404.CAS 

Google Scholar 
61.Giovannetti M, Mosse B. An evaluation of techniques for measuring vesicular arbuscular mycorrhizal infection in roots. New Phytol. 1980; 489–500.62.White TJ, Bruns T, Lee S, Taylor J. Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. In: PCR Protoc. a Guid. to methods Appl. San Diego; 1990. p. 315–22.63.Wahdan SFM, Reitz T, Heintz-Buschart A, Schädler M, Roscher C, Breitkreuz C, et al. Organic agricultural practice enhances arbuscular mycorrhizal symbiosis in correspondence to soil warming and altered precipitation patterns. Environ Microbiol. 2021. https://doi.org/10.1111/1462-2920.15492.64.Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJ, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
65.Weißbecker C, Schnabel B, Heintz-Buschart A. Dadasnake, a Snakemake implementation of DADA2 to process amplicon sequencing data for microbial ecology. Gigascience. 2020. https://doi.org/10.1093/gigascience/giaa135.66.Opik M, Vanatoa A, Vanatoa E, Moora M, Davison J, Kalwij JM, et al. The online database MaarjAM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). New Phytol. 2010;188:223–41.CAS 
PubMed 
Article 

Google Scholar 
67.Katoh K, Asimenos G, Toh H. Multiple alignment of DNA sequences with MAFFT. In: Bioinforma. DNA Seq. Anal. Springer; 2009. p. 39–64.68.Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.CAS 
PubMed 
Article 

Google Scholar 
69.Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
70.Nilsson RH, Larsson KH, Taylor A, Bengtsson-Palme J, Jeppesen TS, Schigel D, et al. The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Res. 2018. https://doi.org/10.1093/nar/gky1022.71.Nguyen NH, Song Z, Bates ST, Branco S, Tedersoo L, Menke J, et al. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016;20:241–8.Article 

Google Scholar 
72.Conway JR, Lex A, Gehlenborg N. UpSetR: an R package for the visualization of intersecting sets and their properties. Bioinformatics. 2017;33:2938–40.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
73.Chytrý M, Tichý L, Holt J, Botta‐Dukát Z. Determination of diagnostic species with statistical fidelity measures. J Veg Sci. 2002;13:79–90.Article 

Google Scholar  More

In total 48 larval specimens of H. tenuicornis were obtained and analysed. We identified 30 larvae of the first instar, 14 of the second instar and 4 of the third instar. Two larvae and one adult specimen of H. tenuicornis were used for molecular phylogenetic placement of the genus within the subfamily Silphinae. The phylogenetic tree was obtained using Bayesian analysis from the concatenated partial 16S (434 bp) and COI (609 bp) sequences (Fig. 1).Figure 1Phylogenetic tree based on Bayesian analysis. Numbers above branches show the posterior probability and bootstrap values (BI)/maximal parsimony (PAUP)/Maximum likelihood (MEGA). Scaphidium quadrimaculatum Olivier, 1790 and Aleochara curtula (Goeze, 1777) (both Staphylinidae) were selected as outgroups.Full size imageSpecies identification based on genetic distancesThe calculated p-distances between concatenated sequences of 16S and COI of larval and adult specimens of H. tenuicornis were between 0.0029 and 0.0078 (the mean calculated p-distance within Heterotemna specimens was 0.01). Conversely, the distance between different species of Silpha was shown to be higher (mean calculated p-distance within the Silpha species was 0.08), thus the larval specimens were confirmed as belonging to the same species as the adult specimen, H. tenuicornis (SM1).Phylogenetic analysesThe Bayesian analysis (posterior probability 99), maximum parsimony bootstrap (84) and maximum likelihood bootstrap (93) strongly supported a clade of the genera Silpha, Heterotemna, Ablattaria and Phosphuga, suggesting close relationships of these genera with Heterotemna inside the genus Silpha, which makes the genus Silpha paraphyletic. The position of H. tenuicornis as a sister lineage to S. tristis Illiger, 1798 was strongly supported by the Bayesian analysis (97) but not strongly supported by the other analyses. The results confirmed the monophyly of the genera Thanatophilus Leach, 1815, Necrodes Leach, 1815, and Oiceoptoma Leach, 1815 within the subfamily Silphinae (Fig. 1).MorphometryThe two commonly used measurements for instar identification, head width and width of protergum , are applicable in the case of H. tenuicornis (Fig. 2c, d) as these two measurements do not overlap between the instars and show significant differences. More specifically, the following measurements were very different between instars; head width (F statistic = 231 on 2, df = 45, p value  More

Study sitesWe sampled ponds in the following five sites representing different regional-scale permafrost integrity: Toolik, Alaska, USA; Qeqertarsuaq, Disko Island, Greenland, Denmark; Whapmagoostui-Kuujjuarapik, Nunavik, Quebec, Canada; Abisko, Sweden and Khanymey, Western Siberia, Russia (Online-only Table 1). The aim was to include representatives of different stages of permafrost thaw in order to understand whether responses can be generalized across different geographic and environmental conditions.The sampling site in Alaska is located in a continuous permafrost area, mostly dominated by moss-tundra characterized by tussock-sedge Eriophorum vaginatum and Carex bigelowii, and dwarf-shrub Betula nana and Salix pulchra15. The average depth of the active layer in 2017 was ~50 cm16. Records of surface air temperature from 1989 to 2014 showed no significant warming trend, and there was no significant increase in the mean maximum thickness of the active layer or maximum thaw depth17.The sampling site in Greenland is located in the Blæsedalen Valley, south of Disko Island, and is characterized as a discontinuous permafrost area. From 1991 to 2011, Hollensen et al.18 observed an increase of the mean annual air temperatures of 0.2 °C per year in the area, while Hansen et al.19 highlighted that sea ice cover reduced 50% from 1991 to 2004. Soil temperatures recorded by the Arctic Station from the active layer of the coarse marine stratified sediments also showed an increase over the years18. The sampling site is comprised of wet sedge tundra, and the dominating species are Carex rariflora, Carex aquatilis, Eriophorum angustifolium, Equisetum arvense, Salix arctophila, Tomentypnum nitens and Aulacomnium turgidum20.The Canadian site is located within a sporadic permafrost zone, in a palsa bog, in the valley of Great Whale river, close to the river mouth to Hudson Bay. The vegetation consists of a coastal forest tundra, dominated by the species Carex sp. and Sphagnum sp.21 Since the mid-1990s, there has been a significant increase in the surface air temperature of the region for spring and fall, which has been correlated to a decline of sea ice coverage in Hudson Bay22. This area has experienced an accelerated thawing of the permafrost over the past decades, resulting in the collapse of palsas and the emergence of thermokarst ponds as well as significant peat accumulation21,23. In this specific site, thermokarst ponds at different development stage can be found, from recently emerging to older, mature thermokarstic waterbodies. The stage of the ponds was estimated based on the distance between the pond and the edge of the closest palsa, as well as based on satellite images14. The edges of the emerging ponds reached a maximum of 1 m from the closest palsa and were less than 0.5 m deep, whereas the edges of the developing ponds had a maximum distance of 2–3 m to the closest palsa and were ~1 m deep. Mature ponds were identified based on satellite images and were up to 60 years old.The Swedish site is located in a discontinuous permafrost zone at the Stordalen palsa mire, on an area of collapsed peatland affected by active thermokarst. The region has experienced an increase in mean annual air temperature and active layer thickness since the 1980s, which has been followed by a shift to wetter conditions24. The vegetation found on the surface of the palsa depressions of Stordalen mire is dominated by sedges (Eriophorum vaginatum, Carex sp.) and mosses (Sphagnum sp.)24,25.The Russian site is located in a discontinuous permafrost area in Western Siberia Lowland, near Khanymey village. The sampling site is a flat frozen palsa bog with a peat depth no more than 2 m, and is affected by active thermokarst, resulting in the emergence of thermokarst ponds26,27. The vegetation is dominated by lichens (Cladonia sp.), schrubs (Ledum palustre, Betula nana, Vaccinium vitis-idaea, Andromeda polifolia, Rubus chamaemorus) and mosses (Sphagnum sp.)28.Sample collectionAt all sites, water from the depth of 10 cm was collected from 12 ponds, totaling 60 ponds for the full dataset. Unfiltered water samples were collected for total P analysis. For analyzing Fe, various dissolved anions and cations, DOC concentrations, and perform optical and mass spectrometry analyses on DOM, water was filtered through GF/F glass fiber filters (0.7 μm, 47 mm, Whatman plc, Maidstone, United Kingdom). Moreover, water samples were collected in order to measure GHG (CO2 and CH4) concentrations. Water, detritus and sediment samples were also collected from ponds for fungal community analyses. Water samples were collected and filtered sequentially first through 5 µm Durapore membrane filter (Millipore, Burlington, Massachusetts, USA) and then through a 0.22 µm Sterivex filter (Millipore) to capture fungal cells of different sizes. The samples were filtered until clogging or up to a maximum of 3.5 liters (filtered volume ranging from 0.1 l to 3.5 l). Surface sediments were sampled from each of the ponds, with the exception of the Canadian site, where only one emerging and three developing ponds were sampled for sediments. From the sites in Alaska, Greenland, and Sweden, also detritus samples (dead plant material) were collected. The detritus was washed in the lab using tap water, followed by overnight incubation in 50 ml tap water to induce sporulation. The use of tap water may have added fungal spores to the samples, which should be kept in mind when using the detritus data. After the incubation, the water was filtered through a 5 μm pore size filter and the filter was stored at −20 °C.All the samples for DNA extraction were transported to the laboratory frozen, with the exception of the Alaskan samples, which were freeze dried prior to transportation. The samples transported frozen were freeze dried prior to DNA extraction to ensure similar treatment of all samples. The samples for nutrient and carbon measurements were transported frozen with the exception of samples for DOC and fluorescence analyses, which were transported cooled.Chemical analysesAll chemical, optical and mass spectrometry results are provided in OSF29. DOC quantification was carried out using a carbon analyzer (TOC-L + TNM-L, Shimadzu, Kyoto, Japan). Accuracy was assessed using EDTA at 11.6 mg C/l as a quality control (results were within + − 5%) and the standard calibration range was of 2–50 mg C/l. Fe(II) and Fe(III) were determined by using the ferrozine method30, but instead of reducing Fe(III) with hydroxylamine hydrochloride, ascorbic acid was used31. Absorbance was measured at 562 nm on a spectrophotometer (UV/Vis Spectrometer Lambda 40, Perkin Elmer, Waltham, Massachusetts, USA). The samples were diluted with milli-Q water if needed. The concentration of total P was determined using persulfate digestion32. The anion NO3− was measured on a Metrohm IC system (883 Basic IC Plus and 919 Autosampler Plus; Riverview, Florida, USA). NO3− were separated with a Metrosep A Supp 5 analytical column (250 × 4.0 mm) which was fit with a Metrosep A Supp 4/5 guard column at a flow rate of 0.7 ml/min, using a carbonate eluent (3.2 mM Na2CO3 + 1.0 mM NaHCO3). SO4 was analyzed using Metrohm IC system (883 Basic IC Plus and 919 Autosampler Plus, Riverview), NH4+ spectrophotometrically as described by Solórzano33, and NO2− and DN as in Greenberg et al.34.For the gas analyses, samples from Alaska and Canada were taken as previously described in Kankaala et al.35, except that room air was used instead of N2 for extracting the gas from the water. Shortly, 30 ml of water was taken into 50 ml syringes, which were warmed to room temperature prior to extraction of the gas. To each syringes 0.5 ml of HNO3 and 10 ml of room air was added and the syringes were shaken for 1 min. Finally, the volumes of liquid and gas phases were recorded and the gas was transferred into glass vials that had been flushed with N2 and vacuumed. For Greenland, Sweden and Russia 5 ml of water was taken for the gas samples with a syringe and immediately transferred to 20 ml glass vials filled with N and with 150 µL H2PO4 to preserve the sample. All gas samples were measured using gas chromatography (Clarus 500, Perkin Elmer, Polyimide Uncoated capillary column 5 m x 0.32 mm, TCD and FID detector respectively).Optical analysesIn order to characterize DOM, we recorded the absorbance of DOM using a UV-visible Cary 100 (Agilent Technologies, Santa Clara, California, USA) or a LAMBDA 40 UV/VIS (PerkinElmer) spectrophotometer, depending on sample origin. SUVA254 is a proxy of aromaticity and the relative proportion of terrestrial versus algal carbon sources in DOM36 and was determined from DOC normalized absorbance at 254 nm after applying a corrective factor based on iron concentration37. S289 enlights the importance of fulvic and humic acids related to algal production38 and were determined for the intervals 279–299 nm by performing regression calculations using SciLab v 5.5.2.39We also recorded fluorescence intensity on a Cary Eclipse spectrofluorometer (Agilent Technologies), across the excitation waveband from 250–450 nm (10 nm increments) and emission waveband of 300–560 nm (2 nm increments), or on a SPEX FluoroMax-2 spectrofluorometer (HORIBA, Kyoto, Japan), across the excitation waveband from 250–445 nm (5 nm increments) and emission waveband of 300–600 nm (4 nm increments), depending on sample origin. Based on the fluorometric scans, we constructed excitation-emission matrices (EEMs) after correction for Raman and Raleigh scattering and inner filter effect40. We calculated the FI as the ratio of fluorescence emission intensities at 450 nm and 500 nm at the excitation wavelength of 370 nm to investigate the origin of fulvic acids41. Higher values (~1.8) indicate microbial derived DOM (autochthonous), whereas lower values (~1.2) indicate terrestrial derived DOM (allochthonous), from plant or soil42. HIX is a proxy of the humic content of DOM and was calculated as the sum of intensity under the emission spectra 435–480 nm divided by the peak intensity under the emission spectra 300–445 nm, at an excitation of 250 nm. Higher values of HIX indicate more complex, higher molecular weight, condensed aromatic compounds43,44. BIX emphasizes the relative freshness of the bulk DOM and was calculated as the ratio of emission at 380 nm divided by the emission intensity maximum observed between 420 and 436 nm at an excitation wavelength of 310 nm45. High values ( >1) are related to higher proportion of more recently derived DOM, predominantly originated from autochthonous production, while lower values (0.6–0.7) indicate lower production and older DOM42,44.High resolution mass spectrometry50 ml water samples were collected from each of the ponds and were filtered with a Whatman GF/F filter for mass spectrometry analyses. For each sample, 1.5 ml of water was dried completely with a vacuum drier, and was then re-dissolved in 100 µL 20% acetonitrile, 80% water with three added compounds as internal standards (Hippuric acid, glycyrrhizic acid and capsaicin, all at 400 ppb v/v). Samples were filtered to an autosampler vials and injected at 50 µL onto the column. In order not to overload the detectors, some of the higher concentration samples were injected at a lower volume, to give a maximum of 20 µg carbon loaded.High-performance liquid chromatography – high resolution mass spectrometry (ESI-HRMS) was conducted as described in Patriarca et al.46 using a C18-Evo column (100 × 2.1 mm, 2.6 µm; Phenomenex, Torrance, California, USA). The ESI-HRMS data was averaged from 2–17 min to allow formula assignment to a single mass list. Formulas considered had masses 150–800 m/z, 4–50 carbon (C) atoms, 4–100 hydrogen (H) atoms, 1–40 oxygen (O) atoms, 0–1 nitrogen (N) atoms and 0–1 13 C atoms. Formulas were only considered if they had an even number of electrons, H/C 0.3–2.2 and O/C ≤ 1. The data are presented as a number of assigned formulas and weighted average O/C ratio, H/C ratio and m/z.The analysis was run in two batches (36 and 24 samples per run, respectively) and to the latter run, three samples of Suwannee River fulvic acid (SRFA, reference material) were added. At the moment of the run, the DOC concentration of these samples was unknown, so 50 µL was injected. From high resolution mass spectrometry, average H/C and a number of assigned formulas were obtained. The H/C can be used as a proxy of DOM aliphatic content; higher H/C values (  > 1) indicate more saturated (aliphatic) compounds, whereas values lower than 1 indicate more unsaturated, aromatic molecules47.DNA extraction, ITS2 amplification and sequencingAll samples for molecular analyses (water and detritus filters and sediments) were extracted using DNeasy PowerSoil® kit (Qiagen, Hilden, Germany), following the manufacturer’s recommendations for low input DNA. Extracts were eluted in 100 µl of Milli-Q water and DNA concentrations were measured with Qubit dsDNA HS kit. The fungal ribosomal internal transcribed spacer 2 (ITS2) sequences were amplified using a modified ITS3 Mix2 forward primer from Tedersoo48, named ITS3-mkmix2 CAWCGATGAAGAACGCAG, and a reverse primer ITS4 (equimolar mix of cwmix1 TCCTCCGCTTAyTgATAtGc and cwmix2 TCCTCCGCTTAtTrATAtGc)14. Each sample received a unique combination of primers containing identification tags generated by Barcrawl49. All tags had a minimum base difference of 3 and a length of 8 nucleotides. Both forward and reverse primer tags were extended by two terminal bases (CA) at the ligation site to avoid bias during ligation of sequencing adaptors, and the forward primer tag also had a linker base (T) added to it50. The list of primers and tags is found in Supplementary Table S1. PCR reactions were performed on a final volume of 50 µl, with an input amount of DNA ranging from 0.07 ng to 10 ng, 0.25 µM of each primer, 200 µM of dNTPs, 1U of Phusion™ High-Fidelity DNA Polymerase (Thermo Fisher Scientific, Waltham, Massachusetts, USA), 1X PhusionTM HF Buffer (1X buffer provides 1.5 mM MgCl2, Thermo Fisher Scientifics) and 0.015 mg of BSA. PCR conditions consisted of an initial denaturation cycle at 95 °C for 3 min, followed by 21–35 cycles for amplification (95 °C for 30 sec, 57 °C for 30 sec and 72 °C for 30 sec), and final extension at 72 °C for 10 min. In order to reduce PCR bias, all samples (in duplicates) were first submitted to 21 amplification cycles. In case of insufficient yield, the number of cycles was increased up to 35 cycles (see the records on the number of cycles for each of the samples in Supplementary Table S2).The PCR products were purified with Sera-MagTM beads (GE Healthcare Life Sciences, Marlborough, Massachusetts, USA), visualized on a 1.5% agarose gel and quantified using Qubit dsDNA HS kit. The purified PCR products were randomly allocated into three DNA pools (20 ng of each sample), which were purified with E.Z.N.A.® Cycle-Pure kit (Omega Bio-Tek, Norcross, Georgia, USA). Nine of the samples (4 water, 1 sediment and 4 detritus) were left out of the pools because of too little PCR product, giving a total of 203 samples for sequencing (Online-only Table 1). Negative PCR controls were added to each pool, as well as a mock community sample containing 10 different fragment sizes from the ITS2 region of a chimera of Heterobasidium irregular and Lophium mytilinum, ranging from 142 to 591 bases, as described by Castaño et al.51. The size distribution and quality of all the pools were verified with BioAnalyzer DNA 7500 (Agilent Technologies), and purity was assessed by spectrophotometry (OD 260:280 and 260:230 ratios) using NanoDrop (Thermo Fisher Scientific). The libraries were sequenced at Science for Life Laboratory (Uppsala University, Sweden), on a Pacific Biosciences Sequel instrument II, using 1 SMRT cell per pool. This PacBio technology allows the generation of highly accurate reads ( >99% accuracy) which are produced based on a consensus sequence after a circularization step.Quality filtering of reads, clustering and taxonomy identification of clustersThe sequencing resulted in a total of 1071489 sequences, ranging from 397 to 9184 sequences per sample (average on 2551 sequences per sample). The raw sequences were filtered for quality and clustered using the SCATA pipeline (https://scata.mykopat.slu.se/, accessed on May 19th, 2020). For quality filtering, sequences from each pool were screened for the primers and tags, requiring a minimum of 90% match for the primers and a 100% match for the tags. Reads shorter than 100 bp were removed, as well as reads with a mean quality lower than 20, or containing any bases with a quality lower than 7. After this filtering, 582234 sequences were retained in the data. The sequences were clustered at the species level by single-linkage clustering at a clustering distance of 1.5%, with penalties of 1 for mismatch, 0 for gap open, 1 for gap extension, and 0 for end gaps. Homopolymers were collapsed to 3 and unique genotypes across all pools were removed. For a preliminary taxonomy affiliation of the clusters, hereafter called OTUs (Operational Taxonomic Units), sequences from the UNITE + INSD dataset for Fungi52 database were included in the clustering process. After the clustering, the data included 518128 sequences, divided among 8218 OTUs. For taxonomical annotation, all OTUs with a minimum of ten total reads in the full dataset were included, retaining 3108 OTUs and 498414 sequences in the taxonomical analysis. More

1.Jaramillo-Villa, U., Maldonado-Ocampo, J. A. & Escobar, F. Altitudinal variation in fish assemblage diversity in streams of the central Andes of Colombia. J. Fish Biol. 76, 2401–2417 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
2.Mercado-silva, N., Lyons, J., Díaz-Pardo, E., Navarrete, S. & Gutiérrez-Hernández, A. Environmental factors associated with fish assemblage patterns in a high gradient river of the Gulf of Mexico slope. Revista Mexicana de Biodiversidad 83, 117–128 (2012).Article 

Google Scholar 
3.Cheng, D. et al. Quantifying the distribution and diversity of fish species along elevational gradients in the Weihe River Basin, Northwest China. Sustainability 11, 6177 (2019).Article 

Google Scholar 
4.Lorion, C. M., Kennedy, B. P. & Braatne, J. H. Altitudinal gradients in stream fish diversity and the prevalence of diadromy in the Sixaola River basin, Costa Rica. Environ. Biol. Fishes 91, 487–499 (2011).Article 

Google Scholar 
5.Li, J. et al. Spatial and temporal variation of fish assemblages and their associations to habitat variables in a mountain stream of north Tiaoxi River, China. Environ. Biol. Fishes 93, 403–417 (2012).Article 

Google Scholar 
6.Súarez, Y. R. et al. Patterns of species richness and composition of fish assemblages in streams of the Ivinhema River basin, Upper Paraná River. Acta Limnol. Bras. 23, 177–188 (2011).Article 

Google Scholar 
7.Vieira, T. B. & Tejerina-Garro, F. L. Relationships between environmental conditions and fish assemblages in tropical Savanna headwater streams. Sci. Rep. 10, 1–12 (2020).Article 
CAS 

Google Scholar 
8.Pokharel, K. K., Basnet, K. B., Majupuria, T. C. & Baniya, C. B. Correlations between fish assemblage structure and environmental variables of the Seti Gandaki River Basin, Nepal. J. Freshw. Ecol. 33, 31–43 (2018).CAS 
Article 

Google Scholar 
9.Carvajal-Quintero, J. D. et al. Variation in freshwater fish assemblages along a regional elevation gradient in the northern Andes, Colombia. Ecol. Evol. 5, 2608–2620 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
10.Li, J. et al. Climate and history explain the species richness peak at mid-elevation for Schizothorax fishes (Cypriniformes: Cyprinidae) distributed in the Tibetan Plateau and its adjacent regions. Glob. Ecol. Biogeogr. 18, 264–272 (2009).Article 

Google Scholar 
11.Fu, C., Wu, J., Chen, J., Wu, Q. & Lei, G. Freshwater fish biodiversity in the Yangtze River basin of China: Patterns, threats and conservation. Biodivers. Conserv. 12, 1649–1685 (2003).Article 

Google Scholar 
12.Orrego, R., Adams, S. M., Barra, R., Chiang, G. & Gavilan, J. F. Patterns of fish community composition along a river affected by agricultural and urban disturbance in south-central Chile. Hydrobiologia 620, 35–46 (2009).Article 

Google Scholar 
13.Nyanti, L. et al. Acidification tolerance of Barbonymus schwanenfeldii (Bleeker, 1854) and Oreochromis niloticus (Linnaeus, 1758)—Implication of fish size. AACL Bioflux 10, 746–753 (2017).
Google Scholar 
14.Nyanti, L. et al. Effects of water temperature, dissolved oxygen and total suspended solids on juvenile Barbonymus schwanenfeldii (Bleeker, 1854) and Oreochromis niloticus (Linnaeus, 1758). AACL Bioflux 11, 394–406 (2018).
Google Scholar 
15.Ling, T. Y. et al. Assessment of the water and sediment quality of tropical forest streams in upper reaches of the Baleh River, Sarawak, Malaysia, subjected to logging activities. J. Chem. 2016, 1–13 (2016).CAS 

Google Scholar 
16.Davies, P. & Nelson, M. Relationships between riparian buffer widths and the effects of logging on stream habitat, invertebrate community composition and fish abundance. Mar. Freshw. Res. 45, 1289–1305 (1994).Article 

Google Scholar 
17.Ikhwanuddin, M., Amal, M., Shohaimi, S., Hasan, H. & Jamil, N. Environmental influences on fish assemblages of the Upper Sungai Pelus, Kuala Kangsar, Perak, Malaysia. Sains Malaysiana 45, 1487–1495 (2016).CAS 

Google Scholar 
18.Zainuddin, Z., Jamal, P. & Akbar, I. Modeling the effect of dam construction and operation towards downstream water quality of Sg. Tawau and Batang Baleh. World J. Appl. Environ. Chem. 1, 57–66 (2012).
Google Scholar 
19.Nyanti, L., Ling, T. & Muan, T. Water quality of Bakun Hydroelectric Dam Reservoir, Sarawak, Malaysia, during the construction of Murum Dam. ESTEEM Acad. J. 11, 81–88 (2015).
Google Scholar 
20.Ling, T. Y. et al. Changes in water and sediment quality of a river being impounded and differences among functional zones of the new large tropical hydroelectric reservoir. Pol. J. Environ. Stud. 28, 4271–4285 (2019).CAS 
Article 

Google Scholar 
21.Osman, N. B., Othman, H. T., Karim, R. A. & Mazlan, M. A. F. Biomass in Malaysia: Forestry-based residues. Int. J. Biomass Renew. 3, 7–14 (2014).
Google Scholar 
22.Inger, R. F. & Chin, P. K. Freshwater Fish of North Borneo (Natural History Publications, 2002).
Google Scholar 
23.Mohsin, A. K. M. & Ambak, M. A. Freshwater fishes of Peninsular Malaysia (Universiti Pertanian Malaysia, 1983).
Google Scholar 
24.Kottelat, M. The fishes of the inland waters of Southeast Asia: A catalogue and core bibliography of the fishes known to occur in freshwaters, mangroves and estuaries. Raffles Bull. Zool. 27, 1–663 (2013).
Google Scholar 
25.Kottelat, M. Conspectus Cobitidum: An inventory of the loaches of the world (Teleostei: Cypriniformes: Cobitoidei). Raffles Bull. Zool. 26, 1–199 (2012).
Google Scholar 
26.Kottelat, M. & Tan, H. H. A synopsis of the genus Lobocheilos in Java, Sumatra and Borneo, with descriptions of six new species (Teleostei: Cyprinidae). Ichthyol. Explor. Freshw. 19, 27–58 (2008).
Google Scholar 
27.Froese, R. & Pauly, D. FishBase. World Wide Web electronic publication. https://www.fishbase.se/search.php (2019).28.van der Laan, R., Fricke, R. & Eschmeyer, W. N. Eschmeyer’s Catalog of Fishes: Classification. http://www.calacademy.org/scientists/catalog-of-fishes-classification/ (2020).29.Shannon, C. E. & Weaver, W. The Mathematical Theory of Communication (The University of Illinois Press, 1964).MATH 

Google Scholar 
30.Margalef, R. Perspectives in Ecological Theory (University of Chicago Press, 1968).
Google Scholar 
31.Pielou, E. C. Species diversity and pattern diversity in the study of ecological succession. J. Theor. Biol. 10, 370–383 (1966).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
32.Ter Braak, C. J. F. Canonical correspondence analysis: A new eigenvector technique for multivariate direct gradient analysis. Ecology 67, 1167–1179 (1986).Article 

Google Scholar 
33.Ter Braak, C. J. F. & Verdonschot, P. F. M. Canonical correspondence analysis and related multivariate methods in aquatic ecology. Aquat. Sci. 57, 255–289 (1995).Article 

Google Scholar 
34.Ward-Campbell, B. M. S., Beamish, F. W. H. & Kongchaiya, C. Morphological characteristics in relation to diet in five coexisting Thai fish species. J. Fish Biol. 67, 1266–1279 (2005).Article 

Google Scholar 
35.Beamish, F. W. H., Sa-ardrit, P. & Tongnunui, S. Habitat characteristics of the cyprinidae in small rivers in Central Thailand. Environ. Biol. Fishes 76, 237–253 (2006).Article 

Google Scholar 
36.Muchlisin, Z. A. & Siti Azizah, M. N. Diversity and distribution of freshwater fishes in Aceh water, northern Sumatra, Indonesia. Int. J. Zool. Res. 5, 62–79 (2009).Article 

Google Scholar 
37.Rashid, Z. A., Asmuni, M. & Amal, M. N. A. Fish diversity of Tembeling and Pahang rivers, Pahang, Malaysia. Check List 11, 1–6 (2015).Article 

Google Scholar 
38.Suvarnaraksha, A., Lek, S., Lek-Ang, S. & Jutagate, T. Fish diversity and assemblage patterns along the longitudinal gradient of a tropical river in the Indo-Burma hotspot region (Ping-Wang River Basin, Thailand). Hydrobiologia 694, 153–169 (2012).CAS 
Article 

Google Scholar 
39.Kottelat, M. Conspectus cobitidum: An inventory of the loaches of the world (Teleostei: Cypriniformes: Cobitoidei). Raffles Bull. Zool. 26, 1–199 (2012).
Google Scholar 
40.Tan, H. H. The Borneo suckers. Revision of the Torrent Loaches of Borneo (Balitoridae: Gastromyzon, Neogastromyzon) (Natural History Publications, 2006).
Google Scholar 
41.Beamish, F. W. H., Sa-Ardrit, P. & Cheevaporn, V. Habitat and abundance of Balitoridae in small rivers of central Thailand. J. Fish Biol. 72, 2467–2484 (2008).Article 

Google Scholar 
42.Ahmad, A., Nek, S. A. R. T. & Ambak, M. A. Preliminary study on fish diversity of Ulu Tungud, Meliau range, Sandakan, Sabah. J. Sustain. Sci. Manag. 1, 21–26 (2006).
Google Scholar 
43.Odum, E. P. & Barret, G. W. Fundamental of Ecology (Cengage Learning, Inc, 2004).
Google Scholar 
44.Magurran, A. E. Ecological Diversity and Its Measurement (Princeton University Press, 1988).Book 

Google Scholar 
45.Au, D. W. T. et al. Chronic effects of suspended solids on gill structure, osmoregulation, growth, and triiodothyronine in juvenile green grouper Epinephelus coioides. Mar. Ecol. Prog. Ser. 266, 255–264 (2004).ADS 
CAS 
Article 

Google Scholar 
46.Kimbell, H. S. & Morrell, L. J. Turbidity influences individual and group level responses to predation in guppies, Poecilia reticulata. Anim. Behav. 103, 179–185 (2015).Article 

Google Scholar 
47.Li, W. et al. Effects of turbidity and light intensity on foraging success of juvenile mandarin fish Siniperca chuatsi (Basilewsky). Environ. Biol. Fishes 96, 995–1002 (2013).Article 

Google Scholar 
48.Kukula, K. & Bylak, A. Synergistic impacts of sediment generation and hydrotechnical structures related to forestry on stream fish communities. Sci. Total Environ. 737, 139751 (2020).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
49.Krause, K. P., Wu, C. L., Chu, M. L. & Knouft, J. H. Fish assemblage–environment relationships suggest differential trophic responses to heavy metal contamination. Freshw. Biol. 64, 632–642 (2019).CAS 
Article 

Google Scholar 
50.Askeyev, A. et al. River fish assemblages along an elevation gradient in the eastern extremity of Europe. Environ. Biol. Fishes 100, 585–596 (2017).Article 

Google Scholar 
51.Zamani Faradonbe, M. & Eagderi, S. Fish assemblages as influenced by environmental factors in Taleghan River (the Caspian Sea basin, Alborz Province, Iran). Caspian J. Environ. Sci. 13, 363–371 (2015).
Google Scholar 
52.Bolner, K. C. S., Copatti, C. E., Rosso, F. L., Loro, V. L. & Baldisserotto, B. Water pH and metabolic parameters in silver catfish (Rhamdia quelen). Biochem. Syst. Ecol. 56, 202–208 (2014).CAS 
Article 

Google Scholar 
53.Abbink, W. et al. The effect of temperature and pH on the growth and physiological response of juvenile yellowtail kingfish Seriola lalandi in recirculating aquaculture systems. Aquaculture 330–333, 130–135 (2012).Article 
CAS 

Google Scholar 
54.Paller, V. G. V., Corpuz, M. N. C. & Ocampo, P. P. Diversity and distribution of freshwater fish assemblages in Tayabas River, Quezon (Philippines). Philip. J. Sci. 142, 55–67 (2013).
Google Scholar 
55.Jeppesen, R. et al. Effects of hypoxia on fish survival and oyster growth in a highly eutrophic estuary. Estuaries Coasts 41, 89–98 (2018).CAS 
Article 

Google Scholar 
56.Rosso, J. J. & Quirós, R. Patterns in fish species composition and assemblage structure in the upper Salado river lakes, Pampa Plain, Argentina. Neotrop. Ichthyol. 8, 135–144 (2010).Article 

Google Scholar 
57.Batzer, D. P., Jackson, C. R. & Mosner, M. Influences of riparian logging on plants and invertebrates in small, depressional wetlands of georgia, U.S.A.. Hydrobiologia 441, 123–132 (2000).Article 

Google Scholar 
58.Cheimonopoulou, M. T., Bobori, D. C., Theocharopoulos, I. & Lazaridou, M. Assessing ecological water quality with macroinvertebrates and fish: A case study from a small mediterranean river. Environ. Manag. 47, 279–290 (2011).ADS 
Article 

Google Scholar 
59.Roberts, T. R. The Freshwater Fishes of Western Borneo (Kalimantan Barat, Indonesia) (California Academy of Science, 1989).
Google Scholar 
60.Tan, H. H. & Leh, C. U. M. Three new species of Gastromyzon (Teleostei: Balitoridae) from southern Sarawak. Zootaxa 19, 1–19 (2006).
Google Scholar 
61.Tan, H. H. & Martin-Smith, K. M. Two new species of Gastromyzon (Teleostei: Balitoridae) from the Kuamut headwaters, Kinabatangan basin, Sabah, Malaysia. Raffles Bull. Zool. 46, 361–371 (1998).
Google Scholar  More

Active ingredients detected in bee collected pollenAll 188 pollen samples had at least 12 active ingredients detected in each sample, with a maximum of 31 AIs and an average of 22.0 ± 0.3 per sample. Over both years, 80 of the 259 screened pesticide active ingredients were detected in the pollen. These included 28 fungicides, 26 insecticides, 21 herbicides, two miticides, and one rodenticide. We also detected one synthetic antioxidant and one pesticide synergist (Table S1). We detected approximately twice as many AIs in pollen collected by honey bees (68 AIs) in 2019 than in pollen collected by bumble bees (32). All AIs detected in pollen from bumble bees were also collected by honey bees in either 2018 or 2019. Honey bee collected pollen also had significantly more AIs on average detected at each site (35.0 ± 0.9 S.E. AIs per site) compared to bumble bees (18.6 ± 0.6) in 2019 (R2m = 0.73; X2 = 68.2, df = 1, p  More

1.Zielińska, K. M., Kiedrzynski, M., Grzyl, A. & Rewicz, A. Forest roadsides harbour less competitive habitats for a relict mountain plant (Pulsatilla vernalis) in lowlands. Sci. Rep. https://doi.org/10.1038/srep31913 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
2.Zielińska, K. M., Kiedrzynski, M., Grzyl, A. & Tomczyk, P. P. Anthropogenic sites maintain the last individuals during the rapid decline of the lowland refugium of the alpine-arctic plant Pulsatilla vernalis (L.) Mill. Pak. J. Bot. 50, 211–215 (2018).3.Grzyl, A. & Ronikier, M. Pulsatilla vernalis (Ranunculaceae) in the Polish lowlands: Current population resources of a declining species. Pol. Bot. J. 56, 185–194 (2011).
Google Scholar 
4.Åström, S. & Stridh, B. The present status of Pulsatilla vernalis in Sweden. Sven. Bot. Tidskr. 97, 117–126 (2003).
Google Scholar 
5.Chappuis, E. Pulsatilla vernalis. The IUCN Red List of Threatened Species 2014: e.T55730086A55730098. (2014). https://doi.org/10.2305/IUCN.UK.2014-1.RLTS.T55730086A55730098.en. Downloaded on 02 December 2020 >.6.Ronikier, M. et al. Phylogeography of Pulsatilla vernalis (L.) Mill. (Ranunculaceae): Chloroplast DNA reveals two evolutionary lineages across central Europe and Scandinavia. J. Biogeogr. 35, 1650–1664. https://doi.org/10.1111/j.1365-2699.2008.01907.x (2008).7.Kiedrzyński, M., Zielińska, K. M., Kiedrzyńska, E. & Rewicz, A. Refugial debate: On small sites according to their function and capacity. Evol. Ecol. 31, 815–827. https://doi.org/10.1007/s10682-017-9913-4 (2017).Article 

Google Scholar 
8.Betz, C., Scheuerer, M. & Reisch, C. Population reinforcement—A glimmer of hope for the conservation of the highly endangered Spring Pasque flower (Pulsatilla vernalis). Biol. Conserv. 168, 161–167. https://doi.org/10.1016/j.biocon.2013.10.004 (2013).Article 

Google Scholar 
9.Nawrocka-Grześkowiak, U. & Frydel, K. Spring pasque-flower (Pulsatilla vernalis (L.) Miller) localities in the Kaliska Forest District. Zarządzanie Ochroną Przyrody w Lasach 6, 77–84 (2012).10.Gutterman, Y. In Seeds: The Ecology of Regeneration in Plant Communities (ed M. Fenner) 59–84 (CAB International, 2000).11.Luzuriaga, A. L., Escudero, A. & Perez-Garcia, F. Environmental maternal effects on seed morphology and germination in Sinapis arvensis (Cruciferae). Weed Res. 46, 163–174. https://doi.org/10.1111/j.1365-3180.2006.00496.x (2006).Article 

Google Scholar 
12.Rao, N. K., Dulloo, M. E. & Engels, J. M. M. A review of factors that influence the production of quality seed for long-term conservation in genebanks. Genet. Resour. Crop Evol. 64, 1061–1074. https://doi.org/10.1007/s10722-016-0425-9 (2017).CAS 
Article 

Google Scholar 
13.Doniak, M., Barciszewska, M. Z., Kaźmierczak, J. & Kaźmierczak, A. The crucial elements of the ‘last step’ of programmed cell death induced by kinetin in root cortex of V. faba ssp. minor seedlings. Plant Cell Rep. 33, 2063–2076. https://doi.org/10.1007/s00299-014-1681-9 (2014).14.Doniak, M., Byczkowska, A. & Kaźmierczak, A. Kinetin-induced programmed death of cortex cells is mediated by ethylene and calcium ions in roots of Vicia faba ssp minor. Plant Growth Regul. 78, 335–343. https://doi.org/10.1007/s10725-015-0096-0 (2016).CAS 
Article 

Google Scholar 
15.Doniak, M., Kaźmierczak, A., Byczkowska, A. & Glińska, S. Reactive oxygen species and sugars may be the messengers in kinetin-induced death of field bean root cortex cells. Biol. Plant. 61, 178–186. https://doi.org/10.1007/s10535-016-0654-y (2017).CAS 
Article 

Google Scholar 
16.Tudela-Isanta, M. et al. Habitat-related seed germination traits in alpine habitats. Ecol. Evol. 8, 150–161. https://doi.org/10.1002/ece3.3539 (2018).Article 
PubMed 

Google Scholar 
17.Baskin, J. M. & Baskin, C. C. A classification system for seed dormancy. Seed Sci. Res. 14, 1–16. https://doi.org/10.1079/ssr2003150 (2004).ADS 
Article 

Google Scholar 
18.Finch-Savage, W. E. & Leubner-Metzger, G. Seed dormancy and the control of germination. New Phytol. 171, 501–523. https://doi.org/10.1111/j.1469-8137.2006.01787.x (2006).CAS 
Article 
PubMed 

Google Scholar 
19.Latrasse, D., Benhamed, M., Bergounioux, C., Raynaud, C. & Delarue, M. Plant programmed cell death from a chromatin point of view. J. Exp. Bot. 20, 5887–5900 (2016).Article 

Google Scholar 
20.Baskin, J. M., Baskin, C. C. & Li, X. Taxonomy, anatomy and evolution of physical dormancy in seeds. Plant Species Biol. 15, 139–152. https://doi.org/10.1046/j.1442-1984.2000.00034.x (2000).Article 

Google Scholar 
21.Grzyl, A. Biology and Ecology of Isolated Populations of Pulsatilla vernalis (L.) Mill. on the Eastern Limits of its RANGE in Poland. (PhD thesis. University of Lodz, Department of Geobotany and Plant Ecology, 2012).22.Grzyl, A., Kiedrzynski, M., Zielinska, K. M. & Rewicz, A. The relationship between climatic conditions and generative reproduction of a lowland population of Pulsatilla vernalis: The last breath of a relict plant or a fluctuating cycle of regeneration?. Plant Ecol. 215, 457–466. https://doi.org/10.1007/s11258-014-0316-0 (2014).Article 

Google Scholar 
23.Oostermeijer, J. G. B., Vaneijck, M. W. & Dennijs, J. C. M. Offspring fitness in relation to population size and genetic variation in the rare perennial plant species Gentiana pneumonanthe (Gentianaceae). Oecologia 97, 289–296. https://doi.org/10.1007/bf00317317 (1994).ADS 
CAS 
Article 
PubMed 

Google Scholar 
24.Ouborg, N. J. & Vantreuren, R. Variation in fitness-related characters among small and large populations of Salvia pratensis. J. Ecol. 83, 369–380. https://doi.org/10.2307/2261591 (1995).Article 

Google Scholar 
25.Fischer, M. & Matthies, D. RAPD variation in relation to population size and plant fitness in the rare Gentianella germanica (Gentianaceae). Am. J. Bot. 85, 811–819. https://doi.org/10.2307/2446416 (1998).CAS 
Article 
PubMed 

Google Scholar 
26.Frankham, R., Ballou, J. D. & Briscoe, D. A. Introduction to Conservation Genetics. (Cambridge University Press, 2002).27.Hensen, I., Oberprieler, C. & Wesche, K. Genetic structure, population size, and seed production of Pulsatilla vulgaris Mill. (Ranunculaceae) in Central Germany. Flora 200, 3–14. https://doi.org/10.1016/j.flora.2004.05.001 (2005).28.Jakobsson, A. & Eriksson, O. A comparative study of seed number, seed size, seedling size and recruitment in grassland plants. Oikos 88, 494–502. https://doi.org/10.1034/j.1600-0706.2000.880304.x (2000).Article 

Google Scholar 
29.Melser, C. & Klinkhamer, P. G. L. Selective seed abortion increases offspring survival in Cynoglossum officinale (Boraginaceae). Am. J. Bot. 88, 1033–1040. https://doi.org/10.2307/2657085 (2001).CAS 
Article 
PubMed 

Google Scholar 
30.Meyer, K. M., Soldaat, L. L., Auge, H. & Thulke, H. H. Adaptive and selective seed abortion reveals complex conditional decision making in plants. Am. Nat. 183, 376–383. https://doi.org/10.1086/675063 (2014).Article 
PubMed 

Google Scholar 
31.Bochenková, M., Hejcman, M. & Karlík, P. Effect of plant community on recruitment of Pulsatilla pratensis in dry grassland. Sci. Agric. Bohem. 2012, 127–133. https://doi.org/10.7160/sab.2012.430402 (2012).Article 

Google Scholar 
32.Ghazoul, J. & Satake, A. Nonviable seed set enhances plant fitness: The sacrificial sibling hypothesis. Ecology 90, 369–377. https://doi.org/10.1890/07-1436.1 (2009).Article 
PubMed 

Google Scholar 
33.Laitinen, P. The Effects of Forest Fires on the Persistence of Pulsatilla vernalis (L.) Mill. edn, (Ms. thesis, University of Jyväskylä, Faculty of Mathematics and Science, Department of Biological and Environmental Science, 2008) [in Finnish with an English abstract].34.Skalická, R., Karlík, P., Hejcman, M. & Bochenková, M. In 17th Symposium of the European Grassland Federation. 388–390.35.Arathi, H. S., Ganeshaiah, K. N., Shaanker, R. U. & Hedge, S. G. Seed abortion in Pongamia pinnata (Fabaceae). Am. J. Bot. 86, 659–662. https://doi.org/10.2307/2656574 (1999).CAS 
Article 
PubMed 

Google Scholar 
36.Brookes, R. H., Jesson, L. K. & Burd, M. A test of simultaneous resource and pollen limitation in Stylidium armeria. New Phytol. 179, 557–565. https://doi.org/10.1111/j.1469-8137.2008.02453.x (2008).Article 
PubMed 

Google Scholar 
37.Yang, C. F., Sun, S. G. & Guo, Y. H. Resource limitation and pollen source (self and outcross) affecting seed production in two louseworts, Pedicularis siphonantha and P. longiflora (Orobanchaceae). Bot. J. Linn. Soc. 147, 83–89. https://doi.org/10.1111/j.1095-8339.2005.00363.x (2005).38.Cendán, C., Sampedro, L. & Zas, R. The maternal environment determines the timing of germination in Pinus pinaster. Environ. Exp. Bot. 94, 66–72. https://doi.org/10.1016/j.envexpbot.2011.11.022 (2013).Article 

Google Scholar 
39.Li, R. et al. Effects of cultivar and maternal environment on seed quality in Vicia sativa. Front. Plant Sci. 8. https://doi.org/10.3389/fpls.2017.01411 (2017).40.Valencia-Diaz, S. & Montaña, C. Temporal variability in the maternal environment and its effect on seed size and seed quality in Flourensia cernua DC. (Asteraceae). J. Arid Environ. 63, 686–695. https://doi.org/10.1016/j.jaridenv.2005.03.024 (2005).41.Chinnusamy, V., Gong, Z. Z. & Zhu, J. K. Abscisic acid-mediated epigenetic processes in plant development and stress responses. J. Integr. Plant Biol. 50, 1187–1195. https://doi.org/10.1111/j.1744-7909.2008.00727.x (2008).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
42.Butuzova, O. G. Peculiarities of seed formation in Pulsatilla vulgaris and Helleborus niger (Ranunculaceae) with embryo postdevelopment. Botanicheskii Zhurnal (St. Petersburg) 103, 313—330 (2018) [in Russian].43.Duncan, C., Schultz, N., Lewandrowski, W., Good, M. K. & Cook, S. Lower dormancy with rapid germination is an important strategy for seeds in an arid zone with unpredictable rainfall. PLoS ONE 14. https://doi.org/10.1371/journal.pone.0218421 (2019).44.Gremer, J. R., Kimball, S. & Venable, D. L. Within and among year germination in Sonoran Desert winter annuals: bet hedging and predictive germination in a variable environment. Ecol. Lett. 19, 1209–1218. https://doi.org/10.1111/ele.12655 (2016).Article 
PubMed 

Google Scholar 
45.Venable, D. L. Bet hedging in a guild of desert annuals. Ecology 88, 1086–1090. https://doi.org/10.1890/06-1495 (2007).Article 
PubMed 

Google Scholar 
46.Evans, M. E. K. & Dennehy, J. J. Germ banking: Bet-hedging and variable release from egg and seed dormancy. Q. R. Biol. 80, 431–451. https://doi.org/10.1086/498282 (2005).Article 

Google Scholar 
47.Goldberg, R. B., de Paiva, G. & Yadegari, R. Plant embryogenesis – zygote to seed. Science 266, 605–614. https://doi.org/10.1126/science.266.5185.605 (1994).ADS 
CAS 
Article 
PubMed 

Google Scholar 
48.Lester, R. N. & Kang, J. H. Embryo and endosperm function and failure in Solanum species and hybrids. Ann. Bot. 82, 445–453. https://doi.org/10.1006/anbo.1998.0695 (1998).Article 

Google Scholar 
49.Lopes, M. A. & Larkins, B. A. Endosperm origin, development, and function. Plant Cell 5, 1383–1399. https://doi.org/10.1105/tpc.5.10.1383 (1993).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
50.Yan, D. W., Duermeyer, L., Leoveanu, C. & Nambara, E. The functions of the endosperm during seed germination. Plant Cell Physiol. 55, 1521–1533. https://doi.org/10.1093/pcp/pcu089 (2014).CAS 
Article 
PubMed 

Google Scholar 
51.Willis, C. G. et al. The evolution of seed dormancy: Environmental cues, evolutionary hubs, and diversification of the seed plants. New Phytol. 203, 300–309. https://doi.org/10.1111/nph.12782 (2014).Article 
PubMed 

Google Scholar 
52.Poisot, T., Bever, J. D., Nemri, A., Thrall, P. H. & Hochberg, M. E. A conceptual framework for the evolution of ecological specialisation. Ecol. Lett. 14, 841–851. https://doi.org/10.1111/j.1461-0248.2011.01645.x (2011).Article 
PubMed 
PubMed Central 

Google Scholar 
53.Pfeifer, E., Holderegger, R., Matthies, D. & Rutishauser, R. Investigation on the population biology of a flagship species of dry meadows: Pulsatilla vulgaris Mill. in north-eastern Switzerland. Bot. Helvet. 112, 153–171 (2002).54.Gargiulo, R. et al. Conservation of the threatened species, Pulsatilla vulgaris Mill. (Pasqueflower), is aided by reproductive system and polyploidy. J. Hered. 110, 618–628. https://doi.org/10.1093/jhered/esz035 (2019).55.Seglias, A. E., Williams, E., Bilge, A. & Kramer, A. T. Phylogeny and source climate impact seed dormancy and germination of restoration-relevant forb species. PLoS ONE 13. https://doi.org/10.1371/journal.pone.0191931 (2018).56.Byczkowska, A., Kunikowska, A. & Kaźmierczak, A. Determination of ACC-induced cell-programmed death in roots of Vicia faba ssp. minor seedlings by acridine orange and ethidium bromide staining. Protoplasma 250, 121–128. https://doi.org/10.1007/s00709-012-0383-9 (2013).57.Dray, S. & Dufour, A. B. The ade4 package: Implementing the duality diagram for ecologists. J. Stat. Softw. 22, 1–20. https://doi.org/10.18637/jss.v022.i04 (2007).58.Kassambara, A. & Mundt, F. factoextra: Extract and Visualize the Results of Multivariate Data Analyses. R package version 1.0.7. https://CRAN.R-project.org/package=factoextra (2020).59.Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer, 2016).Book 

Google Scholar 
60.Oksanen, J. et al. vegan: Community Ecology Package. R package version 2.5-6. https://CRAN.R-project.org/package=vegan (2019).61.Fox, F. & Weisberg, S. An {R} Companion to Applied Regression, Third Edition. (Sage, 2019). https://socialsciences.mcmaster.ca/jfox/Books/Companion/. More

1.Schipper, J. et al. The status of the world’s land and marine mammals: Diversity, threat, and knowledge. Science 322, 225–230 (2008).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
2.Cushman, S. A. Effects of habitat loss and fragmentation on amphibians: A review and prospectus. Biol. Conserv. 128, 231–240 (2006).Article 

Google Scholar 
3.Díaz, S. et al. Pervasive human-driven decline of life on Earth points to the need for transformative change. Science 366, 6471 (2019).Article 
CAS 

Google Scholar 
4.Geary, W. L., Nimmo, D. G., Doherty, T. S., Ritchie, E. G. & Tulloch, A. I. T. Threat webs: Reframing the co-occurrence and interactions of threats to biodiversity. J. Appl. Ecol. 56, 1992–1997 (2019).
Google Scholar 
5.Halpern, B. S. et al. A global map of human impact on marine ecosystems. Science 319, 948–952 (2008).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
6.Johnson, C. N. et al. Biodiversity losses and conservation responses in the Anthropocene. Science 356, 270–275 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
7.Di Prisco, G. et al. A mutualistic symbiosis between a parasitic mite and a pathogenic virus undermines honey bee immunity and Health. Proc. Natl. Acad. Sci. USA. 113, 3203–3208 (2016).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
8.Grant, E. H. C. et al. Identifying management-relevant research priorities for responding to disease-associated amphibian declines. Glob. Ecol. Conserv. 16, 00441 (2018).
Google Scholar 
9.Schindler, D. W. The cumulative effects of climate warming and other human stresses on Canadian freshwaters in the new millennium. Can. J. Fish. Aquat. Sci. 58, 18–29 (2001).Article 

Google Scholar 
10.Cumulative Effects in Wildlife Management: Impact Mitigation. https://doi.org/10.1017/CBO9781107415324.004 (CRC Press, 2011). 11.Frid, A. & Dill, L. Human-caused disturbance stimuli as a form of predation risk. Conserv. Ecol. 6, 11 (2002).
Google Scholar 
12.Rolland, R. M., Hunt, K. E., Kraus, S. D. & Wasser, S. K. Assessing reproductive status of right whales (Eubalaena glacialis) using fecal hormone metabolites. Gen. Comp. Endocrinol. 142, 308–317 (2005).CAS 
PubMed 
Article 

Google Scholar 
13.Sheriff, M. J., Dantzer, B., Delehanty, B., Palme, R. & Boonstra, R. Measuring stress in wildlife: Techniques for quantifying glucocorticoids. Oecologia 166, 869–887 (2011).ADS 
PubMed 
Article 

Google Scholar 
14.Crain, C. M., Kroeker, K. & Halpern, B. S. Interactive and cumulative effects of multiple human stressors in marine systems. Ecol. Lett. 11, 1304–1315 (2008).PubMed 
Article 

Google Scholar 
15.Beal, A., Rodriguez-Casariego, J., Rivera-Casas, C., Suarez-Ulloa, V. & Eirin-Lopez, J. M. Environmental epigenomics and its applications in marine organisms. in Population Genomics: Marine Organisms (eds. Oleksiak, M. F. & Rajora, O. P.) 325–359. https://doi.org/10.1007/13836_2018_28 (Springer, 2018). 16.Eirin-Lopez, J. M. & Putnam, H. M. Marine environmental epigenetics. Ann. Rev. Mar. Sci. 11, 335–368 (2019).PubMed 
Article 

Google Scholar 
17.Laird, P. W. The power and the promise of DNA methylation markers. Nat. Rev. Cancer 3, 253–266 (2003).CAS 
PubMed 
Article 

Google Scholar 
18.Bird, A. P. CpG-rich islands and the function of DNA methylation. Nature 321, 209–213 (1986).ADS 
CAS 
PubMed 
Article 

Google Scholar 
19.Cedar, H. & Bergman, Y. Linking DNA methylation and histone modification: patterns and paradigms. Nat. Rev. Genet. 10, 295–304 (2009).CAS 
PubMed 
Article 

Google Scholar 
20.Matosin, N., Cruceanu, C. & Binder, E. B. Preclinical and clinical evidence of DNA methylation changes in response to trauma and chronic stress. Chronic Stress 1, 247054701771076 (2017).Article 

Google Scholar 
21.Radtke, K. M. et al. Transgenerational impact of intimate partner violence on methylation in the promoter of the glucocorticoid receptor. Transl. Psychiatry 1, e21–e26 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
22.Mueller, B. R. & Bale, T. L. Sex-specific programming of offspring emotionality after stress early in pregnancy. J. Neurosci. 28, 9055–9065 (2008).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
23.Elliott, E., Ezra-Nevo, G., Regev, L., Neufeld-Cohen, A. & Chen, A. Resilience to social stress coincides with functional DNA methylation of the Crf gene in adult mice. Nat. Neurosci. 13, 1351–1353 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
24.Boersma, G. J. et al. Prenatal stress decreases Bdnf expression and increases methylation of Bdnf exon IV in rats. Epigenetics 9, 437–447 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
25.Turecki, G. & Meaney, M. J. Effects of the social environment and stress on glucocorticoid receptor gene methylation: A systematic review. Biol. Psychiatry 79, 87–96 (2016).CAS 
PubMed 
Article 

Google Scholar 
26.Sterrenburg, L. et al. Chronic stress induces sex-specific alterations in methylation and expression of corticotropin-releasing factor gene in the rat. PLoS ONE 6, 1–14 (2011).Article 
CAS 

Google Scholar 
27.Reeder, D. A. M. & Kramer, K. M. Stress in free-ranging mammals: Integrating physiology, ecology, and natural history. J. Mammal. 86, 225–235 (2005).Article 

Google Scholar 
28.Jeanneteau, F. D. et al. BDNF and glucocorticoids regulate corticotrophin-releasing hormone (CRH) homeostasis in the hypothalamus. Proc. Natl. Acad. Sci. USA. 109, 1305–1310 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
29.Smith, S. M. & Vale, W. W. The role of the hypothalamic-pituitary-adrenal axis in neuroendocrine responses to stress. Dialogues Clin. Neurosci. 8, 383–395 (2006).PubMed 
PubMed Central 
Article 

Google Scholar 
30.Turner, J. D. & Muller, C. P. Structure of the glucocorticoid receptor (NR3C1) gene 5′ untranslated region: Identification, and tissue distribution of multiple new human exon 1. J. Mol. Endocrinol. 35, 283–292 (2005).CAS 
PubMed 
Article 

Google Scholar 
31.Bakusic, J., Schaufeli, W., Claes, S. & Godderis, L. Stress, burnout and depression: A systematic review on DNA methylation mechanisms. J. Psychosom. Res. 92, 34–44 (2017).PubMed 
Article 

Google Scholar 
32.Center for Whale Research. Population. https://www.whaleresearch.com. Accessed 11 Jan 2021 (2020).33.Fisheries and Oceans Canada. Recovery Strategy for the Northern and Southern Resident Killer Whales (Orcinus orca) in Canada [Proposed]. Species at Risk Act Recovery Strategy Series, Fisheries & Oceans Canada, Ottawa, x + 84 pp.(2018).34.DFO. Population Status Update for the Northern Resident Killer Whale (Orcinus orca) in 2018. DFO Can. Sci. Advis. Sec. Sci. Resp. 2019/025. (2019).35.Bigg, M. A., Olesiuk, P. F., Ellis, G. M., Ford, J. K. B. & Balcomb, K. C. Social organization and genealogy of resident killer whales (Orcinus orca) in the coastal waters of British Columbia and Washington State. Reports Int. Whal. Comm. 12, 383–405 (1990).
Google Scholar 
36.Ford, J. K. B. & Ellis, G. M. Selective foraging by fish-eating killer whales Orcinus orca in British Columbia. Mar. Ecol. Prog. Ser. 316, 185–199 (2006).ADS 
Article 

Google Scholar 
37.Chen, I.-H. et al. Selection of reference genes for RT-qPCR studies in blood of beluga whales (Delphinapterus leucas). PeerJ 4, e1810 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
38.Horvath, S. & Raj, K. DNA methylation-based biomarkers and the epigenetic clock theory of ageing. Nat. Rev. Genet. 19, 371–384 (2018).CAS 
PubMed 
Article 

Google Scholar 
39.Hoelzel, A. R., Dahlheim, M. E. & Stern, S. J. Low genetic variation among killer whales (Orcinus orca) in the eastern north Pacific and genetic differentiation between foraging specialists. J. Hered. 89, 121–128 (1998).CAS 
PubMed 
Article 

Google Scholar 
40.Yao, M., Stenzel-Poore, M. & Denver, R. J. Structural and functional conservation of vertebrate corticotropin- releasing factor genes: Evidence for a critical role for a conserved cyclic AMP response element. Endocrinology 148, 2518–2531 (2007).CAS 
PubMed 
Article 

Google Scholar 
41.Aguiniga, L. M., Yang, W., Yaggie, R. E., Schaeffer, A. J. & Klumpp, D. J. Acyloxyacyl hydrolase modulates depressive-like behaviors through aryl hydrocarbon receptor. Am. J. Physiol. Regul. Integr. Comp. Physiol. 317, R289–R300 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
42.Hankinson, O. The aryl hydrocarbon receptor complex. Annu. Rev. Pharmacol. Toxicol. 35, 307–340 (1995).CAS 
PubMed 
Article 

Google Scholar 
43.Lundin, J. I. et al. Pre-oil spill baseline profiling for contaminants in Southern Resident killer whale fecal samples indicates possible exposure to vessel exhaust. Mar. Pollut. Bull. 136, 448–453 (2018).CAS 
PubMed 
Article 

Google Scholar 
44.MacDonald, L. H. Evaluating and managing cumulative effects: Process and constraints. Environ. Manag. 26, 299–315 (2000).CAS 
Article 

Google Scholar 
45.National Academies of Sciences Engineering and Medicine. Approaches to Understanding the Cumulative Effects of Stressors on Marine Mammals. https://doi.org/10.17226/23479 (National Academies Press, 2017). 46.Barrett-Lennard, L. G., Smith, T. G. & Ellis, G. M. A cetacean biopsy system using lightweight pneumatic darts, and its effect on the behavior of killer whales. Mar. Mammal Sci. 12, 14–27 (1996).Article 

Google Scholar 
47.Sambrook, J., Fritsch, E. F. & Maniatis, H. Molecular Cloning: A Laboratory Manual (Cold Springs Harbor Laboratory Press, 1989).
Google Scholar 
48.Illumina. 16S Metagenomic Sequencing Library Preparation. Illumina.com 1–28 (2013).49.Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.j 17, 10–12 (2011).Article 

Google Scholar 
50.Bolger, A. M., Lohse, M. & Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
51.Magoč, T. & Salzberg, S. L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27, 2957–2963 (2011).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
52.Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
53.Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar  More

The location, diameter, height and circumference of the coral were measured (Table 1, Fig. 2). The Porites was brown to cream in colour and hemispherical in shape (Fig. 2). It was identified as either Porites lutea (Hump or Pore coral) or P. lobata (Lobe coral)14.The primary habitat on the Porites was live coral (70%), followed by sponge, live coral rock and a small amount of macroalgae (Table 2). No recently dead coral, coral rubble or sand was recorded (Table 2). We observed competition between the Porites and other species of coral and invertebrate including encrusting sponge, plating and branching Acropora spp., Montipora, Chlorodesmis, soft coral and zoanthids (Table 2, Figs. 3, 4).Table 2 Reef Health Impact Survey (RHIS) of habitat and species categories on Porites sp.Full size tableFigure 3Detail of the sub-habitats and competitive interactions Porites sp. and boring sponge Cliona viridis (left) and live coral Porites sp. and Montipora sp. (right) along interspecific contact zones.Full size imageFigure 4Detail of Reef Health Impact Survey (RHIS) of Porites.Full size imageThe boring sponge, Cliona viridis, is abundant on the Great Barrier Reef15. It is a common bioeroding species advancing laterally at around 1 cm and to a depth of 1.2 cm per annum15. Abundance of Cliona viridis is often correlated to substrate availability and water energy with the greatest abundance often on the windward side of bommies15. This correlates to our observations as the large proportion of the substrate estimated to cover the bommie (15%) was on the windward side. The sponge’s advances will likely continue to compromise the colony size and health.We recorded marine debris at the base of the Porites. The debris was 2–3 m of rope that appeared to have been wrapped around the base of an adjacent coral. Adjacent to the bommie were three concrete blocks.How big is the Porites coral at Goolboodi compared to other big corals in the GBR, and the world? Potts et al.6 reported a very large, rounded Porites colony, 6.9 m in diameter which is 3.1 m smaller than this study. Lough et al.16 reported coral cores from colonies between 1.6–8.0 m in height with the largest corals of 6.0 m at Havannah, North Molle and Masthead Islands, 7.5 m at Abraham Reef and 8.0 m at Sanctuary Reef. Recognising the limitations of published data, the Porites coral at Goolboodi is the largest diameter coral that has been measured, and the 6th tallest in the GBR. It is unknown if the other corals are still alive or dead.Other comparatively large massive Porites have previously been located throughout the Pacific. These have included multiple bommies measuring more than 10 m4 and one exceptionally large colony observed measuring 17 m × 12 m in American Samoa17. Additionally, large Porites sp. bommies have been observed at Green Island, 30 km east of Taiwan18 as well as an 11 m diameter Porites australiensis at Sesoko Island, Okinawa, Japan19.How old is this massive Porites? In discussions with the Australian Institute of Marine Science (AIMS), there is a robust, linear relationship ( > 80% variance explained) between Porites average linear extension rate and average annual sea surface temperature (SST)20,21 that provides an estimate of colony age from its height. Using average annual SST at 18.5S, 146.5E of 26.12C (from HadiSST data set), the estimated linear extension rate is determined by (2.97 × 26.12) − 65.46 = 1.21 cm/year. Given the colony height of 5.1–5.3 m, this gives an estimated age of 421–438 years. This is well before European exploration and settlement of Australia. AIMS has investigated over 328 colonies of massive Porites corals from 69 reefs along the GBR and has aged them as being from 10–436 years21. AIMS has not investigated this coral (pers. comm Neal Cantin). Based on limitations of published data, the Porites coral at Goolboodi is one of the oldest corals on the GBR.Why is the Porites partially dead on top and living on the side? The proportion of live coral tissue on a colony reflects the cumulative, integrated effect of both beneficial and adverse environmental factors. Substantial portions of coral tissue can die from exposure to sun at low tides or warm water without lethal consequences to the colony as a whole10. Partial mortality of large bommies provides available real estate for opportunistic, fast growing sessile organisms. In this instance, multiple species of tabulate and branching Acropora sp., encrusting Montipora sp. and encrusting sponges are among the benthic organisms to have colonised 30% of the coral bommies’ surface area. Intraspecific competition is also evident from the skeletal barriers created along contact zones22 (Fig. 3). There was no observation of disease or coral bleaching.The Porites is located in a relatively remote, rarely visited and highly protected Marine National Park (green) zone. Its location had not been previously reported and there is no existing database for significant corals in Australia or globally. Cataloguing the location of massive and long-lived corals can have multiple benefits. Scientific benefits include geochemical and isotopic analyses in coral skeletal cores which can help understand century-scale changes in oceanographic events and can be used to verify climate models. Social and economic benefits can include diving tourism, citizen science23 culture and stewardship. Perhaps the Significant Trees Register, which was designed by the National Trust24 to protect and celebrate Australia’s heritage could be considered as a model. There are risks of cataloguing the location of massive corals. It could be damaged by direct and indirect human uses including anchoring, research and pollution.Indigenous languages are an integral part of Indigenous culture, spirituality, and connection to country. We consulted Manbarra Traditional Owners about protocol and an appropriate cultural name for the Porites and they considered: Big (Muga), Home (Wanga), Coral reef (Muugar), Coral (Dhambi), Old (Anki, Gurgu), Old man (Gulula) and Old person (Gurgurbu)25. The recommendation by Manbarra Traditional Owners is that the Porites is named as Muga dhambi (Big coral). The feedback from the process of consultation was very positive with acknowledgement of the respect that the scientists have demonstrated to acknowledge Traditional Owners in this way.The large Porites coral at Goolboodi (Orpheus) Island is unusually rare and resilient. It has survived coral bleaching, invasive species, cyclones, severely low tides and human activities for almost 500 years. In an attempt to contextualise the resilience of these individual Porites we have reviewed major historic disturbances such as coral bleaching which has occurred since at least 1575 and potentially 99 bleaching events in the GBR over the past 400 plus years26. Other indicators such as high-density ‘stress bands’ were recorded from 1877 and are significantly more frequent in the late twentieth and early twenty-first centuries in accordance with rising temperatures from anthropogenic global warming27. In addition there have been an average of 1–2 tropical cyclones per decade (40–80 in total) that have potentially impacted the coral adjacent to Goolboodi Island28,29; 46 tropical cyclones impacted the area between Ingham and Townsville from 1858 to 200830. The cumulative impact of almost 100 bleaching events and up to 80 major cyclones over a period of four centuries, plus declining nearshore water quality contextualise the high resilience of this Porites coral. Looking to the future there is real concern for corals in the GBR due to many impacts including climate change, declining water quality, overfishing and coastal development31,32. This field note provides important geospatial, environmental, and cultural information of a rare coral that can be monitored, appreciated, potentially restored and hopefully inspire future generations to care more for our reefs and culture. More