Ecology

Sylvan, J. How to protect a coral reef: the public trust doctrine and the law of the sea recommended citation. Sustain. Dev. Law Policy 7, 12 (2006).
Google Scholar 
LaJeunesse, T. C. et al. Systematic revision of symbiodiniaceae highlights the antiquity and diversity of coral endosymbionts. Curr. Biol. 28, 2570–2580.e6 (2018).CAS 
PubMed 
Article 

Google Scholar 
Kopp, C. et al. Highly dynamic cellular-level response of symbiotic coral to a sudden increase in environmental nitrogen. mBio 4, e00052–13 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Muscatine, L. The role of symbiotic algae in carbon and energy flux in reef corals. Coral Reef. 25, 75–87 (1990).
Google Scholar 
Dubinsky, Z. & Stambler, N. Coral reefs: an ecosystem in transition. (Springer, 2011).Wiedenmann, J. et al. Nutrient enrichment can increase the susceptibility of reef corals to bleaching. https://doi.org/10.1038/NCLIMATE1661 (2012).Suggett, D. J., Warner, M. E. & Leggat, W. Symbiotic dinoflagellate functional diversity mediates coral survival under ecological crisis. Trends Ecol. Evolution 32, 735–745 (2017).Article 

Google Scholar 
Morris, L. A., Voolstra, C. R., Quigley, K. M., Bourne, D. G. & Bay, L. K. Nutrient availability and metabolism affect the stability of coral–symbiodiniaceae symbioses. Trends Microbiol. 27, 678–689 (2019).CAS 
PubMed 
Article 

Google Scholar 
Lehnert, E. M. et al. Extensive differences in gene expression between symbiotic and aposymbiotic cnidarians. G3 (Bethesda) 4, 277–95 (2014).CAS 
Article 

Google Scholar 
Dubinsky, Z. & Berman-Frank, I. Uncoupling primary production from population growth in photosynthesizing organisms in aquatic ecosystems. in. Aquat. Sci. 63, 4–17 (2001).CAS 
Article 

Google Scholar 
Burriesci, M. S., Raab, T. K. & Pringle, J. R. Evidence that glucose is the major transferred metabolite in dinoflagellate–cnidarian symbiosis. J. Exp. Biol. 215, 3467–3477 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Davy, S. K., Allemand, D. & Weis, V. M. Cell biology of cnidarian-dinoflagellate symbiosis. Microbiol. Mol. Biol. Rev. 76, 229–61 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rädecker, N., Pogoreutz, C., Voolstra, C. R., Wiedenmann, J. & Wild, C. Nitrogen cycling in corals: the key to understanding holobiont functioning? Trends Microbiol. 23, 490–497 (2015).PubMed 
Article 
CAS 

Google Scholar 
Cui, G. et al. Host-dependent nitrogen recycling as a mechanism of symbiont control in Aiptasia. PLOS Genet. 15, e1008189 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Rädecker, N. et al. Heat stress destabilizes symbiotic nutrient cycling in corals. Proc. Natl Acad. Sci. USA 118, https://doi.org/10.1073/pnas.2022653118 (2021).Weis, V. M. Cellular mechanisms of Cnidarian bleaching: stress causes the collapse of symbiosis. J. Exp. Biol. 211, 3059–3066 (2008).CAS 
PubMed 
Article 

Google Scholar 
Wooldridge, S. A. Breakdown of the coral-algae symbiosis: towards formalising a linkage between warm-water bleaching thresholds and the growth rate of the intracellular zooxanthellae. Biogeosciences Discuss. 9, 8111–8139 (2012).
Google Scholar 
Cziesielski, M. J., Schmidt‐Roach, S. & Aranda, M. The past, present, and future of coral heat stress studies. Ecol. Evol. https://doi.org/10.1002/ece3.5576 (2019).Leggat, W. et al. Differential responses of the coral host and their algal symbiont to thermal stress. PLoS ONE 6, e26687 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pinzón, J. H. et al. Whole transcriptome analysis reveals changes in expression of immune-related genes during and after bleaching in a reef-building coral. R. Soc. Open Sci. 2, 140214 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ziegler, M., Seneca, F. O., Yum, L. K., Palumbi, S. R. & Voolstra, C. R. Bacterial community dynamics are linked to patterns of coral heat tolerance. Nat. Commun. 8, 14213 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Bang, C. et al. Metaorganisms in extreme environments: do microbes play a role in organismal adaptation? Zoology 127, 1–19 (2018).PubMed 
Article 

Google Scholar 
Berkelmans, R. & van Oppen, M. J. H. The role of zooxanthellae in the thermal tolerance of corals: a “nugget of hope” for coral reefs in an era of climate change. Proc. Biol. Sci./R. Soc. 273, 2305–12 (2006).
Google Scholar 
Sampayo, E. M., Ridgway, T., Bongaerts, P. & Hoegh-Guldberg, O. Bleaching susceptibility and mortality of corals are determined by fine-scale differences in symbiont type. Proc. Natl Acad. Sci. 105, 10444–10449 (2008).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Howells, E. J. et al. Coral thermal tolerance shaped by local adaptation of photosymbionts. Nat. Clim. Change https://doi.org/10.1038/nclimate1330 (2011).Cziesielski, M. J. et al. Multi-omics analysis of thermal stress response in a zooxanthellate cnidarian reveals the importance of associating with thermotolerant symbionts. Proc. Biol. Sci. 285, 20172654 (2018).PubMed 
PubMed Central 

Google Scholar 
Baker, A. C., Starger, C. J., McClanahan, T. R. & Glynn, P. W. Corals’ adaptive response to climate change. Nature 430, 741–741 (2004).CAS 
PubMed 
Article 

Google Scholar 
Thornhill, D. J., LaJeunesse, T. C., Kemp, D. W., Fitt, W. K. & Schmidt, G. W. Multi-year, seasonal genotypic surveys of coral-algal symbioses reveal prevalent stability or post-bleaching reversion. Mar. Biol. 148, 711–722 (2006).Article 

Google Scholar 
Palumbi, S. R., Barshis, D. J., Traylor-Knowles, N. & Bay, R. A. Mechanisms of reef coral resistance to environmental stress,making its relative ability to acclimate or adapt extremely important to the to future climate change. Science 344, 895–898 (2014).CAS 
PubMed 
Article 

Google Scholar 
Herrera, M. et al. Temperature transcends partner specificity in the symbiosis establishment of a cnidarian. ISME J. 15, 141–153 (2021).CAS 
PubMed 
Article 

Google Scholar 
Howells, E. J. et al. Corals in the hottest reefs in the world exhibit symbiont fidelity not flexibility. Mol. Ecol. 29, 899–911 (2020).CAS 
PubMed 
Article 

Google Scholar 
Hume, B. C. C., Mejia-Restrepo, A., Voolstra, C. R. & Berumen, M. L. Fine-scale delineation of Symbiodiniaceae genotypes on a previously bleached central Red Sea reef system demonstrates a prevalence of coral host-specific associations. Coral Reefs 1–19 https://doi.org/10.1007/s00338-020-01917-7 (2020).Perez, S. F., Cook, C. B. & Brooks, W. R. The role of symbiotic dinoflagellates in the temperature-induced bleaching response of the subtropical sea anemone Aiptasia pallida. J. Exp. Mar. Biol. Ecol. 256, 1–14 (2001).CAS 
PubMed 
Article 

Google Scholar 
Mieog, J. C. et al. The roles and interactions of symbiont, host and environment in defining coral fitness. PLoS ONE 4, e6364 (2009).Cantin, N. E., van Oppen, M. J. H., Willis, B. L., Mieog, J. C. & Negri, A. P. Juvenile corals can acquire more carbon from high-performance algal symbionts. Coral Reefs 28, 405–414 (2009).Article 

Google Scholar 
Herrera, M. et al. Unfamiliar partnerships limit cnidarian holobiont acclimation to warming. Glob. Change Biol. 26, 5539–5553 (2020).Article 

Google Scholar 
LaJeunesse, T. et al. Closely related Symbiodinium spp. differ in relative dominance in coral reef host communities across environmental, latitudinal and biogeographic gradients. Mar. Ecol. Prog. Ser. 284, 147–161 (2004).Article 

Google Scholar 
Parkinson, J. E. & Baums, I. B. The extended phenotypes of marine symbioses: ecological and evolutionary consequences of intraspecific genetic diversity in coral-algal associations. Front. Microbiol. 5, 445 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
Coffroth, M. A., Poland, D. M., Petrou, E. L., Brazeau, D. A. & Holmberg, J. C. Environmental symbiont acquisition may not be the solution to warming seas for reef-building corals. PLoS ONE 5, e13258 (2010).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Bellantuono, A. J., Granados-Cifuentes, C., Miller, D. J., Hoegh-Guldberg, O. & Rodriguez-Lanetty, M. Coral thermal tolerance: tuning gene expression to resist thermal stress. PLoS ONE 7, e50685 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sunagawa, S. et al. Generation and analysis of transcriptomic resources for a model system on the rise: the sea anemone Aiptasia pallida and its dinoflagellate endosymbiont. BMC Genomics 10, 258 (2009).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Baumgarten, S. et al. The genome of Aiptasia, a sea anemone model for coral symbiosis. Proc. Natl Acad. Sci. 112, 201513318 (2015).
Google Scholar 
Matthews, J. L. et al. Menthol-induced bleaching rapidly and effectively provides experimental aposymbiotic sea anemones (Aiptasia sp.) for symbiosis investigations. J. Exp. Biol. jeb.128934 https://doi.org/10.1242/JEB.128934 (2015).Kenkel, C. D. et al. Evidence for a host role in thermotolerance divergence between populations of the mustard hill coral (Porites astreoides) from different reef environments. Mol. Ecol. 22, 4335–4348 (2013).CAS 
PubMed 
Article 

Google Scholar 
Polato, N. R., Altman, N. S. & Baums, I. B. Variation in the transcriptional response of threatened coral larvae to elevated temperatures. Mol. Ecol. 22, 1366–1382 (2013).CAS 
PubMed 
Article 

Google Scholar 
DeSalvo, M., Sunagawa, S., Voolstra, C. R. & Medina, M. Transcriptomic resonses to heat stress and bleaching in the elkhorn coral Acropora palmata. Mar. Ecol. Prog. Ser. 402, 97–113 (2010).CAS 
Article 

Google Scholar 
Maor-Landaw, K. & Levy, O. Gene expression profiles during short-term heat stress; branching vs. massive Scleractinian corals of the Red Sea. PeerJ 4, e1814 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Yamamoto, K. et al. Control of the heat stress-induced alternative splicing of a subset of genes by hnRNP K. Genes Cells 21, 1006–1014 (2016).CAS 
PubMed 
Article 

Google Scholar 
Seneca, F. O. & Palumbi, S. R. The role of transcriptome resilience in resistance of corals to bleaching. Mol. Ecol. 24, 1467–1484 (2015).PubMed 
Article 

Google Scholar 
Meyer, E. & Weis, V. M. Study of cnidarian-algal symbiosis in the “omics” age. Biol. Bull. 223, 44–65 (2012).CAS 
PubMed 
Article 

Google Scholar 
Oakley, C. A. et al. Thermal shock induces host proteostasis disruption and endoplasmic reticulum stress in the model symbiotic Cnidarian Aiptasia. J. Proteome Res. 16, 2121–2134 (2017).CAS 
PubMed 
Article 

Google Scholar 
Robbart, M. L., Peckol, P., Scordilis, S. P., Curran, H. A. & Brown-Saracino, J. Population recovery and differential heat shock protein expression for the corals Agaricia agaricites and A-tenuifolia in Belize. Mar. Ecol. Prog. Ser. 283, 151–160 (2004).Article 

Google Scholar 
Barshis, D. J. et al. Genomic basis for coral resilience to climate change. Proc. Natl Acad. Sci. 110, 1387–1392 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Traylor-Knowles, N., Rose, N. H. & Palumbi, S. R. The cell specificity of gene expression in the response to heat stress in corals. J. Exp. Biol. 220, 1837–1845 (2017).PubMed 

Google Scholar 
Benchimol, S. p53-dependent pathways of apoptosis. Cell Death Differ. 8, 1049–1051 (2001).CAS 
PubMed 
Article 

Google Scholar 
Moya, A. et al. Functional conservation of the apoptotic machinery from coral to man: The diverse and complex Bcl-2 and caspase repertoires of Acropora millepora. BMC Genomics 17, 62 (2016).Elmore, S. Apoptosis: a review of programmed cell death. Toxicol. Pathol. 35, 495–516 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Karim, W., Nakaema, S. & Hidaka, M. Temperature effects on the growth rates and photosynthetic activities of symbiodinium cells. J. Mar. Sci. Eng. 3, 368–381 (2015).Article 

Google Scholar 
Cunning, R. & Baker, A. C. Excess algal symbionts increase the susceptibility of reef corals to bleaching. Nat. Clim. Change 3, 259–262 (2013).Article 

