Ecology

Worden AZ, Follows MJ, Giovannoni SJ, Wilken S, Zimmerman AE, Keeling PJ. Rethinking the marine carbon cycle: Factoring in the multifarious lifestyles of microbes. Science. 2015;347:1257594–1257594.PubMed 
Article 
CAS 

Google Scholar 
Guillou L, Viprey M, Chambouvet A, Welsh RM, Kirkham AR, Massana R, et al. Widespread occurrence and genetic diversity of marine parasitoids belonging to Syndiniales (Alveolata). Environ Microbiol. 2008;10:3349–65.CAS 
PubMed 
Article 

Google Scholar 
Jephcott TG, Alves-de-Souza C, Gleason FH, van Ogtrop FF, Sime-Ngando T, Karpov SA, et al. Ecological impacts of parasitic chytrids, syndiniales and perkinsids on populations of marine photosynthetic dinoflagellates. Fungal Ecol. 2016;19:47–58.Article 

Google Scholar 
Alacid E, Reñé A, Garcés E. New Insights into the Parasitoid Parvilucifera sinerae Life Cycle: The Development and Kinetics of Infection of a Bloom-forming Dinoflagellate Host. Protist. 2015;166:677–99.PubMed 
Article 

Google Scholar 
Not F, Gausling R, Azam F, Heidelberg JF, Worden AZ. Vertical distribution of picoeukaryotic diversity in the Sargasso Sea. Environ Microbiol. 2007;9:1233–52.CAS 
PubMed 
Article 

Google Scholar 
Lima-Mendez G, Faust K, Henry N, Decelle J, Colin S, Carcillo F, et al. Determinants of community structure in the global plankton interactome. Science 2015;348:1262073.de Vargas C, Audic S, Henry N, Decelle J, Mahe F, Logares R, et al. Eukaryotic plankton diversity in the sunlit ocean. Science. 2015;348:1261605.PubMed 
Article 
CAS 

Google Scholar 
Siano R, Alves-De-Souza C, Foulon E, Bendif M, Simon E, Guillou NL. et al. Distribution and host diversity of Amoebophryidae parasites across oligotrophic waters of the Mediterranean Sea. Biogeosciences. 2011;8:267–78.Article 

Google Scholar 
Coats DW. Parasitic life styles of marine dinoflagellates. J Eukaryot Microbiol. 1999;46:402–9.Article 

Google Scholar 
Farhat S, Le P, Kayal E, Noel B, Bigeard E, Corre E, et al. Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp. BMC Biol. 2021;19:1–21.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gornik SG, Febrimarsa, Cassin AM, MacRae JI, Ramaprasad A, Rchiad Z, et al. Endosymbiosis undone by stepwise elimination of the plastid in a parasitic dinoflagellate. Proc Natl Acad Sci USA. 2015;112:5767–72.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
John U, Lu Y, Wohlrab S, Groth M, Janouškovec J, Kohli GS, et al. An aerobic eukaryotic parasite with functional mitochondria that likely lacks a mitochondrial genome. Sci Adv. 2019;5:1–12.Article 
CAS 

Google Scholar 
McFadden GI, Reith ME, Munholland J, Lang-Unnasch N. Plastid in human parasites. Nature. 1996;381:482.CAS 
PubMed 
Article 

Google Scholar 
Sibbald SJ, Archibald JM. Genomic Insights into Plastid Evolution. Genome Biol Evol. 2020;12:978–90.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Not F, Probert I, Ribeiro CG, Crenn K, Guillou L, Jeanthon C, et al. Photosymbiosis in Marine Pelagic Environments. In: M. S. Cretoiu (ed). The Marine Microbiome. 2016. New York, NY: Springer International Publishing, pp 305–32.Cai R, Kayal E, Alves-de-Souza C, Bigeard E, Corre E, Jeanthon C, et al. Cryptic species in the parasitic Amoebophrya species complex revealed by a polyphasic approach. Sci Rep. 2020;10:1–11.Article 
CAS 

Google Scholar 
Coats DW, Park MG. Parasitism of photosynthetic dinoflagellates by three strains of Amoebophrya (Dinophyta): Parasite survival, infectivity, generation time, and host specificity. J Phycol. 2002;38:520–8.Article 

Google Scholar 
Kayal E, Alves-de-Souza C, Farhat S, Velo-Suarez L, Monjol J, Szymczak J, et al. Dinoflagellate Host Chloroplasts and Mitochondria Remain Functional During Amoebophrya Infection. Front Microbiol. 2020;11:1–11.Article 

Google Scholar 
Miller JJ, Delwiche CF, Coats DW. Ultrastructure of Amoebophrya sp. and its Changes during the Course of Infection. Protist. 2012;163:720–45.PubMed 
Article 

Google Scholar 
Decelle J, Stryhanyuk H, Gallet B, Veronesi G, Schmidt M, Balzano S, et al. Algal Remodeling in a Ubiquitous Planktonic Photosymbiosis. Curr Biol. 2019;29:968–978.e4.CAS 
PubMed 
Article 

Google Scholar 
Uwizeye C, Mars Brisbin M, Gallet B, Chevalier F, LeKieffre C, Schieber NL, et al. Cytoklepty in the plankton: A host strategy to optimize the bioenergetic machinery of endosymbiotic algae. Proc Natl Acad Sci USA 2021;118.Uwizeye C, Decelle J, Jouneau P, Flori S, Gallet B, Keck J, et al. Morphological bases of phytoplankton energy management and physiological responses unveiled by 3D subcellular imaging. Nat Commun. 2021;12:1049.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lowe DG. Distinctive Image Features from Scale-Invariant Keypoints. Int J Comput Vis. 2004;60:91–110.Article 

Google Scholar 
Hennies J, Lleti JMS, Schieber NL, Templin RM, Steyer AM, Schwab Y. AMST: alignment to median smoothed template for focused ion beam scanning electron microscopy image stacks. Sci Rep. 2020;10:1–10.Article 
CAS 

Google Scholar 
Kikinis R, Pieper SD, Vosburgh KG. 3D Slicer: A Platform for Subject-Specific Image Analysis, Visualization, and Clinical Support. Intraoperative Imaging and Image-Guided Therapy. 2014. Springer New York, New York, NY, pp 277–89.Polerecky L, Adam B, Milucka J, Musat N, Vagner T, Kuypers MMM. Look@NanoSIMS – a tool for the analysis of nanoSIMS data in environmental microbiology. Environ Microbiol. 2012;14:1009–23.CAS 
PubMed 
Article 

Google Scholar 
Jia B, Zhu XF, Pu ZJ, Duan YX, Hao LJ, Zhang J, et al. Integrative view of the diversity and evolution of SWEET and semiSWEET sugar transporters. Front Plant Sci. 2017;8:1–18.
Google Scholar 
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20:1160–6.CAS 
PubMed 
Article 

Google Scholar 
Castresana J. Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis. Mol Biol Evol. 2000;17:540–52.CAS 
PubMed 
Article 

Google Scholar 
Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–3.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512.CAS 
PubMed 
Article 

Google Scholar 
Farhat S, Florent I, Noel B, Kayal E, Da Silva C, Bigeard E, et al. Comparative time-scale gene expression analysis highlights the infection processes of two amoebophrya strains. Front Microbiol. 2018;9:1–19.Article 

Google Scholar 
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinforma. 2011;12:323.CAS 
Article 

Google Scholar 
Cachon J. Contribution à l’étude des péridiniens parasites. Cytologie, cycles évolutifs Ann Sci Nat Zool. 1964;6:1–158.
Google Scholar 
Van Dooren GG, Marti M, Tonkin CJ, Stimmler LM, Cowman AF, McFadden GI. Development of the endoplasmic reticulum, mitochondrion and apicoplast during the asexual life cycle of Plasmodium falciparum. Mol Microbiol. 2005;57:405–19.PubMed 
Article 
CAS 

Google Scholar 
Tyler KM, Matthews KR, Gull K. Anisomorphic cell division by African trypanosomes. Protist. 2001;152:367–78.CAS 
PubMed 
Article 

Google Scholar 
Jakob M, Hoffmann A, Amodeo S, Peitsch C, Zuber B, Ochsenreiter T. Mitochondrial growth during the cell cycle of Trypanosoma brucei bloodstream forms. Sci Rep. 2016;6:1–13.Article 
CAS 

Google Scholar 
Hughes L, Borrett S, Towers K, Starborg T, Vaughan S. Patterns of organelle ontogeny through a cell cycle revealed by whole-cell reconstructions using 3D electron microscopy. J Cell Sci. 2017;130:637–47.CAS 
PubMed 

Google Scholar 
Ovciarikova J, Lemgruber L, Stilger KL, Sullivan WJ, Sheiner L. Mitochondrial behaviour throughout the lytic cycle of Toxoplasma gondii. Sci Rep. 2017;7:1–13.Article 
CAS 

Google Scholar 
Nishi M, Hu K, Murray JM, Roos DS. Organellar dynamics during the cell cycle of Toxoplasma gondii. J Cell Sci. 2008;121:1559–68.CAS 
PubMed 
Article 

Google Scholar 
Long M, Marie D, Szymczak J, Toullec J, Bigeard E, Sourisseau M, et al. Dinophyceae can use exudates as weapons against the parasite Amoebophrya sp. (Syndiniales). ISME Commun. 2021;1:34.Article 

Google Scholar 
Harris E, Yoshida K, Cardelli J, Bush J. Rab11-like GTPase associates with and regulates the structure and function of the contractile vacuole system in. Dictyostelium J Cell Sci. 2001;114:3035–45.CAS 
PubMed 
Article 

Google Scholar 
Marshansky V, Rubinstein JL, Grüber G. Eukaryotic V-ATPase: Novel structural findings and functional insights. Biochim Biophys Acta – Bioenerg. 2014;1837:857–79.CAS 
Article 

Google Scholar 
Cox D, Lee DJ, Dale BM, Calafat J, Greenberg S. A Rab11-containing rapidly recycling compartment in macrophages that promotes phagocytosis. Proc Natl Acad Sci USA. 2000;97:680–5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Vines JH, King JS. The endocytic pathways of Dictyostelium discoideum. Int J Dev Biol. 2019;63:461–71.CAS 
PubMed 
Article 

Google Scholar 
Decelle J, Veronesi G, LeKieffre C, Gallet B, Chevalier F, Stryhanyuk H, et al. Subcellular architecture and metabolic connection in the planktonic photosymbiosis between Collodaria (radiolarians) and their microalgae. Environ Microbiol. 2021;23:6569–86.CAS 
PubMed 
Article 

Google Scholar 
Sigee DC, Kearns LP. Levels of dinoflagellate chromosome-bound metals in conditions of low external ion availability: An X-ray microanalytical study. Tissue Cell. 1981;13:441–51.CAS 
PubMed 
Article 

Google Scholar 
Pinchuk GE, Ammons C, Culley DE, Li SMW, McLean JS, Romine MF, et al. Utilization of DNA as a sole source of phosphorus, carbon, and energy by Shewanella spp.: Ecological and physiological implications for dissimilatory metal reduction. Appl Environ Microbiol. 2008;74:1198–208.CAS 
PubMed 
Article 

Google Scholar 
Caffaro CE, Boothroyd JC. Evidence for Host Cells as the Major Contributor of Lipids in the Intravacuolar Network of Toxoplasma-Infected Cells. Eukaryot Cell. 2011;10:1095–9.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lopez J, Bittame A, Massera C, Vasseur V, Effantin G, Valat A, et al. Intravacuolar membranes regulate CD8 T cell recognition of membrane-bound Toxoplasma gondii protective antigen. Cell Rep. 2015;13:2273–86.CAS 
PubMed 
Article 

Google Scholar 
Pszenny V, Ehrenman K, Romano JD, Kennard A, Schultz A, Roos DS, et al. A Lipolytic Lecithin:Cholesterol Acyltransferase Secreted by Toxoplasma Facilitates Parasite Replication and Egress. J Biol Chem. 2016;291:3725–46.CAS 
PubMed 
Article 

Google Scholar 
Nolan SJ, Romano JD, Coppens I. Host lipid droplets: An important source of lipids salvaged by the intracellular parasite Toxoplasma gondii. PLoS Pathog. 2017;13:e1006362.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Fox BA, Guevara RB, Rommereim LM, Falla A, Bellini V, Pètre G, et al. Toxoplasma gondii parasitophorous vacuole membrane-associated dense granule proteins orchestrate chronic infection and GRA12 underpins resistance to host gamma interferon. MBio. 2019; 10:e00589-19.Freeman Rosenzweig ES, Xu B, Kuhn Cuellar L, Martinez-Sanchez A, Schaffer M, Strauss M, et al. The eukaryotic CO2-concentrating organelle is liquid-like and exhibits dynamic reorganization. Cell. 2017;171:148–162.e19.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Zeeman SC, Kossmann J, Smith AM. Starch: Its metabolism, evolution, and biotechnological modification in plants. Annu Rev Plant Biol. 2010;61:209–34.CAS 
PubMed 
Article 

Google Scholar 
Qureshi AA, Suades A, Matsuoka R, Brock J, McComas SE, Nji E, et al. The molecular basis for sugar import in malaria parasites. Nature. 2020;578:321–5.CAS 
PubMed 
Article 

Google Scholar 
Blume M, Rodriguez-Contreras D, Landfear S, Fleige T, Soldati-Favre D, Lucius R, et al. Host-derived glucose and its transporter in the obligate intracellular pathogen Toxoplasma gondii are dispensable by glutaminolysis. Proc Natl Acad Sci USA. 2009;106:12998–3003.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Chen LQ, Qu XQ, Hou BH, Sosso D, Osorio S, Fernie AR, et al. Sucrose efflux mediated by SWEET proteins as a key step for phloem transport. Science (80-). 2012;335:207–11.CAS 
Article 

Google Scholar 
Latorraca NR, Fastman NM, Venkatakrishnan AJ, Frommer WB, Dror RO, Feng L. Mechanism of substrate translocation in an alternating access transporter. Cell. 2017;169:96–107.e12.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Amiar S, Katris NJ, Berry L, Dass S, Duley S, Arnold CS, et al. Division and adaptation to host environment of apicomplexan parasites depend on apicoplast lipid metabolic plasticity and host organelle remodeling. Cell Rep. 2020;30:3778–3792.e9.CAS 
PubMed 
Article 

Google Scholar 
Hu X, Binns D, Reese ML. The coccidian parasites Toxoplasma and Neospora dysregulate mammalian lipid droplet biogenesis. J Biol Chem. 2017;292:11009–20.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gomes AF, Magalhães KG, Rodrigues RM, de Carvalho L, Molinaro R, Bozza PT, et al. Toxoplasma gondii-skeletal muscle cells interaction increases lipid droplet biogenesis and positively modulates the production of IL-12, IFN-g and PGE2. Parasit Vectors. 2014;7:47.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Jacot D, Waller RF, Soldati-Favre D, MacPherson DA, MacRae JI. Apicomplexan energy metabolism: carbon source promiscuity and the quiescence hyperbole. Trends Parasitol. 2016;32:56–70.CAS 
PubMed 
Article 

Google Scholar 
Salcedo-Sora JE, Caamano-Gutierrez E, Ward SA, Biagini GA. The proliferating cell hypothesis: a metabolic framework for Plasmodium growth and development. Trends Parasitol. 2014;30:170–5.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Muñoz-Gómez SA, Wideman JG, Roger AJ, Slamovits CH, Agashe D. The origin of mitochondrial cristae from alphaproteobacteria. Mol Biol Evol. 2017;34:943–56.PubMed 

Google Scholar 
Evers F, Cabrera-Orefice A, Elurbe DM, Kea-te Lindert M, Boltryk SD, Voss TS, et al. Composition and stage dynamics of mitochondrial complexes in Plasmodium falciparum. Nat Commun. 2021;12:1–17.Article 
CAS 

Google Scholar 
Krungkrai J, Prapunwattana P, Krungkrai SR. Ultrastructure and function of mitochondria in gametocytic stage of Plasmodium falciparum. Parasite. 2000;7:19–26.CAS 
PubMed 
Article 

Google Scholar 
Krungkrai J. The multiple roles of the mitochondrion of the malarial parasite. Parasitology 2004;129:511–524. https://doi.org/10.1017/S0031182004005888.Lee JW. Protonic capacitor: elucidating the biological significance of mitochondrial cristae formation. Sci Rep. 2020;10:1–14.Article 
CAS 

Google Scholar 
Stephan T, Brüser C, Deckers M, Steyer AM, Balzarotti F, Barbot M, et al. MICOS assembly controls mitochondrial inner membrane remodeling and crista junction redistribution to mediate cristae formation. EMBO J. 2020;39:1–24.Article 
CAS 

Google Scholar 
Pánek T, Eliáš M, Vancová M, Lukeš J, Hashimi H. Returning to the Fold for Lessons in Mitochondrial Crista Diversity and Evolution. Curr Biol. 2020;30:R575–R588.PubMed 
Article 
CAS 

Google Scholar 
Wideman JG, Muñoz-Gómez SA. The evolution of ERMIONE in mitochondrial biogenesis and lipid homeostasis: An evolutionary view from comparative cell biology. Biochim Biophys Acta – Mol Cell Biol Lipids. 2016. Elsevier B.V., 1861: 900-12.Mühleip A, Kock Flygaard R, Ovciarikova J, Lacombe A, Fernandes P, Sheiner L, et al. ATP synthase hexamer assemblies shape cristae of Toxoplasma mitochondria. Nat Commun. 2021;12:120.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Flegontov P, Michálek J, Janouškovec J, Lai DH, Jirků M, Hajdušková E, et al. Divergent mitochondrial respiratory chains in phototrophic relatives of apicomplexan parasites. Mol Biol Evol. 2015;32:1115–31.CAS 
PubMed 
Article 

Google Scholar 
Painter HJ, Morrisey JM, Mather MW, Vaidya AB. Specific role of mitochondrial electron transport in blood-stage Plasmodium falciparum. Nature. 2007;446:88–91.CAS 
PubMed 
Article 

Google Scholar 
Nishida T, Hatama S, Ishikawa Y, Kadota K. Intranuclear coccidiosis in a calf. J Vet Med Sci. 2009;71:1109–13.PubMed 
Article 

Google Scholar 
Pecka Z. Life cycle and ultrastructure of Eimeria stigmosa, the intranuclear coccidian of the goose (Anser anser domesticus). Folia Parasitol (Praha). 1992;39:105–14.CAS 

Google Scholar 
Voleman L, Doležal P. Mitochondrial dynamics in parasitic protists. PLoS Pathog. 2019;15:e1008008.Bílý T, Sheikh S, Mallet A, Bastin P, Pérez‐Morga D, Lukeš J, et al. Ultrastructural changes of the mitochondrion during the life cycle of Trypanosoma brucei. J Eukaryot Microbiol. 2021;68:e12846.Elliott DA, McIntosh MT, Hosgood HD, Chen S, Zhang G, Baevova P, et al. Four distinct pathways of hemoglobin uptake in the malaria parasite Plasmodium falciparum. Proc Natl Acad Sci USA. 2008;105:2463–8.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

Ottenburghs, J. et al. A history of hybrids? Genomic patterns of introgression in the True Geese. BMC Evol. Biol. 17, 14 (2017).Article 

Google Scholar 
Baiz, M. D., Tucker, P. K. & Cortés-Ortiz, L. Multiple forms of selection shape reproductive isolation in a primate hybrid zone. Mol. Ecol. 28, 1056–1069 (2019).CAS 
PubMed 
Article 

Google Scholar 
Slager, D. L. et al. Cryptic and extensive hybridization between ancient lineages of American crows. Mol. Ecol. 29, 956–969 (2020).CAS 
PubMed 
Article 

Google Scholar 
Grant, P. R. & Grant, B. R. Introgressive hybridization and natural selection in Darwin’s finches. Biol. J. Linnean Soc. 117, 812–822 (2016).Article 

Google Scholar 
Pauquet, G., Salzburger, W. & Egger, B. The puzzling phylogeography of the haplochromine cichlid fish Astatotilapia burtoni. Ecol. Evol. 8, 5637–5648 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
Schluter, D. Ecological character displacement in adaptive radiation. Am. Nat. 156, S4–S16 (2000).Article 

Google Scholar 
Song, Y. et al. Adaptive introgression of anticoagulant rodent poison resistance by hybridization between old world mice. Curr. Biol. 21, 1296–1301 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Anderson, R. P., Peterson, A. T. & Gómez-Laverde, M. Using niche-based GIS modeling to test geographic predictions of competitive exclusion and competitive release in South American pocket mice. Oikos 98, 3–16 (2002).Article 

Google Scholar 
Gramlich, S., Wagner, N. D. & Horandl, E. RAD-seq reveals genetic structure of the F-2-generation of natural willow hybrids (Salix L.) and a great potential for interspecific introgression. BMC Plant Biol. 18, 12 (2018).Article 