Google Scholar 
Rehman, A. U. et al. Symbiodinium sp. cells produce light-induced intra- and extracellular singlet oxygen, which mediates photodamage of the photosynthetic apparatus and has the potential to interact with the animal host in coral symbiosis. N. Phytologist 212, 472–484 (2016).CAS 
Article 

Google Scholar 
Lesser, K. B. & Garcia, F. A. Association between polycystic ovary syndrome and glucose intolerance during pregnancy. J. Matern. Fetal Med. 6, 303–307 (1997).CAS 
PubMed 
Article 

Google Scholar 
Dunn, S. R., Schnitzler, C. E. & Weis, V. M. Apoptosis and autophagy as mechanisms of dinoflagellate symbiont release during cnidarian bleaching: every which way you lose. Proc. R. Soc. Lond. B: Biol. Sci. 274, 3079–3085 (2007).
Google Scholar 
DeSalvo, M. K. et al. Coral host transcriptomic states are correlated with Symbiodinium genotypes. Mol. Ecol. 19, 1174–1186 (2010).CAS 
PubMed 
Article 

Google Scholar 
Levin, R. A. et al. Engineering strategies to decode and enhance the genomes of coral symbionts. Front. Microbiol. https://doi.org/10.3389/fmicb.2017.01220 (2017).Yuyama, I., Ishikawa, M., Nozawa, M., Yoshida, M. & Ikeo, K. Transcriptomic changes with increasing algal symbiont reveal the detailed process underlying establishment of coral-algal symbiosis. Sci. Rep. 8, 16802 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Sproles, A. E. et al. Sub-cellular imaging shows reduced photosynthetic carbon and increased nitrogen assimilation by the non-native endosymbiont Durusdinium trenchii in the model cnidarian Aiptasia. Environ. Microbiol. 22, 3741–3753 (2020).CAS 
PubMed 
Article 

Google Scholar 
Rädecker, N. et al. Using Aiptasia as a model to study metabolic interactions in Cnidarian-Symbiodinium symbioses. Front. Physiol. 9, 214 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Falkowski, P. G., Dubinsky, Z., Muscatine, L. & McCloskey, L. Population control in symbiotic corals. BioScience 43, 606–611 (1993).Article 

Google Scholar 
Wang & Douglas. Nitrogen recycling or nitrogen conservation in an alga-invertebrate symbiosis? J. Exp. Biol. 201, 2445–53 (1998).Loram, J. E., Trapido-Rosenthal, H. G. & Douglas, A. E. Functional significance of genetically different symbiotic algae Symbiodinium in a coral reef symbiosis. Mol. Ecol. 16, 4849–4857 (2007).CAS 
PubMed 
Article 

Google Scholar 
Karako-Lampert, S. et al. Transcriptome analysis of the scleractinian coral Stylophora pistillata. PLoS One 9, e88615 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Hillyer, K. E., Tumanov, S., Villas-Bôas, S. & Davy, S. K. Metabolite profiling of symbiont and host during thermal stress and bleaching in a model cnidarian-dinoflagellate symbiosis. J. Exp. Biol. 219, 516–27 (2016).PubMed 

Google Scholar 
Bertucci, A., Forêt, S., Ball, E. E. & Miller, D. J. Transcriptomic differences between day and night in Acropora millepora provide new insights into metabolite exchange and light-enhanced calcification in corals. Mol. Ecol. 24, 4489–4504 (2015).CAS 
PubMed 
Article 

Google Scholar 
Matthews, J. L. et al. Optimal nutrient exchange and immune responses operate in partner specificity in the cnidarian-dinoflagellate symbiosis. Proc. Natl Acad. Sci. 114, 13194–13199 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lin, M.-F., Takahashi, S., Forêt, S., Davy, S. K. & Miller, D. J. Transcriptomic analyses highlight the likely metabolic consequences of colonization of a cnidarian host by native or non-native Symbiodinium species. Biol. Open 8, bio038281 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Medrano, E., Merselis, D. G., Bellantuono, A. J. & Rodriguez-Lanetty, M. Proteomic Basis of Symbiosis: A Heterologous Partner Fails to Duplicate Homologous Colonization in a Novel Cnidarian– Symbiodiniaceae Mutualism. Front. Microbiol. 10, 1153 (2019).Schoepf, V., Stat, M., Falter, J. L. & McCulloch, M. T. Limits to the thermal tolerance of corals adapted to a highly fluctuating, naturally extreme temperature environment. Sci. Rep. 5, 17639 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Xiang, T., Hambleton, E. A., DeNofrio, J. C., Pringle, J. R. & Grossman, A. R. Isolation of clonal axenic strains of the symbiotic dinoflagellate Symbiodinium and their growth and host specificity1. J. Phycol. 49, 447–458 (2013).CAS 
PubMed 
Article 

Google Scholar 
Bray, N. L., Pimentel, H., Melsted, P. & Pachter, L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol. 34, 525–527 (2016).CAS 
PubMed 
Article 

Google Scholar 
Pimentel, H., Bray, N. L., Puente, S., Melsted, P. & Pachter, L. Differential analysis of RNA-seq incorporating quantification uncertainty. Nat. Methods 14, 687–690 (2017).CAS 
PubMed 
Article 

Google Scholar  More

Shi, M. et al. Redefining the invertebrate RNA virosphere. Nature 540, 539–543 (2016).CAS 
PubMed 
Article 

Google Scholar 
Zhang, Y.-Z., Shi, M. & Holmes, E. C. Using metagenomics to characterize an expanding virosphere. Cell 172, 1168–1172 (2018).CAS 
PubMed 
Article 

Google Scholar 
Li, C.-X. et al. Unprecedented genomic diversity of RNA viruses in arthropods reveals the ancestry of negative-sense RNA viruses. eLife 4, e05378 (2015).PubMed Central 
Article 
CAS 

Google Scholar 
Starr, E. P., Nuccio, E. E., Pett-Ridge, J., Banfield, J. F. & Firestone, M. K. Metatranscriptomic reconstruction reveals RNA viruses with the potential to shape carbon cycling in soil. Proc. Natl Acad. Sci. USA 116, 25900–25908 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wolf, Y. I. et al. Doubling of the known set of RNA viruses by metagenomic analysis of an aquatic virome. Nat. Microbiol. 5, 1262–1270 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zayed, A. A. et al. Cryptic and abundant marine viruses at the evolutionary origins of Earth’s RNA virome. Science 376, 156–162 (2022).CAS 
PubMed 
Article 

Google Scholar 
Simmonds, P. et al. Virus taxonomy in the age of metagenomics. Nat. Rev. Microbiol. 15, 161–168 (2017).CAS 
PubMed 
Article 

Google Scholar 
Trubl, G., Hyman, P., Roux, S. & Abedon, S. T. Coming-of-age characterization of soil viruses: a user’s guide to virus isolation, detection within metagenomes, and viromics. Soil Syst. 4, 23 (2020).CAS 
Article 

Google Scholar 
Jin, M. et al. Diversities and potential biogeochemical impacts of mangrove soil viruses. Microbiome 7, 58 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
Trubl, G. et al. Soil viruses are underexplored players in ecosystem carbon processing. mSystems 3, e00076-18 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Steward, G. F. et al. Are we missing half of the viruses in the ocean? ISME J. 7, 672–679 (2013).CAS 
PubMed 
Article 

Google Scholar 
Paul, K. I., Scott Black, A. & Conyers, M. K. in Advances in Agronomy. Sparks, D.L., Vol. 78 187–214 (Elsevier, 2003).Urayama, S., Takaki, Y. & Nunoura, T. FLDS: a comprehensive dsRNA sequencing method for intracellular RNA virus surveillance. Microbes Environ. 31, 33–40 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Armbrust, E. V. The life of diatoms in the world’s oceans. Nature 459, 185–192 (2009).CAS 
PubMed 
Article 

Google Scholar 
Wu, W., Jin, Y., Bai, F. & Jin, S. in Molecular Medical Microbiology. Tang, Y.W., Liu, D., Schwartzman, J., Sussman, M., Poxton, I., 753–767 (Elsevier, 2015).Cooney, S., O’Brien, S., Iversen, C. & Fanning, S. in Encyclopedia of Food Safety. Motarjemi, Y., 433–441 (Elsevier, 2014).Geoghegan, J. L. et al. Hidden diversity and evolution of viruses in market fish. Virus Evol. 4, vey031 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Lauber, C. et al. Deciphering the origin and evolution of hepatitis B viruses by means of a family of non-enveloped fish viruses. Cell Host Microbe 22, 387–399.e6 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Shi, M., Zhang, Y.-Z. & Holmes, E. C. Meta-transcriptomics and the evolutionary biology of RNA viruses. Virus Res. 243, 83–90 (2018).CAS 
PubMed 
Article 

Google Scholar 
Turnbull, O. M. H. et al. Meta-transcriptomic identification of divergent Amnoonviridae in Fish. Viruses 12, 1254 (2020).CAS 
PubMed Central 
Article 

Google Scholar 
Bauermann, F. V., Hause, B., Buysse, A. R., Joshi, L. R. & Diel, D. G. Identification and genetic characterization of a porcine hepe-astrovirus (bastrovirus) in the United States. Arch. Virol. 164, 2321–2326 (2019).CAS 
PubMed 
Article 

Google Scholar 
Oude Munnink, B. B. et al. A novel astrovirus-like RNA virus detected in human stool. Virus Evol. 2, vew005 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Williamson, K. E. et al. Estimates of viral abundance in soils are strongly influenced by extraction and enumeration methods. Biol. Fertil. Soils 49, 857–869 (2013).Article 

Google Scholar 
Wang, C., Liu, D. & Bai, E. Decreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition. Soil Biol. Biochem. 120, 126–133 (2018).CAS 
Article 

Google Scholar 
Wang, Q. et al. Effects of nitrogen and phosphorus inputs on soil bacterial abundance, diversity, and community composition in Chinese fir plantations. Front. Microbiol. 9, 1543 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Payne, S. in Viruses. Payne, S., 219–226 (Elsevier, 2017).Hillman, B. I. & Cai, G. The family Narnaviridae. Adv. Virus Res. 86, 149–176 (2013).PubMed 
Article 

Google Scholar 
Wolf, Y. I. et al. Origins and evolution of the global RNA virome. mBio 9, e02329-18 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
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 
Li, D., Liu, C.-M., Luo, R., Sadakane, K. & Lam, T.-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31, 1674–1676 (2015).CAS 
PubMed 
Article 

Google Scholar 
Wu, F. et al. A new coronavirus associated with human respiratory disease in China. Nature 579, 265–269 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Buchfink, B., Xie, C. & Huson, D. H. Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60 (2015).CAS 
PubMed 
Article 

Google Scholar 
Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Capella-Gutierrez, S., Silla-Martinez, J. M. & Gabaldon, T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25, 1972–1973 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nguyen, L.-T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).CAS 
PubMed 
Article 

Google Scholar 
Paradis, E. & Schliep, K. ape 5.0: an environment for modern phylogenetics and evolutionary analysis in R. Bioinformatics 35, 526–528 (2019).CAS 
PubMed 
Article 

Google Scholar 
Yu, G., Smith, D. K., Zhu, H., Guan, Y. & Lam, T. T. ggtree: an R package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol. 8, 28–36 (2017).Article 

Google Scholar 
Langmead, B. & Salzberg, S. L. Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nayfach, S. et al. CheckV assesses the quality and completeness of metagenome-assembled viral genomes. Nat. Biotechnol. 39, 578–585 (2021).CAS 
PubMed 
Article 

Google Scholar 
Almagro Armenteros, J. J. et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat. Biotechnol. 37, 420–423 (2019).CAS 
PubMed 
Article 

Google Scholar 
Krogh, A., Larsson, B., von Heijne, G. & Sonnhammer, E. L. L. Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. J. Mol. Biol. 305, 567–580 (2001).CAS 
PubMed 
Article 

Google Scholar 
Gupta, R., Jung, E. & Brunak, S. NetNGlyc 1.0 Server (2017). DTU Health Tech. http://www.cbs.dtu.dk/services/NetNGlyc/Mirdita, M. et al. Uniclust databases of clustered and deeply annotated protein sequences and alignments. Nucleic Acids Res. 45, D170–D176 (2017).CAS 
PubMed 
Article 

Google Scholar 
Remmert, M., Biegert, A., Hauser, A. & Söding, J. HHblits: lightning-fast iterative protein sequence searching by HMM-HMM alignment. Nat. Methods 9, 173–175 (2012).CAS 
Article 

Google Scholar 
Quast, C. et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596 (2012).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Li, H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics 27, 2987–2993 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lagkouvardos, I., Fischer, S., Kumar, N. & Clavel, T. Rhea: a transparent and modular R pipeline for microbial profiling based on 16S rRNA gene amplicons. PeerJ 5, e2836 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
McLeod, A., Xu, C. & Lai, Y. Package ‘bestglm’. CRAN. (2020).Fu, L., Niu, B., Zhu, Z., Wu, S. & Li, W. CD-HIT: accelerated for clustering the next-generation sequencing data. Bioinformatics 28, 3150–3152 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