Google Scholar 
Mavárez, J. et al. Speciation by hybridization in Heliconius butterflies. Nature 441, 868–871 (2006).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Cahill, J. A. et al. Genomic evidence of geographically widespread effect of gene flow from polar bears into brown bears. Mol. Ecol. 24, 1205–1217 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Djogbénou, L. et al. Evidence of introgression of the ace-1(R) mutation and of the ace-1 duplication in West African Anopheles gambiae s. s. PLoS ONE 3, e2172 (2008).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Enciso-Romero, J. et al. Evolution of novel mimicry rings facilitated by adaptive introgression in tropical butterflies. Mol. Ecol. 26, 5160–5172 (2017).PubMed 
Article 

Google Scholar 
Dasmahapatra, K. K. et al. Butterfly genome reveals promiscuous exchange of mimicry adaptations among species. Nature 487, 94–98 (2012).ADS 
CAS 
PubMed Central 
Article 

Google Scholar 
Latch, E. K., Harveson, L. A., King, J. S., Hobson, M. D. & Rhodes, J. R. Assessing hybridization in wildlife populations using molecular markers: a case study in wild turkeys. J. Wildl. Manag. 70, 485–492 (2006).Article 

Google Scholar 
Oliveira, R., Godinho, R., Randi, E. & Alves, P. C. Hybridization versus conservation: are domestic cats threatening the genetic integrity of wildcats (Felis silvestris silvestris) in Iberian Peninsula?. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 363, 2953–2961 (2008).Article 

Google Scholar 
Nichols, P. et al. Secondary contact seeds phenotypic novelty in cichlid fishes. Proc. R. Soc. B Biol. Sci. 282, 8 (2015).
Google Scholar 
Yang, W. Z. et al. Genomic evidence for asymmetric introgression by sexual selection in the common wall lizard. Mol. Ecol. 27, 4213–4224 (2018).CAS 
PubMed 
Article 

Google Scholar 
Boratyński, Z. et al. Introgression of mitochondrial DNA among Myodes voles: consequences for energetics?. BMC Evol. Biol. 11, 355 (2011).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Mondal, M. et al. Genomic analysis of Andamanese provides insights into ancient human migration into Asia and adaptation. Nat. Genet. 48, 1066–1070 (2016).CAS 
PubMed 
Article 

Google Scholar 
Melo-Ferreira, J., Seixas, F. A., Cheng, E., Mills, L. S. & Alves, P. C. The hidden history of the snowshoe hare, Lepus americanus: extensive mitochondrial DNA introgression inferred from multilocus genetic variation. Mol. Ecol. 23, 4617–4630 (2014).CAS 
PubMed 
Article 

Google Scholar 
Mims, M. C., Hulsey, C. D., Fitzpatrick, B. M. & Streelman, J. T. Geography disentangles introgression from ancestral polymorphism in Lake Malawi cichlids. Mol. Ecol. 19, 940–951 (2010).PubMed 
Article 

Google Scholar 
Salazar, C. et al. Genetic evidence for hybrid trait speciation in heliconius butterflies. PLoS Genet. 6, e1000930 (2010).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Naisbit, R. E., Jiggins, C. D. & Mallet, J. Mimicry: developmental genes that contribute to speciation. Evol. Dev. 5, 269–280 (2003).CAS 
PubMed 
Article 

Google Scholar 
Zhang, W., Dasmahapatra, K. K., Mallet, J., Moreira, G. R. P. & Kronforst, M. R. Genome-wide introgression among distantly related Heliconius butterfly species. Genome Biol. 17, 15 (2016).Article 
CAS 

Google Scholar 
Zhang, W., Kunte, K. & Kronforst, M. R. Genome-wide characterization of adaptation and speciation in tiger swallowtail butterflies using De Novo transcriptome assemblies. Genome Biol. Evol. 5, 1233–1245 (2013).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Jones, M. R. et al. Adaptive introgression underlies polymorphic seasonal camouflage in snowshoe hares. Science (New York, NY). 360, 1355–1358 (2018).ADS 
CAS 
Article 

Google Scholar 
Melville, J. Competition and character displacement in two species of scincid lizards. Ecol. Lett. 5, 386–393 (2002).Article 

Google Scholar 
Pfennig, D. W. & Pfennig, K. S. Character displacement and the origins of diversity. Am. Nat. 176, S26–S44 (2010).PubMed 
PubMed Central 
MATH 
Article 

Google Scholar 
Kooyers, N. J., James, B. & Blackman, B. K. Competition drives trait evolution and character displacement between Mimulus species along an environmental gradient. Evol. Int. J. Org. Evol. 71, 1205–1221 (2017).CAS 
Article 

Google Scholar 
Adams, D. C. & Rohlf, F. J. Ecological character displacement in Plethodon: Biomechanical differences found from a geometric morphometric study. Proc. Natl. Acad. Sci. 97, 4106–4111 (2000).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Grant, P. R. & Grant, B. R. Evolution of character displacement in Darwin’s finches. Science (New York, NY) 313, 224–226 (2006).ADS 
CAS 
Article 

Google Scholar 
Pfennig, D. W. & Murphy, P. J. Character displacement in polyphenic tadpoles. Evol. Int. J. Org. Evol. 54, 1738–1749 (2000).CAS 
Article 

Google Scholar 
Jones, G. Acoustic signals and speciation: the roles of natural and sexual selection in the evolution of cryptic species. Adv. Study Behav. 26, 317–354 (1997).Article 

Google Scholar 
Marsteller, S., Adams, D. C., Collyer, M. L. & Condon, M. Six cryptic species on a single species of host plant: morphometric evidence for possible reproductive character displacement. Ecol. Entomol. 34, 66–73 (2009).Article 

Google Scholar 
Tene Fossog, B. et al. Habitat segregation and ecological character displacement in cryptic African malaria mosquitoes. Evol. Appl. 8, 326–345 (2015).PubMed 
PubMed Central 
Article 

Google Scholar 
Ibáñez, C., García-Mudarra, J. L., Ruedi, M., Stadelmann, B. & Juste, J. The Iberian contribution to cryptic diversity in European bats. Acta Chiropterol. 8, 277–297 (2006).Article 

Google Scholar 
Juste, J. et al. Mitochondrial phylogeography of the long-eared bats (Plecotus) in the Mediterranean Palaearctic and Atlantic Islands. Mol. Phylogenet. Evol. 31, 1114–1126 (2004).CAS 
PubMed 
Article 

Google Scholar 
Schreber J. Die Säugthiere in Abbildungen nach der Natur, mit Beschreibungen. Erlangen – Expedition des Schreber’schen säugthier- und des Esper’schen Schmetterlingswerkes. Ernst Mayr Library of the MCZ, 1774–1855 (Harvard, 1774).Temminck, C. J. Monographies de Mammologie, ou description de quelques genres de Mammifères, dont les espèces ont été observes dans les différens Musées de l’Europe, Vol. 2, No. 302, 26–70 (G. Dufour et Ed. D’Ocagne, 1840).Centeno-Cuadros, A. et al. Comparative phylogeography and asymmetric hybridization between cryptic bat species. J. Zool. Syst. Evol. Res. 57, 1004–1018 (2019).Article 

Google Scholar 
Santos, H. et al. Shaping of bat cryptic distribution in Iberia. Biol. J. Linnean Soc. 112, 150–162 (2014).Article 

Google Scholar 
Novella-Fernandez, R. et al. Broad-scale patterns of geographic avoidance between species emerge in the absence of fine-scale mechanisms of coexistence. Divers. Distrib. 27, 1606–1618 (2021).Article 

Google Scholar 
Neubaum, M. A., Douglas, M. R., Douglas, M. E. & O’Shea, T. J. Molecular ecology of the big brown bat (Eptesicus fuscus): genetic and natural history variation in a hybrid zone. J. Mammal. 88, 1230–1238 (2007).Article 

Google Scholar 
Worthington-Wilmer, J. & Barratt, E. A non-lethal method of tissue sampling for genetic studies of chiropterans. Bat Res. News 37(1), 1–4 (1996).
Google Scholar 
Illumination, I.C.o. A colour appearance model for colour management systems: CIECAM02. Technical Report No CIE 159, 2004 (2004).Maroco, J. Análise estatística com utilização do SPSS. 3ª edição. Edições Silabo (2010).Wickham, H. et al. Welcome to the tidyverse. J. Open Source Softw. 4(43), 1686. https://doi.org/10.21105/joss.01686 (2019).ADS 
Article 

Google Scholar 
Karatzoglou, A., Smola, A., Hornik, K. & Zeileis, A. Kernlab: an S4 package for kernel methods in R. J. Stat. Softw. 11(9), 1–20 (2004).Article 

Google Scholar 
Meyer, D., Dimitriadou, E., Hornik, K., Weingessel, A. & Leisch, F. e1071: Misc Functions of the Department of Statistics, Probability Theory Group (Formerly: E1071), TU Wien. https://cran.r-project.org/web/packages/e1071/index.html (2012).Redgwell, R. D., Szewczak, J. M., Jones, G. & Parsons, S. Classification of echolocation calls from 14 species of bat by support vector machines and ensembles of neural networks. Algorithms 2, 907–924 (2009).Article 

Google Scholar 
Ochoa-López, S. et al. Ontogenetic changes in the targets of natural selection in three plant defenses. New Phytol. 226, 1480–1491 (2020).PubMed 
Article 
CAS 

Google Scholar 
Grant, P. R. & Grant, B. R. Phenotypic and genetics effects of hybridization in Darwin’s finches. Evol. Int. J. Org. Evol. 48, 297–316 (1994).Article 

Google Scholar 
Abzhanov, A., Protas, M., Grant, B. R., Grant, P. R. & Tabin, C. J. Bmp4 and morphological variation of beaks in Darwin’s finches. Science (New York, NY). 305, 1462–1465 (2004).ADS 
CAS 
Article 

Google Scholar 
von Holdt, B. M., Kays, R., Pollinger, J. P. & Wayne, R. K. Admixture mapping identifies introgressed genomic regions in North American canids. Mol. Ecol. 25, 2443–2453 (2016).Article 

Google Scholar 
Santana, S. E., Strait, S. & Dumont, E. R. The better to eat you with: functional correlates of tooth structure in bats. Funct. Ecol. 25, 839–847 (2011).Article 

Google Scholar 
Kalcounis, M. C. & Brigham, R. M. Intraspecific variation in wing loading affects habitat use by little brown bats (Myotis lucifugus). Can. J. Zool. 73, 89–95 (1995).Article 

Google Scholar 
Muijres, F. T., Johansson, L. C., Winter, Y. & Anders, H. Comparative aerodynamic performance of flapping flight in two bat species using time-resolved wake visualization. J. R. Soc. Interface 8, 1418–1428 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
Bradley, B. J. & Mundy, N. I. The primate palette: the evolution of primate coloration. Evol. Anthropolo. Issues News Rev. 17, 97–111 (2008).Article 

Google Scholar 
Müller, B. & Peichl, L. Retinal cone photoreceptors in microchiropteran bats. Investig. Ophthalmol. Vis. Sci. 46, 2259–2259 (2005).
Google Scholar 
Winter, Y., López, J. & von Helversen, O. Ultraviolet vision in a bat. Nature 425, 612–614 (2003).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Caro, T. The adaptive significance of coloration in mammals. Bioscience 55, 125–136 (2005).Article 

Google Scholar 
Chaverri, G., Ancillotto, L. & Russo, D. Social communication in bats. Biol. Rev. 93, 1938–1954 (2018).PubMed 
Article 

Google Scholar 
Dietz, C., Von Helversen, O. & Nill, D. Bats of Britain, Europe and Northwest Africa 320–333 (A&C Black Publishers Ltd, 2009).
Google Scholar 
Martinoli, A., Mazzamuto, M.V. & Spada, M. Serotine Eptesicus serotinus (Schreber, 1774). In Handbook of the Mammals of Europe, 1–17 (2020).Dinger, G. Winternachweise von Breitflügelfledermaus (Eptesicus serotinus) in Kirchen. Nyctalus (N.F.) 7, 614–616 (1991).
Google Scholar 
Kowalski, K. & Rzebik-Kowalska, B. Mammals of algeria (1991).Novella-Fernandez, R. et al. Trophic resource partitioning drives fine-scale coexistence in cryptic bat species. Ecol. Evol. 10(24), 14122–14136 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Galván, I., Vargas-Mena, J. C. & Rodríguez-Herrera, B. Tent-roosting may have driven the evolution of yellow skin coloration in Stenodermatinae bats. J. Zool. Syst. Evol. Res. 58, 519–527 (2020).Article 

Google Scholar 
Wang, Z. L., Zhang, D. Y. & Wang, G. Does spatial structure facilitate coexistence of identical competitors. Ecol. Model. 181, 17–23 (2005).Article 

Google Scholar 
Anderson, T. M. et al. Molecular and evolutionary history of melanism in North American gray wolves. Science (New York, NY). 323, 1339–1343 (2009).ADS 
CAS 
PubMed Central 
Article 

Google Scholar 
Mingo-Casas, P. et al. First cases of European bat lyssavirus type 1 in Iberian serotine bats: implications for the molecular epidemiology of bat rabies in Europe. PLoS Negl. Trop. Dis. 12(4), e0006290 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Vázquez-Moron, S., Juste, J., Ibáñez, C., Berciano, J. M. & Echevarria, J. E. Phylogeny of European bat Lyssavirus 1 in Eptesicus isabellinus bats, Spain. Emerg. Infect. Dis. 17, 520–523 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
Burgarella, C. et al. Detection of hybrids in nature: application to oaks (Quercus suber and Q. ilex). Heredity 102, 442–452 (2009).CAS 
PubMed 
Article 

Google Scholar 
Abrams, P. A. Character displacement and niche shift analyzed using consumer-resource models of competition. Theor. Popul. Biol. 29, 107–160 (1986).MathSciNet 
CAS 
PubMed 
MATH 
Article 

Google Scholar  More

Dickie, M., Serrouya, R., McNay, R. S. & Boutin, S. Faster and farther: wolf movement on linear features and implications for hunting behaviour. J. Appl. Ecol. 54, 253–263 (2017).Article 

Google Scholar 
Owen-Smith, N., Fryxell, J. M. & Merrill, E. H. Foraging theory upscaled: The behavioural ecology of herbivore movement. Philos. Trans. R. Soc. B Biol. Sci. 365, 2267–2278. https://doi.org/10.1098/rstb.2010.0095 (2010).CAS 
Article 

Google Scholar 
Holling, C. S. The functional response of predators to prey density and its role in mimicry and population regulation. Mem. Entomol. Soc. Can. 97, 5–60 (1965).Article 

Google Scholar 
Holling, C. The components of predation as revealed by a study of small-mammal predation of the European pine sawfly (1959).Dickie, M., McNay, S. R., Sutherland, G. D., Cody, M. & Avgar, T. Corridors or risk? Movement along, and use of, linear features varies predictably among large mammal predator and prey species. J. Anim. Ecol. https://doi.org/10.1111/1365-2656.13130 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
DeCesare, N. J. Separating spatial search and efficiency rates as components of predation risk. Proc. Biol. Sci. 279, 4626–4633. https://doi.org/10.1098/rspb.2012.1698 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Muhly, T. B., Semeniuk, C., Massolo, A., Hickman, L. & Musiani, M. Human activity helps prey win the predator-prey space race. PLoS ONE 6, e17050. https://doi.org/10.1371/journal.pone.0017050 (2011).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Fleming, P. A. & Bateman, P. W. Novel predation opportunities in anthropogenic landscapes. Anim. Behav. 138, 145–155. https://doi.org/10.1016/j.anbehav.2018.02.011 (2018).Article 

Google Scholar 
Whittington, J. et al. Caribou encounters with wolves increase near roads and trails: A time-to-event approach. J. Appl. Ecol. 48, 1535–1542. https://doi.org/10.1111/j.1365-2664.2011.02043.x (2011).Article 

Google Scholar 
Larivière, S. & Messier, F. Effect of density and nearest neighbours on simulated waterfowl nests: Can predators recognize high-density nesting patches?. Oikos 83, 12–20. https://doi.org/10.2307/3546541 (1998).Article 

Google Scholar 
Taitt, M. J. & Krebs, C. J. Predation, cover, and food manipulations during a spring decline of Microtus townsendii. J. Anim. Ecol. 52, 837–848. https://doi.org/10.2307/4458 (1983).Article 

Google Scholar 
Fisher, J. T. & Wilkinson, L. The response of mammals to forest fire and timber harvest in the North American boreal forest. Mammal. Rev. 35, 51–81 (2005).Article 

Google Scholar 
Fisher, J. T. & Burton, A. C. Wildlife winners and losers in an oil sands landscape. Front. Ecol. Environ. 16, 323–328. https://doi.org/10.1002/fee.1807 (2018).Article 

Google Scholar 
Francis, A. L., Procter, C., Kuzyk, G. & Fisher, J. T. Female Moose Prioritize Forage Over Mortality Risk in Harvested Landscapes. J. Wildl. Manag. (2021).Hebblewhite, M., Munro, R. H. & Merrill, E. H. Trophic consequences of postfire logging in a wolf–ungulate system. For. Ecol. Manag. 257, 1053–1062. https://doi.org/10.1016/j.foreco.2008.11.009 (2009).Article 

Google Scholar 
Pulliam, H. R. Sources, sinks, and population regulation. Am. Nat. 132, 652–661 (1988).Article 

Google Scholar 
Battin, J. When good animals love bad habitats: Ecological traps and the conservation of animal populations. Conserv. Biol. 18, 1482–1491 (2004).Article 

Google Scholar 
Nielsen, S. E., Stenhouse, G. B. & Boyce, M. S. A habitat-based framework for grizzly bear conservation in Alberta. Biol. Conserv. 130, 217–229 (2006).Article 

Google Scholar 
Bentz, B. et al. Salt Lake City 42 (University of Utah Press, 2005).
Google Scholar 
Carroll, A. L., Taylor, S. W., Régnière, J. & Safranyik, L. in Mountain pine beetle symposium: challenges and solutions. 223–232 (Natural Resources Canada, Canadian Forest Service, Pacific Forestry Centre).Lindenmayer, D. B. & Noss, R. F. Salvage logging, ecosystem processes, and biodiversity conservation. Conserv. Biol. 20, 949–958. https://doi.org/10.1111/j.1523-1739.2006.00497.x (2006).CAS 
Article 
PubMed 

Google Scholar 
Leverkus, A. B., Lindenmayer, D. B., Thorn, S. & Gustafsson, L. Salvage logging in the world’s forests: Interactions between natural disturbance and logging need recognition. Glob. Ecol. Biogeogr. 27, 1140–1154. https://doi.org/10.1111/geb.12772 (2018).Article 

Google Scholar 
Kuzyk, G. et al. Moose population dynamics during 20 years of declining harvest in British Columbia. Alces 54, 101–119 (2018).
Google Scholar 
Kuzyk, G. W. Provincial population and harvest estimates of moose in British Columbia. Alces J. Devot. Biol. Manag. Moose 52, 1–11 (2016).Procter, C. et al. Factors affecting moose population declines in British Columbia. 2020 Progress Report: February 2012-May 2020. B.C. Ministry of Forests, Lands, Natural Resource Operations and Rural Development, Victoria, B.C., Wildlife Working Report No. WR-128. Pp. 89. https://www2.gov.bc.ca/gov/content/environment/plants-animals-ecosystems/wildlife/wildlife-conservation/moose/moose-conservation/moose-research. (2020).Wittmer, H. U., Sinclair, A. R. E. & McLellan, B. N. The role of predation in the decline and extirpation of woodland caribou. Oecologia 144, 257–267. https://doi.org/10.1007/s00442-005-0055-y (2005).ADS 
Article 
PubMed 

Google Scholar 
Latham, A. D. M., Latham, M. C., Boyce, M. S. & Boutin, S. Movement responses by wolves to industrial linear features and their effect on woodland caribou in northeastern Alberta. Ecol. Appl. 21, 2854–2865 (2011).Article 

Google Scholar 
James, A. R. C. & Stuart-Smith, A. K. Distribution of caribou and wolves in relation to linear corridors. J. Wildl. Manag. 64, 154–159. https://doi.org/10.2307/3802985 (2000).Article 

Google Scholar 
DeMars, C. A. & Boutin, S. Nowhere to hide: Effects of linear features on predator–prey dynamics in a large mammal system. J. Anim. Ecol. 87, 274–284. https://doi.org/10.1111/1365-2656.12760 (2018).Article 
PubMed 

Google Scholar 
McKenzie, H. W., Merrill, E. H., Spiteri, R. J. & Lewis, M. A. How linear features alter predator movement and the functional response. Interface Focus 2, 205–216. https://doi.org/10.1098/rsfs.2011.0086 (2012).Article 
PubMed 
PubMed Central 