Seebens, H. et al. Projecting the continental accumulation of alien species through to 2050. Glob. Change Biol. 27(5), 970–982 (2021).ADS 

Google Scholar 
Bellard, C., Cassey, P. & Blackburn, T. M. Alien species as a driver of recent extinctions. Biol. Lett. 12(2), 20150623 (2016).PubMed 
PubMed Central 

Google Scholar 
Walsh, J. R., Carpenter, S. R. & Vander Zanden, M. J. Invasive species triggers a massive loss of ecosystem services through a trophic cascade. Proc. Natl. Acad. Sci. 113(15), 4081–4085 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Linders, T. E. W. et al. Direct and indirect effects of invasive species: Biodiversity loss is a major mechanism by which an invasive tree affects ecosystem functioning. J. Ecol. 107(6), 2660–2672 (2019).
Google Scholar 
Diagne, C., Ballesteros-Mejia, L., Bodey, T., Cuthbert, R., Fantle-Lepczyk, J., Angulo, E., Dobigny, G., & Courchamp, F. Economic costs of invasive rodents worldwide: The tip of the iceberg (2021).Schaffner, F., Medlock, J. M. & Van Bortel, A. W. Public health significance of invasive mosquitoes in Europe. Clin. Microbiol. Infect. 19(8), 685–692 (2013).CAS 
PubMed 

Google Scholar 
Schaffner, U. et al. Biological weed control to relieve millions from Ambrosia allergies in Europe. Nat. Commun. 11(1), 1–7 (2020).
Google Scholar 
Shackleton, R. T., Shackleton, C. M. & Kull, C. A. The role of invasive alien species in shaping local livelihoods and human well-being: A review. J. Environ. Manag. 229, 145–157 (2019).
Google Scholar 
Clavero, M. & García-Berthou, E. Invasive species are a leading cause of animal extinctions. Trends Ecol. Evol. 20(3), 110 (2005).PubMed 

Google Scholar 
Crystal-Ornelas, R. & Lockwood, J. L. The ‘known unknowns’ of invasive species impact measurement. Biol. Invasions 22(4), 1513–1525 (2020).
Google Scholar 
Florencio, M., Lobo, J. M. & Bini, L. M. Biases in global effects of exotic species on local invertebrates: A systematic review. Biol. Invasions 21(10), 3043–3061 (2019).
Google Scholar 
Measey, J. et al. Why have a pet amphibian? Insights from YouTube. Front. Ecol. Evol. 7, 52 (2019).
Google Scholar 
Ossiboff, R. J. et al. Differentiating Batrachochytrium dendrobatidis and B. salamandrivorans in amphibian chytridiomycosis using RNAScope in situ hybridization. Front. Vet. Sci. 6, 304 (2019).PubMed 
PubMed Central 

Google Scholar 
Kraus, F. Alien Reptiles and Amphibians: A Scientific Compendium and Analysis, vol. 4. (Springer Science & Business Media, 2009).Kraus, F. Impacts from invasive reptiles and amphibians. Annu. Rev. Ecol. Evol. Syst. 46, 75–97 (2015).
Google Scholar 
Ramsay, N. F., Ng, P. K. A., O’Riordan, R. M., & Chou, L. M. The red-eared slider (Trachemys scripta elegans) in Asia: A review. Biological invaders in inland waters: Profiles, distribution, and threats 161–174 (2007).Lindsay, M. K., Zhang, Y., Forstner, M. R. & Hahn, D. Effects of the freshwater turtle Trachemys scripta elegans on ecosystem functioning: An approach in experimental ponds. Amphibia-Reptilia 34(1), 75–84 (2013).
Google Scholar 
Phillips, B. L. & Shine, R. An invasive species induces rapid adaptive change in a native predator: Cane toads and black snakes in Australia. Proc. R. Soc. B Biol. Sci. 273(1593), 1545–1550 (2006).
Google Scholar 
Shanmuganathan, T. et al. Biological control of the cane toad in Australia: A review. Anim. Conserv. 13, 16–23 (2010).
Google Scholar 
Smart, A. S., Tingley, R. & Phillips, B. L. Estimating the benefit of quarantine: Eradicating invasive cane toads from islands. NeoBiota 60, 117 (2020).
Google Scholar 
Reaser, J. K. et al. Ecological and socioeconomic impacts of invasive alien species in island ecosystems. Environ. Conserv. 34, 98–111 (2007).
Google Scholar 
Fritts, T. H. Economic costs of electrical system instability and power outages caused by snakes on the island of Guam. Int. Biodeterior. Biodegrad. 49(2–3), 93–100 (2002).
Google Scholar 
Rodda, G. H., Fritts, T. H. & Chiszar, D. The disappearance of Guam’s wildlife. Bioscience 47(9), 565–574 (1997).
Google Scholar 
Kraus, F. Reptiles and amphibians. In Encyclopedia of Biological Invasions 590–594. (University of California Press, 2011).Kraus, F. Global trends in alien reptiles and amphibians. Aliens Invasive Species Bull. 28, 13–18 (2009).
Google Scholar 
Capinha, C., Marcolin, F. & Reino, L. Human-induced globalization of insular herpetofaunas. Glob. Ecol. Biogeogr. 29(8), 1328–1349 (2020).
Google Scholar 
Reed, R. N. & Kraus, F. Invasive reptiles and amphibians: Global perspectives and local solutions. Anim. Conserv. 13, 3–4 (2010).
Google Scholar 
Wasserman, R. J., Dick, J. T., Welch, R. J., Dalu, T. & Magellan, K. Site and species selection for religious release of non-native fauna. Conserv. Biol. 33(4), 969–971 (2019).PubMed 

Google Scholar 
Li, X., Liu, X., Kraus, F., Tingley, R. & Li, Y. Risk of biological invasions is concentrated in biodiversity hotspots. Front. Ecol. Environ. 14(8), 411–417 (2016).
Google Scholar 
Bellard, C. & Jeschke, J. M. A spatial mismatch between invader impacts and research publications. Conserv. Biol. 30(1), 230–232 (2016).CAS 
PubMed 

Google Scholar 
Diagne, C. et al. InvaCost, a public database of the economic costs of biological invasions worldwide. Sci. Data 7(1), 1–12 (2020).
Google Scholar 
Diagne, C. et al. High and rising economic costs of biological invasions worldwide. Nature 592(7855), 571–576 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Cuthbert, R. N., Diagne, C., Haubrock, P. J., Turbelin, A. J., & Courchamp, F. Are the “100 of the world’s worst” invasive species also the costliest? Biol. Invasions 1–10 (2021).Cuthbert, R. N. et al. Global economic costs of aquatic invasive alien species. Sci. Total Environ. 775, 145238 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Haubrock, P. J. et al. Biological invasions in Singapore and Southeast Asia: Data gaps fail to mask potentially massive economic costs. NeoBiota 67, 131–152 (2021).
Google Scholar 
Van Wilgen, N. J., Gillespie, M. S., Richardson, D. M. & Measey, J. A taxonomically and geographically constrained information base limits non-native reptile and amphibian risk assessment: A systematic review. PeerJ 6, e5850 (2018).PubMed 
PubMed Central 

Google Scholar 
Cuthbert, R. et al. Economic costs of biological invasions in the United Kingdom. Neobiota 67, 299–328 (2021).
Google Scholar 
Heringer, G. et al. The economic costs of biological invasions in Central and South America: A first regional assessment. NeoBiota 67, 401 (2021).
Google Scholar 
Paini, D. R. et al. Global threat to agriculture from invasive species. Proc. Natl. Acad. Sci. 113(27), 7575–7579 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Capinha, C. et al. Diversity, biogeography and the global flows of alien amphibians and reptiles. Divers. Distrib. 23(11), 1313–1322 (2017).
Google Scholar 
Kumschick, S. et al. How repeatable is the Environmental Impact Classification of Alien Taxa (EICAT)? Comparing independent global impact assessments of amphibians. Ecol. Evol. 7(8), 2661–2670 (2017).PubMed 
PubMed Central 

Google Scholar 
Diagne, C., Catford, J. A., Essl, F., Nuñez, M. A. & Courchamp, F. What are the economic costs of biological invasions? A complex topic requiring international and interdisciplinary expertise. NeoBiota 63, 25 (2020).
Google Scholar 
Diagne, C. et al. The economic costs of biological invasions in Africa: A growing but neglected threat?. NeoBiota 67, 11–51 (2021).
Google Scholar 
Bradshaw, C. J. et al. Detailed assessment of the reported economic costs of invasive species in Australia. NeoBiota 67, 511–550 (2021).
Google Scholar 
Dorcas, M. E. et al. Severe mammal declines coincide with proliferation of invasive Burmese pythons in Everglades National Park. Proc. Natl. Acad. Sci. 109(7), 2418–2422 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Mazzotti, F. J. et al. Implications of the 2013 Python Challenge® for ecology and management of Python molorus bivittatus (Burmese Python) in Florida. Southeast. Nat. 15(sp8), 63–74 (2016).
Google Scholar 
Smith, B. J. et al. Betrayal: Radio-tagged Burmese pythons reveal locations of conspecifics in Everglades National Park. Biol. Invasions 18(11), 3239–3250 (2016).
Google Scholar 
Easteal, S. The history of introductions of Bufo marinus (Amphibia: Anura); A natural experiment in evolution. Biol. J. Lin. Soc. 16(2), 93–113 (1981).
Google Scholar 
Haubrock, P. J., Bernery, C., Cuthbert, R. N., Liu, C., Kourantidou, M., Leroy, B., Turbelin, A., Kramer, A. M., Verbrugge, L., Diagne, C., Courchamp, F., & Gozlan, R. E. What is the recorded economic cost of alien invasive fishes worldwide? (2021).Angulo, E., Hoffmann, B., Ballesteros-Mejia, L., Taheri, A., Balzani, P., Renault, D., Cordonnier, M., Bellard, C., Diagne, C., Ahmed, D. A., Watari, Y., & Courchamp, F. Economic costs of invasive alien ants worldwide. (2021).Kouba, A., Oficialdegui, F., Cuthbert, R., Kourantidou, M., Tricarico, E., Leroy, B., Gozlan, R., Courchamp, F., & Haubrock, P. Feeling the pinch: Global economic costs of crayfish invasions and comparison with other aquatic crustaceans (2021).Dufresnes, C. et al. Cryptic invasion of Italian pool frogs (Pelophylax bergeri) across Western Europe unraveled by multilocus phylogeography. Biol. Invasions 19(5), 1407–1420 (2017).
Google Scholar 
Kumschick, S. et al. Impact assessment with different scoring tools: How well do alien amphibian assessments match?. NeoBiota 33, 53 (2017).
Google Scholar 
Crystal-Ornelas, R. et al. Economic costs of biological invasions within North America. NeoBiota 67, 485 (2021).
Google Scholar 
Angulo, E. et al. Non-English languages enrich scientific knowledge: The example of economic costs of biological invasions. Sci. Total Environ. 775, 144441 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Taylor, R., & Edwards, G. A review of the impact and control of cane toads in Australia with recommendations for future research and management approaches. A report to the Vertebrate Pests Committee from the National Cane Toad Taskforce (2005).Burnett, K., Pongkijvorasin, S. & Roumasset, J. Species invasion as catastrophe: The case of the brown tree snake. Environ. Resour. Econ. 51(2), 241–254 (2012).
Google Scholar 
Haubrock, P. J., Cuthbert, R. N., Ricciardi, A., Diagne, C., & Courchamp, F. Massive economic costs of invasive bivalves in freshwater ecosystems (2021).Pimentel, D., Zuniga, R. & Morrison, D. Update on the environmental and economic costs associated with alien-invasive species in the United States. Ecol. Econ. 52(3), 273–288 (2005).
Google Scholar 
Fantle-Lepczyk, J. E. et al. Economic costs of biological invasions in the United States. bioRxiv 89, 89 (2021).
Google Scholar 
European Environment Agency. The impacts of invasive alien species in Europe. Publications Office of the European Union (2013).Measey, J. et al. Invasive amphibians in southern Africa: A review of invasion pathways. Bothalia-Afr. Biodivers. Conserv. 47(2), 1–12 (2017).
Google Scholar 
Anton, A., Geraldi, N. R., Ricciardi, A. & Dick, J. T. Global determinants of prey naiveté to exotic predators. Proc. R. Soc. B 287(1928), 20192978 (2020).PubMed 
PubMed Central 

Google Scholar 
Rico-Sánchez, A. E. et al. Economic costs of invasive alien species in Mexico. NeoBiota 67, 459–483 (2021).
Google Scholar 
McNeely, J. Invasive species: A costly catastrophe for native biodiversity. Land Use Water Resour. Res. 1(1732-2016-140260) (2001).Sax, D. F. & Gaines, S. D. Species invasions and extinction: The future of native biodiversity on islands. Proc. Natl. Acad. Sci. 105, 11490–11497 (2008).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Christie, M., Fazey, I., Cooper, R., Hyde, T. & Kenter, J. O. An evaluation of monetary and non-monetary techniques for assessing the importance of biodiversity and ecosystem services to people in countries with developing economies. Ecol. Econ. 83, 67–78 (2012).
Google Scholar 
Essl, F. et al. Socioeconomic legacy yields an invasion debt. Proc. Natl. Acad. Sci. 108(1), 203–207 (2011).ADS 
CAS 
PubMed 

Google Scholar 
Seebens, H. et al. No saturation in the accumulation of alien species worldwide. Nat. Commun. 8(1), 1–9 (2017).
Google Scholar 
Ahmed, D. A., Hudgins, E. J., Cuthbert, R. N., Kourantidou, M., Diagne, C., Haubrock, P. J., et al. Managing biological invasions: The cost of inaction. Biol. Invasions. 1–20. (2022).Leung, B. et al. An ounce of prevention or a pound of cure: Bioeconomic risk analysis of invasive species. Proc. R. Soc. Lond. Ser. B Biol. Sci. 269(1508), 2407–2413 (2002).
Google Scholar 
Haubrock, P. J. et al. Geographic and taxonomic trends of rising biological invasion costs. Sci. Total Environ. 817, 152948 (2022).ADS 
CAS 
PubMed 

Google Scholar 
Kopecký, O., Kalous, L. & Patoka, J. Establishment risk from pet-trade freshwater turtles in the European Union. Knowl. Manag. Aquat. Ecosyst. 410, 02 (2013).
Google Scholar 
Mohanty, N. P. & Measey, J. The global pet trade in amphibians: Species traits, taxonomic bias, and future directions. Biodivers. Conserv. 28(14), 3915–3923 (2019).
Google Scholar 
Altherr, S. & Lameter, K. The rush for the rare: Reptiles and amphibians in the European pet trade. Animals 10, 2085 (2020).PubMed Central 

Google Scholar 
Cuthbert, R. N. et al. Biological invasion costs reveal insufficient proactive management worldwide. Sci. Total Environ. 819, 153404 (2022).ADS 
CAS 
PubMed 

Google Scholar 
Ricciardi, A. Invasive species. In Ecological Systems 161–178. (Springer, 2013).Leroy, B., Kramer, A. M., Vaissière, A. C., Courchamp, F., & Diagne, C. Analysing global economic costs of invasive alien species with the invacost R package. bioRxiv (2020). More

Tedersoo L, Bahram M, Põlme S, Kõljalg U, Yorou NS, Wijesundera R, et al. Fungal biogeography. Global diversity and geography of soil fungi. Science (80-). 2014;346:1256688.Article 
CAS 

Google Scholar 
Oliverio AM, Geisen S, Delgado Baquerizo M, Maestre FT, Turner BL, Fierer N. The global-scale distributions of soil protists and their contributions to belowground systems. Sci Adv. 2020;6:eaax8787.Article 
CAS 

Google Scholar 
Delgado Baquerizo M, Oliverio AM, Brewer TE, Benavent-Gonzalez A, Eldridge DJ, Bardgett RD, et al. A global atlas of the dominant bacteria found in soil. Science. 2018;359:320–5.CAS 
PubMed 
Article 

Google Scholar 
Bates ST, Clemente JC, Flores GE, Walters WA, Parfrey LW, Knight R, et al. Global biogeography of highly diverse protistan communities in soil. ISME J. 2012;7:652–9.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Crowther TW, van den Hoogen J, Wan J, Mayes MA, Keiser AD, Mo L, et al. The global soil community and its influence on biogeochemistry. Science. 2019;365:eaav0550.CAS 
PubMed 
Article 

Google Scholar 
Xiong W, Jousset A, Li R, Delgado-Baquerizo M, Bahram M, Logares R, et al. A global overview of the trophic structure within microbiomes across ecosystems. Environ Int. 2021;151:106438.PubMed 
Article 

Google Scholar 
Singh BK, Bardgett RD, Smith P, Reay DS. Microorganisms and climate change: terrestrial feedbacks and mitigation options. Nat Rev Microbiol. 2010;8:779–90.CAS 
PubMed 
Article 

Google Scholar 
Nowicka B, Kruk J. Powered by light: Phototrophy and photosynthesis in prokaryotes and its evolution. Microbiol Res. 2016;186-7:99–118.Article 
CAS 

Google Scholar 
Hamard S, Céréghino R, Barret M, Sytiuk A, Lara E, Dorrepaal E, et al. Contribution of microbial photosynthesis to peatland carbon uptake along a latitudinal gradient. J Ecol. 2021;109:3424–41.CAS 
Article 

Google Scholar 
Seppey CVW, Singer D, Dumack K, Fournier B, Belbahri LL, Mitchell EAD, et al. Distribution patterns of soil microbial eukaryotes suggests widespread algivory by phagotrophic protists as an alternative pathway for nutrient cycling. Soil Biol Biochem. 2017;112:68–76.CAS 
Article 

Google Scholar 
Schmidt O, Dyckmans J, Schrader S. Photoautotrophic microorganisms as a carbon source for temperate soil invertebrates. Biol Lett. 2016;12:20150646.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Halvorson HM, Barry JR, Lodato MB, Findlay RH, Francoeur SN, Kuehn KA. Periphytic algae decouple fungal activity from leaf litter decomposition via negative priming. Funct Ecol. 2019;33:188–201.PubMed 
Article 

Google Scholar 
Wyatt KH, Turetsky MR. Algae alleviate carbon limitation of heterotrophic bacteria in a boreal peatland. J Ecol. 2015;103:1165–71.CAS 
Article 

Google Scholar 
Elbert W, Weber B, Burrows S, Steinkamp J, Büdel B, Andreae MO, et al. Contribution of cryptogamic covers to the global cycles of carbon and nitrogen. Nat Geosci. 2012;5:459–62.CAS 
Article 

Google Scholar 
Jassey VEJ, Walcker R, Kardol P, Geisen S, Heger T, Lamentowicz M, et al. Contribution of soil algae to the global carbon cycle. New Phytol. 2022;234:64–76.CAS 
PubMed 
Article 

Google Scholar 
Tahon G, Tytgat B, Willems A. Diversity of phototrophic genes suggests multiple bacteria may be able to exploit sunlight in exposed soils from the Sør Rondane Mountains, East Antarctica. Front Microbiol. 2016;7:2026.PubMed 
PubMed Central 
Article 

Google Scholar 
Maier S, Tamm A, Wu D, Caesar J, Grube M, Weber B. Photoautotrophic organisms control microbial abundance, diversity, and physiology in different types of biological soil crusts. ISME J. 2018;12:1032–46.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Büdel B. Ecology and diversity of rock-inhabiting cyanobacteria in tropical regions. Eur J Phycol. 1999;34:361–70.Article 

Google Scholar 
Hamard S, Küttim M, Céréghino R, Jassey VEJ. Peatland microhabitat heterogeneity drives phototrophic microbes distribution and photosynthetic activity. Environ Microbiol. 2021;23:6811–27.CAS 
PubMed 
Article 

Google Scholar 
Cano-Díaz C, Maestre FT, Eldridge DJ, Singh BK, Bardgett RD, Fierer N, et al. Contrasting environmental preferences of photosynthetic and non-photosynthetic soil cyanobacteria across the globe. Glob Ecol Biogeogr. 2020;29:2025–38.Article 

Google Scholar 
Rodriguez-Caballero E, Belnap J, Büdel B, Crutzen PJ, Andreae MO, Pöschl U, et al. Dryland photoautotrophic soil surface communities endangered by global change. Nat Geosci. 2018;11:185–9.CAS 
Article 

Google Scholar 
Pointing SB, Belnap J. Microbial colonization and controls in dryland systems. Nat Rev Microbiol. 2012;10:551–62.CAS 
PubMed 
Article 

Google Scholar 
Bates ST, Clemente JC, Flores GE, Walters WA, Parfrey LW, Knight R, et al. Global biogeography of highly diverse protistan communities in soil. ISME J. 2013;7:652–9.CAS 
PubMed 
Article 

Google Scholar 
Küttim L, Küttim M, Puusepp L, Sugita S. The effects of ecotope, microtopography and environmental variables on diatom assemblages in hemiboreal bogs in Northern Europe. Hydrobiologia. 2017;792:137–49.Article 
CAS 

Google Scholar 
Mahé F, de Vargas C, Bass D, Czech L, Stamatakis A, Lara E, et al. Parasites dominate hyperdiverse soil protist communities in Neotropical rainforests. Nat Ecol Evol. 2017;1:91.PubMed 
Article 

Google Scholar 
Lindo Z, Gonzalez A. The Bryosphere: An Integral and Influential Component of the Earth’s Biosphere. Ecosystems. 2010;13:612–27.Article 

Google Scholar 
Sporn SG, Bos MM, Kessler M, Gradstein SR. Vertical distribution of epiphytic bryophytes in an Indonesian rainforest. Biodivers Conserv. 2010;19:745–60.Article 

Google Scholar 
Cornelissen JHC, Lang SI, Soudzilovskaia NA, During HJ. Comparative cryptogam ecology: a review of bryophyte and lichen traits that drive biogeochemistry. Ann Bot. 2007;99:987–1001.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Van Breemen N. How Sphagnum bogs down other plants. Trends Ecol Evol. 1995;10:270–5.PubMed 
Article 

Google Scholar 
Jonsson M, Kardol P, Gundale MJ, Bansal S, Nilsson M-C, Metcalfe DB, et al. Direct and Indirect Drivers of Moss Community Structure, Function, and Associated Microfauna Across a Successional Gradient. Ecosystems. 2014;18:1–16.
Google Scholar 
Bragina A, Berg C, Cardinale M, Shcherbakov A, Chebotar V, Berg G. Sphagnum mosses harbour highly specific bacterial diversity during their whole lifecycle. ISME J. 2012;6:802–13.CAS 
PubMed 
Article 

Google Scholar 
Bay G, Nahar N, Oubre M, Whitehouse MJ, Wardle DA, Zackrisson O, et al. Boreal feather mosses secrete chemical signals to gain nitrogen. New Phytol. 2013;200:54–60.CAS 
PubMed 
Article 

Google Scholar 
Kip N, van Winden JF, Pan Y, Bodrossy L, Reichart G-J, Smolders AJP, et al. Global prevalence of methane oxidation by symbiotic bacteria in peat-moss ecosystems. Nat Geosci. 2010;3:617–21.CAS 
Article 

Google Scholar 
Lindo Z, Nilsson M-C, Gundale MJ. Bryophyte-cyanobacteria associations as regulators of the northern latitude carbon balance in response to global change. Glob Chang Biol. 2013;19:2022–35.PubMed 
Article 

Google Scholar 
Jassey VEJ, Shimano S, Dupuy C, Toussaint M-L, Gilbert D. Characterizing the feeding habits of the testate amoebae Hyalosphenia papilio and Nebela tincta along a narrow ‘fen-bog’ gradient using digestive vacuole content and 13C and 15N isotopic analyses. Protist. 2012;163:451–64.PubMed 
Article 

Google Scholar 
Raanan H, Oren N, Treves H, Keren N, Ohad I, Berkowicz SM, et al. Towards clarifying what distinguishes cyanobacteria able to resurrect after desiccation from those that cannot: The photosynthetic aspect. Biochim Biophys Acta – Bioenerg. 2016;1857:715–22.CAS 
Article 

Google Scholar 
Puente-Sánchez F, Arce-Rodríguez A, Oggerin M, García-Villadangos M, Moreno-Paz M, Blanco Y, et al. Viable cyanobacteria in the deep continental subsurface. Proc Natl Acad Sci USA. 2018;115:10702–7.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Küttim M, Küttim L, Ilomets M, Laine AM. Controls of Sphagnum growth and the role of winter. Ecol Res. 2020;35:219–34.Article 
CAS 

Google Scholar 
Jassey VEJ, Chiapusio G, Mitchell EAD, Binet P, Toussaint M-L, Gilbert D. Fine-scale horizontal and vertical micro-distribution patterns of testate amoebae along a narrow Fen/Bog gradient. Microb Ecol. 2011;61:374–85.PubMed 
Article 

Google Scholar 
Wilken S, Huisman J, Naus-Wiezer S, Van Donk E. Mixotrophic organisms become more heterotrophic with rising temperature. Ecol Lett. 2012;16:225–33.PubMed 
Article 