Google Scholar 
Houle, M., Fortin, D., Dussault, C., Courtois, R. & Ouellet, J.-P. Cumulative effects of forestry on habitat use by gray wolf (Canis lupus) in the boreal forest. Landsc. Ecol. 25, 419–433. https://doi.org/10.1007/s10980-009-9420-2 (2010).Article 

Google Scholar 
Kuzyk, G. W., Kneteman, J. & Schmiegelow, F. K. Winter habitat use by wolves, Canis lupus, in relation to forest harvesting in west-central Alberta. Can. Field Nat. 118, 368–375 (2004).Article 

Google Scholar 
Mumma, M. A. et al. Regional moose (Alces alces) responses to forestry cutblocks are driven by landscape-scale patterns of vegetation composition and regrowth. For. Ecol. Manag. 481, 118763 (2021).Article 

Google Scholar 
Scheideman, M. Use and selection at two spatial scales by female moose (Alces alces) across central British Columbia following a mountain pine beetle outbreak MSc thesis, University of Northern British Columbia (2018).Alfaro, R. I., van Akker, L. & Hawkes, B. Characteristics of forest legacies following two mountain pine beetle outbreaks in British Columbia Canada. Can. J. For. Res. 45, 1387–1396 (2015).Article 

Google Scholar 
Dhar, A., Parrott, L. & Hawkins, C. D. B. Aftermath of mountain pine beetle outbreak in British Columbia: Stand dynamics, management response and ecosystem resilience. Forests 7, 171 (2016).Article 

Google Scholar 
Shackelford, N., Standish, R. J., Ripple, W. & Starzomski, B. M. Threats to biodiversity from cumulative human impacts in one of North America’s last wildlife frontiers. Conserv. Biol. 32, 672–684 (2018).Article 

Google Scholar 
Corbett, L. J., Withey, P., Lantz, V. A. & Ochuodho, T. O. The economic impact of the mountain pine beetle infestation in British Columbia: Provincial estimates from a CGE analysis. For. Int. J. For. Res. 89, 100–105. https://doi.org/10.1093/forestry/cpv042 (2015).Latham, A. D. M. Wolf ecology and caribou-primary prey-wolf spatial relationships in low productivity peatland complexes in northeastern Alberta PhD thesis, University of Alberta, (2009).Person, D. K. & Russell, A. L. Reproduction and den site selection by wolves in a disturbed landscape. Northw. Sci. 83, 211–224. https://doi.org/10.3955/046.083.0305 (2009).Article 

Google Scholar 
Gillingham, M. Documentation for using Find Points Cluster Identification Program (Version 2) (University of Northern British Columbia, 2009).
Google Scholar 
Avgar, T., Potts, J. R., Lewis, M. A. & Boyce, M. S. Integrated step selection analysis: Bridging the gap between resource selection and animal movement. Methods Ecol. Evol. 7, 619–630. https://doi.org/10.1111/2041-210X.12528 (2016).Article 

Google Scholar 
Signer, J., Fieberg, J. & Avgar, T. Animal movement tools (amt): R package for managing tracking data and conducting habitat selection analyses. Ecol. Evol. 9, 880–890 (2019).Article 

Google Scholar 
Thurfjell, H., Ciuti, S. & Boyce, M. S. Applications of step-selection functions in ecology and conservation. Mov. Ecol. 2, 4. https://doi.org/10.1186/2051-3933-2-4 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Benson, J. F. & Patterson, B. R. Spatial overlap, proximity, and habitat use of individual wolves within the same packs. Wildl. Soc. Bull. (2011-) 39, 31–40 (2015).Fieberg, J., Matthiopoulos, J., Hebblewhite, M., Boyce, M. S. & Frair, J. L. Correlation and studies of habitat selection: problem, red herring or opportunity?. Philos. Trans. R. Soc. B Biol. Sci. 365, 2233–2244 (2010).Article 

Google Scholar 
Ladle, A. et al. Grizzly bear response to spatio-temporal variability in human recreational activity. J. Appl. Ecol. 56, 375–386. https://doi.org/10.1111/1365-2664.13277 (2019).Article 

Google Scholar 
Kohl, M. T. et al. Diel predator activity drives a dynamic landscape of fear. Ecol. Monogr. 88, 638–652 (2018).Article 

Google Scholar 
Scrafford, M. A., Avgar, T., Heeres, R. & Boyce, M. S. Roads elicit negative movement and habitat-selection responses by wolverines (Gulo gulo luscus). Behav. Ecol. 29, 534–542. https://doi.org/10.1093/beheco/arx182 (2018).Article 

Google Scholar 
Prokopenko, C. M., Boyce, M. S. & Avgar, T. Characterizing wildlife behavioural responses to roads using integrated step selection analysis. J. Appl. Ecol. 54, 470–479. https://doi.org/10.1111/1365-2664.12768 (2017).Article 

Google Scholar 
Avgar, T., Lele, S. R., Keim, J. L. & Boyce, M. S. Relative selection strength: Quantifying effect size in habitat- and step-selection inference. Ecol. Evol. 7, 5322–5330. https://doi.org/10.1002/ece3.3122 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Boyce, M. S., Vernier, P. R., Nielsen, S. E. & Schmiegelow, F. K. A. Evaluating resource selection functions. Ecol. Model. 157, 281–300. https://doi.org/10.1016/S0304-3800(02)00200-4 (2002).Article 

Google Scholar 
Visscher, D. R. & Merrill, E. H. Temporal dynamics of forage succession for elk at two scales: Implications of forest management. For. Ecol. Manag. 257, 96–106. https://doi.org/10.1016/j.foreco.2008.08.018 (2009).Article 

Google Scholar 
Stelfox, J. G., Lynch, G. M. & McGillis, J. R. Effects of clearcut logging on wild ungulates in the Central Albertan foothills. For. Chron. 52, 65–70. https://doi.org/10.5558/tfc52065-2 (1976).Article 

Google Scholar 
Gagné, C., Mainguy, J. & Fortin, D. The impact of forest harvesting on caribou–moose–wolf interactions decreases along a latitudinal gradient. Biol. Conserv. 197, 215–222. https://doi.org/10.1016/j.biocon.2016.03.015 (2016).Article 

Google Scholar 
Potvin, F., Breton, L. & Courtois, R. Response of beaver, moose, and snowshoe hare to clear-cutting in a Quebec boreal forest: a reassessment 10 years after cut. Can. J. For. Res. 35, 151–160 (2005).Article 

Google Scholar 
Rempel, R. S., Elkie, P. C., Rodgers, A. R. & Gluck, M. J. Timber-management and natural-disturbance effects on moose habitat: landscape evaluation. J. Wildl. Manag. 61, 517–524. https://doi.org/10.2307/3802610 (1997).Article 

Google Scholar 
Kunkel, K. E. & Pletscher, D. H. Habitat factors affecting vulnerability of moose to predation by wolves in southeastern British Columbia. Can. J. Zool. 78, 150–157. https://doi.org/10.1139/z99-181 (2000).Article 

Google Scholar 
Mech, L. D. & Boitani, L. Wolves: behavior, ecology, and conservation. (University of Chicago Press, 2007).Charnov, E. L. Optimal foraging, the marginal value theorem. (1976).Hebblewhite, M. & Merrill, E. H. Trade-offs between predation risk and forage differ between migrant strategies in a migratory ungulate. Ecology 90, 3445–3454. https://doi.org/10.1890/08-2090.1 (2009).Article 
PubMed 

Google Scholar 
Lendrum, P. E., Anderson Jr, C. R., Long, R. A., Kie, J. G. & Bowyer, R. T. Habitat selection by mule deer during migration: effects of landscape structure and natural-gas development. Ecosphere 3, art82. https://doi.org/10.1890/ES12-00165.1 (2012).Mumma, M. & Gillingham, M. Determining factors that affect survival of moose in Central British Columbia. Technical report to the Habitat Conservation Trust Foundation for Grant Agreement CAT19-0-522 (1 April 2017 through 31 March 2019). 56 (2019).Roffler, G. H., Gregovich, D. P. & Larson, K. R. Resource selection by coastal wolves reveals the seasonal importance of seral forest and suitable prey habitat. For. Ecol. Manag. 409, 190–201. https://doi.org/10.1016/j.foreco.2017.11.025 (2018).Article 

Google Scholar 
Lesmerises, F., Dussault, C. & St-Laurent, M.-H. Wolf habitat selection is shaped by human activities in a highly managed boreal forest. For. Ecol. Manag. 276, 125–131. https://doi.org/10.1016/j.foreco.2012.03.025 (2012).Article 

Google Scholar 
Muhly, T. B. et al. Functional response of wolves to human development across boreal North America. Ecol. Evol. 9, 10801–10815. https://doi.org/10.1002/ece3.5600 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Mladenoff, D. J., Sickley, T. A. & Wydeven, A. P. Predicting gray wolf landscape recolonization: logistic regression models vs. new field data. Ecol. Appl. 9, 37–44. https://doi.org/10.1890/1051-0761(1999)009[0037:PGWLRL]2.0.CO;2 (1999).Rogala, J. K. et al. Human activity differentially redistributes large mammals in the Canadian Rockies National Parks. Ecol. Soc. 16 (2011).Robertson, B. A. & Hutto, R. L. A framework for understanding ecological traps and an evaluation of existing evidence. Ecology 87, 1075–1085. https://doi.org/10.1890/0012-9658(2006)87[1075:AFFUET]2.0.CO;2 (2006).Article 
PubMed 

Google Scholar 
Finnegan, L. et al. Natural regeneration on seismic lines influences movement behaviour of wolves and grizzly bears. PLoS ONE 13, e0195480. https://doi.org/10.1371/journal.pone.0195480 (2018).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Dickie, M., Serrouya, R., DeMars, C., Cranston, J. & Boutin, S. Evaluating functional recovery of habitat for threatened woodland caribou. Ecosphere 8, e01936. https://doi.org/10.1002/ecs2.1936 (2017).Article 

Google Scholar  More

Gattuso JP, Frankignoulle M, Wollast R. Carbon and carbonate metabolism in coastal aquatic ecosystems. Annu Rev Ecol Syst. 1998;29:405–34.Article 

Google Scholar 
Arnosti C, Wietz M, Brinkhoff T, Hehemann JH, Probandt D, Zeugner L, et al. The biogeochemistry of marine polysaccharides: sources, inventories, and bacterial drivers of the carbohydrate cycle. Ann Rev Mar Sci. 2021;13:81–108.CAS 
PubMed 
Article 

Google Scholar 
Moran MA, Kujawinski EB, Stubbins A, Fatland R, Aluwihare LI, Buchan A, et al. Deciphering ocean carbon in a changing world. Proc Natl Acad Sci USA. 2016;113:3143–51.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Poretsky RS, Sun S, Mou X, Moran MA. Transporter genes expressed by coastal bacterioplankton in response to dissolved organic carbon. Environ Microbiol. 2010;12:616–27.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Azam F. Microbial control of oceanic carbon flux: The plot thickens. Science.1998;280:694–96.CAS 
Article 

Google Scholar 
Buchan A, LeCleir GR, Gulvik CA, González JM. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat Rev Microbiol. 2014;12:686–98.CAS 
PubMed 
Article 

Google Scholar 
Bunse C, Bertos-Fortis M, Sassenhagen I, Sildever S, Sjoqvist C, Godhe A, et al. Spatio-temporal interdependence of bacteria and phytoplankton during a baltic sea spring bloom. Front Microbiol. 2016;7:517.PubMed 
PubMed Central 
Article 

Google Scholar 
Cram JA, Chow CE, Sachdeva R, Needham DM, Parada AE, Steele JA, et al. Seasonal and interannual variability of the marine bacterioplankton community throughout the water column over ten years. ISME J. 2015;9:563–80.PubMed 
Article 

Google Scholar 
Seymour JR, Amin SA, Raina JB, Stocker R. Zooming in on the phycosphere: the ecological interface for phytoplankton-bacteria relationships. Nat Microbiol. 2017;2:17065.CAS 
PubMed 
Article 

Google Scholar 
Taylor JD, Cottingham SD, Billinge J, Cunliffe M. Seasonal microbial community dynamics correlate with phytoplankton-derived polysaccharides in surface coastal waters. ISME J. 2014;8:245–8.CAS 
PubMed 
Article 

Google Scholar 
Teeling H, Fuchs BM, Becher D, Klockow C, Gardebrecht A, Bennke CM, et al. Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Science.2012;336:608–11.CAS 
PubMed 
Article 

Google Scholar 
Hernando-Morales V, Varela MM, Needham DM, Cram J, Fuhrman JA, Teira E. Vertical and seasonal patterns control bacterioplankton communities at two horizontally coherent coastal upwelling sites off galicia (NW Spain). Microb Ecol. 2018;76:866–84.Needham DM, Fuhrman JA. Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom. Nat Microbiol. 2016;1:16005.Tada Y, Taniguchi A, Nagao I, Miki T, Uematsu M, Tsuda A, et al. Differing growth responses of major phylogenetic groups of marine bacteria to natural phytoplankton blooms in the western North Pacific Ocean. Appl Environ Microbiol. 2011;77:4055–65.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Williams TJ, Wilkins D, Long E, Evans F, DeMaere MZ, Raftery MJ, et al. The role of planktonic Flavobacteria in processing algal organic matter in coastal East Antarctica revealed using metagenomics and metaproteomics. Environ Microbiol. 2013;15:1302–17.CAS 
PubMed 
Article 

Google Scholar 
Ottesen EA, Young CR, Eppley JM, Ryan JP, Chavez FP, Scholin CA, et al. Pattern and synchrony of gene expression among sympatric marine microbial populations. Proc Natl Acad Sci USA. 2013;110:E488–97.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Ottesen EA, Young CR, Gifford SM, Eppley JM, Marin R 3rd, Schuster SC, et al. Ocean microbes. Multispecies diel transcriptional oscillations in open ocean heterotrophic bacterial assemblages. Science. 2014;345:207–12.CAS 
PubMed 
Article 

Google Scholar 
Beier S, Rivers AR, Moran MA, Obernosterer I. The transcriptional response of prokaryotes to phytoplankton-derived dissolved organic matter in seawater. Environ Microbiol. 2015;17:3466–80.CAS 
PubMed 
Article 

Google Scholar 
Sarmento H, Gasol JM. Use of phytoplankton-derived dissolved organic carbon by different types of bacterioplankton. Environ Microbiol. 2012;14:2348–60.CAS 
PubMed 
Article 

Google Scholar 
Shi Y, McCarren J, DeLong EF. Transcriptional responses of surface water marine microbial assemblages to deep-sea water amendment. Environ Microbiol. 2012;14:191–206.CAS 
PubMed 
Article 

Google Scholar 
Vorobev A, Sharma S, Yu M, Lee J, Washington BJ, Whitman WB, et al. Identifying labile DOM components in a coastal ocean through depleted bacterial transcripts and chemical signals. Environ Microbiol. 2018;20:3012–30.CAS 
PubMed 
Article 

Google Scholar 
Moran MA, Belas R, Schell MA, González JM, Sun F, Sun S, et al. Ecological genomics of marine Roseobacters. Appl Environ Microbiol. 2007;73:4559–69.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Nowinski B, Moran MA. Niche dimensions of a marine bacterium are identified using invasion studies in coastal seawater. Nat Microbiol. 2021;6:524–32.CAS 
PubMed 
Article 

Google Scholar 
Rinta-Kanto JM, Sun S, Sharma S, Kiene RP, Moran MA. Bacterial community transcription patterns during a marine phytoplankton bloom. Environ Microbiol. 2012;14:228–39.CAS 
PubMed 
Article 

Google Scholar 
Sharma AK, Becker JW, Ottesen EA, Bryant JA, Duhamel S, Karl DM, et al. Distinct dissolved organic matter sources induce rapid transcriptional responses in coexisting populations of Prochlorococcus, Pelagibacter and the OM60 clade. Environ Microbiol. 2013;16:2815–30.PubMed 
Article 
CAS 

Google Scholar 
Kieft B, Li Z, Bryson S, Hettich RL, Pan C, Mayali X, et al. Phytoplankton exudates and lysates support distinct microbial consortia with specialized metabolic and ecophysiological traits. Proc Natl Acad Sci USA. 2021;118:e2101178118.McCarren J, Becker JW, Repeta DJ, Shi Y, Young CR, Malmstrom RR, et al. Microbial community transcriptomes reveal microbes and metabolic pathways associated with dissolved organic matter turnover in the sea. Proc Natl Acad Sci USA. 2010;107:16420–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Sosa OA, Gifford SM, Repeta DJ, DeLong EF. High molecular weight dissolved organic matter enrichment selects for methylotrophs in dilution to extinction cultures. ISME J. 2015;9:2725–39.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pontiller B, Martínez-García S, Lundin D, Pinhassi J. Labile dissolved organic matter compound characteristics select for divergence in marine bacterial activity and transcription. Front Microbiol. 2020;11:588778.PubMed 
PubMed Central 
Article 

Google Scholar 
Bryson S, Li Z, Chavez F, Weber PK, Pett-Ridge J, Hettich RL, et al. Phylogenetically conserved resource partitioning in the coastal microbial loop. ISME J. 2017;11:2781–92.PubMed 
PubMed Central 
Article 

Google Scholar 
Joglar V, Prieto A, Barber-Lluch E, Hernández-Ruiz M, Fernández E, Teira E. Spatial and temporal variability in the response of phytoplankton and prokaryotes to B-vitamin amendments in an upwelling system. Biogeosciences.2020;17:2807–23.Article 

Google Scholar 
Martínez-García S, Fernández E, Álvarez-Salgado XA, González J, Lønborg C, Marañón E, et al. Differential responses of phytoplankton and heterotrophic bacteria to organic and inorganic nutrient additions in coastal waters off the NW Iberian Peninsula. Mar Ecol Prog Ser. 2010;416:17–33.Article 
CAS 

Google Scholar 
Welschmeyer NA. Fluorometric analysis of chlorophyll a in the presence of chlorophyll b and pheopigments. Limnol Oceanogr. 1994;39:1985–92.CAS 
Article 

Google Scholar 
Parsons TR, Maita Y, Lalli CM. Fluorometric determination of chlorophylls. A manual of chemical & biological methods for seawater analysis: Oxford, UK: Pergamon Press; 1984. p. 107–09.Hansen H, Grasshoff K. Automated chemical analysis. Methods of seawater analysis. 2nd ed: Verlag Chemie, Weinheim; 1983. p. 347–95.Álvarez-Salgado XA, Miller AEJ. Simultaneous determination of dissolved organic carbon and total dissolved nitrogen in seawater by high temperature catalytic oxidation: conditions for precise shipboard measurements. Mar Chem. 1998;62:325–33.Article 

Google Scholar 
Calvo-Díaz A, Morán XAG. Seasonal dynamics of picoplankton in shelf waters of the southern Bay of Biscay. Aquat Micro Ecol. 2006;42:159–74.Article 

Google Scholar 
Smith DC, Farooq A. A simple, economical method for measuring bacterial protein synthesis rates in seawater using 3H-leucine. Mar Micro Food Webs. 1992;6:107–14.
Google Scholar 
Parada AE, Needham DM, Fuhrman JA. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol. 2016;18:1403–14.CAS 
PubMed 
Article 

Google Scholar 
Logares R, Sunagawa S, Salazar G, Cornejo-Castillo FM, Ferrera I, Sarmento H, et al. Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities. Environ Microbiol. 2014;16:2659–71.CAS 
PubMed 
Article 

Google Scholar 
Straub D, Blackwell N, Langarica-Fuentes A, Peltzer A, Nahnsen S, Kleindienst S. Interpretations of environmental microbial community studies are biased by the selected 16S rRNA (Gene) amplicon sequencing pipeline. Front Microbiol. 2020;11:1–18.Article 

Google Scholar 
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.CAS 
PubMed 
Article 

Google Scholar 
Poretsky RS, Gifford S, Rinta-Kanto J, Vila-Costa M, Moran MA. Analyzing gene expression from marine microbial communities using environmental transcriptomics. J Vis Exp. 2009;24:1–6.
Google Scholar 
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. 2011;17:10–12.Article 

Google Scholar 
Joshi NA, Fass JN. Sickle: a sliding-window, adaptive, quality-based trimming tool for FastQ files. 2011. https://github.com/najoshi/sickle.Del Fabbro C, Scalabrin S, Morgante M, Giorgi FM. An extensive evaluation of read trimming effects on Illumina NGS data analysis. PLoS One. 2013;8:1–13.
Google Scholar 
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics.2015;31:1674–6.CAS 
PubMed 
Article 

Google Scholar 
Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:119.Article 
CAS 

Google Scholar 
Rognes T, Flouri T, Nichols B, Quince C, Mahe F. VSEARCH: a versatile open source tool for metagenomics. PeerJ.2016;4:e2584.PubMed 
PubMed Central 
Article 

Google Scholar 
Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.CAS 
PubMed 
Article 

Google Scholar 
Huson DH, Beier S, Flade I, Górska A, El-Hadidi M, Mitra S, et al. MEGAN community edition – interactive exploration and analysis of large-scale microbiome sequencing data. PLoS Comput Biol. 2016;12:1–12.Article 
CAS 