Google Scholar 
Jassey VEJ, Signarbieux C. Effects of climate warming on Sphagnumphotosynthesis in peatlands depend on peat moisture and species‐specific anatomical traits. Glob Chang Biol. 2019;182:12–65.
Google Scholar 
McDonald D, Price MN, Goodrich J, Nawrocki EP, Desantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2011;6:610–8. 2012 63PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Lozupone CA, Stombaugh J, Gonzalez A, Ackermann G, Wendel D, Vázquez-Baeza Y, et al. Meta-analyses of studies of the human microbiota. Genome Res. 2013;23:1704–14.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pawluczyk M, Weiss J, Links MG, Egaña Aranguren M, Wilkinson MD, Egea-Cortines M. Quantitative evaluation of bias in PCR amplification and next-generation sequencing derived from metabarcoding samples. Anal Bioanal Chem. 2015;407:1841–8.CAS 
PubMed 
Article 

Google Scholar 
Ramirez KS, Knight CG, de Hollander M, Brearley FQ, Constantinides B, Cotton A, et al. Detecting macroecological patterns in bacterial communities across independent studies of global soils. Nat Microbiol. 2018;3:189–96.CAS 
PubMed 
Article 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2019.Opelt K, Berg C, Schönmann S, Eberl L, Berg G. High specificity but contrasting biodiversity of Sphagnum-associated bacterial and plant communities in bog ecosystems independent of the geographical region. ISME J. 2007;1:502–16.CAS 
PubMed 
Article 

Google Scholar 
Hamard S, Robroek BJM, Allard P-M, Signarbieux C, Zhou S, Saesong T, et al. Effects of Sphagnum Leachate on Competitive Sphagnum Microbiome Depend on Species and Time. Front Microbiol. 2019;10:3317.Article 

Google Scholar 
Holland-Moritz H, Stuart J, Lewis LR, Miller S, Mack MC, McDaniel SF, et al. Novel bacterial lineages associated with boreal moss species. Environ Microbiol. 2018;20:2625–38.CAS 
PubMed 
Article 

Google Scholar 
Singer D, Metz S, Unrein F, Shimano S, Mazei Y, Mitchell EAD, et al. Contrasted Micro-Eukaryotic Diversity Associated with Sphagnum Mosses in Tropical, Subtropical and Temperate Climatic Zones. Microb Ecol. 2019;78:714–24.CAS 
PubMed 
Article 

Google Scholar 
Holland-Moritz H, Stuart JEM, Lewis LR, Miller SN, Mack MC, Ponciano JM, et al. The bacterial communities of Alaskan mosses and their contributions to N2-fixation. Microbiome. 2021;9:1–14.Article 
CAS 

Google Scholar 
Righetti D, Vogt M, Gruber N, Psomas A, Zimmermann NE. Global pattern of phytoplankton diversity driven by temperature and environmental variability. Sci Adv. 2019;5:eaau6253.PubMed 
PubMed Central 
Article 

Google Scholar 
Amend AS, Cobian GM, Laruson AJ, Remple K, Tucker SJ, Poff KE, et al. Phytobiomes are compositionally nested from the ground up. PeerJ. 2019;2019:e6609.Article 

Google Scholar 
Dedysh SN, Pankratov TA, Belova SE, Kulichevskaya IS, Liesack W. Phylogenetic analysis and in situ identification of Bacteria community composition in an acidic Sphagnum peat bog. Appl Environ Microbiol. 2006;72:2110–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Robroek BJM, Martí M, Svensson BH, Dumont MG, Veraart AJ, Jassey VEJ. Rewiring of peatland plant–microbe networks outpaces species turnover. Oikos. 2021;303:605–15.
Google Scholar 
Holland-Moritz H, Stuart J, Lewis LR, Miller S, Mack MC, Mcdaniel SF, et al. Novel bacterial lineages associated with boreal moss species. Environ Microbiol. 2018;20:2625–38.CAS 
PubMed 
Article 

Google Scholar 
Sytiuk A, Céréghino R, Hamard S, Delarue F, Guittet A, Barel JM, et al. Predicting the structure and functions of peatland microbial communities from Sphagnum phylogeny, anatomical and morphological traits and metabolites. J Ecol. 2021;1365-2745:13728.
Google Scholar 
Rudolph H, Samland J. Occurrence and metabolism of sphagnum acid in the cell walls of bryophytes. Phytochemistry. 1985;24:745–9.CAS 
Article 

Google Scholar 
Chiapusio G, Jassey VEJ, Bellvert F, Comte G, Weston LA, Delarue F, et al. Sphagnum species modulate their phenolic profiles and mycorrhizal colonization of surrounding Andromeda polifolia along peatland microhabitats. J Chem Ecol. 2018;27:1–12.
Google Scholar 
Rasmussen S, Wolff C, Rudolph H. Compartmentalization of phenolic constituents in sphagnum. Phytochemistry. 1995;38:35–39.CAS 
Article 

Google Scholar 
Sytiuk A, Céréghino R, Hamard S, Delarue F, Dorrepaal E, Küttim M, et al. Biochemical traits enhance the trait concept in Sphagnum ecology. Oikos 2022;00:00.Hájek T, Ballance S, Limpens J, Zijlstra M, Verhoeven JTA. Cell-wall polysaccharides play an important role in decay resistance of Sphagnum and actively depressed decomposition in vitro. Biogeochemistry. 2011;103:45–57.Article 
CAS 

Google Scholar 
Bengtsson F, Rydin Hå, Hájek T. Biochemical determinants of litter quality in 15 species of Sphagnum. Plant Soil. 2018;425:161–76.CAS 
Article 

Google Scholar 
Fudyma JD, Lyon J, AminiTabrizi R, Gieschen H, Chu RK, Hoyt DW, et al. Untargeted metabolomic profiling of Sphagnum fallax reveals novel antimicrobial metabolites. Plant Direct. 2019;3:e00179–17.Article 

Google Scholar 
He L, Mazza Rodrigues JL, Soudzilovskaia NA, Barceló M, Olsson PA, Song C, et al. Global biogeography of fungal and bacterial biomass carbon in topsoil. Soil Biol Biochem. 2020;151:108024.CAS 
Article 

Google Scholar 
Hanson CA. Beyond biogeographic patterns: processes shaping the microbial landscape. Nat Rev Microbiol. 2012;10:497–506.CAS 
PubMed 
Article 

Google Scholar 
Waddington JM, Morris PJ, Kettridge N, Granath G, Thompson DK, Moore PA. Hydrological feedbacks in northern peatlands. Ecohydrology. 2015;8:113–27.Article 

Google Scholar 
Reczuga MK, Lamentowicz M, Mulot M, Mitchell EAD, Buttler A, Chojnicki B, et al. Predator–prey mass ratio drives microbial activity under dry conditions in Sphagnum peatlands. Ecol Evol. 2018;8:5752–64.PubMed 
PubMed Central 
Article 

Google Scholar 
Ritchie RJ. Consistent sets of spectrophotometric chlorophyll equations for acetone, methanol and ethanol solvents. Photosynth Res. 2006;89:27–41.CAS 
PubMed 
Article 

Google Scholar 
Perrine Z, Negi S, Sayre RT. Optimization of photosynthetic light energy utilization by microalgae. Algal Res. 2012;1:134–42.Article 

Google Scholar 
Carey CC, Ibelings BW, Hoffmann EP, Hamilton DP, Brookes JD. Eco-physiological adaptations that favour freshwater cyanobacteria in a changing climate. Water Res. 2012;46:1394–407.CAS 
PubMed 
Article 

Google Scholar 
Gorbunov MY, Falkowski PG. Using chlorophyll fluorescence kinetics to determine photosynthesis in aquatic ecosystems. Limnol Ocean. 2020;66:1–13.Article 
CAS 

Google Scholar 
MacIntyre HL, Kana TM, Anning T, Geider RJ. Photoacclimation of photosynthesis irradiance response curves and photosynthetic pigments in microalgae and cyanobacteria. J Phycol. 2002;38:17–38.Article 

Google Scholar 
Grote EE, Belnap J, Housman DC, Sparks JP. Carbon exchange in biological soil crust communities under differential temperatures and soil water contents: implications for global change. Glob Chang Biol. 2010;16:2763–74.Article 

Google Scholar 
Robarts RD, Zohary T. Temperature effects on photosynthetic capacity, respiration, and growth rates of bloom‐forming cyanobacteria. New Zealand Journal of Marine and Freshwater Research. 1987;21:391–9.CAS 
Article 

Google Scholar 
Le Quéré C, Andrew RM, Friedlingstein P, Sitch S, Pongratz J, Manning AC, et al. Global Carbon Budget 2017. Earth Syst Sci Data. 2018;10:405–48.Article 

Google Scholar  More

Joos, F., Plattner, G. K., Stocker, T. F., Marchal, O. & Schmittner, A. Global warming and marine carbon cycle feedbacks on future atmospheric CO2. Science 284(5413), 464–467 (1999).ADS 
CAS 
Article 

Google Scholar 
Smith, K. L. et al. Climate, carbon cycling, and deep-ocean ecosystems. Proc. Nat. Acad. Sci USA 106, 19211–19218 (2009).ADS 
CAS 
Article 

Google Scholar 
Ramirez-Llodra, E. et al. Man and the last great wilderness: Human impact on the deep sea. PLoS ONE 6, e22588 (2011).ADS 
CAS 
Article 

Google Scholar 
Pham, C. K. et al. Marine litter distribution and density in European Seas, from the shelves to deep basins. PLoS ONE 9, e95839 (2014).ADS 
Article 

Google Scholar 
Angel, M. What is the deep sea? In Deep-sea fishes (eds Randall, D. & Farrell, A.) 1–41 (Academic Publishing, 1997).
Google Scholar 
Smith, C. R., De Leo, F. C., Bernardino, A. F., Sweetman, A. K. & Arbizu, P. M. Abyssal food limitation, ecosystem structure and climate change. Trends Ecol. Evol. 23, 518–528 (2008).Article 

Google Scholar 
Thurber, A. R. et al. Ecosystem function and services provided by the deep sea. Biogeosciences 11, 3941–3963 (2014).ADS 
Article 

Google Scholar 
Solan, M. et al. Extinction and ecosystem function in the marine benthos. Science 306(5699), 1177–1180 (2004).ADS 
CAS 
Article 

Google Scholar 
Danise, S., Twitchett, R. J., Little, C. T. & Clemence, M. E. The impact of global warming and anoxia on marine benthic community dynamics: An example from the Toarcian (Early Jurassic). PLoS ONE 8(2), e56255 (2013).ADS 
CAS 
Article 

Google Scholar 
Nomaki, H. et al. In situ experimental evidences for responses of abyssal benthic biota to shifts in phytodetritus compositions linked to global climate change. Glob. Chang. Biol. 27, 6139–6155 (2021).Article 

Google Scholar 
Viehman, H. A. & Zydlewski, G. B. Fish interactions with a commercial-scale tidal energy device in the natural environment. Estuaries Coast 38(1), 241–252 (2015).Article 

Google Scholar 
Danovaro, R. et al. Implementing and innovating marine monitoring approaches for assessing marine environmental status. Front. Mar. Sci. 3, 213 (2016).Article 

Google Scholar 
Mizuno, K. et al. An efficient coral survey method based on a large-scale 3-D structure model obtained by Speedy Sea Scanner and U-Net segmentation. Sci. Rep. 10(1), 12416. https://doi.org/10.1038/s41598-020-69400-5 (2020).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Eleftheriou, A., & Moore, D. C. (2013). Macrofauna techniques. Methods for the study of marine benthos, 175–251.Solan, M. et al. In situ quantification of bioturbation using time lapse fluorescent sediment profile imaging (f SPI), luminophore tracers and model simulation. Mar. Ecol. Prog. Ser. 271, 1–12 (2004).ADS 
Article 

Google Scholar 
Hale, R. et al. High-resolution computed tomography reconstructions of invertebrate burrow systems. Sci. Data 2(1), 1–5 (2015).Article 

Google Scholar 
Plets, R. M. et al. The use of a high-resolution 3D Chirp sub-bottom profiler for the reconstruction of the shallow water archaeological site of the Grace Dieu (1439), River Hamble, UK. J. Archaeol. Sci. 36(2), 408–418 (2009).Article 

Google Scholar 
Mizuno, K. et al. Automatic non-destructive three-dimensional acoustic coring system for in situ detection of aquatic plant root under the water bottom. Case Stud. Nondestruct. Test. Evaluat. 5, 1–8 (2016).CAS 
Article 

Google Scholar 
Suganuma, H., Mizuno, K. & Asada, A. Application of wavelet shrinkage to acoustic imaging of buried asari clams using high-frequency ultrasound. J. Appl. Phys. 57(7S1), 07LG08 (2018).Article 

Google Scholar 
Dorgan, K. M. et al. Impacts of simulated infaunal activities on acoustic wave propagation in marine sediments. J. Acoust. Soc. Am. 147(2), 812–823 (2020).ADS 
Article 

Google Scholar 
Mizuno, K., Cristini, P., Komatitsch, D. & Capdeville, Y. Numerical and experimental study of wave propagation in water-saturated granular media using effective method theories and a full-wave numerical simulation. IEEE J. Ocean. Eng. 45(3), 772–785 (2020).ADS 
Article 

Google Scholar 
Schulze, I. et al. Laboratory measurements to image endobenthos and bioturbation with a high-frequency 3D seismic lander. Geosciences 11(12), 508 (2021).ADS 
Article 

Google Scholar 
Hashimoto, J. et al. Deep-sea communities dominated by the giant clam, Calyptogena soyoae, along the slope foot of Hatsushima Island, Sagami Bay, central Japan. Palaeogeogr. Palaeoclimatol. Palaeoecol. 71(12), 179–192 (1989).Article 

Google Scholar 
Fujikura, K., Hashimoto, J. & Okutani, T. Estimated population densities of megafauna in two chemosynthesisbased communities: A cold seep in Sagami Bay and a hydrothermal vent in the Okinawa Trough. Benthos. Res. 57(1), 21–30 (2002).Article 

Google Scholar 
Childress, J. J. & Girguis, P. R. The metabolic demands of endosymbiotic chemoautotrophic metabolism on host physiological capacities. J. Exp. Biol. 214(2), 312–325 (2011).CAS 
Article 

Google Scholar 
Okuba, K. (2021). Basic study on sonar system development for exploring infaunal bivalves. Master thesis, GSFS, The University of Tokyo (in Japanese).Stoll, R. D. & Bryan, G. M. Wave attenuation in saturated sediments. The J. Acoust. Soc. Am. 47(5B), 1440–1447 (1970).ADS 
Article 

Google Scholar 
Schwartz, L. & Plona, T. J. Ultrasonic propagation in close-packed disordered suspensions. J. Appl. Phys. 55(11), 3971–3977 (1984).ADS 
Article 

Google Scholar 
Seike, K., Shirai, K. & Murakami-Sugihara, N. Using tsunami deposits to determine the maximum depth of benthic burrowing. PLoS ONE 12(8), e0182753. https://doi.org/10.1371/journal.pone.0182753 (2017).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar  More

Our research involved work with animal subjects (unionid mussels and Round Goby fishes) and was conducted following relevant regulations and standard procedures. The field collections were carried out under Pennsylvania Fish and Boat Commission permits (# 2018-01-0136 and 2019-01-0026). The experimental protocols were approved by Penn State University’s Institutional Animal Care and Use Committee (IACUC# 201646941 and 201646962). All new DNA sequencing data are made publicly available in GenBank (with accession numbers provided in Table 1) and a BioProject (# PRJNA813547) of the National Center for Biotechnology Information40.Propensity of Round Goby to consume unionid mussels in a controlled lab settingStream table setupWe conducted lab experiments to observe the potential predation of juvenile freshwater mussels by the Round Goby, following standard research protocols for work with animal subjects (IACUC# 201646962, Penn State University). We constructed four artificial stream tables in an aquatic laboratory, each measuring 3 × 2 m and featuring two run and two pool sections (each 0.63 × 0.56 × 0.46 m). Water flow was produced using eight Homsay 920 GPH submersible water pumps, which pumped water from a central reservoir tub into each table at the start of each run section. The water flow direction was clockwise for stream tables 1 and 3, and counterclockwise for stream tables 2 and 4. Water pumped into the stream tables exited via two drains located medially of each run section, where it flowed back to the central reservoir tub. Each stream table was filled with a 6 mm layer of substrate consisting of a mixture of sand, gravel (4–6 mm), and crushed stone (size 2B, with an average size of ~ 19 mm). The day before each experiment, field technicians traveled to local streams and collected macroinvertebrates using one minute D-frame kick net samples for each of the four stream tables. The macroinvertebrates and associated substrate were transported back to the facility and were introduced into each stream table system.Preferential feeding experimentsBefore each experiment commenced, juvenile Plain Pocketbook mussels (Lampsillis cardium) were introduced into each stream table (with 165 mussel specimens for experiment 1 and 100 mussels for experiments 2 and 3). This  widespread and abundant species is not imperiled in Pennsylvania, and mussels were provided for this study by the White Sulphur Springs National Fish Hatchery located in southeast West Virgina. The mussels were allowed to acclimate in the stream tables for 2 h before commencing each experiment. Ten Round Gobies were introduced into each stream system (stream tables 1 and 2 for experiment 1, and all tables for experiments 2 and 3). The total length (from nose tip to caudal tip) of each fish was measured prior to introduction and after the termination of experiments 2 and 3. Experiment 1 was conducted for 3 weeks, while experiments 2 and 3 were conducted for 8 days. During these experiments, Round Gobies were allowed to exist in the systems and feed preferentially, on the mussels and macroinvertebrates, for the allotted time before each investigation concluded. We acknowledge that in these experiments, the mussel abundances are higher and macroinvertebrate densities lower and less rich than commonly occur in the natural stream environment. Further, the Round Goby fish densities used are much higher than currently in the French Creek watershed, though are comparable to what is currently seen in parts of the Great Lakes basin. Nonetheless, the experiment scenarios allowed us to observe if Round Gobies would consume the mussels when given the choice to feed on a variety of food items.Evaluation of unionids consumed by fishRound Gobies were removed from the stream tables upon completion of each experiment. They were euthanized using  > 250 mg/L buffered (pH ~ 7) tricaine-S (MS222) solution. The fish were submerged for 10 min beyond the cessation of opercular movement to ensure proper euthanasia, and tissues were collected after we confirmed complete euthanasia—compliant with AVMA guidelines and approved by the IACUC protocol. The Round Gobies were placed in a 10% solution of formalin for preservation, and after 2 weeks, they were rinsed with clean water and were placed in 70% ethanol for long-term storage. After fish were removed from the system, the water was drained, and the substrate was sifted to recover the remaining mussels. Mussels were counted, and live individuals were returned to holding tanks for use in subsequent experiments. To further assess whether Round Gobies had consumed mussels during the investigation, Round Gobies were x-rayed using a Bruker Skyscan 1176 micro-CT scanner. After that, the stomachs of each fish were excised, and the contents examined using a Leica CME dissection scope to confirm the identity of Plain Pocketbook mussels. Contents posterior to the stomach were not analyzed because they could not be reliably counted and identified.DNA metabarcoding to identify mussel species consumed by Round Goby in a stream settingFish and mussel sample acquisitionWe collected 39 Round Gobies directly from streams within the French Creek watershed—their newly invaded natural stream habitats—in June 2018. We aimed to quantify which species, if any, of unionid mussels they consumed. Fish collection locations included LeBoeuf Creek at Moore Road and 100 m below the confluence of French Creek and LeBoeuf Creek. The unionid mussel populations and the environmental field settings at these locations are detailed by Clark et al.19. A team of field technicians collected fish by kick seining (3 m × 1 m × 9.5 mm nylon mesh) while moving downstream. Seining was the sampling method of choice compared to electrofishing to avoid possible regurgitation of food items prior to excision of fishes’ stomachs. The stream reaches sampled at each location were between 100 and 200 m in length and included riffle, run, and pool habitats. In addition to fish samples, unionid mussel samples from French Creek were also collected for analysis (under Pennsylvania Fish & Boat Commission collectors permits # 2018-01-0136 and 2019-01-0026). Following standard research protocols (under IACUC# 201646941, Penn State University), the Round Gobies collected were euthanized using buffered Tricaine-S (MS222) solution; and stomachs were excised using sterilized utensils before being placed in sterilized tubes filled with 97% ethanol. After excision of stomachs, fishes were placed in a 10% formalin solution for preservation. After 2 weeks, fishes were rinsed with clean water and transferred to 70% ethanol for storage. The stomach samples were immediately placed in ethanol and on ice in the field. Samples were stored in a freezer before being shipped to the US Geological Survey’s Eastern Ecological Science Center for various molecular ecology analyses. Once the fish and mussel samples arrived at this lab, they were recorded and stored at four °C until analysis.Primer developmentSpecific primers targeting a moderately conserved region of the mitochondrial COI gene for 25 species of unionids inhabiting French Creek were designed. Previously a PCR-based amplification method utilizing restriction enzyme digests was used to identify genetic fingerprints of 25 unionid species inhabiting French Creek41. Here, we designed a new degenerate PCR primer set modified with sequencing overhangs to facilitate compatibility with a MiSeq amplicon sequencing method previously designed for 16S Amplicon sequencing. We targeted the locus of the mitochondrial COI gene of unionids known to inhabit the Atlantic Slope Drainage. Consensus sequences were derived using Multalin analysis and a tiling method to identify conserved primer binding regions flanking an ~ 300 bp region of the COI gene. This gene was targeted in part due to the availability of partial or complete sequences representing these target species in the NCBI reference database40. Cytochrome oxidase sequences were downloaded for the 25 unionid mussel species of interest. However, a COI sequence for the Rabbitsfoot (Theliderma cylindrica) mussel was absent from the NCBI database, which required us to sequence this region for an in-house reference (which is described later in the paper). We designed a degenerate primer cocktail specific to all mussel species of interest that amplified a ~ 289 bp product, with forward and reverse primers used for the amplification of unionid specific COI presented as supplemental information (see Table S-230. We evaluated the suitability of the primers using samples from field identified mussels. For primer optimization, PCR was run across a gradient of annealing temperatures to determine suitability. In addition, we used Round Goby DNA as a template to evaluate specificity. In addition to Round Goby stomach samples, mussel samples of several species collected from French Creek were included as positive controls.DNA extraction from tissue samplesFollowing the manufacturer’s protocols, tissue samples (including fish stomach and mussel tissue) were extracted with the Zymo Research ZymoBIOMICS 96 MagBead DNA Kit (San Diego, CA). Random samples of DNA extracts were analyzed on an Agilent 2100 Bioanalyzer using a high-sensitivity assay kit. Fragments in the target amplicon range were apparent (albeit not known to be of mussel origin). All samples were stored at − 20 °C until PCR was performed. DNA from both the T. cylindrica and L. complanata samples were analyzed for DNA quality.Rolling circle amplification of mitochondrial genomesTo acquire COI sequences for T. cylindrica and L. complanata, we subjected archived DNA samples to rolling circle amplification (RCA) followed by amplicon sequencing on the MiSeq. In short, 2 µl of DNA template was added to 2 µl Equiphi29 DNA polymerase reaction buffer containing 1 µl of Exonuclease-resistant random primers (ThermoFisher). Samples were denatured by heating to 95 °C for 3 min followed immediately by cooling on ice for more than 5 min. A volume of 5 µl was added to an RCA master mix containing 1.5 µl of 10 × Equiphi29 DNA polymerase reaction buffer, 0.2 µl of 100 mM dithiothreitol, 8 µl of 2.5 mM dNTPs, 1 µl of Eqiphi29 DNA polymerase (10U) and 4.3 µl of nuclease-free water. The samples were heated to 45 °C for 3 h and then 65 °C for 10 min. Samples were then placed in ice and then frozen at − 20 °C. All RCA products were normalized to 0.2 ng/µl in 10 mM Tris–HCl, pH 8.5. Normalized RCA product was utilized as a template for an Illumina Nextera XT library preparation. Sequencing libraries were prepared following the Nextera XT Library Preparation Reference Guide (CT# 15031942 v01) using the Nextera XT Library Preparation Kit (Illumina, San Diego, CA). Final libraries were analyzed for size and quality using the Agilent BioAnalyzer with the accompanying DNA 1000 Kit (Agilent, Santa Clara, CA). Libraries were quantified using the Qubit H.S. Assay Kit (Invitrogen, Carlsbad, CA) and normalized to 4 nM using 10 mM Tris, pH 8.5. Libraries were pooled and run on the Illumina MiSeq at a concentration of 10 pM with a 5% PhiX spike with run parameters of 1 × 150. Bioinformatic processing of this data is outlined below.Amplification of the cytochrome oxidase 1 geneExtracted genomic DNA was used as template for end-point PCR. Samples evaluated were from mussels and round gobies (see supporting Table S-330). The ~ 289 bp COI region was amplified with the mussel primers as follows. The amplification reaction contained 0.15 µM of each primer, 1 µL of the initial amplification product, and Promega Go Taq Green Master Mix following manufacturer recommendations for a 25 µL reaction. The thermocycler program consisted of an initial denaturing step of 95 °C for 3 min, followed by 30 cycles of 30 s at 95 °C, 30 s at 52 °C, and 1 min at 72 °C. Products were subjected to a final extension of 72 °C for 5 min then held until collection at 12 °C. An appropriately sized amplification product was confirmed for each reaction by electrophoresis of 5 µL of the reaction product through a 1.5% I.D. N.A. agarose gel (FMC Bioproducts) at 100 V for 45 min. PCR products were cleaned with the Qiagen Qiaquick PCR purification kit (Valencia, CA) and quantified using the Qubit dsDNA H.S. Assay Kit (Thermofisher Scientific, Grand Island, NY). Samples were diluted in 10 mM Tris buffer (pH 8.5) to a final concentration of 5 ng/µL.Generation of mock mussel samplesTo better understand and minimize sources of error or bias in the taxonomic assignment, we created a mock extraction by mixing sequences from known mussel taxa at defined concentrations. For each mussel, approximately 25-mg of tissue was extracted with the ZymoBIOMICS 96 MagBead DNA Kit (San Diego, CA) following the manufacturer’s protocol. The COI sequence was amplified from each species using the same primer-protocol combination described above. A total of 5 PCR products were mixed at equal concentration (mass/volume) to generate the mock sample (“Mock” hereafter). To confirm the identity of these inputs, each COI region was amplified and sequenced on the Illumina MiSeq during the same run as the Mock and samples.Sequencing library preparation and quality assessmentNext-generation sequencing was performed on the Illumina MiSeq platform to observe species-specific sequences and determine the diet of the Round Goby. Inclusion of the overhangs on the amplification primers allowed us to utilize the Illumina 16S Metagenomic Sequencing Library Preparation protocol42. Amplicon libraries were prepared following the same manufacturer’s protocol. All samples were indexed using the Illumina Nextera XT multiplex library indices. DNA read size spectra were determined with the Agilent 2100 Bioanalyzer using the Agilent DNA 1000 Kit (Santa Clara, Calif.). Libraries were quantified with the Qubit dsDNA H.S. Assay Kit (ThermoFisher Scientific, Grand Island, N.Y.) and normalized to 4 nM (nM) using 10 mM (mM) Tris (hydroxymethyl) aminomethane buffer pH 8.5. A final concentration of 10 picomolar library with a 6.5% PhiX control spike was created with the combined pool of all indexed libraries. All bioinformatic operations were completed on CLC Genomic Workbench v20 (Qiagen, Valencia, Calif.).Read filtering, trimming, and RNAseq metabarcoding assemblyFASTQ files from the sequencing runs were imported as paired-end reads into CLC Genomics Workbench v20.0.4 (Qiagen Bioinformatics, Redwood City, Calif.) for initial filtering of exogenous sequence adaptors and poor-quality base calls. The trimmed overlapping paired-end reads were mapped to the 25 target unionid sequences specific for the species of interest. Several mapping iterations were run using different levels of stringency. We utilized + 2/− 3 match-mismatch scoring and set the length fraction to 0.90. Analyses were iterated using different similarity fractions ranging from 0.90 to 0.99. Reads were annotated, and relative abundance was determined using a curated reference library (see supporting Datasets S-1 and S-230). More

Hughes, W. O. H., Oldroyd, B. P., Beekman, M. & Ratnieks, F. L. W. Ancestral monogamy shows kin selection is key to the evolution of eusociality. Science 320, 1213–1216 (2008).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Ross, K. G. & Matthews, R. W. The Social Biology of Wasps (Cornell University Press, 1991).Book 

Google Scholar 
Itô, Y. & Higashi, S. Spring behaviour of Ropalidia plebeiana (Hymenoptera: Vespidae) within a huge aggregation of nests. Appl. Entomol. Zool. 22, 519–527 (1987).Article 

Google Scholar 
Saito, F. & Kojima, J.-I. Colony cycle in the south-eastern coastal populations of Ropalidia plebeiana, the only Ropalidia wasp occurring in temperate Australia. Entomol. Sci. 8, 263–275 (2005).Article 

Google Scholar 
Richards, O. W. The Australian social wasps (Hymenoptera: Vespidae). Aust. J. Zool. Suppl. Ser. 26, 1–132 (1978).Article 

Google Scholar 
Makino, S., Yamane, S., Itô, Y. & Spradbery, J. P. Process of comb division of reused nests in the Australian paper wasp Ropalidia plebeinana (Hymenoptera, Vespidae). Ins. Soc. 41, 411–422 (1994).Article 

Google Scholar 
Goodnight, K. F. & Queller, D. C. Computer software for performing likelihood tests of pedigree relationship using genetic markers. Mol. Ecol. 8, 1231–1234 (1999).Article 

Google Scholar 
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).CAS 
PubMed 
Article 