Google Scholar 
Zhang H, Yohe T, Huang L, Entwistle S, Wu P, Yang Z, et al. dbCAN2: A meta server for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2018;46:W95–101.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Gloor GB, Wu JR, Pawlowsky-Glahn V, Egozcue JJ. It’s all relative: analyzing microbiome data as compositions. Ann Epidemiol. 2016;26:322–9.PubMed 
Article 

Google Scholar 
R Core Team. R: A language and environment for statistical computing. 4.1.0 ed. Vienna, Austria: R Foundation for Statistical Computing; 2021.Cottrell MT, Kirchman DL. Transcriptional control in marine copiotrophic and oligotrophic bacteria with streamlined genomes. Appl Environ Microbiol. 2016;82:6010–18.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Giovannoni SJ. SAR11 bacteria: The most abundant plankton in the oceans. Ann Rev Mar Sci. 2017;9:231–55.PubMed 
Article 

Google Scholar 
Newton RJ, Griffin LE, Bowles KM, Meile C, Gifford S, Givens CE, et al. Genome characteristics of a generalist marine bacterial lineage. ISME J. 2010;4:784–98.CAS 
PubMed 
Article 

Google Scholar 
Allers E, Gomez-Consarnau L, Pinhassi J, Gasol JM, Simek K, Pernthaler J. Response of Alteromonadaceae and Rhodobacteriaceae to glucose and phosphorus manipulation in marine mesocosms. Environ Microbiol. 2007;9:2417–29.CAS 
PubMed 
Article 

Google Scholar 
Alonso C, Pernthaler J. Roseobacter and SAR11 dominate microbial glucose uptake in coastal North Sea waters. Environ Microbiol. 2006;8:2022–30.CAS 
PubMed 
Article 

Google Scholar 
Noell SE, Giovannoni SJ. SAR11 bacteria have a high affinity and multifunctional glycine betaine transporter. Environ Microbiol. 2019;21:2559–75.CAS 
PubMed 
Article 

Google Scholar 
Sowell SM, Wilhelm LJ, Norbeck AD, Lipton MS, Nicora CD, Barofsky DF, et al. Transport functions dominate the SAR11 metaproteome at low-nutrient extremes in the Sargasso Sea. ISME J. 2009;3:93–105.CAS 
PubMed 
Article 

Google Scholar 
Chenard C, Wijaya W, Vaulot D, Lopes Dos Santos A, Martin P, Kaur A, et al. Temporal and spatial dynamics of bacteria, archaea and protists in equatorial coastal waters. Sci Rep. 2019;9:16390.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Pajares S, Varona-Cordero F, Hernández-Becerril DU. Spatial distribution patterns of bacterioplankton in the oxygen minimum zone of the tropical mexican pacific. Micro Ecol. 2020;80:519–36.CAS 
Article 

Google Scholar 
Signori CN, Pellizari VH, Enrich-Prast A, Sievert SM. Spatiotemporal dynamics of marine bacterial and archaeal communities in surface waters off the northern Antarctic Peninsula. Deep-Sea Res Pt Ii. 2018;149:150–60.Article 

Google Scholar 
Ling SK, Xia J, Liu Y, Chen GJ, Du ZJ. Agarilytica rhodophyticola gen. nov., sp. nov., isolated from Gracilaria blodgettii. Int J Syst Evol Microbiol. 2017;67:3778–83.CAS 
PubMed 
Article 

Google Scholar 
Cheng H, Zhang S, Huo YY, Jiang XW, Zhang XQ, Pan J, et al. Gilvimarinus polysaccharolyticus sp. nov., an agar-digesting bacterium isolated from seaweed, and emended description of the genus Gilvimarinus. Int J Syst Evol Microbiol. 2015;65:562–69.CAS 
PubMed 
Article 

Google Scholar 
Malmstrom RR, Kiene RP, Cottrell MT, Kirchman DL. Contribution of SAR11 bacteria to dissolved dimethylsulfoniopropionate and amino acid uptake in the North Atlantic ocean. Appl Environ Microbiol. 2004;70:4129–35.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kirchman DL. Growth rates of microbes in the oceans. Ann Rev Mar Sci. 2016;8:285–309.PubMed 
Article 

Google Scholar 
Durham BP, Dearth SP, Sharma S, Amin SA, Smith CB, Campagna SR, et al. Recognition cascade and metabolite transfer in a marine bacteria-phytoplankton model system. Environ Microbiol. 2017;19:3500–13.CAS 
PubMed 
Article 

Google Scholar 
Ferrer-González FX, Widner B, Holderman NR, Glushka J, Edison AS, Kujawinski EB, et al. Resource partitioning of phytoplankton metabolites that support bacterial heterotrophy. ISME J. 2021;15:762–73.PubMed 
Article 
CAS 

Google Scholar 
Landa M, Burns AS, Roth SJ, Moran MA. Bacterial transcriptome remodeling during sequential co-culture with a marine dinoflagellate and diatom. ISME J. 2017;11:2677–90.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pedler BE, Aluwihare LI, Azam F. Single bacterial strain capable of significant contribution to carbon cycling in the surface ocean. Proc Natl Acad Sci USA. 2014;111:7202–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Koch H, Dürwald A, Schweder T, Noriega-Ortega B, Vidal-Melgosa S, Hehemann JH, et al. Biphasic cellular adaptations and ecological implications of Alteromonas macleodii degrading a mixture of algal polysaccharides. ISME J. 2019;13:92–103.CAS 
PubMed 
Article 

Google Scholar 
López-Pérez M, Gonzaga A, Martin-Cuadrado AB, Onyshchenko O, Ghavidel A, Ghai R, et al. Genomes of surface isolates of Alteromonas macleodii: the life of a widespread marine opportunistic copiotroph. Sci Rep. 2012;2:696.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Bergauer K, Fernandez-Guerra A, Garcia JAL, Sprenger RR, Stepanauskas R, Pachiadaki MG, et al. Organic matter processing by microbial communities throughout the Atlantic water column as revealed by metaproteomics. Proc Natl Acad Sci USA. 2018;115:E400–8.CAS 
PubMed 
Article 

Google Scholar 
Hou S, López-Pérez M, Pfreundt U, Belkin N, Stüber K, Huettel B, et al. Benefit from decline: the primary transcriptome of Alteromonas macleodii str. Te101 during Trichodesmium demise. ISME J. 2018;12:981–96.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Pinhassi J, Sala MM, Havskum H, Peters F, Guadayol O, Malits A, et al. Changes in bacterioplankton composition under different phytoplankton regimens. Appl Environ Microbiol. 2004;70:6753–66.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Riemann L, Steward GF, Azam F. Dynamics of bacterial community composition and activity during a mesocosm diatom bloom. Appl Environ Microbiol. 2000;66:578–87.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Cottrell MT, Kirchman DL. Natural assemblages of marine Proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl Environ Microbiol. 2000;66:1692–7.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Krüger K, Chafee M, Ben Francis T, Glavina Del Rio T, Becher D, Schweder T, et al. In marine Bacteroidetes the bulk of glycan degradation during algae blooms is mediated by few clades using a restricted set of genes. ISME J. 2019;13:2800–16.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Ben Hania W, Joseph M, Bunk B, Spröer C, Klenk HP, Fardeau ML, et al. Characterization of the first cultured representative of a Bacteroidetes clade specialized on the scavenging of cyanobacteria. Environ Microbiol. 2017;19:1134–48.PubMed 
Article 
CAS 

Google Scholar 
Fernández-Gómez B, Richter M, Schüler M, Pinhassi J, Acinas SG, González JM, et al. Ecology of marine Bacteroidetes: a comparative genomics approach. ISME J. 2013;7:1026–37.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Orsi WD, Smith JM, Liu S, Liu Z, Sakamoto CM, Wilken S, et al. Diverse, uncultivated bacteria and archaea underlying the cycling of dissolved protein in the ocean. ISME J. 2016;10:2158–73.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Xing P, Hahnke RL, Unfried F, Markert S, Huang S, Barbeyron T, et al. Niches of two polysaccharide-degrading Polaribacter isolates from the North Sea during a spring diatom bloom. ISME J. 2015;9:1410–22.CAS 
PubMed 
Article 

Google Scholar 
Sichert A, Corzett CH, Schechter MS, Unfried F, Markert S, Becher D, et al. Verrucomicrobia use hundreds of enzymes to digest the algal polysaccharide fucoidan. Nat Microbiol. 2020;5:1026–39.Ivanova AA, Naumoff DG, Miroshnikov KK, Liesack W, Dedysh SN. Comparative genomics of four Isosphaeraceae Planctomycetes: a common pool of plasmids and glycoside hydrolase genes shared by Paludisphaera borealis PX4(T), Isosphaera pallida IS1B(T), Singulisphaera acidiphila DSM 18658(T), and Strain SH-PL62. Front Microbiol. 2017;8:412.PubMed 
PubMed Central 
Article 

Google Scholar 
Vidal-Melgosa S, Sichert A, Francis TB, Bartosik D, Niggemann J, Wichels A, et al. Diatom fucan polysaccharide precipitates carbon during algal blooms. Nat Commun. 2021;12:1150.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Le Costaouëc T, Unamunzaga C, Mantecon L, Helbert W. New structural insights into the cell-wall polysaccharide of the diatom Phaeodactylum tricornutum. Algal Res. 2017;26:172–9.Article 

Google Scholar 
Francis TB, Bartosik D, Sura T, Sichert A, Hehemann JH, Markert S, et al. Changing expression patterns of TonB-dependent transporters suggest shifts in polysaccharide consumption over the course of a spring phytoplankton bloom. ISME J. 2021;15:2336–50.Bode A, Estévez MG, Varela M, Vilar JA. Annual trend patterns of phytoplankton species abundance belie homogeneous taxonomical group responses to climate in the NE Atlantic upwelling. Mar Environ Res. 2015;110:81–91.CAS 
PubMed 
Article 

Google Scholar 
Cermeño P, Marañón E, Pérez V, Serret P, Fernández E, Castro CG. Phytoplankton size structure and primary production in a highly dynamic coastal ecosystem (Ría de Vigo, NW-Spain): Seasonal and short-time scale variability. Estuar Coast Shelf Sci. 2006;67:251–66.Article 

Google Scholar 
Nogueira E, Pérez FF, Rı́os AF. Modelling thermohaline properties in an estuarine upwelling ecosystem (Rı́a de Vigo: NW Spain) using box-jenkins transfer function models. Estuar Coast Shelf Sci. 1997;44:685–702.Article 

Google Scholar 
Broullón E, López-Mozos M, Reguera B, Chouciño P, Doval MD, Fernández-Castro B, et al. Thin layers of phytoplankton and harmful algae events in a coastal upwelling system. Prog Oceanogr. 2020;189:102449.Fraga F. Upwelling off the Galician Coast, Northwest Spain. In: Richards FA, editor. Coastal Upwelling. Washington: American Geophysical Union; 1981. p. 176–82.Nogueira E, Figueiras FG. The microplankton succession in the Ría de Vigo revisited: species assemblages and the role of weather-induced, hydrodynamic variability. J Mar Syst. 2005;54:139–55.Article 

Google Scholar 
Pitcher GC, Walker DR, Mitchellinnes BA, Moloney CL. Short-term variability during an anchor station study in the southern Benguela Upwelling System – phytoplankton dynamics. Prog Oceanogr. 1991;28:39–64.Article 

Google Scholar 
Smayda TJ, Trainer VL. Dinoflagellate blooms in upwelling systems: Seeding, variability, and contrasts with diatom bloom behaviour. Prog Oceanogr. 2010;85:92–107.Article 

Google Scholar 
Wilkerson FP, Lassiter AM, Dugdale RC, Marchi A, Hogue VE. The phytoplankton bloom response to wind events and upwelled nutrients during the CoOP WEST study. Deep-Sea Res Pt Ii. 2006;53:3023–48.Article 

Google Scholar 
Reintjes G, Arnosti C, Fuchs B, Amann R. Selfish, sharing and scavenging bacteria in the Atlantic Ocean: a biogeographical study of bacterial substrate utilisation. ISME J. 2019;13:1119–32.CAS 
PubMed 
Article 

Google Scholar 
Mayali X, Weber PK, Pett-Ridge J. Taxon-specific C/N relative use efficiency for amino acids in an estuarine community. FEMS Microbiol Ecol. 2013;83:402–12.CAS 
PubMed 
Article 

Google Scholar  More

Newsome, S. D., Clementz, M. T. & Koch, P. L. Using stable isotope biogeochemistry to study marine mammal ecology. Mar. Mamm. Sci. 26, 509–572. https://doi.org/10.1111/j.1748-7692.2009.00354.x (2010).CAS 
Article 

Google Scholar 
Layman, C. A. et al. Applying stable isotopes to examine food-web structure: An overview of analytical tools. Biol. Rev. Camb. Philos. Soc. 87, 545–562. https://doi.org/10.1111/j.1469-185X.2011.00208.x (2011).Article 
PubMed 

Google Scholar 
Larsen, T. et al. Tracing carbon sources through aquatic and terrestrial food webs using amino acid stable isotope fingerprinting. PLoS ONE 8, e73441. https://doi.org/10.1371/journal.pone.0073441 (2013).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Post, D. M. Using stable isotopes to estimate trophic position: Models, methods and assumptions. Ecology 83, 703–718 (2002).Article 

Google Scholar 
Inger, R. & Bearhop, S. Applications of stable isotope analyses to avian ecology. Ibis 150, 447–461 (2008).Article 

Google Scholar 
McCutchan, J. H., Lewis, W. M., Kendall, C. & McGrath, C. C. Variation in trophic shift for stable isotope ratios of carbon, nitrogen, and sulfur. Oikos 102, 378–390 (2003).CAS 
Article 

Google Scholar 
Olive, P. J. W., Pinnegar, J. K., Polunin, N. V. C., Richards, G. & Welch, R. Isotope trophic-step fractionation: A dynamic equilibrium model. J. Anim. Ecol. 72, 608–617 (2003).Article 

Google Scholar 
McMahon, K. W., Polito, M. J., Abel, S., McCarthy, M. D. & Thorrold, S. R. Carbon and nitrogen isotope fractionation of amino acids in an avian marine predator, the gentoo penguin (Pygoscelis papua). Ecol. Evol. 5, 1278–1290. https://doi.org/10.1002/ece3.1437 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Webb, E. C. et al. Compound-specific amino acid isotopic proxies for distinguishing between terrestrial and aquatic resource consumption. Archaeol. Anthropol. Sci. 10, 1–18. https://doi.org/10.1007/s12520-015-0309-5 (2016).Article 

Google Scholar 
Whiteman, J. P., Kim, S. L., McMahon, K. W., Koch, P. L. & Newsome, S. D. Amino acid isotope discrimination factors for a carnivore: Physiological insights from leopard sharks and their diet. Oecologia 188, 977–989. https://doi.org/10.1007/s00442-018-4276-2 (2018).ADS 
Article 
PubMed 

Google Scholar 
Rogers, M., Bare, R., Gray, A., Scott-Moelder, T. & Heintz, R. Assessment of two feeds on survival, proximate composition, and amino acid carbon isotope discrimination in hatchery-reared Chinook salmon. Fisher. Res. https://doi.org/10.1016/j.fishres.2019.06.001 (2019).Article 

Google Scholar 
Wang, Y. V., Wan, A. H. L., Krogdahl, A., Johnson, M. & Larsen, T. (13)C values of glycolytic amino acids as indicators of carbohydrate utilization in carnivorous fish. PeerJ 7, e7701. https://doi.org/10.7717/peerj.7701 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
McMahon, K. W., Fogel, M. L., Elsdon, T. S. & Thorrold, S. R. Carbon isotope fractionation of amino acids in fish muscle reflects biosynthesis and isotopic routing from dietary protein. J. Anim. Ecol. 79, 1132–1141. https://doi.org/10.1111/j.1365-2656.2010.01722.x (2010).Article 
PubMed 

Google Scholar 
McMahon, K. W., Thorrold, S. R., Houghton, L. A. & Berumen, M. L. Tracing carbon flow through coral reef food webs using a compound-specific stable isotope approach. Oecologia 180, 809–821. https://doi.org/10.1007/s00442-015-3475-3 (2016).ADS 
Article 
PubMed 

Google Scholar 
Wang, Y. V. et al. Know your fish: A novel compound-specific isotope approach for tracing wild and farmed salmon. Food Chem 256, 380–389. https://doi.org/10.1016/j.foodchem.2018.02.095 (2018).CAS 
Article 
PubMed 

Google Scholar 
Jim, S., Jones, V., Ambrose, S. H. & Evershed, R. P. Quantifying dietary macronutrient sources of carbon for bone collagen biosynthesis using natural abundance stable carbon isotope analysis. Br J. Nutr. 95, 1055–1062. https://doi.org/10.1079/bjn20051685 (2006).CAS 
Article 
PubMed 

Google Scholar 
Newsome, S. D., Fogel, M. L., Kelly, L. & del Rio, C. M. Contributions of direct incorporation from diet and microbial amino acids to protein synthesis in Nile tilapia. Funct. Ecol. 25, 1051–1062. https://doi.org/10.1111/j.1365-2435.2011.01866.x (2011).Article 

Google Scholar 
Griffiths, H. Applications of stable isotope technology in physiological ecology. Funct. Ecol. 5, 254–269 (1991).Article 

Google Scholar 
Lorrain, A. et al. Differential δ13C and δ15N signatures among scallop tissues: Implications for ecology and physiology. J. Exp. Mar. Biol. Ecol. 275, 47–61 (2002).CAS 
Article 

Google Scholar 
Li, P., Mai, K., Trushenski, J. & Wu, G. New developments in fish amino acid nutrition: Towards functional and environmentally oriented aquafeeds. Amino Acids 37, 43–53. https://doi.org/10.1007/s00726-008-0171-1 (2009).CAS 
Article 
PubMed 

Google Scholar 
Boecklen, W. J., Yarnes, C. T., Cook, B. A. & James, A. C. On the use of stable isotopes in trophic ecology. Annu. Rev. Ecol. Evol. Syst. 42, 411–440. https://doi.org/10.1146/annurev-ecolsys-102209-144726 (2011).Article 

Google Scholar 
Perga, M. E. & Gerdeaux, D. “Are fish what they eat” all year round?. Oecologia 144, 598–606. https://doi.org/10.1007/s00442-005-0069-5 (2005).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Sponheimer, M. et al. Turnover of stable carbon isotopes in the muscle, liver, and breath CO2 of alpacas (Lama pacos). Rapid Commun. Mass Spectrom. 20, 1395–1399. https://doi.org/10.1002/rcm.2454 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Logan, J. M. & Lutcavage, M. E. Stable isotope dynamics in elasmobranch fishes. Hydrobiologia 644, 231–244. https://doi.org/10.1007/s10750-010-0120-3 (2010).CAS 
Article 

Google Scholar 
Madigan, D. J. et al. Tissue turnover rates and isotopic trophic discrimination factors in the endothermic teleost, pacific bluefin tuna (Thunnus orientalis). PLoS ONE 7, e49220. https://doi.org/10.1371/journal.pone.0049220 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Skinner, M. M., Cross, B. K. & Moore, B. C. Estimating in situ isotopic turnover in Rainbow Trout (Oncorhynchus mykiss) muscle and liver tissue. J. Freshw. Ecol. 32, 209–217. https://doi.org/10.1080/02705060.2016.1259127 (2016).CAS 
Article 

Google Scholar 
Kaushik, S. J. & Seiliez, I. Protein and amino acid nutrition and metabolism in fish: Current knowledge and future needs. Aquac. Res. 41, 322–332. https://doi.org/10.1111/j.1365-2109.2009.02174.x (2010).CAS 
Article 

Google Scholar 
Hou, Y., Hu, S., Li, X., He, W. & Wu, G. Amino Acid Metabolism in the Liver: Nutritional and Physiological Significance. Vol. 1265 (2020).Gannes, L. Z., O’Brien, D. M. & Del Rio, C. M. Stable isotopes in animal ecology: Assumptions, caveats and a call for more laboratory experiments. Ecology 78, 1271–1276 (1997).Article 

Google Scholar 
Martinez del Rio, C. M., Wolf, N., Carleton, S. A. & Gannes, L. Z. Isotopic ecology ten years after a call for more laboratory experiments. Biol. Rev. Camb. Philos Soc. 84, 91–111. https://doi.org/10.1111/j.1469-185X.2008.00064.x (2009).Article 

Google Scholar 
Hendry, A. P., Peichel, C. L., Boughman, J. W., Matthews, B. & Nosil, P. Stickleback research: The now and the next. Evol. Ecol. Res. 15, 111–141 (2013).
Google Scholar 
Fang, B., Merila, J., Ribeiro, F., Alexandre, C. M. & Momigliano, P. Worldwide phylogeny of three-spined sticklebacks. Mol Phylogenet Evol 127, 613–625. https://doi.org/10.1016/j.ympev.2018.06.008 (2018).Article 
PubMed 