Google Scholar 
Crochet, P.-A. Genetic structure of avian populations—Allozymes revisited. Mol. Ecol. 9, 1463–1469 (2000).CAS 
PubMed 
Article 

Google Scholar 
Meirmans, P. G. & Hedrick, P. W. Assessing population structure: FST and related measures. Mol. Ecol. Res. 11, 5–18 (2011).Article 

Google Scholar 
Hedrick, P. W. A standardized genetic differentiation measure. Evolution 59, 1633–1638 (2005).CAS 
PubMed 
Article 

Google Scholar 
Itô, Y., Yamane, S. & Spradbery, J. P. Population consequences of huge nesting aggregations of Ropalidia plebeiana (Hymenoptera: Vespidae). Res. Popul. Ecol. 30, 279–295 (1988).Article 

Google Scholar 
Boomsma, J. J. & d’Ettorre, P. Nice to kin and nasty to non-kin: revisiting Hamilton’s early insights on eusociality. Biol. Lett. 9, 20130444 (2013).PubMed 
PubMed Central 
Article 

Google Scholar 
Hannonen, M. & Sundström, L. Worker nepotism among polygynous ants. Nature 421, 910 (2003).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Parsons, P. J., Grinsted, L. & Field, J. Partner choice correlates with fine scale kin structuring in the paper wasp Polistes dominula. PLoS ONE 14, e0221701 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Strassmann, J. E., Queller, D. C., Solis, C. R. & Hughes, C. R. Relatedness and queen number in the Neotropical wasp, Parachartergus colobopterus. Anim. Behav. 42, 461–470 (1991).Article 

Google Scholar 
Leadbeater, E., Carruthers, J. M., Green, J. P., Rosser, N. S. & Field, J. Nest inheritance is the missing source of direct fitness in a primitively eusocial insect. Science 333, 874–876 (2011).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Field, J. & Leadbeater, E. Cooperation between non-relatives in a primitively eusocial paper wasp, Polistes dominula. Philos. Trans. R. Soc. B 371, 20150093 (2016).Article 

Google Scholar 
Bhadra, A. & Gadagkar, R. We know that the wasps “know”: cryptic successors to the queen in Ropalidia marginata. Biol. Lett. 4, 634–637 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
Bang, A. & Gadagkar, R. Reproductive queue without overt conflict in the primitively eusocial wasp Ropalidia marginata. Proc. Nat. Acad. Sci. USA 109, 14494–14499 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Frank, S. A. Hierarchical selection theory and sex ratios. II. On applying the theory, and a test with fig wasps. Evolution 39, 949–964 (1985).PubMed 
Article 

Google Scholar 
Silk, J. B. & Brown, G. R. Local resource competition and local resource enhancement shape primate birth sex ratios. Proc. R. Soc. Lond. B. 275, 1761–1765 (2008).
Google Scholar 
Schwarz, M. P. Local resource enhancement and sex ratios in a primitively social bee. Nature 331, 346–348 (1988).ADS 
Article 

Google Scholar 
Cronin, A. L. & Schwarz, M. P. Sex ratios, local fitness enhancement and eusociality in the allodapine bee Exoneura richardsoni. Evol. Ecol. 11, 567–577 (1997).Article 

Google Scholar 
Schwarz, M. P., Bull, N. J. & Hogendoorn, K. Evolution of sociality in the allodapine bees: A review of sex allocation, ecology and evolution. Ins. Soc. 45, 349–368 (1998).Article 

Google Scholar 
Gamboa, G. J., Wacker, T. L., Duffy, K. G., Dobson, S. W. & Fishwild, T. G. Defence against intraspecific usurpation by paper wasp cofoundresses (Polistes fuscatus, Hymenoptera: Vespidae). Can J. Zool. 70, 2369–2372 (1992).Article 

Google Scholar 
Katada, S. & Iwahashi, O. Characteristics of usurped colonies in the subtropical paper wasp, Ropalidia fasciata (Hymenoptera: Vespidae). Ins. Soc. 43, 247–253 (1996).Article 

Google Scholar 
Yamane, S. Ecological factors influencing the colony cycle of Polistes wasps. in Natural History and Evolution of Paper-Wasps (Turillazzi, S. & West-Eberhard, M. J. eds.). 75–97. (Oxford University Press, 1996).Clouse, R. Some effects of group size on the output of beginning nests of Mischocyttarus mexicanus (Hymenoptera: Vespidae). Flor. Entomol. 84, 418–425 (2001).Article 

Google Scholar 
Strassmann, J. E. Female-biased sex ratios in social insects lacking morphological castes. Evolution 38, 256–266 (1984).PubMed 

Google Scholar 
Suzuki, T. Production schedule of males and reproductive females, investment sex ratios, and worker-queen conflict in paper wasps. Am. Nat. 128, 366–378 (1986).Article 

Google Scholar 
Tsuchida, K. & Suzuki, T. Conflict over sex ratio and male production in paper wasps. Ann. Zool. Fenn. 43, 468–480 (2006).
Google Scholar 
Ohtsuki, H. & Tsuji, K. Adaptive reproduction schedule as a cause of worker policing in social Hymenoptera: A dynamic game analysis. Am. Nat. 173, 747–758 (2009).PubMed 
Article 

Google Scholar 
Walsh, P. S., Metzger, D. A. & Higuchi, R. Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. Biotechniques 10, 506–513 (1991).CAS 
PubMed 

Google Scholar 
Bassam, B. J., Caetano-Anolles, G. & Gresshoff, P. M. Fast and sensitive silver staining of DNA in polyacrylamide gels. Anal. Biochem. 196, 80–83 (1991).CAS 
PubMed 
Article 

Google Scholar 
van Oosterhout, C., Hutchinson, W. F., Wills, D. P. M. & Shipley, P. MICRO-CHECHER: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes. 4, 535–538 (2004).Article 
CAS 

Google Scholar 
Peakall, R. & Smouse, P. E. GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes. 6, 288–295 (2006).Article 

Google Scholar 
Raymond, M. & Rousset, F. GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. J. Hered. 86, 248–249 (1995).Article 

Google Scholar 
Meirmans, P. G. GenoDive version 3.0: Easy-to-use software for the analysis of genetic data of diploids and polyploids. Mol. Ecol. Res. 20, 1126–1131 (2020).CAS 
Article 

Google Scholar 
Michalakis, Y. & Excoffier, L. A generic estimation of population subdivision using distances between alleles with special reference for microsatellite loci. Genetics 142, 1061–1064 (1996).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Earl, D. & vonHoldt, B. STRUCTURE HARVESTER: A website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv. Genet. Res. 4, 359–361 (2012).Article 

Google Scholar 
Kopelman, N. M., Mayzel, J., Jakobsson, M., Rosenberg, N. A. & Mayrose, I. Clumpak: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Res. 15, 1179–1191 (2015).CAS 
Article 

Google Scholar 
Goodnight, K. F. Relatedness 4.2c Release (Rice University, 1996).
Google Scholar 
Queller, D. C. A method for detecting kin discrimination within natural colonies of social insects. Anim. Behav. 47, 569–576 (1994).Article 

Google Scholar  More