Google Scholar 
Kume, M. & Kitano, J. Genetic and stable isotope analyses of threespine stickleback from the Bering and Chukchi seas. Ichthyol. Res. 64, 478–480. https://doi.org/10.1007/s10228-017-0580-9 (2017).Article 

Google Scholar 
Reimchen, T. E., Ingram, T. & Hansen, S. C. Assessing niche differences of sex, armour and asymmetry phenotypes using stable isotope analyses in Haida Gwaii sticklebacks. Behaviour 145, 561–577 (2008).Article 

Google Scholar 
Pinnegar, J. Unusual stable isotope fractionation patterns observed for fish host–parasite trophic relationships. J. Fish Biol. 59, 494–503. https://doi.org/10.1006/jfbi.2001.1660 (2001).Article 

Google Scholar 
Power, M. & Klein, G. M. Fish host-cestode parasite stable isotope enrichment patterns in marine, estuarine and freshwater fishes from northern Canada. Isotopes Environ. Health Stud. 40, 257–266 (2004).CAS 
Article 

Google Scholar 
Li, X., Zheng, S. & Wu, G. Nutrition and metabolism of glutamate and glutamine in fish. Amino Acids 52, 671–691. https://doi.org/10.1007/s00726-020-02851-2 (2020).CAS 
Article 
PubMed 

Google Scholar 
Vander Zanden, M. J., Clayton, M. K., Moody, E. K., Solomon, C. T. & Weidel, B. C. Stable isotope turnover and half-life in animal tissues: A literature synthesis. PLoS ONE 10, e0116182. https://doi.org/10.1371/journal.pone.0116182 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Newsome, S. D., del Rio, C. M., Bearhop, S. & Phillips, D. L. A niche for isotopic ecology. Front. Ecol. Environ. 5, 429–436. https://doi.org/10.1890/060150.01 (2007).Article 

Google Scholar 
Voigt, C. C., Rex, K., Michener, R. H. & Speakman, J. R. Nutrient routing in omnivorous animals tracked by stable carbon isotopes in tissue and exhaled breath. Oecologia 157, 31–40. https://doi.org/10.1007/s00442-008-1057-3 (2008).ADS 
Article 
PubMed 

Google Scholar 
Tieszen, L. L., Boutton, T. W., Tesdahl, K. G. & Slade, N. A. Fractionation and turnover of stable carbon isotopes in animal tissues: Implications for δ13C analysis of diet. Oecologia 57, 21–37 (1983).Article 

Google Scholar 
Cerling, T. E. et al. Determining biological tissue turnover using stable isotopes: The reaction progress variable. Oecologia 151, 175–189. https://doi.org/10.1007/s00442-006-0571-4 (2007).ADS 
Article 
PubMed 

Google Scholar 
Martínez del Rio, C. & Carleton, S. A. How fast and how faithful: The dynamics of isotopic incorporation into animal tissues: Fig. 1. J. Mammal. 93, 353–359. https://doi.org/10.1644/11-mamm-s-165.1 (2012).Article 

Google Scholar 
McCullagh, J. S., Juchelka, D. & Hedges, R. E. Analysis of amino acid 13C abundance from human and faunal bone collagen using liquid chromatography/isotope ratio mass spectrometry. Rapid Commun. Mass Spectrom. 20, 2761–2768. https://doi.org/10.1002/rcm.2651 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Raghavan, M., McCullagh, J. S., Lynnerup, N. & Hedges, R. E. Amino acid δ13C analysis of hair proteins and bone collagen using liquid chromatography/isotope ratio mass spectrometry: Paleodietary implications from intra-individual comparisons. Rapid Commun. Mass Spectrom. 24, 541–548. https://doi.org/10.1002/rcm.4398 (2010).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Newsome, S. D., Wolf, N., Peters, J. & Fogel, M. L. Amino acid δ13C analysis shows flexibility in the routing of dietary protein and lipids to the tissue of an omnivore. Integr. Comp. Biol. 54, 890–902. https://doi.org/10.1093/icb/icu106 (2014).CAS 
Article 
PubMed 

Google Scholar 
Walton, M. J. & Cowey, C. B. Aspects of intermediary metabolism in salmonid fish. Comp. Biochem. Physiol. 73B, 59–79 (1982).CAS 

Google Scholar 
Fernandes, R., Nadeau, M.-J. & Grootes, P. M. Macronutrient-based model for dietary carbon routing in bone collagen and bioapatite. Archaeol. Anthropol. Sci. 4, 291–301. https://doi.org/10.1007/s12520-012-0102-7 (2012).Article 

Google Scholar 
Ohkouchi, N., Ogawa, N. O., Chikaraishi, Y., Tanaka, H. & Wada, E. Biochemical and physiological bases for the use of carbon and nitrogen isotopes in environmental and ecological studies. Prog. Earth Planet Sci. 2, 1–17. https://doi.org/10.1186/s40645-015-0032-y (2015).ADS 
Article 

Google Scholar 
Wu, G. & Morris, M. Arginine metabolism: Nitric oxide and beyond. Biochem. J. 336, 1–17 (1998).CAS 
Article 

Google Scholar 
Metges, C. C., Petzke, K. J. & Henning, U. Gas chromatography/combustion/isotope ratio mass spectrometric comparison of N-acetyl- and N-pivaloyl amino acid esters to measure 15N isotopic abundances in physiological samples : A pilot study on amino acid synthesis in the upper gastro-intestinal tract of minipigs. J. Mass Spectrom. 31, 367–376 (1996).ADS 
CAS 
Article 

Google Scholar 
Dunn, P. J., Honch, N. V. & Evershed, R. P. Comparison of liquid chromatography-isotope ratio mass spectrometry (LC/IRMS) and gas chromatography-combustion-isotope ratio mass spectrometry (GC/C/IRMS) for the determination of collagen amino acid δ13C values for palaeodietary and palaeoecological reconstruction. Rapid Commun. Mass Spectrom. 25, 2995–3011. https://doi.org/10.1002/rcm.5174 (2011).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Ayayee, P. A., Jones, S. C. & Sabree, Z. L. Can (13)C stable isotope analysis uncover essential amino acid provisioning by termite-associated gut microbes?. PeerJ 3, e1218. https://doi.org/10.7717/peerj.1218 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Ayayee, P. A., Larsen, T. & Sabree, Z. Symbiotic essential amino acids provisioning in the American cockroach, Periplaneta americana (Linnaeus) under various dietary conditions. PeerJ 4, e2046. https://doi.org/10.7717/peerj.2046 (2016).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Larsen, T. et al. The dominant detritus-feeding invertebrate in Arctic peat soils derives its essential amino acids from gut symbionts. J. Anim. Ecol. 85, 1275–1285. https://doi.org/10.1111/1365-2656.12563 (2016).Article 
PubMed 

Google Scholar 
Romero-Romero, S., Miller, E. C., Black, J. A., Popp, B. N. & Drazen, J. C. Abyssal deposit feeders are secondary consumers of detritus and rely on nutrition derived from microbial communities in their guts. Sci. Rep. 11, 12594. https://doi.org/10.1038/s41598-021-91927-4 (2021).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
McCullagh, J. S. Mixed-mode chromatography/isotope ratio mass spectrometry. Rapid Commun. Mass Spectrom. 24, 483–494. https://doi.org/10.1002/rcm.4322 (2010).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Tsai, Y. et al. Histamine contents of fermented fish products in Taiwan and isolation of histamine-forming bacteria. Food Chem. 98, 64–70. https://doi.org/10.1016/j.foodchem.2005.04.036 (2006).CAS 
Article 

Google Scholar 
Landete, J. M., De Las Rivas, B., Marcobal, A. & Munoz, R. Updated molecular knowledge about histamine biosynthesis by bacteria. Crit. Rev. Food Sci. Nutr. 48, 697–714. https://doi.org/10.1080/10408390701639041 (2008).CAS 
Article 
PubMed 

Google Scholar 
Kanki, M., Yoda, T., Tsukamoto, T. & Baba, E. Histidine decarboxylases and their role in accumulation of histamine in tuna and dried saury. Appl. Environ. Microbiol. 73, 1467–1473. https://doi.org/10.1128/AEM.01907-06 (2007).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Fernandez-Salguero, J. & Mackie, I. M. Histidine metabolism in mackerel (Scomber scombrus). Studies on histidine decarboxylase activity and histamine formation during storage of flesh and liver under sterile and non-sterile conditions. J. Fd Technol. 14, 131–139 (1979).CAS 
Article 

Google Scholar 
Sánchez-Muros, M.-J., Barroso, F. G. & Manzano-Agugliaro, F. Insect meal as renewable source of food for animal feeding: A review. J. Clean. Prod. 65, 16–27. https://doi.org/10.1016/j.jclepro.2013.11.068 (2014).CAS 
Article 

Google Scholar 
Khan, M. A. Histidine requirement of cultivable fish species: A review. Oceanogr Fish Open Access J. 8, 1–7. https://doi.org/10.19080/ofoaj.2018.08.555746 (2018).Article 

Google Scholar 
Hatch, K. A. in Comparative Physiology of Fasting, Starvation, and Food Limitation Ch. Chapter 20, 337–364 (2012).Bertinetto, C., Engel, J. & Jansen, J. ANOVA simultaneous component analysis: A tutorial review. Anal. Chim. Acta X 6, 100061. https://doi.org/10.1016/j.acax.2020.100061 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Nogales-Mérida, S. et al. Insect meals in fish nutrition. Rev. Aquac. 11, 1080–1103. https://doi.org/10.1111/raq.12281 (2018).Article 

Google Scholar 
Thongprajukaew, K., Pettawee, S., Muangthong, S., Saekhow, S. & Phromkunthong, W. Freeze-dried forms of mosquito larvae for feeding of Siamese fighting fish (Betta splendens Regan, 1910). Aquac. Res. 50, 296–303. https://doi.org/10.1111/are.13897 (2018).CAS 
Article 

Google Scholar 
Jackson, G. P., An, Y., Konstantynova, K. I. & Rashaid, A. H. Biometrics from the carbon isotope ratio analysis of amino acids in human hair. Sci. Justice 55, 43–50. https://doi.org/10.1016/j.scijus.2014.07.002 (2015).Article 
PubMed 

Google Scholar 
Werner, R. A. & Brand, W. A. Referencing strategies and techniques in stable isotope ratio analysis. Rapid. Commun. Mass Spectrom. 15, 501–519. https://doi.org/10.1002/rcm.258 (2001).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Marks, R. G. H., Jochmann, M. A., Brand, W. A. & Schmidt, T. C. How to couple LC-IRMS with HRMS─a proof-of-concept study. Anal. Chem. 94, 2981–2987 (2022).CAS 
Article 

Google Scholar 
Lynch, A. H., McCullagh, J. S. & Hedges, R. E. Liquid chromatography/isotope ratio mass spectrometry measurement of δ13C of amino acids in plant proteins. Rapid Commun. Mass Spectrom. 25, 2981–2988. https://doi.org/10.1002/rcm.5142 (2011).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Falco, F., Stincone, P., Cammarata, M. & Brandelli, A. Amino acids as the main energy source in fish tissues. Aquac. Fish Stud. 3, 1–11 (2020).
Google Scholar  More

Ma, Z., Xue, L. & Du, S. Theory and Technology of High-Efficient Use of Water and Fertilizer for Watermelon and Melon in Gravel-Mulched Field 59–62 (Science Press, Beijing, 2018).
Google Scholar 
Qiu, Y., Xie, Z., Wang, Y., Malhi, S. S. & Ren, J. Long-term effects of gravel—Sand mulch on soil organic carbon and nitrogen in the Loess Plateau of northwestern China. J. Arid. Land 7, 46–53. https://doi.org/10.1007/s40333-014-0076-7 (2015).Article 

Google Scholar 
Zhang, K., Zhang, W., Tan, L., An, Z. & Zhang, H. Effects of gravel mulch on aeolian transport: A field wind tunnel simulation. J. Arid. Land 7, 296–303. https://doi.org/10.1007/s40333-015-0121-1 (2015).Article 

Google Scholar 
Yamanaka, T., Inoue, M. & Kaihotsu, I. Effects of gravel mulch on water vapor transfer above and below the soil surface. Agric. Water Manag. 67, 145–155. https://doi.org/10.1016/j.agwat.2004.01.002 (2004).Article 

Google Scholar 
Wang, J., Xie, Z., Guo, Z. & Wang, Y. Simulating the effect of gravel-sand mulched field degradation on soil temperature and evaporation. J. Desert Res. 30, 6 (2010).
Google Scholar 
Kaseke, K. F. et al. The effects of desert pavements (gravel mulch) on soil micro-hydrology. Pure Appl. Geophys. 169, 873–880. https://doi.org/10.1007/s00024-011-0367-2 (2012).ADS 
Article 

Google Scholar 
Inagaki, M. N. How does a stone mulch increase transpiration and grain yield in wheat under soil water deficit stress?. Cereal Res. Commun. 40, 486–493 (2012).CAS 
Article 

Google Scholar 
Abdelfattah, M. A. Pedogenesis, land management and soil classification in hyper-arid environments: Results and implications from a case study in the United Arab Emirates. Soil Use Manag. 29, 279–294 (2013).Article 

Google Scholar 
Lightfoot, D. The cultural ecology of Puebloan pebble-mulch gardens. Hum. Ecol. 21, 115–143. https://doi.org/10.1007/BF00889356 (1993).Article 

Google Scholar 
Graf, A., Kuttler, W. & Werner, J. Mulching as a means of exploiting dew for arid agriculture?. Atmos. Res. 87, 369–376 (2008).Article 

Google Scholar 
Shao Ping, D. U., Ma, Z. M. & Xue, L. Distribution characteristics of soil aggregates and their associated organic carbon in gravel-mulched land with different cultivation years. Ying Yong Sheng Tai Xue Bao 28, 1619–1625 (2017).
Google Scholar 
Zhong-Ming, M. A., Shao-Ping, D. U. & Xue, L. Influences of sand-mulching years on soil temperature, water content, and growth and water use efficiency of watermelon. J. Desert Res. 33, 1433–1439 (2013).
Google Scholar 
Pang, L. et al. Effect of different gravel mulched years on soil microflora and physicochemical properties in gravel-sand mulched field. Agric. Res. Arid Areas (2017).Hao, H. et al. Effects of gravel-sand mulching on soil bacterial community and metabolic capability in the semi-arid Loess Plateau, China. World J. Microbiol. Biotechnol. 33, 209 (2017).PubMed 
Article 
CAS 

Google Scholar 
Gregorich, E. G., Carter, M. R., Angers, D. A., Monreal, C. & Ellert, B. H. Towards a minimum data set to assess soil organic matter quality in agricultural soils. Can. J. Soil Sci. 74, 367–385 (1994).CAS 
Article 

Google Scholar 
Jandl, R. et al. Current status, uncertainty and future needs in soil organic carbon monitoring. Sci. Total Environ. 468–469, 376–383 (2014).ADS 
PubMed 
Article 
CAS 

Google Scholar 
Lal, R. Enhancing crop yields in the developing countries through restoration of the soil organic carbon pool in agricultural lands. Land Degrad. Dev. 17, 197–209 (2010).Article 

Google Scholar 
Razafimbelo, T. M. et al. Aggregate associated-C and physical protection in a tropical clayey soil under Malagasy conventional and no-tillage systems. Soil Tillage Res. 98(2), 140–149 (2007).Article 

Google Scholar 
Hongbing, Z. et al. Effect of long-term tillage on soil aggregates and aggregate-associated carbon in black soil of Northeast China. PLoS ONE 13, e0199523 (2018).Article 
CAS 

Google Scholar 
Bajracharya, R. M., Lal, R. & Kimble, J. M. Soil organic carbon distribution in aggregates and primary particle fractions as influenced by erosion phases and landscape position. In Soil Processes & the Carbon Cycle (eds Lal, R. et al.) (CRC Press, 1998).
Google Scholar 
Sekaran, U., Sagar, K. L. & Kumar, S. Soil aggregates, aggregate-associated carbon and nitrogen, and water retention as influenced by short and long-term no-till systems. Soil Tillage Res. 208, 104885 (2020).Article 

Google Scholar 
Tang, X., Liu, S., Liu, J. & Zhou, G. Effects of vegetation restoration and slope positions on soil aggregation and soil carbon accumulation on heavily eroded tropical land of Southern China. J. Soils Sediments 10, 505–513 (2010).CAS 
Article 

Google Scholar 
Wang, Y. et al. 23-Year manure and fertilizer application increases soil organic carbon sequestration of a rice–barley cropping system. Biol. Fertil. Soils 51(5), 583–591 (2015).Article 

Google Scholar 
Zhou, H., Fang, H., Mooney, S. J. & Peng, X. Effects of long-term inorganic and organic fertilizations on the soil micro and macro structures of rice paddies. Geoderma 266, 66–74 (2016).ADS 
CAS 
Article 

Google Scholar 
Congreves, K. A., Hooker, D. C., Hayes, A., Verhallen, E. A. & Eerd, L. V. Interaction of long-term nitrogen fertilizer application, crop rotation, and tillage system on soil carbon and nitrogen dynamics. Plant Soil 410, 113–127 (2017).CAS 
Article 

Google Scholar 
Tang, H., Xiao, X., Chao, L., Ke, W. & Pan, X. Impact of long-term fertilization practices on the soil aggregation and humic substances under double-cropped rice fields. Environ. Sci. Pollut. Res. Int. 25, 11034–11044 (2018).CAS 
PubMed 
Article 

Google Scholar 
Rong, Y., Su, Y. Z., Wang, T. & Qin, Y. Effect of chemical and organic fertilization on soil carbon and nitrogen accumulation in a newly cultivated farmland. J. Integr. Agric. 15, 658–666 (2019).
Google Scholar 
Chen, Y. et al. Rotation and organic fertilizers stabilize soil water-stable aggregates and their associated carbon and nitrogen in flue-cured tobacco production. J. Soil Sci. Plant Nutr. 20, 192–205 (2020).Article 
CAS 

Google Scholar 
Li, T. et al. Contrasting impacts of manure and inorganic fertilizer applications for nine years on soil organic carbon and its labile fractions in bulk soil and soil aggregates. CATENA 194, 104739. https://doi.org/10.1016/j.catena.2020.104739 (2020).CAS 
Article 

Google Scholar 
Hati, K. M. et al. 50 Years of continuous no-tillage, stubble retention and nitrogen fertilization enhanced macro-aggregate formation and stabilisation in a Vertisol. Soil Tillage Res. 214, 105163. https://doi.org/10.1016/j.still.2021.105163 (2021).Article 

Google Scholar 
Ma, P. et al. Macroaggregation is promoted more effectively by organic than inorganic fertilizers in farmland ecosystems of China—A meta-analysis. Soil Tillage Res. 221, 105394. https://doi.org/10.1016/j.still.2022.105394 (2022).Article 

Google Scholar 
Lu, R. Chemical Analysis Methods of Soil and Agriculture (China Agricultural Science and Technology Press, 2000).
Google Scholar 
Dorodnikov, M., Blagodatskaya, E., Blagodatsky, S., Marhan, S. & Kuzyakov, Y. Stimulation of microbial extracellular enzyme activities by elevated CO2 depends on soil aggregate size. Glob. Change Biol. 15, 1603–1614 (2010).ADS 
Article 

Google Scholar 
Kemper, W. D. & Rosenau, R. C. Aggregate stability and size distribution. In Methods of Soil Analysis: Part 1 Physical and Mineralogical Methods. Agronomy Monograph No. 9, ASA and SSSA 2nd edn (ed. Klute, A.) 425–442 (Wiley, 1986).
Google Scholar 
Jones, J. B. Laboratory Guide for Conducting Soil Tests and Plant Analysis (Nurse Educ, 2001).Book 

Google Scholar 
Blair, G., Lefroy, R. & Lisle, L. Soil carbon fractions based on their degree of oxidation, and the development of a carbon management index for agricultural systems. Aust. J. Agric. Res. 46, 393–406 (1995).Article 

Google Scholar 
Li, W., Zheng, Z., Li, T. & Liu, M. Distribution characteristics of soil aggregates and its organic carbon in different tea plantation age. Acta Ecol. Sin. 34, 6326–6336 (2014).
Google Scholar 
Matos, E. S., Freese, D., Böhm, C., Quinkenstein, A. & Hüttl, R. Organic matter dynamics in reclaimed lignite mine soils under Robinia pseudoacacia L. plantations of different ages in Germany. Commun. Soil Sci. Plant Anal. 43, 745–755 (2012).CAS 
Article 

Google Scholar 
Zádorová, T., Jakšík, O., Kodešová, R. & Penížek, V. Influence of terrain attributes and soil properties on soil aggregate stability. Soil Water Res. 6, 111–119 (2011).Article 