By pooling the results of the entire known range analysis of 14 GPS-collared elephants living in the Kinabatangan, our study suggests that this populations range covers at least 628 km2 (Table 3). Nine different locations were identified as hot spots, representing 266.9 km2 or 43% of this range, suggesting that just under half is highly used and/or frequented (Fig. 1). We found that the size of individual’s hot spots was positively related to the size of the entire range, meaning the larger the entire range the larger the summed area of an elephants hot spots. On average, hot spots represented a relatively small percent of an animal’s entire range (ranging from 4 to 20%, averaging 12%, Table 3). However, time spent within these hot spots ranged from 10 to 60% (averaging 34% across elephants, Table 5), with time spent in hot spots being related to the overall size of the hot spots (the larger the hot spot the more time elephants spent in them).Identifying the location of these hot spots is essential in designing appropriate management practices in collaboration with land users and identifying the best location for elephant corridors. In the last 25 years, forest cover in the Lower Kinabatangan has been drastically reduced and fragmented46, eroding the biodiversity value of this landscape. Today, this region has little remaining forests, and what is left is insufficient for sustaining the local elephant population10. Moreover, forests are highly fragmented along the Kinabatangan River, with a number of bottlenecks constraining elephant movements9. The situation in this landscape is getting worse because of further land clearances for agriculture, namely oil palm; as well as for the highly controversial Sukau Bridge and new road/highway that is planned for the region.Our analyses revealed a highly significant difference between the average proportions of protected area, unprotected forest, and oil palm estate extents within the elephant’s entire range; and a substantive, but not significant, difference across these land use/land cover types within hot spots (Table SI 4). At the individual level, there was a highly significant negative relationship between the proportion of protected areas and oil palm estates both within the elephant’s entire range and within the hot spots.At the pooled level, we found that around 45% of the entire known range and hot spots were within forested environments (280.44 km2 and 120.29 km2 respectively). Our results showed strong fidelity of certain elephants to these forested habitats. Our k-means cluster analysis found that within elephant entire ranges and hot spots, two out of the three cluster groups had high or very high usage of forests. Both cluster 1, for the entire range, and cluster 1 for hot spots extents, had five females that on average used forest environments 90% of their time, with protected areas being used 64% and 59%, and unprotected forested being used on average 26% and 31%, respectively (Table 7).Individuals in cluster 2, for the entire range analysis, on average, spent 73% of their time in forests (57% of this in protected areas and 16% in unprotected forests; Table 7). For the hot spot analysis, the individuals in cluster 2 spent on average 65% of their time in forests (52% of this in the unprotected forests and 13% in protected forests; Table 7). Elephants within these clusters were all females. Our results suggest that forest may be of particular importance for females as they had forest as their dominant land cover type within their entire range, hot spot extents and time spent analyses (Fig. 3, Table 5). Several studies have shown that adult females influence and guide the movement patterns and habitat utilization for their family group and that females in family units tend to inhabit less risky areas, such as within natural forest habitat60,61,62.However, the unprotected forest is at risk. We identified about 8% (or 49 km2) of forest identified within the pooled entire known range were not protected, with half potentially being on state land, and the remaining half on land titles of various types (Table SI 4). For the pooled hot spot areas, unprotected forest was proportionally higher, comprising of 11% (or 29 km2) of the total extent, with 54% being potentially on State land and 46% on land titles (Table SI 4). Protecting these forests would be an essential and efficient way to secure key elephant habitat since all collared individuals were using these forest fragments in their entire range (averaging 11%, and ranging from 8 to 18%), and hot spot extents (averaging 20%, and ranging from 0 to 53%) (Table SI 4, Fig. 3). On average, 24% of time was spent in unprotected forests within hot spots, though this varied widely from 0% (for the male elephant known as Gading) to 61% (for the female matriarch named Jasmine) (Table 5). In fact, five females had large proportions of their hot spot extents (24–53%) in unprotected forests, spending substantial periods of their time (33–61%) within these threatened areas.Our findings show that unprotected forests around the villages of Bilit and Sukau, were of particular significance (Figs. 1, 2). These unprotected forests largely consist of lowland dry forest, seasonally flooded swamp forest, and swamp forest, which are considered important habitats for elephants for feeding, resting and moving47,63. Within these forests, and along the forest margins and river banks there are also natural open grasslands that consist of Phragmites karka and Dinochloa scabrida that provide essential forage, mainly in the riparian areas for elephants9,21,23. Forested environments are also considered to be important in providing natural refugee from human activities and disturbance. For example, elephants have been documented to form significantly larger group sizes, as well as engaging in significantly more social interactions, in natural forest habitat compared to, for example, oil palm landscapes63. Adult females, generally, avoid areas considered unsafe for their respective social units, are more selective in the resources they use, and require regular access to water because of the presence of young64,65,66. This may be why our results, strongly suggest that forest habitats seem to be most important for adult females.Another significant issue faced by these elephants is the threat from the controversial planned Sukau bridge and road/highway that is set out in the Sabah Structure Plan, an overarching policy document for the State58. Currently, a new road/highway is under construction on the northern bank of the village of Sukau, and this has already cleared areas of unprotected forest. This public road could link to a potential new bridge that would cross over the Kinabatangan River, cutting through unprotected forest and a protected area (Lower Kinabatangan Wildlife Sanctuary), before going through oil palm estates then through another protected area to the south and through the Tabin elephant population range. For the Kinabatangan, creating a public highway will cut the elephant population range into two parts (Figs. 2, 3). All collared elephants use this area, as it is a key bottleneck and the only alternative option to pass around Sukau village9. We found that nine elephants have hot spots that intersect or meet up with the current road (which will be up-graded and get considerably busier) and/or the planned road/highway alignment (Figs. SI 1 and 2). For these elephants, we calculated that they spent from 2 to 44% (average 14%) of their time within these hot spots (Table 4). Our statistical analyses suggest that if the road/highway goes ahead it will have a significant impact on the elephants’ behaviour with respect to time spent in the hot spots. Indeed, this infrastructure project could have dire consequences for these elephants and their family groups, by disrupting their ranging patterns and segmenting the entire elephant range into two (Figs. 2, 4). The existing road in Batu Putih has already proven to be an impassable barrier for this elephant population, as no elephants have been observed crossing this road since the early 2000s14. For elephants that do try and cross, vehicle collisions may become a significant threat to elephants and drivers alike67, and potentially increasing human–elephant conflict in the nearby villages, as well as in plantations14,68,69, exacerbating an already difficult situation for this small and fragmented population.Results from the pooled analysis show that about 53% of the entire known population range is within oil palm estates; and 51% for the pooled hot spots (Fig. 3, Table SI 4). Our k-means clustering analysis grouped 6 elephants into cluster 3 that on average spent 57% of time in oil palm estates; and 7 elephants into cluster 2 within the hot spot analysis that on average spent 73% of their time in oil palm estates (Table 6). All the males, were clustered within these groups (Table 5). In fact, the three collared males were amongst the highest users of oil palm estates (Fig. 3, Table SI4, and 5). This could be related to a ‘‘high risk, high gain’’ strategy, often adopted by males to increase body size and enhance reproductive success32,33,60. However, it is interesting to see that three females (Ita, Ratu and Koyah) and their respective social units, also seemed to have high levels of oil palm use, while other individuals had zero or very little use of oil palm (e.g. Aqeela, Jasmin, Sandi, Kasih; Table SI 4, Fig. 3). Differential choices may result from differences in individual knowledge and experience with people during past encounters, for example70,71.We identified that collared elephants were ranging in 11 known oil palm estates, with the five most regularly used being Melangking Oil Palm Plantation (with 12 elephants entire range overlapping with this estate and six hot spots), IOI Corporation (with 11 overlapping entire ranges, and eight hot spots), Genting Plantations (14 and seven, respectively), Sime Darby Plantation (five and two, respectively), and Karangan Agriculture (8 and 2, respectively) (Table 6; Fig. 4). Presence of bottlenecks and barriers (e.g. electric fences) may explain hot spot occurrences in these estates, as well as feeding opportunities, management strategies of specific estates, and historical and seasonal ranges.Linear features like major highways, electric fences and drainage ditches hamper elephant movements within the Lower Kinabatangan9. A previous study identified 20 bottlenecks in the Lower Kinabatangan with the two main ones (of 9 km and 6.5 km in length) found around the village of Sukau9. In addition, the unplanned and chaotic erection of electric fences by large estates and smallholdings has disrupted significantly elephant movement patterns and resulted in artificial hot spots for certain individuals (e.g. Liun, Ita, Gading and Sejati)35,72. Electric fences have widely been used to mitigate human–elephant conflicts. The establishment of fences rarely consider the traditional elephant routes nor the location of existing fences in neighbouring estates. If elephants manage to enter such areas, they often become trapped and experience difficulties in returning to nearby forests, exacerbating conflicts with people35.Certain estates such as Melangking Oil Palm Plantation have allowed elephants to roam freely in their estate (Muhammad Al-Shafieq, personal communication). Since 2017, this plantation has shown a drastic reduction in damages to their oil palms following the removal of their permanent electric fences surrounding their entire estate. Instead, this plantation is using a temporary electric fencing regime around newly planted palm areas. Concurrently, they now do not push elephants out of their estate, which can explain why Melangking Oil Palm Plantation is a significant hotspot in the region.Another reason why elephant ranges incorporate oil palm estates is to move between forest patches that are becoming completely isolated following forest conversion, as is the case close to Sukau (Fig. SI1 and SI2; Fig. 1). Unlike other elephant species that increase their speed of movement rates in highly disturbed areas27,30,66, the Bornean elephant has been observed doing the opposite, which may explain some of the hot spots within oil palm estates. This movement strategy may allow for better vigilance as seen on a few occasions when elephants spent 2–5 days in the Bukit Melapi-Yu Kwang Corridor, near the village of Sukau, before leaving the area (Othman, personal observation).Hot spots in the oil palm landscape can also be explained by feeding opportunities, since elephants feed on palm shoots, leaves and hearts73. Elephants are known to eat the shoots of newly planted oil palms, often killing the palms and causing significant economic damages35. Since 2010, many estates located in the Lower Kinabatangan have started a new palm rotation. Palms are replanted every 25 years. A new rotation includes land clearing, bole and root mass removal, and the shredding or chipping of felled palms. Elephants are attracted to the shredded palm hearts since it gives them easy access to one of their favourite food72. This particular behaviour does not cause economic damage, and some estate managers allow the elephants to stay and forage in the chipping areas. This was documented for several collared elephants, whose hot spots and time spent were particularly high within oil palm (e.g. Gading and Sandy, two males; and Ratu and Ita, two females). Once the shredded palms have dried, however, elephants will leave these areas and move elsewhere. Within oil palm estates, some elephants have been found to travel more directly and rapidly suggesting ‘exploratory’ behaviour, which could be associated with searching for young palms or areas of palm felling and chipping of palm hearts15.Lastly, elephants may still be using their historical range that used to be covered with forest before conversion to oil palm. Other factors potentially explaining the relatively high use of oil palm estates include seasonal variations of ranging patterns. Indeed events of drought or floods limit the access to various parts of the floodplain and will tend to confine the animals in some areas9,63.In Sabah the state authorities have recorded at least 200 elephant deaths from the year 2010 to 2021 and most of these have occurred on, or near, oil palm estates14,74,75,76. Deaths from non-natural causes are largely due to poisoning (both accidental and intentional), gunshot wounds, poaching for tusks and other body parts, and snares35. Stopping killing and enabling a safe coexistence between people and elephants within multiple-use landscapes that are dominated by oil palm is one of the key strategies developed in the Bornean Elephant Action Plan for Sabah (2020–2029), which was endorsed by the State14. Based on our results in Lower Kinabatangan, a series of recommendations are proposed.This study underscores the importance of remaining forested areas for the Lower Kinabatangan elephant population. Full protection of all forest fragments left in the Lower Kinabatangan is urgently needed. Several official mechanisms are available to fulfil this request that has been proposed for the past 20 years by various organizations46.The current network of forests available in the Lower Kinabatangan is too small and fragmented to sustain a viable elephant population. Forest corridors must be created across the landscape through reforestation exercises, whilst concurrently undertaking enrichment planting of native understory forage within forested areas as this may minimize the need for elephants to search for easily accessible food in high-risk oil palm landscapes21,22,23.Current governmental plans to build a road bridge and public road/highway linking the southern bank of the Kinabatangan River to Tabin Wildlife Reserve to the south will irreversibly impact the Lower Kinabatangan elephant population by cutting the current range into two isolated parts. This will impact the elephants ranging patterns, potentially even fragmenting the already small population into two groups, and potentially leading to elephant deaths by vehicle collisions (which is becoming increasingly common in Peninsular Malaysia), and increase the risk of poaching activities, all resulting in a decrease in the genetic diversity of the, already small and isolated, population14,67.Eventually, the future of the Kinabatangan elephant population resides in improving land use and management practices within oil palm estates currently used by elephants. We recommend that priority should be given at improving elephant movements in oil palm estates by removing unnecessary man-made barriers and only cautiously installing temporary electric fences to protect sensitive areas. For example, the use of electric fences around mature oil palm and areas whereby palms are being removed and chipped could be prohibited, and electric fences permitted solely for protecting oil palm nurseries, new plantings and young oil palms (e.g. up to 7–8 years old), and staff and office quarters. This would greatly allow for landscape permeability for elephants, and other species that need to cross the landscape for their ecological and biological needs14.A handful of guidelines exist to assist oil palm managers and staff in managing elephant populations in their respective estates72,77. However, there is a need for a more comprehensive set of guidelines, which delineate better practices with the aim to increase the protection of people and elephants outside protected areas. Guidelines should specify “do’s” and “don’ts” (based on best available data and knowledge) of actions needed before, during and after elephants visit oil palm estates and smallholdings.Sabah now is in an interesting transition within their palm oil sector. On the 21st October 2015, the Sabah State Cabinet committed to produce 100% certified sustainable palm oil, by 2025, under the Roundtable for Sustainable palm Oil (RSPO) Jurisdictional Certification approach. Under this approach, areas of High Conservation Value and areas identified within the High Carbon Stock Approach need specific management and monitoring, in order to comply with RSPO principles and criteria78,79,80. Sabah government can use this platform to build an integrated landscape level approach to better manage landscapes within known elephant ranges (which is considered a High Conservation Value species) to allow for a safe and permeable movement through the landscape.Eventually, long-term survival of the Bornean elephant will mainly depend on how people and elephants can co-exist. It is our hope that this study illustrates the importance of protecting all forested habitat and effectively managing areas outside of protected areas to allow for long-term elephant coexistence with humans in this landscape. More