Google Scholar 
Sajjadi, S. A. & Mahmoodabadi, M. Aggregate breakdown and surface seal development influenced by rain intensity, slope gradient and soil particle size. Solid Earth 6(1), 311–321 (2015).ADS 
Article 

Google Scholar 
Yang, W. et al. Mechanical properties and soil stability affected by fertilizer treatments for an Ultisol in subtropical China. Plant Soil 363(1), 157–174 (2013).CAS 
Article 

Google Scholar 
Lal, R. Soil health and carbon management. Food Energy Secur. 5, 212–222 (2016).Article 

Google Scholar 
Jagadamma, S., Lal, R., Hoeft, R. G., Nafziger, E. D. & Adee, E. A. Nitrogen fertilization and cropping systems effects on soil organic carbon and total nitrogen pools under chisel-plow tillage in Illinois. Soil Tillage Res. 95, 348–356 (2007).Article 

Google Scholar 
Yang, X. M. et al. Long-term effects of fertilization on soil organic carbon changes in continuous corn of northeast China: RothC model simulations. Environ. Manage. 32, 459–465 (2003).CAS 
PubMed 
Article 

Google Scholar 
Russell, A. E., Laird, D. A., Parkin, T. B. & Mallarino, A. P. Impact of nitrogen fertilization and cropping system on carbon sequestration in Midwestern mollisols. Soil Sci. Soc. Am. J. 69, 413–422 (2005).ADS 
CAS 
Article 

Google Scholar 
Hartwig, N. L. Cover crop and living mulches. Weed Sci. 50, 688–699 (2002).CAS 
Article 

Google Scholar 
Yu, H. et al. Effects of long-term compost and fertilizer application on stability of aggregate-associated organic carbon in an intensively cultivated sandy loam soil. Biol. Fertil. Soils 48, 325–336 (2012).CAS 
Article 

Google Scholar 
Li, C., Yan, L. & Tang, L. The effects of long-term fertilization on the accumulation of organic carbon in the deep soil profile of an oasis farmland. Plant Soil 369, 645–656 (2013).CAS 
Article 

Google Scholar 
Puget, P., Chenu, C. & Balesdent, J. Dynamics of soil organic matter associated with particle-size fractions of water-stable aggregates. Eur. J. Soil Sci. 51, 595–605 (2000).Article 

Google Scholar 
Ashman, M. R., Hallett, P. D. & Brookes, P. C. Are the links between soil aggregate size class, soil organic matter and respiration rate artefacts of the fractionation procedure?. Soil Biol. Biochem. 35, 435–444 (2003).CAS 
Article 

Google Scholar 
Six, J., Elliott, E. T. & Paustian, K. Soil macroaggregate turnover and microaggregate formation: a mechanism for C sequestration under no-tillage agriculture. Soil Biol. Biochem. 32, 2099–2103 (2000).CAS 
Article 

Google Scholar 
Chivenge, P. P., Murwira, H. K., Giller, K. E., Mapfumo, P. & Six, J. Long-term impact of reduced tillage and residue management on soil carbon stabilization: Implications for conservation agriculture on contrasting soils. Soil Tillage Res. 94, 328–337 (2006).Article 

Google Scholar 
Kölbl, A. & Knabner, I. K. Content and composition of free and occluded particulate organic matter in a differently textured arable Cambisol as revealed by solid-state 13C NMR spectroscopy. J. Plant Nutr. Soil Sci. 167, 45–53 (2004).Article 
CAS 

Google Scholar 
Tong, L. et al. Response of organic carbon fractions and microbial community composition of soil aggregates to long-term fertilizations in an intensive greenhouse system. J. Soils Sediments 20, 641–652. https://doi.org/10.1007/s11368-019-02436-x (2020).CAS 
Article 

Google Scholar 
Li, C., Li, Y. & Tang, L. The effects of long-term fertilization on the accumulation of organic carbon in the deep soil profile of an oasis farmland. Plant Soil 369, 645–656. https://doi.org/10.1007/s11104-013-1605-4 (2013).CAS 
Article 

Google Scholar 
Su, Y. Z., Wang, F., Suo, D. R., Zhang, Z. H. & Du, M. W. Long-term effect of fertilizer and manure application on soil-carbon sequestration and soil fertility under the wheat–wheat–maize cropping system in northwest China. Nutr. Cycl. Agroecosyst. 75, 285–295 (2006).CAS 
Article 

Google Scholar 
Du, S. P., Ma, Z. M. & Xue, L. Effect of manure combined with chemical fertilizers on fruit yield, fruit quality and water and nitrogen use efficiency in watermelon grown in gravel-mulched field. J. Fruit Sci. 37, 10. https://doi.org/10.13925/j.cnki.gsxb.20190380 (2020).CAS 
Article 

Google Scholar 
Zhengchao, Z., Zhuoting, G., Zhouping, S. & Fuping, Z. Effects of long-term repeated mineral and organic fertilizer applications on soil organic carbon and total nitrogen in a semi-arid cropland. Eur. J. Agron. 45, 20–26. https://doi.org/10.1016/j.eja.2012.11.002 (2013).CAS 
Article 

Google Scholar 
Lv, W. et al. Effects of organic fertilizers on continuous cropping watermelon growth and soil microflora. Acta Agric. Shanghai 22, 96–98 (2006).
Google Scholar 
Zhong, W. et al. The effects of mineral fertilizer and organic manure on soil microbial community and diversity. Plant Soil 326, 511–522 (2010).CAS 
Article 

Google Scholar  More

Historical changes and current status of the Dongzhaigang mangrove areaBased on the literature and remote sensing data, we calculated the changes in the area of mangrove forests in Dongzhaigang since the 1950s presented in Fig. 2. In the last 60 years, the area of mangrove forests in Dongzhaigang has experienced large fluctuations mainly due to human destruction and protection activities such as mariculture reclamation, cofferdams, and restoration: it decreased from 3416 hm2 in 195617 to 3213 hm2 in 195919,29 and then decreased sharply to 1733 hm2 in 1983 and to 1537 hm2 in 198720,30. Since the establishment of the national nature reserve in 1986, the decline in area of Dongzhaigang mangrove has stopped19, which are now protected and restored owing to the law and regulations that prohibit human activities from destroying the mangrove resource. In 1988, the area was restored to 1809 hm2, and since the 1990s, it has no longer decreased, remaining constant at approximately 1711 hm2 (in the range of 1575–1812 hm2) based on the literature)18,20,31,32,33,34 (Fig. 2). The area of the Dongzhaigang mangrove forest in 2019 was estimated to be 1842 hm2 based on the latest 2 m resolution remote sensing data21. Hence, we wonder how SLR has impacted Dongzhaigang mangrove in the past decades. However, it is very difficult to analyze how SLR has historically impacted the spatial changes in the Dongzhaigang mangrove; the same can be said regarding the influence of human activities, such as destruction before mid-1980s and protection after 1990s. However, the dynamic changes among low plant edges in the intertidal zone can be used to analyze the impact of natural driving forces such as SLR35, based on the latest remote sensing data for the period of 1986–2020. Thus, we analyzed the dynamic changes in low mangrove edges (hereafter, the edges), which are mainly impacted by natural impact drivers, as shown in Fig. 3. The dynamic low mangrove edges represented by 1986, 2000, and 2020 reveal the changes in spatial distribution of Dongzhaigang mangrove. As shown in Fig. 3. Most of the edges along the coast of Dongzhaigang between 1986 and 2020 migrated landward, but not significantly. However, if we look at the changes in detail, some edges such as those in Daoxue, Sanjiang (purple circles in Figs. 3a,b–d,e–g) more clearly retreated landward compared to other places. Besides, some edges of Luodou along the northeastern coast of Dongzhaigang outside the reserve and an unnamed small island (pruple circles in Fig. 3a,h–j) also migrated landward very distinctly. On the contrary, the two smaller shore lines (black circles) in the northern part of Yangfeng and Daxue districts showed seaward expansion (Fig. 3a).Figure 2Changes in the mangrove area in Dongzhaigang from 1956 to 2019. The equation in the upper-right-hand corner of the plot refers to the fitting equation of historical changes in the total area of Dongzhaigang mangrove.Full size imageFigure 3The dynamic changes in low mangrove edges in Dongzhaigang from 1986 to 2020. Maps generated in ArcMap v10.0 (https://www.esri.com/en-us/home).Full size imageVertical rate of sediment accretion in mangrove wetlandsThe vertical rate of sediment accumulation in mangrove wetlands can reflect whether the mangroves can adjust the soil surface elevation change through sediment trapping to adapt to SLR6,11. The vertical sediment accretion rates at two sites of Dongzhaigang mangrove (i.e., Linshi and Daoxue villages in Fig. 1b) can be obtained from historical documents, which are 0.41 cm year−1 at LS and 0.64 cm year−1 at DX, respectively27,28. Since historical data may not be enough to reflect the vertical sediment accretion rates in time and space, we conducted a supplementary investigation on the sediment accumulation rates at site HG in Yanfeng and SJ site in Sanjiang farms, respectively (Fig. 1b), based on the assumption that they can reflect the sediment supplies from main reivers such as Yanfeng West River and Yanzhou River, respectively. Sediment accretion rates measured using 210Pbex specific activity in the cores from sites HG and SJ showed that 210Pbex decayed exponentially with increasing depth, and the R2 values of both cores were approximately 0.80 after curve fitting. This analysis resulted in vertical sediment accretion rates of 0.53 and 0.40 cm year−1 at HG and SJ, respectively (Fig. 4). Therefore, the locations of sediment cores at sites LS, DX, HG, and SJ can basically represent the whole Dongzhaigang mangrove forest area.Figure 4210Pbex activity profiles in selected cores such as from (a) station HG and (b) station SJ.Full size imageRate of relative sea level rise in Dongzhaigang mangroveThe global mean sea level (GMSL) is accelerating due to global warming-induced thermal expansion of the oceans and melting of land-based glaciers and ice caps into the sea36. Between 1901 and 2010, the GMSL rose by 0.19 m9. Coastal China is among the regions that experience the highest levels of SLR23. The rate of RSLR along China’s coast from 1980 to 2019 was 3.4 mm year−1, higher than the global average23. In the future, under the premise of increasing anthropogenic GHG emissions, global sea levels will rise rapidly, and it is projected that the GMSL may rise by 0.84 m (0.61–1.10 m) relative to the current levels by the end of the twenty-first century9. Based on the observations from the tide gauge stations in the Haikou area and model data from the Coupled Model Intercomparison Projection 5 (CMIP5), the rate of RSLR around Dongzhaigang reached 4.6 mm year−1 from 1980 to 2018. This rate is much higher than the global and China’s average values23,25 and will likely accelerate further in the future. Based on the results of the CMIP5 model simulations under different GHG emission scenarios24, the RSLR in coastal Haikou waters, including in Dongzhaigang, is expected to be significant by 2030, 2050, and 2100 for the low, intermediate, and very high GHG emission scenarios RCPs 2.6, 4.5, and 8.5, respectively (Table 1, Fig. 5). Under RCPs 2.6, 4.5, and 8.5, the sea level will rise by 65 (42–90, likely range), 75 (51–102, likely range), and 96 (70–125, likely range) cm by 2100, respectively, with the average RSLR rates of 6.84 (4.42–9.47, likely range), 7.89 (5.37–10.74, likely range), and 10.1 (7.37–13.12, likely range) mm year−1, respectively.Table 1 Estimated coastal relative sea level rise (cm) and its rate (mm year−1) in the Haikou area under different GHG emission scenarios (data from Kopp et al.24).Full size tableFigure 5Historical and future relative sea level changes along coastal Dongzhaigang, Haikou City from 1980 to 2100; the 5–95% uncertainty ranges are shaded for RCPs 2.6, 4.5, and 8.5, respectively.Full size imageImpact of relative sea level rise on Dongzhaigang mangroveMangroves cannot easily adapt to rising sea levels if the rate of GMSL rise exceeds 6.1 mm year−1 ( > 90% probability, very likely), whereas the survival threshold for mangroves is extremely likely to be exceeded ( > 95% probability, extremely likely) when the rate of GMSL exceeds 7.6 mm year−17. Although these values are based on global levels7, they still reflect the threat of SLR to local mangroves. In view of this, we further analyzed the potential impact and risks to Dongzhaigang mangrove from future SLR under different climate scenarios.Based on the predicted future rates of SLR under RCPs 2.6, 4.5, and 8.5 and on the vertical sediment accretion rates of Dongzhaigang mangrove wetlands, the mangroves are likely to be affected by rising sea levels by 2030, 2050, and 2100, respectively (Table 2, Fig. 6). Under the low GHG emission scenario (RCP 2.6), the area of the Dongzhaigang mangrove forest will only experience a small reduction: 16.40% (1.20–16.95%, likely range), 302 hm2 (22–312 hm2, likely range); 16.73% (1.20–17.82%, likely range), 308 hm2 (22–328 hm2, likely range); and 17.60% (1.14–31.02%, likely range), 324 hm2 (21–571 hm2, likely range) by 2030, 2050, and 2100, respectively (Table 2, Fig. 6a). This is because the vertical sediment accretion rate of Dongzhaigang mangrove will remain largely constant with increasing RSLR rate. Moreover, it should be noted that compared with 2030, the increase areas of mangroves inundation caused by SLR will be small by 2050 under three RCPs scenarios (Table 2). In contrast, under the intermediate and very high GHG emission scenarios (RCPs 4.5 and 8.5), Dongzhaigang mangrove is expected to be more significantly affected by SLR. Under RCP 4.5, 26.56% (16.19–40.74%, likely range) or 489 hm2 (298–750 hm2, likely range) of mangrove forest will likely be lost by the end of the century (Table 3, Fig. 6b). Under RCP 8.5, it is projected that 31.99% (18.14–50.73%, likely range) or 589 hm2 (334–934 hm2, likely range) of mangrove forest will be lost by 2100 (Table 2, Fig. 6c). Therefore, under RCPs 4.5 and 8.5, the impact of SLR on mangrove wetlands by 2100 is much higher than that of RCP 2.6, and is likely to result in  > 26% of mangroves being lost, whereas under RCP 2.6, only 17% of mangroves are likely to be lost.Table 2 Area (hm2) and percentage of future mangrove loss in Dongzhaigang under different climate scenarios (RCPs 2.6, 4.5, and 8.5) (likely ranges).Full size tableFigure 6Potential loss of mangrove forests in Dongzhaigang under different climate scenarios (RCPs 2.6, 4.5, and 8.5). Maps generated in ArcMap v10.0 (https://www.esri.com/en-us/home).Full size imageTable 3 Core stations and depths.Full size tableUnder RCP 2.6, the rate of RSLR around Dongzhaigang will reach 0.72 cm year−1 in 2030 and then decrease in 2050 and 2080 to 0.69 and 0.68 cm year−1, respectively (Table 1). However, under RCP 4.5 (8.5), by 2030, 2050, and 2100, the rate of RSLR will reach 0.72 (0.72), 0.73 (0.80), and 0.79 (10.1) cm year−1, respectively. By 2100, some mangroves in the northern part of Tashi village, the eastern part of Yanfeng, the northern part of Daoxue Village, and the northeastern part of the Sanjiang farm will likely be lost owing to SLR, and other coastal wetlands will also be impacted. Since the rate of RSLR around Dongzhaigang is higher than the global average survival threshold for mangroves (i.e., the SLR rate exceeds 7.0 mm year−1), the Dongzhaigang mangrove will be significantly affected by SLR, with a potential loss of 31–32%; however, the survival threshold will not increase (Table 2, Fig. 6). More

Site selectionOur study focused on the mixed-conifer zone of the Sierra Nevada which is dominated by Ponderosa pine (Pinus ponderosa), Jeffrey pine (Pinus jeffreyi), Incense cedar (Calocedrus decurrens), White fir (Abies concolor) and California black oak (Quercus kelloggii). Common shrubs include Deer brush (Caenothus integerrimus), Mountain whitethorn (Caenothus cordulatus), Greenleaf manzanita (Arctostaphylos patula) and Prostrate ceanothus (Caenothus prostrates). Study sites in the Eldorado National Forest near Placerville, CA (38°45′N 120°20′W), were in and near the area burned during the King Fire (Fig. 2).We sampled sites in the mixed-conifer zone between 4000–6000 ft in three burn categories: unburned, low-to-moderate severity and high severity. We selected sites that occurred in similar pre-burn habitat (moderate to dense conifer forest) based on remotely-sensed vegetation class data from the California Wildlife Habitat Relationships program (CWHR)16. No site experienced wildfire or controlled burns in the preceding century23. Burn categories were based on remote sensing relative differenced Normalized Burn Ratio (RdNBR) canopy cover calibrated for the Sierra Nevada24. We validated these burn categories as meaningful and discrete with remote sensing (immediately post fire) and field data (3 years post-fire). Monitoring Trends in Burn Severity (MTBS) maps classify burn severity with Landsat reflectance imagery of pre-fire and post-fire conditions at 30-m resolution25. MTBS assigns pixels a value based on burn severity: 0 = outside the burn boundary; 1 = unburned-low severity within the fire perimeter; 2 = low severity; 3 = moderate severity; 4 = high severity). At each site, we determined the mean MTBS value of pixels in the small-mammal trapping grid (90 × 90 m) and a 225 m buffer (the largest home range diameter for small mammals we captured)26 extended on all sides. Immediately after the fire, the three burn categories differed significantly in tree cover and remote sensing scores of burn severity25 and are still different in tree surveys three years later.We sampled 27 sites, each 4 ha in size. To minimize site-specific influences, we paired sites across treatments to account for elevation, slope, pre-burn vegetation characteristics, pre-burn forest management, ownership, and soil type. Each burn category received nine sites. We blocked sites across burn severities (one site in each burn category per block). Site locations were chosen for accessibility, but all sites were at least 50 m from access roads and at least 200 m apart (sites were >1 km apart on average). We also excluded areas that experienced or were scheduled to experience salvage logging post fire. Occasionally, field conditions at a site location did not align with remotely sensed data classifications, a site was not large enough for homogeneous sampling plots or was dominated by slopes >30 degrees. In these cases, sites were moved to nearby locations that satisfied the site selection criteria.Sampling designIn this study we report the body size, abundance and biomass density of plants and free-living animals using a variety of methods in three different fire severities. Each site consisted of a 200 m × 200 m plot (4 ha) around which sampling methods were organized (Fig. 3). We collected data for organisms with 19 different sampling methods that were scaled to the abundance and body size of targeted organisms. Some methods were not performed at all sites, and some methods were pooled across sites (within treatments) due to low sampling success. Below, we explain each sampling method and detail any variation in its application across sites. To minimize seasonal effects, we sampled all sites in a block at the same time, over 4–5 continuous days, between late June and early September 2017. Weather was consistently hot and dry during this period with no pronounced variation. All animal survey procedures were approved by UC Davis IACUC (protocol number CA-17-451736) and carried out under CDFW permit #SC-3638.Species inclusionTo evaluate the effects of wildfire-burn severity on community structure we assembled a list of organisms and life stages encountered during our sampling of the Eldorado National Forest. Every animal in our list was broken into ontogenetic life stages. Plants were broken into constituent parts (e.g. leaves, roots, seeds) rather than stages. Every organism in our lists had at least one life stage observed by this study in the Eldorado National Forest. Not all life stages in the list were directly observed, many were inferred. These life stages were suspected to occur in the habitat but were not observed (or quantified) because (1) it was impractical (e.g. mycorrhizae), (2) our sampling methods did not capture them (e.g. larval insects in plant tissues) or (3) they could not be identified to species (e.g. arthropod eggs). Stages (not species) were excluded when they did not occur in the terrestrial habitat (e.g. aquatic larvae). Stages (or parts) were lumped when they could not be distinguished in terms of their resources or consumers (e.g. parasitoid wasp eggs, larvae and pupae). Unobserved stages were omitted from a burn category community if they were feeding (e.g. not eggs) but did not have any resources. All life stages, observed or not, received a body size estimate.This approach has several benefits. First, it fills life-cycle gaps without artificially inflating species richness (e.g. having a node labelled “bird eggs”). Second, unobserved life stages may not be major biomass contributors, but they do make important contributions to food-web structure and population dynamics. Including life stages without abundance information is useful because it allows their inclusion in analyses of network structure, informs consumer-resource body-size ratios and allows comparisons to other food datasets organized around body size, but lacking abundance estimates.A comprehensive food web required the inclusion of non-living nodes like detritus which is an important resource for many consumers. Fire creates strong spatial structure in woody detritus, so we collected mass-density information on it and partitioned detritus into types (e.g., woody vs. leaf) and size classes. We did not collect mass-density information on carcasses but treated them similarly by partitioning them into logarithmic size classes to reflect their availability to different consumers (e.g., tick carcasses vs. deer carcasses). We have also included the biotic products honeydew and dung due to their importance to certain consumer groups in the system. We have not done so here but future efforts may wish to include nutrients as resources for plants to capture gradients and competition. Finally, analyses may wish to augment our lists by explicitly including fire as a consumer/herbivore as opposed to an external force shaping treatment effects.Below we describe the methods used to quantify the richness, abundance, and biomass density of these organisms in each of our burn categories. Unless noted, sampling effort did not vary with burn category.Species resolutionMost entries in the species list represent ontogenetic stages or constituent parts. With the exception of a few nodes (i.e. algae, moss, lichen, mycorrhizae, saprophytic bacteria, saprophytic fungi and nematodes) this ecological resolution is consistent throughout the list. Taxonomic information was assigned using the Global Biodiversity Information Facility database27, but taxonomic resolution varies. Vertebrates are identified at the species level. Invertebrates, while distinguished as morphospecies, were not identified below the family level in many cases. Despite not always being identified to the lowest taxonomic category, entries represent life stages of distinct species or groups, and are accompanied by all the taxonomic information we could provide at that level of resolution.Biomass density estimationEach species observed in our sampling methods received a biomass density estimate. Biomass density was estimated as the product of mean individual body mass and density. We estimated body-mass either by weighing individuals directly, or conservatively estimating their mass volumetrically. For the latter, we measured the length (tip of abdomen to tip of head), width (max width of body excluding appendages) and depth (max depth of body excluding appendages) of individual organisms and converted into an approximate volumetric shape (e.g. ellipsoid, cylinder, hemisphere, etc.)28,29. Mass was estimated by multiplying this volume by a tissue density of 1.1 g/mL30,31,32. Body sizes for stages that were not directly observed were inferred from published records and databases. Density estimates were derived from the sampling methods discussed below and varied with burn category. Mean and standard errors were derived from the nine replicate sites within each treatment, unless otherwise stated. Due to the brief sampling window at each site, these estimates should be regarded as point estimates rather than integrated over time.Plant samplingVegetation transectsTo estimate density and biomass of trees, shrubs, and ground cover, we conducted vegetation transects at each site. Transects were located within each 200 m × 200 m plot but were offset from other sampling methods by 10 m to avoid trampling. Transect direction was chosen at random and two transects were run in parallel. Each transect was 50 m in length but width (and distance between transects and sampling area) varied with sampling method as detailed below.Ground coverTo estimate plant ground cover, we surveyed 1 m2 quadrats at 5 m intervals on each 50 m transect. Summed, this gave a pooled ground cover sample area of 20 m2 at each site. Within each quadrat we estimated the percent cover and average height of grasses, forbs, woody litter, soft-loose and soft-rooted litter, shrubs, and trees. We identified small trees, shrubs and dominant forbs (i.e., forbs with the greatest cover at a site) to species. We then used cover area and height to estimate a volume for each species or group. Volumes were later multiplied by taxon-specific measurements of mass-to-volume ratios derived from vegetation box quadrats to obtain biomass estimates for each ground cover species or group.Canopy coverTo estimate canopy cover, above each of the 20 ground cover sampling locations, we identified each tree species overhead and estimated its absolute canopy cover. This gave a pooled canopy cover sample area of 20 m2 at each site. Canopy cover estimates were incorporated into arthropod biomass density estimates from fogging.TreesTo estimate tree density, we identified and measured all trees that exceeded 15 cm diameter-at-breast-height (DBH) within a 15 m wide band extending the length of each 50 m vegetation transect. Small trees (DBH 1 cm. We estimated shrub volume (length × width × height) and biomass density by combining volume and density estimates with taxon-specific measurements of mass-to-volume ratios derived from vegetation box quadrats. We estimated above ground biomass for small trees using species-specific DBH-to-biomass conversions33.

Understory plants
To estimate the volume and cover of understory vegetation we used three-dimensional box quadrats. These box quadrats consisted of a canvas-walled rectangular box with the floor panel removed, that was placed over the understory to be surveyed. The canvas walls prevent arthropods from escaping as they were also used in arthropod sampling. Each box quadrat measured 1 m × 0.5 m × 0.5 m, covering a ground area of 0.25 m2. At each site, we placed box quadrats at three locations representative of the dominant vegetation. This gave a pooled sample area of 0.75 m2 at each site. To generate a volume estimate, within each box we estimated the percent cover and maximum height of each plant species. Leaf litter was treated separately but also quantified. Then, using a Velcroed flap as access (designed to prevent bugs from escaping), we collected and weighed all the above ground plant biomass, separated by species and type (litter, branches, leaves, etc.). A subsample of material (up to 1 kg) from each plant species and type was retained and returned to the station to extract the associated arthropods via Berlese funnels (detailed below).
For each plant species measured in the vegetation box samples, we estimated a mass-to-volume ratio. We estimated the biomass of understory plants by combining our mass-to-volume ratios with volume estimates from ground cover transects. For taxa present in ground cover transects but not vegetation boxes, or with fewer than three measurements of mass-to-volume ratio, we pooled measurements from higher taxonomic levels. For example, if there was only one measurement of species A, but two measurements for a congener, species B, we used measurements for both species A and B to represent the mass-to-volume ratio of each. In this way, we estimated mean mass-to-volume ratios for every species-type encountered in the shrub and cover transects.
Coarse woody debrisTo estimate the mass density of woody detritus we quantified it according to the method detailed in Waddell (2002). Specifically, we used the center of the vegetation transect as a line-intercept for woody detritus. Woody debris was recorded if its longitudinal axis intersected the transect line, its diameter at the point of intersection was ≧ 12.5 cm; it was ≧ 1 m long and it was not decayed to the point of disintegration. For each piece of woody detritus, we measured the diameter at both ends, the length, and the stage of decay. Length was measured only for the portion exceeding 12.5 cm in diameter. We estimated the volume (m3) for each piece of debris:$$V{m}^{3}=left.frac{(pi /8)({D}_{s}^{2}+{D}_{L}^{2})l}{10,000}right)$$
(4)
({D}_{S}^{2}) is small diameter (cm), ({D}_{L}^{2}) is the large diameter (cm) and l is the length (m)34. We converted this estimate to a volume-density (m3 ha−1) estimate:$${m}^{3}h{a}^{-1}=left(frac{pi }{2;{L}_{t}}right)left(frac{V{m}^{3}}{{l}_{i}}right)f$$
(5)
Lt is the combined transect length (100 m), lt is the length of the individual piece, f is a conversion for area (10,000 m3 ha−1)34. Finally, we converted this is a dry-weight biomass density (kg ha−1):$$kg;h{a}^{-1}=({m}^{3}h{a}^{-1})left(frac{1000,kg}{{m}^{3}}right)SpGast D$$
(6)
SpG is a specific gravity estimate and D is the correction for the state of decay34. Specific gravity can be a species-specific estimate. Woody debris is difficult to identify to species so we applied a single specific gravity estimate to all debris weighted by the relative abundance of tree species in our surveys. We used a weighted specific gravity of 0.382, estimated by multiplying the mean specific gravity of Pinaceae (0.372) by their relative abundance (0.947) and adding that to the specific gravity of Quercus (0.566) multiplied by their relative abundance (0.053). We applied the weighted means to the decay corrections in Waddell (2002).Species of uncertain identificationSome plant identifications were difficult, particularly at burned sites, where some individuals were identified to genus or family. To assign species identities to individuals classified a higher taxonomic rank at a focal site, we first randomly selected a proxy site from the nine unburned sites. This proxy approach is possible because the composition of burned and unburned sites was similar before the fire. Proxy sites were randomly selected from among all the unburned sites because the spatial grouping didn’t lend itself to maintaining distinctions between blocks. Using the species abundance distribution of the focal plant taxon at the proxy site, we randomly assigned (with replacement) plant species identities to unidentified individual plants at the focal site. We repeated this random sampling 999 times, then estimated the resulting mean species abundance and bootstrapped standard errors at the treatment level.Invertebrate samplingGiven the diversity of invertebrates in the ecosystem, and the biases associated with each invertebrate sampling method, we employed several methods to quantitatively survey the entire community.Arthropods on understory plantsTo estimate the density of arthropods on understory plants we collected them along with plant material in the box quadrats. Up to 1 kg of plant material of each species was bagged inside box quadrats and transported to the field station for processing. Insects were extracted from the plant material in Berlese funnels that were hung in the shade at ambient temperatures and powered with 60 W frosted-incandescent bulbs. All plant samples were processed in funnels for 72 hours, after which they were checked twice a day (collecting jars changed). After no additional arthropods had been collected for 24 hours samples were removed. Arthropods were placed in ethanol for later identification. After identification, arthropod densities were estimated as the product of the arthropod-to-volume ratio for each plant and the total volume of each plant estimated from transects. We collected 4543 individual arthropods from 170 morphospecies associated with plant material from box quadrats, processed in Berlese funnels.Soil-dwelling arthropodsTo estimate the density of soil-dwelling arthropods we collected soil cores. Alongside each box, we collected four cores, each 10 cm in diameter. Two shallow cores were sunk to a depth of 5 cm, giving a combined (6 replicates total) sample area of 0.0471 m2 at each site. Shallow cores, which targeted small arthropods (e.g. collembola, diplura, acari), were collected and transported back to the field station for processing in Berlese funnels. Soil samples were processed in Berlese funnels in the same manner as plant material. Two deep cores were sunk to a depth of 15 cm, giving a combined (6 replicates total) sample area of 0.0471 m2 at each site. Deep cores targeted large invertebrates (e.g. lumbricidae, myriapoda) and were processed by hand in the field. Any large invertebrates encountered were placed in ethanol for later identification. After identification, densities were estimated as the quotient of counts and area sampled. We collected 690 individual arthropods from 57 morphospecies from soil cores.Arthropods on hard substrateTo estimate the density of arthropods like wasps and spiders on tree trunks, rock and logs we visually surveyed these hard substrates. In three locations at each site, we surveyed all hard substrates within a cylinder 5 m in diameter and 2 m in height. The total sample area at each site was 58.9 m2. All invertebrates larger than 2 cm in length were collected and fixed in ethanol for later identification. After identification, hard substrate densities were estimated as the quotient of counts and sample area. Hard substrate surveys yielded 616 individual arthropods from 74 morphospecies.Understory arthropodsTo estimate the arthropod densities associated with understory shrubs at each site we supplemented box quadrats with sweep net surveys. We used nets with a 38 cm diameter opening (Bioquip) to sweep the same 5 m diameter area as the hard substrate surveys. Three replicates at each site gave a total sample area of 58.9 m2. Net sweeps were performed at a constant pace of 2 arcs per meter. All invertebrates collected were fixed in ethanol for later identification. After identification, arthropod densities were estimated as the quotient of counts and sample area. Sweep net surveys yielded 706 individual arthropods from 168 morphospecies.Canopy arthropodsTo survey arboreal arthropods, we used a thermal fogger to loft insecticidal fog into the tree canopy. While this method can sample a large area, we were hampered by permit and accessibility issues at some sites. As a result, we sampled a representative subset of 12 sites (6 unburned, 3 high severity, and 3 low-to-moderate severity) with canopy fogging. Chemical fogging remains the most widely used method for sampling canopy arthropods35,36,37. It is notoriously difficult to assess canopy arthropods with any other method35. We used an IGEBA TF-35 Thermal Fogger to disperse a pyrethrum-based insecticide (EverGreen Crop Protection EC 60-6) into the forest canopy. During each fog we dispersed 4 L of solution (a 7% concentration of insecticide in a water-based carrier-dispersant) in 10 minutes over two representative trees and one live shrub at each site. Live trees were sampled at unburned sites, dead trees were sampled at high severity sites, and one dead and one live tree were sampled at low-to-moderate severity sites. Fogging was done under low-wind conditions allowing the thermal fogger and temperature gradients to lift the fog into the canopy. White 2.25 m2 tarps were placed beneath the fogged area to collect insects. The number of tarps used varied with the size and shape of the canopy. After an optimal elapsed drop-time of 120 minutes37, arthropods were collected from tarps and placed in ethanol for later identification. To estimate site-level arthropod abundance using canopy fogging, we estimated the density per square meter of tarp for each morphospecies and habitat type (living tree, dead tree, or shrub), then multiplied these densities by the total area covered by each habitat type. The area covered by each fogged habitat type was obtained from canopy cover transects (trees) and shrub transects (shrubs). Mean and standard errors for density estimates were estimated using the 3 (or 6) sites per treatment. Canopy fogging yielded 16,353 individual arthropods from 489 morphospecies.Strong-flying insectsTo estimate the diversity and relative abundance of strong-flying arthropods, we utilized blacklight traps (Miniature Downdraft Blacklight (UV) Trap Model 912). We supplemented our sampling with blacklight traps, because strong flying arthropods (e.g. vespid wasps) may not be adequately sampled by fogging or sweep netting. Black lights were deployed at sites prior to dusk on the third night of bat acoustic surveys (see below). Each site received one trap and black lights were not deployed on rainy or windy nights. Samples were collected the following morning. To preserve their integrity for later identification, lepidoptera samples were removed and frozen, other flying insects were fixed in ethanol. Black light traps collected 5,508 individuals from 279 morphospecies.Black light samples were only used to estimate densities for those species not regularly captured with other methods. Two hundred fifteen morphospecies were only found in black light samples. Black light traps, which can sample large areas, are an efficient means of generating diversity estimates for nocturnal flying insects. However, black light traps generate activity densities rather than absolute densities, over sample areas that vary within and between nights. To convert counts to density estimates for morphospecies collected only in black lights we paired them with an analogous species (based on taxonomy and body size). These density analogues were captured at the same site in both black lights and at least one other method with explicit absolute densities. We estimated the mean relative abundance ratios of these species in black lights across all sites in which they co-occurred. We then multiplied the density of the ecological analog by these abundance ratios to estimate densities of all species that were collected only in the black light traps.Larval biomass densityTo estimate the density of larval Formicidae, lepidoptera and coleoptera, we partitioned them according to the relative abundance of their adults present in the same treatments. Larvae for many hemimetabolous invertebrates can be difficult to identify without molecular methods. When larvae could not be identified to species, they were binned into categories of Formicidae, lepidoptera and coleoptera. We then collected body size information and density estimates for these stages. Larval densities were then partitioned according to the relative abundance of adults within each treatment. Adults were eligible to receive larvae if (1) they belonged to the same taxonomic rank as the larvae, (2) they did not already have any observed larval stages, (3) there were available resources for their larvae in the treatment.Vertebrate samplingSmall mammalsTo estimate the density of mammals less than 2.5 kg we used live traps uniformly distributed over a 90 m × 90 m grid38. We placed 100 traps 10 m apart, alternating between large (7.5 cm × 9 cm × 23 cm) and extra-large Sherman traps (10 cm × 11.5 cm × 38 cm)39. Traps were baited (a mix of oats, peanut butter, bird seed and molasses) and covered with natural materials for insulation. Traps were locked open and pre-baited for three days and then operated for three consecutive nights. Traps were closed during the day due to high daytime temperatures and lack of shade, particularly in the high severity burn areas. Traps were re-opened in the late afternoon. Captured mammals were identified to species (or individuals during recaptures), weighed, measured and fit with ear tags for future identification, then released. In 7,906 trap nights (adjusted for trap failures), we had 988 captures of 11 species of small mammals.Live-trapping and marking of small mammals allowed us to estimate their abundance using spatial capture-recapture (SCR) models40. These models use individual detections in combination with detection locations (here, location of a given live trap within the grid) to estimate density while accounting for imperfect detection and variation in detection due to variability in individual exposure to the trapping grid. When data are collected across several trapping grids, as in the present study, the joint modeling of all data allows estimation of grid-level covariate effects on density41. We used this framework to analyze data from all 27 plots jointly and allowed for density to vary according to burn category of a site (unburned, low-to-moderate severity or high severity), so that, for example, all unburned sites had the same density of a given species. If a species was never caught in a burn category, we fixed its density (and consequently, its biomass) to 0 for all sites in that burn category.Because trapping data for most species was too sparse to fit species-specific models, we grouped ecologically similar species and built models that shared parameters among grouped species. In this manner, we analyzed the joint data of mice and chipmunks: Brush mouse (Peromyscus boylii), Deer mouse (Peromyscus maniculatus), Pinyon mouse (Peromyscus truei) and Western harvest mouse (Reithrodontomys megalotis), Yellow-pine chipmunks (Tamias amoenus), Long-eared chipmunks (Tamias quadrimaculatus) and Shadow chipmunks (Tamias senex).In the mouse model, Brush mouse and Pinyon mouse densities were assumed to have the same response to burn category (but not the same densities), whereas Deer mice were allowed a different response to burn category. This choice was made based on capture frequencies of species in the different burn categories. Pinyon mice and Western harvest mice had a fixed density of 0 in low-to-moderate and high severity burn sites. All species shared the same detection parameters.Even combined, the chipmunk data captures were too sparse to estimate species specific densities. We therefore treated all chipmunks like a single species to estimate overall chipmunk density, then calculated species specific densities by multiplying overall density with the proportion of individuals of a given species in the data. For example, if species A made up 50% of all individuals in the data, we would calculate density for species A by multiplying overall chipmunk density by 0.5. This model entails the assumption that the effect of burn category on density was the same for all chipmunks, which seems reasonable based on capture frequencies.We built species-specific SCR models for California ground squirrels (Otospermophilus beecheyi) and Dusky-footed woodrats (Neotoma fuscipes); for the latter we fixed density to 0 in high severity sites. Northern flying squirrels were only caught 3 times at unburned sites, so we set its density to 0 in low-to-moderate and high severity sites, and used a published density estimate from a Sierra Nevada site42 for unburned sites. Because individuals of Trowbridge’s shrew (Sorex trowbridgii) were never recaptured and had high rates of trap mortality, we estimated their abundance using non-spatial removal models43 for unburned and moderately burned sites and set their abundance to 0 at high severity sites. We transformed abundance to density by dividing it by the 8100 m2 area (90 × 90 m) of the live trapping grid.We fit SCR models in R using the package “secr” ver. 3.1.344, and removal models using the package RMark ver. 2.2.445.BirdsTo estimate the diversity and density of birds, we conducted point count surveys46. Each site had two point-count stations, spaced at least 200 m apart, that were surveyed on the same day. Counts were conducted at each site on three different days per season. To mitigate migration effects all counts were conducted between mid-June and mid-July. Bird surveys did not take place when it was raining, extremely cold (20 kmh). Wind speeds were obtained with handheld anemometers. Point counts began 15 minutes after sunrise and were completed by 10:00 AM, corresponding to when passerine birds are most active. Each survey consisted of one ten-minute count, split into two five-minute periods. Each individual bird was recorded only once over the entire count. Trained observers identified all birds to species by sight or sound and estimated the number of individuals of each species within 100 m of the sampling point. This gave us a sample area of 15,708 m2 at each site. Birds that flew through or were detected outside of the survey area or survey time were documented but not included in richness or density estimates. During 162 point-surveys we observed 1,039 birds from 52 species (107 individuals could not be identified).The repeated point counts allowed us to estimate bird abundance and density using N-mixture models47. These models use repeated counts of individuals to estimate abundance within a sampling unit (in this case, the 100 m radius point count) while accounting for imperfect detection. Because abundance estimates refer to a set area, these can be converted to densities. Because data for some bird species were sparse, we fit data of all species jointly in a community model48. In community models, each species has its own set of parameters, but species-specific parameters are modeled as coming from a common underlying distribution that is shared by the community of species (essentially, a species-level random effect). This constitutes a form of information sharing, which improves parameter estimates for data-sparse species. In our community N-mixture model, we allowed for species-specific detection probabilities, as well as species-specific effects of burn category on abundance. We further included a fixed (across species) effect of the amount of wind during a given survey on detection, as wind can impair auditory detection of birds. We implemented the community N-mixture model in a Bayesian framework using the program JAGS ver. 4.2.049 accessed through R with the package jagsUI ver. 1.5.050. JAGS fits models using Markov chain Monte Carlo (MCMC) algorithms; we ran 100,000 MCMC iterations and discarded the first 50,000 as burn-in. We checked for chain convergence using the Gelman-Rubin statistic51; all parameters had a value 2.5 kg we used camera traps. At each site, we deployed two Reconyx HC600 Hyperfire cameras, which have a no-glow infrared flash, preventing disturbance to wildlife and detection by humans. To reduce false triggers, we set cameras in shaded locations and cleared vegetation in front. Cameras were deployed along landscape features that suggested animal movement: game trails, forest openings, abandoned dirt roads, etc. Cameras were set at least 50 m apart and strapped to trees 40 cm above the ground and 1–2 m away from the edge of a trail or opening. Cameras were operated for 24 hours a day, with 3 consecutive pictures per trigger and no time lapse between triggers. All 54 camera traps were installed by early-July and retrieved in mid-September. Cameras were checked after four to six weeks for SD card and battery replacement. All pictures were reviewed manually for species identification; identified pictures of animals were organized into camera and species-specific folders for post-processing in the R computing environment ver. 3.4.352 using the package camtrapR ver. 0.99.553. Even though some bird and small mammal species were occasionally photographed by camera traps, we excluded their photo-records from further analysis. After adjusting for malfunction, cameras operated for 4,238 trapping nights over which they assembled 12,243 independent records (detections that were at least one hour apart if of the same species at the same camera) comprising 10 species of mammals >2.5 kg.Because camera-trap images do not allow identification or counting of individuals for analysis with SCR or N-mixture models, we used the Royle-Nichols (RN) occupancy modeling framework54 to estimate abundance of medium/large mammal species. Occupancy models55 use repeated species detection/non-detection data from a collection of sampling locations to estimate species occurrence probability while accounting for imperfect detection. The RN model makes use of the fact that the probability of detecting a species at a site increases with the abundance of that species at that site, allowing estimation of site-level abundance from species level detection/non-detection data. We used the R package camtrapR53 to convert raw camera trap data for each species into a binary (detected = 1, not detected = 0) location-by-occasion format. Because of their proximity, we considered both cameras at a given plot to constitute a single sampling location. We defined an occasion as 10 consecutive days of sampling. To account for malfunctioning of some cameras, we calculated effort as the number of days per occasion that each pair of cameras was functional.Because data for some species were sparse, we jointly analyzed data for all species in a community RN model (see Birds for a description of community models); we allowed for species-specific detection probabilities as well as species-specific effects of burn category on species abundance and included a fixed effect (across species) of effort on detection. The RN model returns sampling location level estimates of abundance, but in contrast to bird and reptile surveys, which were explicitly linked to a specific sampled area, the area sampled by a camera trap is not easily defined. Mammals recorded in the relatively small detection zone of a camera use much larger areas. We calculated densities of large mammals by dividing the abundance estimates from the RN model by the average home ranges of the recorded species. We fit the community RN model in a Bayesian framework as described for birds, running 30,000 MCMC iterations and discarding the first 15,000 as burn-in. All chains converged, according to the Gelman-Rubin statistic.To obtain estimates for home range size and individual mass of large mammals, we searched the scientific literature for studies conducted in similar habitat (temperate coniferous forests) and for home ranges, during the summer and fall seasons. When no information from a similar habitat was available, we used studies from forested areas. When no information for target species was available, we used information from congeners. When multiple studies were available, we calculated the mean across studies, weighted by sample size if provided. Similarly, when studies provided information separately for males and females, we calculated a weighted mean. For some large mammal species, we only found information on annual home ranges, and to approximate a summer/fall home range, we multiplied these annual ranges by 0.67.BatsTo survey bat community composition and abundance, we conducted acoustic sampling at each site. We placed a Wildlife Acoustics Songmeter SM4BAT FS echolocation detector at each site for a minimum of three consecutive nights. The detectors were set to record from sunset to sunrise when bats are most active. We attached microphones to t-posts placed in open areas, elevating them 2.5 meters from the ground to minimize signal attenuation from nearby vegetation and canopy. Bat calls were analyzed using SonoBat software to assess likely species for each file. These identifications were vetted by an experienced bat researcher, Ted Weller (USDA Forest Service). One site (B3U1) was not sampled with bat detectors, whereas two sites had one and two additional nights sampled, respectively. During 81 recorder nights, acoustic recorders captured 17,484 calls, of which 7348 were identified to 17 species.Similar to camera traps, bat calls recorded by acoustic recorders do not allow counting or identifying individuals. Therefore, we built a community RN model to estimate abundance for bats. The model structure was identical to that described for medium/large mammals, except that it did not include any covariate on detection probability. Analogous to our medium/large mammal analysis, we used literature information on average home range size to convert bat abundance to bat density. We ran the bat model for 150,000 iterations and discarded the first 75,000 iterations as burn-in. All chains converged according to the Gelman-Rubin statistic.ReptilesTo estimate community composition and density of reptiles, we conducted timed searches within the 90 m × 90 m small mammal trapping grids at each site. This gave us a sample area of 8100 m2 at each site. Reptile surveys were conducted either prior to small mammal trapping or two weeks post to minimize impacts from sampling disturbance. Reptile surveys were conducted between 8:00 am and 10:00 am by teams of two to four, for a total of one person-hour. We performed reptile surveys three times at each site, with at least one week between surveys. During searches, all refugia (i.e., rocks, logs) within the plot were carefully overturned and then replaced. In addition to reptile species, we documented time of day and weather, as well as cover type and cover length if applicable. Our reptile surveys yielded 145 records of two snake species and four lizard species. Two additional snake species were observed incidentally during other sampling methods.Repeated counts of reptiles allowed us to use N-mixture models47 to estimate abundance and density. However, reptile data were extremely sparse and contained few species, precluding the use of a full community model as fit for birds. Instead, we structured reptile data into three groups: snakes (Western terrestrial garter snake (Thamnophis elegans) and Yellow-bellied racer (Coluber constrictor)), Alligator lizards (Elgaria coerulea and E. multicarinata) and Western fence lizards (Sceloporus occidentalis). We then combined data from all groups in a single N-mixture model, allowing for group-specific abundances. Due to sparse data, we only estimated an effect of burn category on abundance for fence lizards. Abundance for other species groups was assumed to be constant across burn categories in the model. To obtain species and burn category specific estimates of abundance, we combined model estimates of abundance with raw counts of individuals. For the two snakes, we calculated species level abundance by burn category by multiplying overall model-estimated snake abundance with the proportion of individuals made up by each species in each burn category. Because some Elgaria sp. observations could not be identified to species level, we did not attempt to calculate species specific abundances (but we note that the majority of observations was of E. coerula). We calculated species/group density by dividing abundance by the size of the sampling unit, in this case, a 90 × 90 m square. We fit the N-mixture model in R using the package unmarked ver. 0.12.256.Species identificationVertebrates were identified in the field or from photographs. Invertebrate specimens were collected and fixed in the field for identification in the lab where they were split into morphospecies using published keys57,58 and consultation with experts at the UC Davis Bohart Museum. Morphospecies were not always identified to lower taxonomic levels but were distinguished from similar species by coarse morphological characters. Prior to density estimation some morphospecies were aggregated based on taxonomy, body size, sample methodology, co-occurrence patterns and the difficulty in distinguishing members of the group. The rationales behind aggregations are documented in the raw invertebrate sampling data.Body size estimationFood webs are commonly organized around body size59,60,61,62,63 and we include length and mass estimates for all life stages in our data set. We assume mean body size does not change substantially across burn severities and use a single estimate for all life stages regardless of treatment. When possible, we measured 10 individuals for each morphospecies and life stage (for invertebrates and small mammals only, as individuals from other taxonomic groups were not handled). Body mass for small mammals was measured directly with a Pesola scale. Body length estimates for invertebrates were derived from the longest measurement. Volumetric estimates were derived by applying the measured length, height and width to the three-dimensional shape that most closely approximates that of the organism (e.g. ellipsoid, rectangular prism, cylinder, hemisphere, cone). We estimated biomass by multiplying these volumetric estimates by a tissue density of 1.1 g/mL30.Body sizes for organismal life stages observed but not captured (e.g. adult birds and large mammals) were estimated from the literature. Life stages for many organisms were not directly observed (e.g. lepidoptera eggs) but are almost certainly present. For these stages we again used the literature to estimate the body sizes. For example, we used a data set of egg sizes from 6,700 insect species to inform estimates for the species in our study64. The species list indicates the estimation method and literature source used to estimate the body size of each individual life stage is included in the species list.LinksEcological networks can be broken into two types: Undirected and directed. Undirected networks include bipartite networks (e.g. plant-pollinator, host-parasite) and social-networks. In undirected networks, interactions (links) and their participants (nodes) are observed at the same time. Links are not inferred in undirected ecological networks (unless false-negatives due to sampling error are taken into account) because they are directly observed (e.g. tick removed from a lizard). Undirected links can be weighted (interaction strength) by counting observations.Directed networks include food webs. Most food web links are inferred because it is not feasible to evaluate all possible consumer-resource interactions in a system through direct observation. The number of possible consumer-resource interactions in a food web is equivalent to the number of nodes squared. There were 3,084 nodes in the Eldorado National Forest, resulting in 9,511,056 potential consumer-resource interactions. In the scope of an ecological study, it is rarely feasible to directly observe most feeding links for most species15, much less an entire web, though molecular methods are bringing this closer to possibility for some guilds (e.g. large mammalian herbivores)65. While often detailed for birds and mammals, published accounts of direct observations of diets are often general (e.g. “Species A eats organic matter”) or entirely lacking for other groups. Restricting link assignments only to those that are directly observed will not create an accurate or unbiased food web.To assign resources to consumers we supplemented our own field observations with published diet records, expert opinion and rule-based filters. In the absence of direct observation, rule-based filters are an important tool for sorting through the huge number of potential feeding interactions. Filters varied with the species and ontogenetic stages to which they were applied but consisted of two main types: Encounter and compatibility. Encounter filters determine the potential for consumers to interact with resources. Encounter filters are based on habitat co-occurrence, forging strategy and diel activity patterns. Applied to links that have passed through encounter filters, compatibility filters determine the outcome of a potential encounter. Compatibility filters are based on consumer diet information, consumer diet breadth, resource palatability, and consumer-resource body size ratios. For inclusion in the food web, a link must pass through both filters.Link assignmentAll species were considered potential resources. For feasibility and reproducibility, species were broken into groups of encounter filters. First, we separated plants, fungi and detritus from other metazoans.Fig. 2Sampling sites map and King Fire perimeter. We sampled 27 sites total in the Eldorado National Forest, California, three years after the King Fire, nine in each burn category: Unburned, moderate severity, and high severity. Features not indicated in legend are typical of topological maps.Full size imageFig. 3Representative map of sampling design. We employed 19 different methods to estimate the richness and biomass density of organisms in the Eldorado National Forest, California, three years after the King Fire. Methods were conducted entirely within or centered within the 200 × 200 m site perimeter. Methods were paired in space and time when useful (e.g. black lights and acoustic bat surveys), and separated when necessary to avoid interference (e.g. small mammal trapping grid and vegetation transects). All methods at a site were conducted over 4–5 consecutive days.Full size imageThe first encounter filter was a loose group consisting of primary producers, non-living resources (e.g. detritus, carcasses), and saprophytes. In this first filter, primary producers were broken into parts (e.g. root, seed, leaves, etc). With few exceptions (i.e. moss, lichens, algae) primary producers in this filter group were evaluated as families or species. Saprophytes were evaluated as spores or adults. These groups served as the encounter filter for primary consumers and fungivores.Next, we partitioned all other metazoans into smaller encounter filters based on phylogeny, behavior, activity, habitat and palatability. A species’ life stage can be a member of multiple resource groups. These groups served as the encounter filter for most predatory or omnivorous organism stages. These encounter filters included:

Flying invertebrates

Nocturnal flying invertebrates

Diurnal flying invertebrates

Non-flying invertebrates

Ground-dwelling invertebrates

Ground-dwelling invertebrates excluding spider eggs

Soft-bodied ground dwelling invertebrates

Invertebrates on plants

Invertebrates on plants excluding spider eggs

Soft-bodied invertebrates on plants

Invertebrates on trees

Reptiles

Birds

Small mammals

Large mammals

Encounter filter exposure was tailored to consumer type. For consistency, consumers were broken into phylogenetic and ontogenetic guilds (e.g. lepidoptera larvae). For example, web-building spiders were exposed to flying insects, but not small-mammals encounter filters. Compatibility filters are specific to species and stage. For example, while both are exposed to the flying-insects encounter filter, adult and juvenile web-building spiders will have different compatibility filters because they have different consumer-resource body size ratios. Phylogenetic and ontogenetic consumer guilds are detailed below, but the decisions for each of the 178,655 links assignments in the food web can be found in and reproduced with the R code accompanying this manuscriptSpidersAs common, generalist insectivores, spiders are important consumers in the network66. To assign them resource links, spiders were broken into three ontogenetic stages: adult, egg, juvenile. Spiders were then separated into 17 consumer guilds by Family. The 87 spider species that could not be identified at the family level were assumed to be web-building spiders. To assign feeding interactions, we applied one or more invertebrate resource group filters (see above) to the members (species and stages) of each spider guild. Next, each link passed through a consumer-resource body size filter before being included in the network. For example, web-building spiders were assigned to capture flying insects and consume insects larger than 20% of their own body-length (as mesh size sets a lower limit of prey size) and less than 200% of their own body length67. Because of the generality of their filter-feeding hunting strategy, web-building spiders accounted for nearly half (32,663) of interactions with spiders as consumers. Additionally, when it was reported in the literature, adult members of a guild were assigned to consume nectar.Herbivorous beetlesTo assign consumer links to non-predatory beetles, their larval and adult stages were separated into Family-level feeding guilds. When information was available, links were evaluated at the species-stage level rather than the Family level. Resource links containing primary producers and non-living resources (detritus, carcasses, dung, etc) were assessed based on direct observation, published reports and expert opinions. When feeding on living plants, beetles were assigned to specific tissue types (e.g. flower, leaves, roots). For example, Cerambycidae_sp_2 was assigned to feed on the stems of Arctostaphylos patula as well as stems of plants from the family Rhamnaceae. This assigned 10,289 consumer links to herbivorous beetles.Larval lepidopteraCaterpillars were the most speciose herbivore group in the network. Larval lepidoptera vary greatly in host plant range and the quantity of information on their feeding habitats and were not lumped into guilds whenever possible. Microlepidoptera were not identified below the order level and were assumed to feed like Gelechiidae caterpillars (many species of which have been documented in Eldorado National Forest). We assigned links to caterpillars with resource filters of host plants and parts based on published diet records, palatability, expert opinion and observed co-occurrence between potential hosts and adults. This assigned 2,434 consumer links to caterpillars.Adult lepidopteraAdult lepidoptera feed primarily or exclusively on nectar from flowering plants. The host ranges for most adult lepidoptera are not well described. To assign nectar links, adult lepidoptera were grouped into families and passed through resource filters based on published records, co-occurrence, and expert opinion. Non-feeding adult moths were removed from link consideration. These filters assigned 5,854 nectar links to adult lepidoptera.Adult hymenoptera as herbivoresMost species of hymenoptera are parasitoids as larvae and free-living as adults (eusocial hymenoptera are a notable exception). Many adult parasitic wasps supplement their protein intake with nectar and pollen. The plant host ranges for solitary adult wasps are not well known. We applied nectar and pollen resource filters to adult hymenoptera based on published records, co-occurrence and expert opinion. These filters assigned 4,880 links to adult hymenoptera.Predatory beetlesPredatory and scavenging beetles are important consumers in terrestrial ecosystems. At the family level, the diets of predatory beetles are well-described relative to other arthropods. We grouped predatory and scavenging beetles into family guilds. We applied resource filters to predatory and scavenging beetles based on published records and expert opinion, habitat overlap, foraging strategy, palatability, and body size. These filters assigned 10,320 resource links to predatory beetles.FliesFlies in the Eldorado National Forest are a speciose group with diverse consumer strategies: herbivores, predators, scavengers, detritivores, parasitoids, and micropredators. Because of their trophic diversity each dipteran family was treated as its own guild, even then there were often large ontogenetic diet shifts across stages within a guild. We applied resource filters to dipterans based on published records and expert opinion of their diet. These filters were based on habitat overlap, palatability, and for predatory stages consumer-resource body size ratios. These filters assigned 9,541 resource links to flies.Hymenopteran parasitoidsParasitoid wasps are a diverse group that can have strong, even regulatory, effects on their host populations. We applied resource filters to hymenopteran parasitoid larvae in a three-step process. First, we used published records and expert opinion to establish their host range and host specificity. Next, we looked for the subset of potential hosts that co-occured with adult parasitoids across the most sites. Finally, we retained the subset of those host species whose adult body sizes were equivalent to or slightly larger than adult parasitoids. These filters assigned 4,695 host links to parasitoid wasps.Other hymenopteraThe remaining hymenoptera were composed of wood wasps and gall wasps, as well as bees, predatory wasps and ants. Host plant links were assigned to wood wasps and gall wasps in the same way as larval lepidoptera. Nectar links were assigned to bees in the same way as adult lepidoptera. Prey links were assigned to solitary and eusocial wasps in the same way as insectivorous gleaning birds. Ants were treated as generalist omnivores whose links were assigned links in a manner similar to all of the above.HemipteraHemiptera were the most abundant consumer group in the Eldorado National Forest. Hemiptera use their proboscis to feed on fluids, either as herbivores or predators. To accommodate this variation in consumer strategies resource filters were assigned to individual hemiptera species based on phylogeny. Host plants were assigned to herbivores based on field observations, published records and expert opinion. Plant fluids (with the exception of sap from soft woods) were not included in the node list, so herbivores were assigned feeding links based on the plant tissues that they pierced (i.e. stem or leaves). Predatory hemiptera were assigned insects on plants passed through consumer-resource body size filters. When more specific dietary information was available resource filters were further refined. These filters assigned 3,177 links to the hemipteraCollembolaFrom soil to canopy, springtails were widely distributed across Eldorado National Forest habitats and are important resources for arthropod secondary consumers. Collembola are herbivores, detritivores and fungivores. Little diet information is available for collembola, which were assigned resource filters based on habitat and phylogeny (order). Habitat was determined by collection method, and resource availability was determined by habitat. These, phylogenetic-habitat-resource filters assigned 140 links to the collembola.Bark licePsocoptera were common on plants in the Eldorado National Forest, feeding on algae, lichen, moss, fungi and detritus. They exhibit little variation in diet and were assigned resource filters as a group at the Order level. This small set of resource filters assigned 120 links to the psocoptera.ThripsThysanoptera were common and abundant consumers in the Eldorado National Forest, whose diets can range from herbivory to facultative predation and predation. To accommodate this variation in consumer strategy, resource filters were assigned to individual thysanoptera species based on phylogeny. Host plants for herbivores were based on field observations, supplemented with published records and expert opinion. Thysanoptera predators were assigned soft-bodied insects on plants passed through consumer-resource body size ratio filters. These resource filters assigned 924 links to thysanoptera.BatsBats are important consumers of arthropods. Arthropod encounter filters were determined by bat foraging strategies reported in the literature. Juvenile bats were treated separately. Aerial hawking bats capture nocturnal flying insects. Gleaning bats capture insects on plants. Pallid bats (Antrozous pallidus) which fly close to the ground, capture arthropods on the ground and understory plants. These arthropod filters then passed through consumer-resource body size ratio filters, which assigned 4,621 links to adult bats.BirdsBirds were the most speciose vertebrate group in the Eldorado National Forest. This diversity is reflected in their diets. Bird resource filters were based on published records and expert opinion. Herbivorous birds that were reported to feed on a type of plant tissue were assumed to feed on all tissues of that type, unless the literature indicated resource specialization. For example, if a bird was reported to feed on a fruit, it was assumed to feed on all fruit in the system. Resource links for predatory and omnivorous birds were passed through consumer-resource body size ratio filters. These filters assigned 56,584 links to birds as consumers.Small mammalsMammals 2.5 kg were comparatively uncommon but important consumers in the Eldorado National Forest. Large mammal resource filters were based on published records and expert opinion. Herbivorous large mammals that were reported to feed on a type of plant tissue were assumed to feed on all tissues of that type, unless the literature indicated resource specialization. Resource links for omnivorous and predatory large mammals were passed through consumer-resource body size ratio filters. These filters assigned 840 links to large mammals as consumers.Juvenile mammalsJuvenile mammals were considered unweaned. Nursing mammals were included in the species list because they are more vulnerable than adults and are potential resources to different consumers as a result. Unweaned mammals “feed” on their mothers, who are considered their only resource link. Nursing filters assigned 38 links to juvenile mammals.ReptilesWhile not diverse or common, reptiles on the mountain slopes of Eldorado National Forest are potentially important consumers for many groups. Reptile resource filters were assigned based on published records and expert opinion. All resource links were passed through consumer-resource body size ratio filters. These filters assigned 1,487 links to reptiles.MitesMites are ubiquitous and important consumers in many habitats but are often overlooked because of their small size and taxonomic difficulty. Mites were treated as consumer groups based on taxonomy. Resources for predatory mites were based on habitat and passed through consumer-resource body size ratios filters. Resource filters assigned 2,971 links to mites.Miscellaneous arthropodsA few arthropods not discussed above were not speciose enough to be dealt with separately here. These remaining miscellaneous arthropods were treated as groups at various taxonomic levels. Encounter and compatability filters for these groups were based on published records. Encounter filters for predatory (centipedes, odonates, pseudoscorpions, solfugids, mantids, opiliones, neuroptera, embioptera, pauropoda, eusocial wasps) and omnivorous (dermaptera) groups were based on habitat and passed through consumer-resource body size filters. Encounter filters were assigned by habitat for detritivores, herbivores and fungivores (archeognatha and isopods). These filters assigned 3,157 links to miscellaneous arthropods as consumers. More