Ecology

Darrah, S. E. et al. Improvements to the Wetland Extent Trends (WET) index as a tool for monitoring natural and human-made wetlands. Ecol. Ind. 99, 294–298 (2019).
Google Scholar 
Gardner, R.C., & Finlayson, C., (eds) Global Wetland Outlook: State of the World’s Wetlands and Their Services to People. (Ramsar Convention Secretariat, 2018).Tao, S. et al. Rapid loss of lakes on the Mongolian Plateau. Proc. Natl. Acad. Sci. 112(7), 2281–2286 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Davidson, N. Wetland losses and the status of wetland-dependent species. in The wetland book II: distribution, description, and conservation 369–381 (Springer, 2018).Xu, Y. et al. Loss of functional connectivity in migration networks induces population decline in migratory birds. Ecol. Appl. 29(7), e10960 (2019).
Google Scholar 
Wang, X. et al. Effects of anthropogenic landscapes on population maintenance of waterbirds. Conserv. Biol. https://doi.org/10.1111/cobi.13808 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Ma, Z., Cai, Y., Li, B. & Chen, J. Managing wetland habitats for waterbirds: An international perspective. Wetlands 30(1), 15–27 (2010).CAS 

Google Scholar 
Rönkä, M. T., Saari, C. L. V., Lehikoinen, E. A., Suomela, J., Häkkilä, K. (eds) Environmental changes and population trends of breeding waterfowl in northern Baltic Sea. in Annales Zoologici Fennici JSTOR (2005).Pöysä, H. & Paasivaara, A. Shifts in fine-scale distribution and breeding success of boreal waterbirds along gradients in ice-out timing and habitat structure. Freshw. Biol. 66(11), 2038–2050 (2021).
Google Scholar 
Borgo, J. S. & Conover, M. R. Visual and olfactory concealment of duck nests: Influence on nest site selection and success. Hum. Wildl. Interact. 10(1), 110–121 (2016).
Google Scholar 
Batbayar, N., Takekawa, J. Y., Natsagdorj, T., Spragens, K. A. & Xiao, X. Site selection and nest survival of the bar-headed goose (Anser indicus) on the Mongolian Plateau. Waterbirds Int. J. Waterbird Biol. 37(4), 381–393 (2014).
Google Scholar 
Syroechkovsky, E. V., Litvin K. E., Gurtovaya E. N. Nesting ecology of Bewick’s Swans on Vaygach Island, Russia. Waterbirds 221–6 (2002).Meijer, T. & Drent, R. Re-examination of the capital and income dichotomy in breeding birds. Ibis 141(3), 399–414 (1999).
Google Scholar 
Rothenbach, C. A. & Kelly, J. P. The parental dilemma under variable predation pressure: Adaptive variation in nest attendance by great egrets. Condor 114(1), 90–99 (2012).
Google Scholar 
Zhang, L., An, B., Shu, M. & Yang, X. Nest-site selection, reproductive ecology and shifts within core-use areas of black-necked cranes at the northern limit of the Tibetan Plateau. PeerJ 5(1), e2939 (2017).PubMed 
PubMed Central 

Google Scholar 
Newton, I. The Migration Ecology of Birds (Elsevier, 2010).
Google Scholar 
Donnelly, J. P. et al. Migration efficiency sustains connectivity across agroecological networks supporting sandhill crane migration. Ecosphere 12(6), e03543 (2021).
Google Scholar 
Barzen, J. A. & Serie, J. R. Nutrient reserve dynamics of breeding canvasbacks. Auk 107(1), 75–85 (1990).
Google Scholar 
Burnham, J. W. et al. Novel foraging by wintering Siberian cranes leucogeranus leucogeranus at China’s Poyang lake indicates broader changes in the ecosystem and raises new challenges for a critically endangered species. Bird Conserv. Int. 27, 204–223 (2017).
Google Scholar 
Krapu, G. L., Iverson, G. C., Reinecke, K. J. & Boise, C. M. Fat deposition and usage by arctic-nesting sandhill cranes during spring. Auk 102(2), 362–368 (1985).
Google Scholar 
Tøttrup, A. P. et al. Avian migrants adjust migration in response to environmental conditions en route. Biol. Lett. 4(6), 685–688 (2008).PubMed 
PubMed Central 

Google Scholar 
Sorte, F. A. L. et al. The role of atmospheric conditions in the seasonal dynamics of North American migration flyways. J. Biogeogr. 41(9), 1685–1696 (2014).
Google Scholar 
Marra, P. P., Francis, C. M., Mulvihill, R. S. & Moore, F. R. The influence of climate on the timing and rate of spring bird migration. Oecologia 142(2), 307–315 (2005).ADS 
PubMed 

Google Scholar 
Studds, C. E. & Marra, P. P. Rainfall-induced changes in food availability modify the spring departure programme of a migratory bird. Proc. R. Soc. B Biol. Sci. 278(1723), 3437–3443 (2011).
Google Scholar 
Saino, N. et al. Ecological conditions during winter predict arrival date at the breeding quarters in a trans-Saharan migratory bird. Ecol. Lett. 7(1), 21–25 (2004).
Google Scholar 
Barzen, JA. Chapter 15 – Ecological Implications of Habitat Use by Reintroduced and Remnant Whooping Crane Populations. In: Nyhus PJ, French JB, Converse SJ, Austin JE, Delap Jh (eds.), Whooping Cranes: Biology and Conservation. (Academic Press, New York, 2019).Norris, D. R., Marra, P. P., Kyser, T. K., Sherry, T. W. & Ratcliffe, L. M. Tropical winter habitat limits reproductive success on the temperate breeding grounds in a migratory bird. Proc. R. Soc. B Biol. Sci. 271(1534), 59–64 (2004).
Google Scholar 
Norris, D. R. & Taylor, C. M. Predicting the consequences of carry-over effects for migratory populations. Biol. Let. 2(1), 148–151 (2006).
Google Scholar 
Newton, I. Population Limitation in Birds (Academic press, 1998).
Google Scholar 
Altwegg, R. & Anderson, M. D. Rainfall in arid zones: Possible effects of climate change on the population ecology of blue cranes. Funct. Ecol. 23(5), 1014–1021 (2009).
Google Scholar 
Layton-Matthews, K., Hansen, B. B., Grøtan, V., Fuglei, E. & Loonen, M. J. J. E. Contrasting consequences of climate change for migratory geese: Predation, density dependence and carryover effects offset benefits of high-arctic warming. Glob. Chang. Biol. 26(2), 642–657 (2020).ADS 
PubMed 

Google Scholar 
Madsen, J. & Fox, A. D. Impacts of hunting disturbance on waterbirds-a review. Wildl. Biol. 1(4), 193–207 (1995).
Google Scholar 
Minias, P. Reproduction and survival in the city: Which fitness components drive urban colonization in a reed-nesting waterbird?. Curr. Zool. 62(2), 79–87 (2016).PubMed 
PubMed Central 

Google Scholar 
Cheng, Y., Fiedler, W., Wikelski, M. & Flack, A. “Closer-to-home” strategy benefits juvenile survival in a long-distance migratory bird. Ecol. Evol. 9, 8945–8952 (2019).PubMed 
PubMed Central 

Google Scholar 
Sargeant, A. B. & Raveling, D. G. Mortality during the breeding season. In Ecology and Management of Breeding Waterfowl (eds Batt, B. D. J. et al.) 396–422 (University of Minnesota Press, 1992).
Google Scholar 
IUCN. The IUCN Red List of Threatened Species. http://www.iucnredlist.org Downloaded on 22 November 2021. Version 2021–2 (2021).Mudrik, E. A. et al. Gene pool homogeneity of western and eastern populations of the white-naped crane antigone vipio in different flyways. Russ. J. Genet. 58(5), 566–575 (2022).CAS 

Google Scholar 
Wang, W., Fraser, J. D. & Chen, J. Distribution and long-term population trends of wintering waterbirds in Poyang lake China. Wetlands 39(S1), 125–135 (2019).
Google Scholar 
Wang, W., Fraser, J. D. & Chen, J. Wintering waterbirds in the middle and lower Yangtze river floodplain: Changes in abundance and distribution. Bird Conserv. Int. 27(2), 167–186 (2017).CAS 

Google Scholar 
Gilbert, M., Buuveibaatar, B., Fine, A. E., Jambal, L. & Strindberg, S. Declining breeding populations of white-naped cranes in eastern Mongolia, a ten-year update. Bird Conserv. Int. 26(04), 490–504 (2016).
Google Scholar 
Mirande, C.M., Harris J.T. Crane Conservation Strategy. in International Crane Foundation (Baraboo, Wisconsin, USA, 2019).Bouchard, F. et al. Paleolimnology of thermokarst lakes: A window into permafrost landscape evolution. Arct. Sci. 3, 91–117 (2016).
Google Scholar 
Fernández-Tizón, M., Emmenegger, T., Perner, J. & Hahn, S. Arthropod biomass increase in spring correlates with NDVI in grassland habitat. Sci. Nat. 107(5), 42 (2020).
Google Scholar 
Dong, J. et al. Mapping paddy rice planting area in northeastern Asia with landsat 8 images, phenology-based algorithm and google earth engine. Remote Sens. Environ. 185, 142–154 (2016).ADS 
PubMed 
PubMed Central 

Google Scholar 
Pekel, J.-F., Cottam, A., Gorelick, N. & Belward, A. S. High-resolution mapping of global surface water and its long-term changes. Nature 540(7633), 418–422 (2016).ADS 
CAS 
PubMed 

Google Scholar 
Hersbach, H. et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 146(730), 1999–2049 (2020).ADS 

Google Scholar 
Bivand, R., Rundel, C. Rgeos: Interface to geometry engine-open source (GEOS). R package version 03–26. (2017).Zuur, A. F., Ieno, E. N. & Elphick, C. S. A protocol for data exploration to avoid common statistical problems. Methods Ecol. Evol. 1(1), 3–14 (2010).
Google Scholar 
Dormann, C. F. et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 36(1), 27–46 (2013).
Google Scholar 
Barton, K., Barton, M. K. Package ‘MuMIn’. R package version. 1(6), (2019).Anderson, D. & Burnham, K. Model selection and multi-model inference. Second NY Springer-Verlag 2004(63), 10 (2020).
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(2), 880–890 (2019).PubMed 
PubMed Central 

Google Scholar 
Batbayar, N. et al. Combining tracking and remote sensing to identify critical year-round site, habitat use and migratory connectivity of a threatened waterbird species. Remote Sens. 13(20), 4049 (2021).ADS 

Google Scholar 
Lang, X. et al. Luan river upper reaches: The important stopover site of the white-naped crane (Grus vipio) western population. Biodivers. Sci. 28(10), 1213 (2020).
Google Scholar 
Jia, Y. et al. Shifting of the migration route of white-naped crane (Antigone vipio) due to wetland loss in China. Remote Sens. 13(15), 2984 (2021).ADS 
CAS 

Google Scholar 
Lloyd, C. T. et al. Global spatio-temporally harmonised datasets for producing high-resolution gridded population distribution datasets. Big Earth Data 3(2), 108–139 (2019).PubMed 
PubMed Central 

Google Scholar 
Kéry, M. & Schaub, M. Bayesian Population Analysis Using WinBUGS: A Hierarchical Perspective (Academic Press, 2011).
Google Scholar 
Lê, S., Josse, J. & Husson, F. FactoMineR: An R package for multivariate analysis. J. Stat. Softw. 25(1), 1–18 (2008).
Google Scholar 
Hastie, T. J. & Tibshirani, R. J. Generalized Additive Model (Routledge, 2017).MATH 

Google Scholar 
Fair, J. M. & Jones, J. Guidelines to the Use of Wild Birds in Research (Ornithological council, 2010).
Google Scholar 
RStudio Team RStudio. Integrated Development Environment for R (Studio, PBC, 2022).
Google Scholar 
Higuchi, H. et al. Using a remote technology in conservation: Satellite tracking white-naped cranes in Russia and Asia. Conserv. Biol. 18(1), 136–147 (2004).
Google Scholar 
Weller, M. W. Wetland Birds: Habitat Resources and Conservation Implications (University Press, Cambridge, 1999).
Google Scholar 
Vandandorj, S., Gantsetseg, B. & Boldgiv, B. Spatial and temporal variability in vegetation cover of Mongolia and its implications. J. Arid Land 7(4), 450–461 (2015).
Google Scholar 
Bradter, U. T. E., Gombobaatar, S., Uuganbayar, C., Grazia, T. E. & Exo, K.-M. Reproductive performance and nest-site selection of white-naped cranes grus vipio in the Ulz river valley, north-eastern Mongolia. Bird Conserv. Int. 15(4), 313–326 (2005).
Google Scholar 
Bradter, U., Gombobaatar, S., Uuganbayar, C., Grazia, T. E. & Exo, K. M. Time budgets and habitat use of white-naped cranes grus vipio in the Ulz river valley, northeastern Mongolia during the breeding season. Bird Conserv. Int. 17(3), 259–271 (2007).
Google Scholar 
Zheng, Y. et al. Dynamic changes and driving factors of wetlands in inner Mongolia plateau China. PLoS ONE 14(8), e0221177 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Linscott, J. A. & Senner, N. R. Beyond refueling: Investigating the diversity of functions of migratory stopover events. Ornithol. Appl. 123, duaa074 (2021).
Google Scholar 
Xu, Y., Kieboom, M., van Lammeren, R. J. A., Si, Y. & de Boer, W. F. Indicators of site loss from a migration network: Anthropogenic factors influence waterfowl movement patterns at stopover sites. Global Ecol. Conserv. 25, e01435 (2021).
Google Scholar 
Piersma, T. et al. Simultaneous declines in summer survival of three shorebird species signals a flyway at risk. J. Appl. Ecol. 53(2), 479–490 (2016).
Google Scholar 
Wilson, S., Gil-Weir, K. C., Clark, R. G., Robertson, G. J. & Bidwell, M. T. Integrated population modeling to assess demographic variation and contributions to population growth for endangered whooping cranes. Biol. Conserv. 197, 1–7 (2016).
Google Scholar 
Wheeler, M. E., Barzen, J. A., Crimmins, S. M. & Deelen, T. R. V. Effects of territorial status and life history on sandhill crane (Antigone canadensis) population dynamics in south-central Wisconsin, USA. Can. J. Zool. 97(2), 112–120 (2019).
Google Scholar 
Gerber, B. D. & Kendall, W. L. Considering transient population dynamics in the conservation of slow life-history species: An application to the sandhill crane. Biol. Conserv. 200, 228–239 (2016).
Google Scholar 
Servanty, S., Converse, S. J. & Bailey, L. L. Demography of a reintroduced population: Moving toward management models for an endangered species, the whooping crane. Ecol. Appl. 24(5), 927–937 (2014).PubMed 

Google Scholar 
Xu, W. et al. Hidden loss of wetlands in China. Curr. Biol. 29(18), 3065–71.e2 (2019).CAS 
PubMed 

Google Scholar 
Niu, Z. et al. Mapping wetland changes in China between 1978 and 2008. Chin. Sci. Bull. 57(22), 2813–2823 (2012).
Google Scholar 
Austin, J., Morrison, K. & Harris, J. Cranes and Agriculture: A Global Guide for Sharing the Landscape (International Crane Foundation, Baraboo, 2018).
Google Scholar 
Lee, S. D., Jabłoński, P. G. & Higuchi, H. Winter foraging of threatened cranes in the demilitarized zone of Korea: Behavioral evidence for the conservation importance of unplowed rice fields. Biol. Conserv. 138(1), 286–289 (2007).
Google Scholar 
Nilsson, L., Bunnefeld, N., Persson, J., Žydelis, R. & Månsson, J. Conservation success or increased crop damage risk? The Natura 2000 network for a thriving migratory and protected bird. Biol. Conserv. 236, 1–7 (2019).
Google Scholar 
Mukherjee, A., Borad, C. K. & Parasharya, B. M. Breeding performance of the Indian sarus crane in the agricultural landscape of western India. Biol. Conserv. 105(2), 263–269 (2002).
Google Scholar 
Gopi Sundar, K. S. Are rice paddies suboptimal breeding habitat for sarus cranes in Uttar Pradesh, India?. Condor 111(4), 611–623 (2009).
Google Scholar 
Pekarsky, S., Schiffner, I., Markin, Y. & Nathan, R. Using movement ecology to evaluate the effectiveness of multiple human-wildlife conflict management practices. Biol. Conserv. 262, 109306 (2021).
Google Scholar 
Hemminger, K., König, H., Månsson, J., Bellingrath-Kimura, S.-D. & Nilsson, L. Winners and losers of land use change: A systematic review of interactions between the world’s crane species (Gruidae) and the agricultural sector. Ecol. Evol. 12(3), e8719 (2022).PubMed 
PubMed Central 

Google Scholar 
Amano, T. Conserving bird species in Japanese farmland: Past achievements and future challenges. Biol. Conserv. 142(9), 1913–1921 (2009).
Google Scholar 
Okuya, K. et al. Isolation and characterization of influenza A viruses from environmental water at an overwintering site of migratory birds in Japan. Adv. Virol. 160(12), 3037–3052 (2015).CAS 

Google Scholar 
Wille, M. & Barr, I. G. Resurgence of avian influenza virus. Science 376(6592), 459–460 (2022).ADS 
CAS 
PubMed 

Google Scholar  More

All animal care and experimental protocols complied with local and federal laws and guidelines and were approved by the appropriate governing body in Berlin, Germany, the Landesamt fur Gesundheit und Soziales (LaGeSo G-0224/20).Experimental breeding and designThe all-female Amazon molly (Poecilia formosa) is a naturally clonal, live-bearing fish species that gives birth to broods of genetically identical offspring. Like all unisexual vertebrates, Amazon mollies are the result of inter-specific hybridization44,45. As such, this ‘frozen hybrid’ has a heterozygous genome from its ancestral P. mexicana mother and P. latipinna father alleviating concerns about reduced genetic variation and the resulting inbreeding depression often associated with artificially selected isogenic animals. Additionally, despite their clonal nature, the Amazon’s genome shows no evidence of increased mutation accumulation, genomic decay or transposable element activity suggesting the genomes of these animals are evolving in similar ways as sexual species46. They reproduce through gynogenesis where the meiotic process is disrupted so that the eggs contain a full maternal genome. The egg must be fused with a sperm from one of their ancestral species to stimulate embryogenesis, but this paternal DNA is not incorporated into the egg. This provides the opportunity to control when reproduction occurs by controlling the females’ access to male sperm donors.We placed adult females, as potential mothers of experimental fish, in individual (5-gallon) breeding tanks with two Atlantic molly (P. mexicana) males for one week to act as sperm donors. Amazon mollies give birth to broods of generally ~8-30 individuals. A brood is born at once (i.e. all individuals are born within minutes of each other) and birth generally happens early in the day close to dawn. These parental fish were lab-bred and themselves sisters, so of the same age and lineage, and were kept at similar social densities and under standardized environmental conditions throughout their lives to further minimize potential variation in maternal experience. Each breeding tank contained an artificial plant as refuge and was checked frequently each day for the presence of offspring, especially during the morning hours when births are most likely. Newborn mollies were always found in the morning and then singly netted by trained animal caretakers, into individual experimental tanks where their behavior was automatically recorded for the next 70 days (see below). Moving the fish from the maternal tank to the experimental tanks was done in a standardized manner (i.e. individual fish were netted and placed into small dishes of water and then placed in the tracking tanks to limit exposure to the air) by the same caretakers to minimize variation in experience among individual fish. Altogether, eight mothers provided offspring that completed the entire 10-week experiment (Supplementary Table 1).Experimental tanks (27 x 27 cm), made of white Perspex, consisted of four equally sized compartments, and were evenly lit from below using 6500K-LEDs. Environmental conditions were highly standardized across tanks: all tanks were on the same 11:13 (L:D) light schedule, water depth was maintained at 10 cm depth, temperature was maintained at 25 ± 1 °C by a room air conditioning system, and fish received a standardized amount of powdered flake fish food (TetraMin™) twice daily. Opaque blinds surrounded the tanks to further limit outside disturbances. All experimental tanks were connected to the same filtration system where water could mix in the sump tank, allowing chemical cues to be shared across all experimental fish. Previous work has shown exposure to just chemical cues of conspecifics is sufficient in preventing the developmental of pathological behavior that could be associated with development in complete isolation14. We initially placed a total of 40 newborn individuals into the tracking tanks. At the end of the 10-week experiment, we were able to achieve complete tracking data on 26 individuals; camera malfunctions prevented data collection on four individuals, two individuals jumped into neighboring tanks causing the loss of data of all four individuals as we could not verify their identity; four newborn individuals escaped through holes in the water outlet of the tanks; and four individuals died as newborns. All results in the manuscript are on these 26 animals, though including data from all 40 (e.g. patterns of individual variation on the first day post birth) did not change the results or their interpretation (see Supplementary Table 2).Behavioral trackingWe developed a custom recording system using Raspberry Pi computers, which are an upcoming low-cost, highly adaptable solution for many applications in the biological sciences25. Specifically, we created a local network of Raspberry Pi 3B + ’s, each connected to a Raspberry Pi camera positioned exactly above an experimental tank, commanded by a lab computer, and connected to the server on the institute network (Supplementary Fig. 1). We programmed the Raspberry Pi’s using pirecorder26 to take timestamped photos every 3 s across the daily light period, each day, for 10 weeks, and store them automatically in dedicated, automatically named folders on the server. Image settings and resolution were thereby optimized to minimize file size while assuring image quality. After the experimental period, we created videos of all the recorded images of each fish of each day. These videos were subsequently tracked with the Biotracker software27, using background subtraction, providing the x, y coordinates of each fish in each frame. We then processed the data, including scaling and converting the coordinates to mm, and, for each frame, computed fish’s swimming speed (cm/s) and distance from the tank walls (cm). We then summarized these variables both on an hourly and daily basis to compute fish’s median swimming speed, inter-quartile range of swimming speeds, activity (proportion of time spent moving >0.5 cm/s), and median border distance. To quantify fish’s body size over time, we randomly selected five photos per week of each compartment, making sure the fish was away from the compartment walls and did not show strong body curvature, and then used ImageJ software to measure total body length (mm) from the tip of the snout to the end of the body. By averaging the measurements of the five images, we acquired one body size measurement per week.Error checkingWe collected up to 924,000 photos on each individual throughout the experimental period resulting in a total of over 24 million data points collected on our experimental animals (N = 26 individuals). To ensure that our tracking software accurately captured the behavior of our fish, we checked for potential tracking errors in two ways. First, we estimated overall error rates. To do this, we selected at random a starting frame from within a day; then we manually checked each of the subsequent 200 frames and identified whether an error was made (fish was not properly located by BioTracker) or not (fish was properly located) by visual inspection of the videos. We estimated the error rate as the number of errors divided by the total number of checked frames. The overall median error rate over the entire observation period was estimated to be 7%. Error rates increased earlier in the observation period when the fish were smaller (Supplementary Note I). As such, as a second step, we manually went through and corrected all frames for the very first day of tracking (i.e. day 1 post-birth) for all fish (~13,200 frames per individual) as this is a critical time period for one of our research questions. This ensured that the resulting behavioral data were completely accurate for this day. This manual correction allowed us the additional opportunity to compare how well our automatically tracked (i.e. not manually corrected) data performed compared to the manually corrected data. We found that the automatically tracked data re-created near identical estimates of among- and within-individual variance components and most importantly the among-individual correlation between the automatically tracked and manually corrected data was over 0.98 for our behavioral variables (Supplementary Note I). This strongly suggests that any errors introduced by our automated tracking software have minimal influence of our behavioral variables at best and do not affect our interpretation of the results.Statistical analysesWe used linear mixed, or hierarchical, models to partition the behavioral variation across different times periods into its among- and within-individual components. Throughout we focused our analysis on the 26 individuals for which we had complete data for the entire 10-week observation period to ensure comparable variation over time and across models.Our first question of interest was to test when individual differences in behavior first appeared over the course of the experiment. We started by investigating behavior on the first day post birth (Fig. 1A, Supplementary Table 2) and then planned to proceed in a day-by-day fashion until significant repeatability in behavior was apparent (Supplementary Table 3). We used hourly median swimming speed (11 observations for each of 26 individuals) as our response variable and included ‘hour’ and ‘total length (TL)’ as fixed effects and ‘individual’ was included as our random effect of interest. Including TL as a covariate allowed us to test whether behavior was related to an offspring’s body size on its first day of life. We set the first hour of the day as 0 and mean-centered TL as this would allow the among- (and within-) individual variance components to be estimated at these values (i.e. the earliest possible moment from when we could record behavior in the fish). We estimated the adjusted repeatability of median swimming speed as the variance attributable to individual identity over the total variance not explained by the fixed effects. We additionally estimated both marginal and conditional R-squared values which estimate the variance explained by the fixed effects only and the variance explained by the fixed and random effects combined, respectively. As our individual experimental fish came from different mothers, we first explored a number of different variance structures including random intercepts and slopes for both individual ID and maternal ID. This allowed us to test whether maternal identity explained variation in individual behavior. However, the most supported model included random intercepts and slopes for individual ID and not for mother ID, indicating that our methods to reduce variation among mothers were successful (Table 1). We used median swimming speed as our behavioral variable of interest throughout the main manuscript, as this behavior was tightly correlated with most of our other behavioral variables (Supplementary Fig. 2); though results using the other behavioral variables yielded the same interpretation (i.e. that significant individuality in (any) behavior was present on the very first day post-birth; Supplementary Table 2).Our second research question was to investigate how individual behavioral variance changed over the course of the entire observation period (70 days). Again, we first explored several different variance structures to test the importance of maternal identity and/or individual identity on behavioral variation. We found support for the inclusion of random slopes at the individual level, but not maternal level (Table 1). This indicates that levels of among- (and within-) individual variation may differ throughout the observation period. To investigate patterns of change in the variance components, we ran a series of models where we centered the observation covariate on different days. Individual intercepts are estimated when all covariates are set to zero, so this allowed us to ‘slice’ the data to estimate the among- and within-individual variance at different time points over the ten weeks. We ran 11 models as we chose to center the data every 7 days (first model was centered on observation 1; 11th model was centered on observation 70). The predicted individual intercepts (best linear unbiased predictors) and estimated variance components from each model are plotted in Fig. 3.We also closely investigated any potential influence of body size and/or growth rate differences on behavioral expression and individual behavioral variation in this entire 10-week data set. First, we estimated the repeatability of both weekly total length and weekly growth rates to determine if individuals consistently differed in these traits. Then, we ran a series of models with median weekly swimming speed as the response variable and included either weekly total length, weekly growth rate, and/or overall growth rate (estimated over the entire 10 weeks), as our fixed effects of interest. Each model also included the random effects of individual intercepts and slopes. Finally, because body size varies both among individuals (some individuals are on average larger than others) and within individuals (as they grow), we also performed within-individual centering of total length. In this fifth model, we included each individual’s average total length and their weekly deviation from their average length as the two fixed effects of interest. Individual identity and slopes were included as random effects. For all models, we estimated the variance explained by the fixed effects (marginal R2) and the fixed and random effects together (conditional R2). These results are reported in Table 2.For our third and final research question, we tested whether early-life behavior predicted later-life behavior. To test this, we estimated the among-individual correlation (including ‘individual ID’ as our random effect) in behavior using multivariate mixed models where the daily median swimming speeds in each week were the response variables (7 observations per week per individual; 10 weeks total; Fig. 4A). Then to investigate how the strength of these correlations may change over development, we used a linear model to test whether the correlation strength was predicted by the interaction between the first week included in the correlation and distance to the next week in the correlation (1, 2, 3, 4 or 5 weeks away in time; Fig. 4B).All models were performed using Markov Chain Monte Carlo estimation with the MCMCglmm package38 in R v3.6.139. We set our models to run 510,000 iterations with a 10,000 burn-in and thinning every 200 iterations. To ensure proper model mixing and convergence, we initially ran 5 independent chains and inspected posterior trace plots of parameter estimates (Supplementary Note II). In a preliminary analysis we tested three different prior settings (Supplementary Note II); results did not change with prior settings so we chose parameter-expanded priors for all models reported here as these are generally considered to be more robust. An R Markdown file with all the results presented here is included in Supplementary Note II.Reporting summaryFurther information on research design is available in the Nature Research Reporting Summary linked to this article. More

Estes, J. et al. Megafaunal impacts on structure and function of ocean ecosystems. Annu. Rev. Environ. Resour. 41, 83–116 (2016).
Google Scholar 
Ferretti, F., Worm, B., Britten, G., Heithaus, M. & Lotze, H. Patterns and ecosystem consequences of shark declines in the ocean. Ecol. Lett. 13, 1055–1071 (2010).PubMed 

Google Scholar 
Heithaus, M. R. et al. Seagrasses in the age of sea turtle conservation and shark overfishing. Front. Mar. Sci. 1, 1–6 (2014).
Google Scholar 
McCauley, D. et al. Marine defaunation: Animal loss in the global ocean. Science 347, 1255641 (2015).PubMed 

Google Scholar 
Pereira, H. et al. Scenarios for global biodiversity in the 21st century. Science 330, 1496–1501 (2010).ADS 
CAS 
PubMed 

Google Scholar 
Oliver, S., Braccini, M., Newman, S. & Harvey, E. S. Global patterns in the bycatch of sharks and rays. Mar. Policy 54, 86–97 (2015).
Google Scholar 
Pacoureau, N. et al. Half a century of global decline in oceanic sharks and rays. Nature 589, 567–574 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Gilman, E. et al. Shark interactions in pelagic longline fisheries. Mar. Policy 32, 1–18 (2008).
Google Scholar 
Worm, B. et al. Global catches, exploitation rates, and rebuilding options for sharks. Mar. Policy 40, 194–204 (2013).
Google Scholar 
Bowlby, H. & Gibson, A. Implications of life history uncertainty when evaluating status in the Northwest Atlantic population of white shark (Carcharodon carcharias). Ecol. Evol. 10, 4990–5000 (2020).PubMed 
PubMed Central 

Google Scholar 
Dulvy, N. et al. Overfishing drives over one third of all sharks and rays toward a global extinction crisis. Curr. Biol. 31, 4773-4787.e8 (2021).CAS 
PubMed 

Google Scholar 
Heino, M., Pauli, B. & Dieckmann, U. Fisheries-induced evolution. Annu. Rev. Ecol. Evol. Syst. 46, 461–480 (2015).
Google Scholar 
Mitchell, J., McLean, D., Collins, S. & Langlois, T. Shark depredation in commercial and recreational fisheries. Rev. Fish Biol. Fish 28, 715–748 (2018).
Google Scholar 
Jaiteh, V. F., Loneragan, N. & Warren, C. The end of shark finning? Impacts of declining catches and fin demand on coastal community livelihoods. Mar. Policy 82, 224–233 (2017).
Google Scholar 
Seidu, I. et al. Fishing for survival: Importance of shark fisheries for the livelihoods of coastal communities in Western Ghana. Fish. Res. 246, 106157 (2022).
Google Scholar 
Gilman, E., Weijerman, M. & Suuronen, P. Ecological data from observer programs underpin ecosystem-based fisheries management. ICES J. Mar. Sci. 74, 1481–1495 (2017).
Google Scholar 
Melnychuk, M. et al. Identifying management actions that promote sustainable fisheries. Nat. Sustain. https://doi.org/10.1038/s41893-020-00668-1 (2021).Article 

Google Scholar 
Musyl, M. & Gilman, E. Meta-analysis of post-release fishing mortality in apex predatory pelagic sharks and white marlin. Fish Fish. 20, 466–500 (2019).
Google Scholar 
Clarke, S. A status snapshot of key shark species in the western and central pacific and potential management options. in WCPFC-SC7-2011/EB-WP-04. Western and Central Pacific Fisheries Commission, Kolonia, Federated States of Micronesia (2011).Dapp, D., Walker, T., Huveneers, C. & Reina, R. Respiratory mode and gear type are important determinants of elasmobranch immediate and post-release mortality. Fish Fish. 17, 507–524 (2016).
Google Scholar 
ICES. Report of the working group on elasmobranch fishes. in ICES CM 2018/ACOM:16. International Council for the Exploration of the Sea, Copenhagen (2018).Dicks, L. et al. A transparent process for “evidence-informed” policy making. Conserv. Lett. 7, 119–125 (2014).
Google Scholar 
Nichols, J., Kendall, W. & Boomer, G. Accumulating evidence in ecology: Once is not enough. Ecol. Evol. 9, 13991–14004 (2019).PubMed 
PubMed Central 

Google Scholar 
Nakagawa, S. et al. Meta-analysis of variation: Ecological and evolutionary applications and beyond. Methods Ecol. Evol. 6, 143–152 (2015).
Google Scholar 
Pfaller, J., Chaloupka, M., Bolten, A. & Bjorndal, K. Phylogeny, biogeography and methodology: A meta-analytic perspective on heterogeneity in adult marine turtle survival rates. Sci. Rep. 8, 5852. https://doi.org/10.1038/s41598-018-24262-w (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Godin, A., Carlson, J. & Burgener, V. The effect of circle hooks on shark catchability and at-vessel mortality rates in longlines fisheries. Bull. Mar. Sci. 88, 469–483 (2012).
Google Scholar 
Reinhardt, J. et al. Catch rate and at-vessel mortality of circle hooks versus J-hooks in pelagic longline fisheries: A global meta-analysis. Fish Fish. 19, 413–430 (2018).
Google Scholar 
Rosa, D., Santos, C. & Coelho, R. Assessing the effects of hook, bait and leader type as potential mitigation measures to reduce bycatch and mortality rates of shortfin mako: A meta-analysis with comparisons for target, bycatch and vulnerable fauna interactions. in ICCAT Collective Volume of Scientifics Papers 76, 247–278 (2020).Santos, C., Rosa, D. & Coelho, R. Hook, bait and leader type effects on surface pelagic longline retention and mortality rates: A meta-analysis with comparisons for target, bycatch and vulnerable fauna interactions. in IOTC-2019-WPEB15-39. Indian Ocean Tuna Commission, Mahe, Seychelles (2019).Santos, C., Rosa, D. & Coelho, R. Progress on a meta-analysis for comparing hook, bait and leader effects on target, bycatch and vulnerable fauna interactions. in Collective Volume of Scientifics Papers ICCAT 77, 182–217 (2020).Condamine, F., Romieu, J. & Guinot, G. Climate cooling and clade competition likely drove the decline of lamniform sharks. PNAS 116, 20584–20590 (2019).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Vehtari, A., Gelman, A., Simpson, D., Carpenter, B. & Bürkner, P. Rank-normalization, folding, and localization: An improved Rhat for assessing convergence of MCMC (with Discussion). Bayesian Anal. 16, 667–718 (2021).MathSciNet 

Google Scholar 
Cinar, O., Nakagawa, S. & Viechtbauer, W. Phylogenetic multilevel meta-analysis: A simulation study on the importance of modelling the phylogeny. Methods Ecol. Evol. 13, 383–395 (2022).
Google Scholar 
Lajeunesse, M. Meta-analysis and the comparative phylogenetic method. Am. Nat. 174, 369–381 (2009).PubMed 

Google Scholar 
Chamberlain, S. et al. Does phylogeny matter? Assessing the impact of phylogenetic information in ecological meta-analysis. Ecol. Lett. 15, 627–636 (2012).PubMed 

Google Scholar 
Burns, J. & Strauss, S. More closely related species are more ecologically similar in an experimental test. Proc. Natl. Acad. Sci. USA 108, 5302–5307 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Cachera, M. & Le Loc’h, F. Assessing the relationships between phylogenetic and functional singularities in sharks (Chondrichthyes). Ecol. Evol. 7, 6292–6303 (2017).PubMed 
PubMed Central 

Google Scholar 
Münkemüller, T., Boucher, F. C., Thuiller, W. & Lavergne, S. Phylogenetic niche conservatism—Common pitfalls and ways forward. Funct. Ecol. 29, 627–639 (2015).PubMed 
PubMed Central 

Google Scholar 
Bazzi, M., Campione, N., Kear, B., Pimiento, C. & Ahlberg, P. Feeding ecology has shaped the evolution of modern sharks. Curr. Biol. 31, 5138–5148 (2021).CAS 
PubMed 

Google Scholar 
Sepulveda, C., Wegner, N., Bernal, D. & Graham, J. The red muscle morphology of the thresher sharks (family Alopiidae). J. Exp. Biol. 208, 4255–4261 (2005).CAS 
PubMed 

Google Scholar 
Wosnick, N. et al. Multispecies thermal dynamics of air-exposed ectothermic sharks and its implications for fisheries conservation. J. Exp. Mar. Biol. Ecol. 513, 1–9 (2019).
Google Scholar 
French, R. et al. High survivorship after catch-and-release fishing suggests physiological resilience in the endothermic shortfin mako shark (Isurus oxyrinchus). Conserv. Physiol. https://doi.org/10.1093/conphys/cov044 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Davis, M. Key principles for understanding fish bycatch discard mortality. Can. J. Fish. Aquat. Sci. 59, 1834–1843 (2002).
Google Scholar 
Massey, Y., Sabarros, P., Rabearisoa, N. & Bach, P. Drivers of at-haulback mortality of sharks caught during pelagic longline fishing experiments. in IOTC-2019-WPEB15-14_Rev1. Indian Ocean Tuna Commission, Mahe, Seychelles (2019).Musyl, M., Moyes, C., Brill, R. & Fragoso, N. Factors influencing mortality estimates in post-release survival studies: Comment on Campana et al. (2009). Mar. Ecol. Prog. Ser. 396, 157–159 (2009).Pimiento, C., Cantalapiedra, J., Shimada, K., Field, D. & Smaers, J. Evolutionary pathways towards gigantism in sharks and rays. Evolution 73, 588–599 (2019).PubMed 

Google Scholar 
Musyl, M. & Gilman, E. Post-release fishing mortality of blue (Prionace glauca) and silky shark (Carcharhinus falciformes) from a Palauan-based commercial longline fishery. Rev. Fish Biol. Fish. 28, 567–658 (2018).
Google Scholar 
Childs, D., Sheldon, B. & Rees, M. The evolution of labile traits in sex- and age-structured populations. J. Anim. Ecol. 85, 329–342 (2016).PubMed 
PubMed Central 

Google Scholar 
Comte, L., Murienne, J. & Grenouillet, G. Species traits and phylogenetic conservatism of climate-induced range shifts in stream fishes. Nat. Commun. 5, 5053 (2014).ADS 
PubMed Central 

Google Scholar 
IUCN. The IUCN Red List of Threatened Species. Version 2021-3. www.iucnredlist.org. ISSN 2307-8235 (International Union for the Conservation of Nature, Gland, Switzerland, 2022).García, V., Lucifora, L. & Ransom, M. The importance of habitat and life history to extinction risk in sharks, skates, rays and chimaeras. Proc. R. Soc. B 275, 83–89 (2008).PubMed 

Google Scholar 
Cortes, E. Perspectives on the intrinsic rate of population growth. Methods Ecol. Evol. 7, 1136–1145 (2016).
Google Scholar 
Ellis, J. et al. A review of capture and post-release mortality of elasmobranchs. J. Fish Biol. 90, 653–722 (2017).CAS 
PubMed 

Google Scholar 
Gallagher, A., Orbesen, E., Hammerschlag, N. & Serafy, J. Vulnerability of oceanic sharks as pelagic longline bycatch. Glob. Ecol. Conserv. 1, 50–59 (2014).
Google Scholar 
Afonso, A., Santiago, R., Hazin, H. & Hazin, F. Shark bycatch and mortality and hook bite-offs in pelagic longlines: Interactions between hook types and leader materials. Fish. Res. 131–133, 9–14 (2012).
Google Scholar 
Gilman, E., Chaloupka, M. & Musyl, M. Effects of pelagic longline hook size on species- and size-selectivity and survival. Rev. Fish Biol. Fish. 28, 417–433 (2018).
Google Scholar 
Epperly, S., Watson, J., Foster, D. & Shah, A. Anatomical hooking location and condition of animals captured with pelagic longlines: The grand banks experiments 2002–2003. Bull. Mar. Sci. 88, 513–527 (2012).
Google Scholar 
Amorim, S., Santos, M., Coelho, R. & Fernandez-Carvalho, J. Effects of 17/0 circle hooks and bait on fish catches in a southern Atlantic swordfish longline fishery. Aquat. Conserv. 25, 518–533 (2014).
Google Scholar 
Coelho, R., Fernandez-Carvalho, J., Lino, P. & Santos, M. An overview of the hooking mortality of elasmobranchs caught in a swordfish pelagic longline fishery in the Atlantic Ocean. Aquat. Living Resour. 25, 311–319 (2012).
Google Scholar 
Gilman, E. et al. A decision support tool for integrated fisheries bycatch management. Rev. Fish Biol. Fish. https://doi.org/10.1007/s11160-021-09693-5 (2022).Article 
PubMed 
PubMed Central 

Google Scholar 
Pascoe, S. et al. Use of incentive-based management systems to limit bycatch and discarding. Int. Rev. Environ. Resour. Econ. 4, 123–161 (2010).
Google Scholar 
Somers, K., Pfeiffer, L., Miller, S. & Morrison, W. Using incentives to reduce bycatch and discarding: Results under the west coast catch share program. Coast. Manag. 46, 1–17 (2019).
Google Scholar 
Abbott, J. & Wilen, J. Regulation of fisheries bycatch with common-pool output quotas. J. Environ. Econ. Manag. 57, 195–204 (2009).MATH 

Google Scholar 
Gilman, E. et al. Increasing the functionalities and accuracy of fisheries electronic monitoring systems. Aquat. Conserv. 29, 901–926 (2019).
Google Scholar 
Watling, J. Fishing observers ‘intimidated and bribed by EU crews’. Quota checks allegedly being compromised aboard Northwest Atlantic fishery boats, as observers report surveillance and theft. The Guardian (2012, accessed 21 July 2022). https://www.theguardian.com/environment/2012/may/18/fishing-inspectors-intimidated-bribed-crews.Clarke, S., Harley, S., Hoyle, S. & Rice, J. Population trends in Pacific oceanic sharks and the utility of regulations on shark finning. Conserv. Biol. 27, 197–209 (2013).PubMed 

Google Scholar 
Tolotti, M. T. et al. Banning is not enough: The complexities of oceanic shark management by tuna regional fisheries management organizations. Glob. Ecol. Conserv. 4, 1–7 (2015).
Google Scholar 
Gilman, E., Chaloupka, M., Merrifield, M., Malsol, N. & Cook, C. Standardized catch and survival rates, and effect of a ban on shark retention, Palau pelagic longline fishery. Aquat. Conserv. 26, 1031–1062 (2016).
Google Scholar 
CITES. Appendices I, II and III. Valid from 22 June 2021. Convention on International Trade in Endangered Species of Wild Fauna and Flora, United Nations Environment Program, Geneva (2021).Ward-Paige, C. A global overview of shark sanctuary regulations and their impact on shark fisheries. Mar. Policy 82, 87–97 (2017).
Google Scholar 
E.U. Regulation (E.U.) No 1380/2013 of the European Parliament and of the Council of 11 December 2013 on the Common Fisheries Policy, amending Council Regulations (EC) No 1954/2003 and (EC) No 1224/2009 and repealing Council Regulations (EC) No 2371/2002 and (EC) No 639/2004 and Council Decision 2004/585/EC. Official Journal of the European Union L354, 22–61 (2013).FAO. International Guidelines on Bycatch Management and Reduction of Discards (Food and Agriculture Organization of the United Nations, Rome, 2011).CCSBT. Resolution to Align CCSBT’s Ecologically Related Species Measures with those of other Tuna RFMOs (Commission for the Conservation of Southern Bluefin Tuna, Deakin West, Australia, 2021).IATTC. Active Resolutions and Recommendations (Inter-American Tropical Tuna Commission, La Jolla, 2022).ICCAT. Compendium. Management Recommendations and Resolutions Adopted by ICCAT for the Conservation of Atlantic Tunas and Tuna-like Species (International Commission for the Conservation of Atlantic Tunas, Madrid, 2021).IOTC. Compendium of Active Conservation and Management Measures for the Indian Ocean Tuna Commission (Indian Ocean Tuna Commission, Mahe, 2021).WCPFC. Conservation and Management Measures and Resolutions of the Western and Central Pacific Fisheries Commission. Compiled 31 August 2021 (Western and Central Pacific Fisheries Commission, Kolonia, Federated States of Micronesia, 2021).Faith, D. Threatened species and the potential loss of phylogenetic diversity: Conservation scenarios based on estimated extinction probabilities and phylogenetic risk analysis. Conserv. Biol. 22, 1461–1470 (2008).PubMed 

Google Scholar 
Dolce, J. & Wilga, C. Evolutionary and ecological relationships of gill slit morphology in extant sharks. Bull. Mus. Comp. 161, 79–109 (2013).
Google Scholar 
MacLeod, N. & Forey, P. Morphology, Shape and Phylogeny (CRC Press, 2002).
Google Scholar 
Haddaway, N., Macura, B., Whaley, P. & Pullin, A. ROSES RepOrting standards for Systematic Evidence Syntheses: pro forma, flow-diagram and descriptive summary of the plan and conduct of environmental systematic reviews and systematic maps. Environ. Evid. https://doi.org/10.1186/s13750-018-0121-7 (2018).Article 

Google Scholar 
Pullin, A., Frampton, G., Livoreil, B. & Petrokofsky, G., Eds. Section 5. Conducting a Search. Key CEE Standards for Conduct and Reporting. In Pullin, A., Frampton, G., Livoreil, B., Petrokofsky, G., Eds. Guidelines and Standards for Evidence Synthesis in Environmental Management. Version 5.0. Collaboration for Environmental Evidence (2020).Pullin, A., Frampton, G., Livoreil, B. & Petrokofsky, G. (eds) Section 3. Planning a CEE Evidence Synthesis. In Pullin, A., Frampton, G., Livoreil, B., Petrokofsky, G. (eds) Guidelines and Standards for Evidence Synthesis in Environmental Management. Version 5.0. Collaboration for Environmental Evidence (2021).Page, M. et al. The PRISMA statement: An updated guideline for reporting systematic reviews. BMJ https://doi.org/10.1136/bmj.n.71 (2020).Article 
PubMed 

Google Scholar 
Tuyl, F., Gerlach, R. & Mengersen, K. Comparison of Bayes-Laplace, Jeffreys, and other priors: The case of zero events. Am. Stat. 62, 40–44 (2008).MathSciNet 

Google Scholar 
Dorai-Raj, S. binom: Binomial confidence intervals for several parameterizations. R package version 1.1-1. https://CRAN.R-project.org/package=binom (2014).van Lissa, C. Small sample meta-analyses: Exploring heterogeneity using MetaForest. Chapter 13. In Small Sample Size Solutions: A Guide for Applied Researchers and Practitioners (eds Van De Schoot, R. & Miočević, M.) 186–202 (Routledge, Oxford, 2020).
Google Scholar 
Wright, M. & Ziegler, A. ranger: A fast implementation of random forests for high dimensional data in C++ and R. J. Stat. Softw. 77, 1–17 (2017).
Google Scholar 
Janitza, S., Celik, E. & Boulesteix, A. A computationally fast variable importance test for random forests for high-dimensional data. Adv. Data Anal. Classif. 12, 885–915 (2018).MathSciNet 
MATH 

Google Scholar 
Mayer, M. missRanger: Fast imputation of missing values. R package version 2.1.3. https://CRAN.R-project.org/package=missRanger (2021).Konstantopoulos, S. Fixed effects and variance components estimation in three-level meta-analysis. Research Synthesis. Methods 2, 61–76 (2011).
Google Scholar 
Amaral, C., Pereira, F., Silva, D., Amorim, A. & de Carvalho, E. The mitogenomic phylogeny of the Elasmobranchii (Chondrichthyes). Mitochondrial DNA A 29, 867–878 (2017).
Google Scholar 
Hara, Y. et al. Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates. Nat. Ecol. Evol. 2, 1761–1771 (2018).PubMed 

Google Scholar 
Naylor, G. et al. A DNA sequence-based approach to the identification of shark and ray species and its implications for global elasmobranch diversity and parasitology. Bull. Am. Mus. Nat. 367, 1–262 (2012).
Google Scholar 
Stein, R. et al. Global priorities for conserving the evolutionary history of sharks, rays and chimaeras. Nat. Ecol. Evol. 2, 288–298 (2018).PubMed 

Google Scholar 
Maddison, D., Swofford, D., Maddison, W. & Cannatella, D. Nexus: An extensible file format for systematic information. Syst. Biol. 46, 590–621 (1997).CAS 
PubMed 

Google Scholar 
Upham, N., Esselstyn, J. & Jetz, W. Inferring the mammal tree: Species-level sets of phylogenies for questions in ecology, evolution, and conservation. PLoS Biol. 17, e3000494 (2019).CAS 
PubMed 
PubMed Central 

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

Google Scholar 
Yu, G. Using ggtree to visualize data on tree-like structures. Curr. Protoc. Bioinform. 69, e96. https://doi.org/10.1002/cpbi.96 (2020).Article 

Google Scholar 
Hadfield, J. & Nakagawa, S. General quantitative genetic methods for comparative biology: Phylogenies, taxonomies and multi-trait models for continuous and categorical characters. J. Evol. Biol 23, 494–508 (2010).CAS 
PubMed 

Google Scholar 
Lin, L. & Chu, H. Meta-analysis of proportions using generalized linear mixed models. Epidemiology 31, 713–717 (2020).PubMed 
PubMed Central 

Google Scholar 
Carpenter, B. et al. Stan: A probabilistic programming language. J. Stat. Softw. 76, 1–32 (2017).
Google Scholar 
Bürkner, P. brms: An R Package for Bayesian multilevel models using Stan. J. Stat. Softw. 81, 1–28 (2017).
Google Scholar 
Günhan, B., Röver, C. & Friede, T. Random-effects meta-analysis of few studies involving rare events. Res. Synth. Methods 11, 74–90 (2020).PubMed 

Google Scholar 
Pappalardo, P. et al. Comparing traditional and Bayesian approaches to ecological meta-analysis. Methods Ecol. Evol. 11, 1286–1295 (2020).
Google Scholar 
Ott, M., Plummer, M. & Roos, M. How vague is vague? How informative is informative? Reference analysis for Bayesian meta-analysis. Stat. Med. 40, 4505–4521 (2021).MathSciNet 
PubMed 
PubMed Central 

Google Scholar 
Wood, S. Generalized Additive Models: An Introduction with R 2nd edn. (Chapman and Hall, 2017).MATH 

Google Scholar 
Kruschke, J. & Liddell, T. The Bayesian new statistics: Hypothesis testing, estimation, meta-analysis, and power analysis from a Bayesian perspective. Psychon. Bull. Rev. 25, 178–206 (2018).PubMed 

Google Scholar 
Kay, M. tidybayes: Tidy data and geoms for Bayesian models. R package version 2.1.1. https://doi.org/10.5281/zenodo.1308151 (2020).Makowski, D., Ben-Shachar, M. & Lüdecke, D. bayestestR: Describing effects and their uncertainty, existence and significance within the Bayesian framework. J. Open Source Softw. 4, 1541 (2019).ADS 

Google Scholar 
Searle, S., Speed, F. & Milliken, G. Population marginal means in the linear model: An alternative to least squares means. Am. Stat. 34, 216–221 (1980).MathSciNet 
MATH 

Google Scholar 
Lenth, R. emmeans: Estimated marginal means, aka least-squares means. R package version 1.5.2-1. https://CRAN.R-project.org/package=emmeans (2020).Pagel, M. Inferring the historical patterns of biological evolution. Nature 401, 877–884 (1999).ADS 
CAS 
PubMed 

Google Scholar 
Münkemüller, T. et al. How to measure and test phylogenetic signal. Methods Ecol. Evol. 3, 743–756 (2012).
Google Scholar 
Vehtari, A., Gelman, A. & Gabry, J. Practical Bayesian model evaluation using leave-one-out cross-validation and WAIC. Stat. Comput. 27, 1413–1432 (2017).MathSciNet 
MATH 

Google Scholar 
Yao, Y., Vehtari, A., Simpson, D. & Gelman, A. Using stacking to average Bayesian predictive distributions (with Discussion). Bayesian Anal. 13, 917–1003 (2018).MathSciNet 
MATH 

Google Scholar 
Gabry, J., Simpson, D., Vehtari, A., Betancourt, M. & Gelman, A. Visualization in Bayesian workflow. J. R. Soc. Ser. A 182, 1–14 (2019).MathSciNet 

Google Scholar 
Lazic, S., Mellor, J., Ashby, M. & Munafo, M. A Bayesian predictive approach for dealing with pseudoreplication. Sci. Rep. 10, 2020. https://doi.org/10.1038/s41598-020-59384-7 (2020).Article 
CAS 

Google Scholar 
Page, M., Sterne, J., Higgins, J. & Egger, M. Investigating and dealing with publication bias and other reporting biases in meta-analyses of health research: A review. Res. Synth. Methods 12, 248–259 (2021).PubMed 

Google Scholar 
Peters, J., Sutton, A., Jones, D., Abrams, K. & Rushton, L. Contour-enhanced meta-analysis funnel plots help distinguish publication bias from other causes of asymmetry. J. Clin. Epidemiol. 61, 991–996 (2008).PubMed 

Google Scholar 
Viechtbauer, W. Conducting meta-analyses in R with the metafor package. J. Stat. Softw. 36, 1–48 (2010).
Google Scholar 
Gasparrini, A., Armstrong, B. & Kenward, M. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat. Med. 31, 3821–3839 (2012).MathSciNet 
CAS 
PubMed 
PubMed Central 

Google Scholar  More

FAO. The State of the World’s Forests 2018—Forest Pathways to Sustainable Development (FAO, 2018).
Google Scholar 
Pan, Y. et al. A large and persistent carbon sink in the world’s forests. Science 333(6045), 988–993. https://doi.org/10.1126/science.1201609 (2011).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Luyssaert, S. et al. Tradeoffs in using European forests to meet climate objectives. Nature 562(7726), 259–262. https://doi.org/10.1038/s41586-018-0577-1 (2018).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Büntgen, U. et al. Limited capacity of tree growth to mitigate the global greenhouse effect under predicted warming. Nat. Commun. https://doi.org/10.1038/s41467-019-10174-4 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Anderegg, W. R. L. et al. Climate-driven risks to the climate mitigation potential of forests. Science 368, 1327. https://doi.org/10.1126/science.aaz7005 (2020).Article 
CAS 

Google Scholar 
Jiang, M. et al. The fate of carbon in a mature forest under carbon dioxide enrichment. Nature 580(7802), 227–231. https://doi.org/10.1038/s41586-020-2128-9 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Naudts, K. et al. Europe’s forest management did not mitigate climate warming. Science 351(6273), 597–599. https://doi.org/10.1126/science.aad7270 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Tong, X. et al. Forest management in southern China generates short term extensive carbon sequestration. Nat. Commun. https://doi.org/10.1038/s41467-019-13798-8 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Yu, K. et al. Effects of stand age on soil respiration in Pinus massoniana plantations in the hilly red soil region of Southern China. CATENA 178, 313–321. https://doi.org/10.1016/j.catena.2019.03.038 (2019).Article 
CAS 

Google Scholar 
Mei, G., Sun, Y. & Sajjad, S. Models for predicting the biomass of Cunninghamia lanceolata trees and stands in southeastern China. PLoS ONE 12, e0169747. https://doi.org/10.1371/journal.pone.0169747 (2017).Article 
CAS 

Google Scholar 
Wu, H. et al. Soil phosphorus bioavailability and recycling increasedwith stand age in Chinese fir plantations. Ecosystems 23, 973–988. https://doi.org/10.1007/s10021-019-00450-1 (2019).Article 

Google Scholar 
State Forestry Administration. General situation of forest resources in China. The 8th National Forest Inventory (State Forestry Administration, 2014).Wang, X. et al. Vegetation carbon storage and density of forest ecosystems in China. Chin. J. Appl. Ecol. 12(1), 13–16 (2001) (in Chinese with English Abstract).ADS 
CAS 

Google Scholar 
Kang, H. et al. Simulating the impact of climate change on the growth of Chinese fir plantations in Fujian province, China. NZ J. For. Sci. 47(1), 20. https://doi.org/10.1186/s40490-017-0102-6 (2017).Article 

Google Scholar 
Lu, Y. et al. A process-based approach to estimate Chinese fir (Cunninghamia lanceolata) distribution and productivity in southern China under climate change. Forests 6, 360–379. https://doi.org/10.3390/f6020360 (2015).Article 

Google Scholar 
Zhang, X. et al. Relative contributions of competition, stand structure, age, and climate factors to tree mortality of Chinese fir plantations: Long-term spacing trials in southern China. For. Ecol. Manag. 465, 118103. https://doi.org/10.1016/j.foreco.2020.118103 (2020).Article 

Google Scholar 
You, R. et al. Variation in wood physical properties and effects of climate for different geographic sources of Chinese fir in subtropical area of China. Sci. Rep. https://doi.org/10.1038/s41598-021-83500-w (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Djomo, A. N., Ibrahima, A., Saborowski, J. & Gravenhorst, G. Allometric equations for biomass estimations in Cameroon and pan moist tropical equations including biomass data from Africa. For. Ecol. Manag. 260(10), 1873–1885. https://doi.org/10.1016/j.foreco.2010.08.034 (2010).Article 

Google Scholar 
Peng, D. et al. Estimating the aboveground biomass for planted forests based on stand age and environmental variables. Remote Sens. 11(19), 2270. https://doi.org/10.3390/rs11192270 (2019).Article 
ADS 

Google Scholar 
Zhou, X. et al. Dynamic allometric scaling of tree biomass and size. Nat. Plants. 7(1), 42–49. https://doi.org/10.1038/s41477-020-00815-8 (2021).Article 
PubMed 

Google Scholar 
Li, L. Study on the tree volume table compilation of Chinese fir in Kaihua Forest Farm (Beijing Forestry University, 2011) http://cdmd.cnki.com.cn/Article/CDMD-10022-1011134655.htm (in Chinese).Wang, J. P. et al. Study on the effect of Chinese fir volume formula on estimating the volume of fir standing trees in different sites. Guizhou For. Technol. 19(1), 26–29 (1991) (in Chinese).
Google Scholar 
Zeng, W. S. et al. Compatible tree volume and aboveground biomass equations for Chinese fir plantation in Guizhou. J. Beijing For. Univ. 33(4), 1–6 (2011) (in Chinese).
Google Scholar 
Xia, Z. S. et al. Construction of tree volume equations for Chinese fir plantation in Guizhou Province, southwestern China. J. Beijing For. Univ. 34(1), 1–5 (2012) (in Chinese).
Google Scholar 
Lin, H. Study on biomass and carbon storage of main coniferous forest in Jiangle state-owned forestry farm. J. Fujian For. Sci. Technol. 45(1), 30–34. https://doi.org/10.13428/j.cnki.fjlk.2018.01.007 (2018) (in Chinese with English Abstract).Article 
ADS 

Google Scholar 
Cai, Z. A study on biomass models of Cunninghamia lanceolata plantation in Fujian. (Beijing Forestry University, 2014), http://cdmd.cnki.com.cn/Article/CDMD-10022-1014327550.htm (in Chinese).Chen, G. et al. Carbon storage in a chronosequence of Chinese fir plantations in southern China. For. Ecol. Manag. 300, 68–76. https://doi.org/10.1016/j.foreco.2012.07.046 (2013).Article 

Google Scholar 
Zhang, G. et al. Biomass Characteristics of dominant tree species (group) at Lingnan forest farm in Anhui province. Scientia Silvae Sinicae. 48(5), 136–140. https://doi.org/10.1007/s11783-011-0280-z (2012) (in Chinese with English abstract).Article 
ADS 
CAS 

Google Scholar 
Shi, W. et al. Biomass model and carbon storage of Chinese fir plantation in Dabieshan Mountains in Anhui. Resour. Environ. Yangtze Basin. 24(5), 758–764. https://doi.org/10.11870/cjlyzyyhj201505007 (2015) (in Chinese with English abstract).Article 

Google Scholar 
Li, H. & Zhao, P. Improving the accuracy of tree-level aboveground biomass equations with height classification at a large regional scale. For. Ecol. Manag. 289, 153–163. https://doi.org/10.1016/j.foreco.2012.10.002 (2013).Article 

Google Scholar 
Zeng, W. & Tang, S. A new general allometric biomass model. Nat. Precedings. https://doi.org/10.1038/npre.2011.6704.1 (2011).Article 

Google Scholar 
Schumacher, F. X. & Hall, F. D. S. Logarithmic expression of timber-tree volume. J. Agric. Res. 47(9), 719–734 (1933).
Google Scholar 
Honer, T. G. A new total cubic foot volume function. For. Chron. 41(4), 476–493. https://doi.org/10.5558/tfc41476-4 (1965).Article 

Google Scholar 
Burkhart, H. E. Cubic-foot volume of loblolly pine to any merchantable top limit. South. J. Appl. For. 2, 7–9. https://doi.org/10.1093/sjaf/1.2.7 (1977).Article 

Google Scholar 
Lee, D., Seo, Y. & Choi, J. Estimation and validation of stem volume equations for Pinus densiflora, Pinus koraiensis, and Larix kaempferi in South Korea. For. Sci. Technol. 13(2), 77–82. https://doi.org/10.1080/21580103.2017.1315963 (2017).Article 

Google Scholar 
Chen, B. H. & Chen, C. Y. A preliminary study on the biomass and productivity of Picea koraiensis forests in the dunes. Scientia Silvae Sinicae 4, 269–278 (1980) (in Chinese).
Google Scholar 
Niklas, K. J. Plant Allometry: The Scaling of Form and Process (University of Chicago Press, 1994).
Google Scholar 
Ketterings, Q. M. et al. Reducing uncertainty in the use of allometric biomass equations for predicting above-ground tree biomass in mixed secondary forests. For. Ecol. Manag. 146, 199–209. https://doi.org/10.1016/S0378-1127(00)00460-6 (2001).Article 

Google Scholar 
Chen, X. G. The biomass and allometric equation of a 20-years-old Cunninghamia lanceolata plantation. Prot. For. Sci. Technol. 4, 28–29, 40. https://doi.org/10.3969/j.issn.1005-5215.2007.04.010.(inChinese) (2007).Article 

Google Scholar 
Wang, X. P. et al. Climatic control of primary forest structure and DBH–height allometry in Northeast China. For. Ecol. Manag. 234, 264–274. https://doi.org/10.1016/j.foreco.2006.07.007 (2006).Article 

Google Scholar 
Peng, C. et al. Developing and evaluating tree height–diameter models at three geographic scales for black spruce in Ontario. N. J. Appl. For. 21(2), 83–92. https://doi.org/10.1093/njaf/21.2.83 (2004).Article 

Google Scholar 
López-Serrano, F. R. et al. Site and weather effects in allometries: A simple approach to climate change effect on pines. For. Ecol. Manag. 215(1–3), 251–270. https://doi.org/10.1016/j.foreco.2005.05.014 (2005).Article 

Google Scholar 
Zhang, C. et al. Developing aboveground biomass equations both compatible with tree volume equations and additive systems for single trees in Poplar plantations in Jiangsu Province, China. Forests 7, 32. https://doi.org/10.3390/f7020032 (2016).Article 

Google Scholar 
Liu, J. C. et al. Comparing non-destructive methods to estimate volume of three tree taxa in Beijing, China. Forests 10, 92. https://doi.org/10.3390/f10020092 (2019).Article 

Google Scholar 
Thangjam, U. et al. Developing tree volume equation for Parkia timoriana grown in home gardens and shifting cultivation areas of North-East India. For. Trees Livelihoods 28(12), 1–13. https://doi.org/10.1080/14728028.2019.1624200 (2019).Article 

Google Scholar 
Dutcă, I. et al. Does slope aspect affect the aboveground tree shape and volume allometry of European Beech (Fagus sylvatica L.) trees?. Forests 13, 1071. https://doi.org/10.3390/f13071071 (2022).Article 

Google Scholar 
Segura, M. & Kanninen, M. Allometric models for tree volume and total aboveground biomass in a tropical humid forest in Costa Rica. Biotropica 37(1), 2–8. https://doi.org/10.2307/30045500 (2005).Article 

Google Scholar 
Wang, X. W. et al. Additive tree biomass equations for Betula platyphylla Suk. plantations in Northeast China. Ann. For. Sci. 75, 60. https://doi.org/10.1007/s13595-018-0738-2 (2018).Article 

Google Scholar 
Niklas, K. J. & Enquist, B. J. Canonical rules for plant organ biomass partitioning and annual allocation. Am. J. Bot. 89(5), 812–819. https://doi.org/10.3732/ajb.89.5.812 (2002).Article 
PubMed 

Google Scholar 
Xiang, W. H. et al. General allometric equations and biomass allocation of Pinus massoniana trees on a regional scale in southern China. Ecol. Res. 26, 697–711. https://doi.org/10.1007/s11284-011-0829-0 (2011).Article 

Google Scholar 
Brown, S. Measuring carbon in forests: Current status and future challenges. Environ. Pollut. 116, 363–372. https://doi.org/10.1016/s0269-7491(01)00212-3 (2002).Article 
CAS 
PubMed 

Google Scholar 
Brassard, B. W. et al. Influence of environmental variability on root dynamics in northern forests. Crit. Rev. Plant Sci. 28, 179–197. https://doi.org/10.1080/07352680902776572 (2009).Article 

Google Scholar 
Montagu, K. D. et al. Developing general allometric relationship for regional estimates of carbon sequestration—An example using Eucalyptus pilularis from seven contrasting sites. For. Ecol. Manag. 204, 113–127. https://doi.org/10.1016/j.foreco.2004.09.003 (2005).Article 

Google Scholar 
Williams, R. J. et al. Allometry for estimating aboveground tree biomass in tropical and subtropical eucalypt woodlands: Towards general predictive equations. Aust. J. Bot. 53, 607–619. https://doi.org/10.1071/BT04149 (2005).Article 

Google Scholar 
Ouimet, R. et al. Estimation of coarse root biomass and nutrient content for sugar maple, jack pine, and black spruce using stem diameter at breast height. Can. J. For. Res. 38, 92–100. https://doi.org/10.1139/x07-134 (2008).Article 

Google Scholar 
Peichl, M. & Arain, M. A. Allometry and partitioning of above-and belowground tree biomass in an age-sequence of white pine forests. For. Ecol. Manag. 253, 68–80. https://doi.org/10.1016/j.foreco.2007.07.003 (2007).Article 

Google Scholar 
Bond-Lamberty, B. et al. Aboveground and below-ground biomass and sapwood area allometric equations for six boreal tree species of northern Manitoba. Can. J. For. Res. 32, 1441–1450. https://doi.org/10.1139/x02-063 (2002).Article 

Google Scholar 
King, J. S. et al. Biomass partitioning in red pine (Pinus resinosa) along a chronosequence in the Upper Peninsula of Michigan. Can. J. For. Res. 37(1), 93–102. https://doi.org/10.1139/x06-217 (2007).Article 

Google Scholar 
Ziania, D. & Mencuccini, M. Aboveground biomass relation-ships for beech (Fagus moesiaca Cz.) trees in Vermio Mountain, northern Greece, and generalised equations for Fagus sp. Ann. For. Sci. 60(5), 439–448. https://doi.org/10.1051/forest:2003036 (2003).Article 

Google Scholar 
Martin, J. G. et al. Aboveground biomass and nitrogen allocation of ten deciduous southern Appalachian tree species. Can. J. For. Res. 28(11), 1648–1659. https://doi.org/10.1139/x98-146 (1998).Article 

Google Scholar 
Wang, C. Biomass allometric equations for 10 co-occurring tree species in Chinese temperate forests. For. Ecol. Manag. 222, 9–16. https://doi.org/10.1016/j.foreco.2005.10.074 (2006).Article 

Google Scholar  More

The Santa Barbara Coastal Long Term Ecological Research program has monitored benthic communities in five kelp forests seasonally since 2008 using fixed transect diver surveys, and moored sensors at each reef have recorded bottom temperatures every 15 min. Blob-associated positive bottom temperature anomalies began in winter 2014 and persisted through autumn 2016 (Fig. 1a)18. Peak temperature anomalies occurred during the summer and autumn of 2014 and 2015 (Fig. 1a), and the average temperature anomaly in autumn 2015 was +3.1 °C, equivalent to an average daily temperature of 19.6 °C. In 2014 and 2015, 91 and 69% of autumn days, respectively, were classified as heatwave days as defined by Hobday et al.20. Seasonal chlorophyll-a concentration, a proxy for phytoplankton abundance, was obtained from satellite imagery at each of the five reefs over the 14-year period. The average chlorophyll-a concentration was anomalously low throughout the warming period, and exceptionally low during the springs of 2014 and 2015 (Fig. 1a), when upwelling-driven nutrient enrichment typically supports dense phytoplankton blooms.Fig. 1: Average seasonal bottom temperature anomaly, chlorophyll-a concentration anomaly, and percent cover and species richness of sessile invertebrates across five sites.The Blob, an anomalous warming period from spring of 2014 to winter of 2016, is highlighted in gray, coincident with (a) positive temperature anomalies (°C; solid line), negative chlorophyll-a anomalies (mg/m3; dashed line), and declines in (b) invertebrate cover (solid line) and species richness (number of unique species/taxa/80 contact points; dashed line). Seasons are denoted by Sp (Spring), Su (Summer), A (Autumn) and W (Winter).Full size imageMean sessile invertebrate cover averaged across all sites declined 71% during the Blob, reaching a 14-year minimum of 7% in autumn of 2015 (Fig. 1b and Supplementary Fig. 1). Species richness declined 69% during the same period (Fig. 1b and Supplementary Fig. 1). The responses of invertebrates to warming were not consistent across time even though the duration and intensity of warming was similar in 2014 and 2015, suggesting that extended periods of elevated seawater temperature were not solely responsible for the most severe loss of invertebrates. For ectotherms, increases in ambient seawater temperature should be met with increases in metabolic rate and food requirements to sustain metabolism21. Because of their sedentary lifestyle, sessile invertebrates cannot actively forage for food or seek spatial refuge from thermal extremes, and limitations in their planktonic food supply can result in metabolic stress over extended periods22,23. Anomalously low chlorophyll-a concentrations during the Blob (Fig. 1a), particularly in the spring of 2015, indicated that food limitation was a likely driver of invertebrate decline. Results from piecewise structural equation modeling (Fig. 2) that incorporated biological interactions with space competitors (understory macroalgae), predators (sea urchins), and foundation species (giant kelp) showed that the severity of warming had both a direct and indirect effect on the sessile invertebrate community. The proportion of heatwave days was a direct negative predictor of sessile invertebrate cover (−0.11) and species richness (−0.21). The proportion of heatwave days was an even stronger negative predictor of chlorophyll-a concentration (−0.26), yielding negative indirect effects on invertebrate cover (−0.07) and species richness (−0.05) due to the positive influence of chlorophyll-a concentration on sessile invertebrate cover (+0.26) and richness (+0.20).Fig. 2: Piecewise structural equation modeling (SEM) for sessile invertebrate cover and species richness.Arrows indicate directionality of effects on (a) invertebrate cover and (b) species richness. Red arrows show negative relationships; black arrows show positive relationships. R2 values are conditional R2. Arrow widths are proportional to effect sizes as measured by standardized regression coefficients (shown next to arrows). ***p  More

Conception of the workflow to demonstrate the microbial associations from co-occurrence networks with microbial cultivationMicrobial co-occurrence networks are composed of nodes and edges, which usually represent microbes and statistically significant associations between microbes, respectively. We hypothesized that the microbial associations could be validated if the topological properties of networks are simplified, and if the microbes representing the nodes can be cultivated. To test this hypothesis, we designed a workflow as shown in Fig. 1. A total of 12,096 wells from 126 96-well plates were inoculated with droplets of series diluted environmental samples, wells from each 96-well plate represented the same combination of given culture condition, sample type (plants, roots, and sediments) and dilution rate (from 10–1 to 10–7). After being cultivated at 30 °C for 10 days, 69 effective (Supplementary Table S3) plates with  > 30% wells showing microbial growth were retained for downstream microbial community analysis. Microbial DNA in each well was extracted, bar-coded, and sequenced for the inference of co-occurrence networks. The wells of plates showing high abundances of target Zotus were targeted for microbial isolations. Lastly, the cultivated microbial isolates were matched to Zotus in the network and used for demonstration of microbial interactions.Figure 1Overview of experimental demonstration of microbial interactions in co-occurrence networks. For detailed description, please refer to the method section.Full size imagePrevalent Zotu pairs in the co-occurrence networksDepending on the microbial density in samples, the 96-well plates harbored different numbers of wells with microbial growth. We obtained 65 96-well plates (6,091 wells) that were effective with microbial growth and data analysis for co-occurrence network reconstruction. After quality control and denoise, we obtained 130 Gbp sequence data. A total of 14,377 Zotus were annotated (Supplementary Table S4). There were 217 ± 94 (average ± standard deviation) prevalent Zotus, i.e., these Zotus appeared at frequencies ≥ 30% of wells in a given 96-well plate.Next, we analyzed Zotus compositions and abundances in each well of the 65 plates. Accordingly, we reconstructed 65 independent microbial co-occurrence networks and further retrieved the robust (Spearman’s |ρ| > 0.6 and P  More

“This new study is inspired by some our earlier theoretical work applied to killer whales that suggests that age-related changes in relatedness are important for the evolution of menopause,” says Samuel Ellis, the first author of the study. “Reproduction can be thought of as a form of generalized harm as the birth of an offspring increases within-group competition for resources. Kinship dynamics — the ways in which local relatedness changes over an individual’s lifetime — are one way that menopause could be favored, because older females are more inclined to cease reproduction to not harm their group mates than younger females. Here we wanted to generalize this concept to both sexes, and to other species without menopause.” More

After massive parallel sequencing of the nucleic acids obtained from fur seal scats, a wide variety of invertebrate and vertebrate viral hosts assignations with low nucleotidic and amino-acidic identities were obtained, most of them corresponding to animal species not described before in Antarctica. These results make us reconsider the use of closed RefSeq databases for viral discovery, especially because the studied area was a remote geographical area where a high number of new viral species is expected to occur22.After repeating the analysis of the contigs obtained using BLASTn, a high number of miss-assignments was observed, corresponding almost entirely to contigs newly assigned as unclassified Eukaryotic Circular Rep-Encoding Single-Stranded DNA (CRESS-DNA) viral sequences. CRESS viruses have been detected ubiquitously in many different animals without any recognised role in the development of any disease23,24,25,26.These results are in accordance with the recent reporting of CRESS sequences also being ubiquitous in a wide variety of environments and at high proportions, including Antarctica, where they have been described to represent more than 50% of sequences obtained from glacier waters27.Viral-host distributionVirome studies in other Arctocephalus species from subantarctic and South American regions revealed a 5% of viral sequences with predominance of bacteriophages followed by viruses from the Parvoviridae family28. The methodology here applied provided an increase of 12–25% viral reads when probe-based Target Enrichment Sequencing (TES) was applied, that in comparison with Untargeted Viral Metagenomics (UVM) approaches conducted in these type of samples28 could be considered an optimal result.Most of the viral species detected in feces corresponded to unknown viruses, 83.59% from the total of sequences, followed by viruses that infect invertebrates, 8.75%, bacteriophages, 4.46%, and vertebrate viruses, 3.11% (Fig. 1).Figure 1Host distribution of viral assignations sequenced from fecal (A) and serum (B) samples collected from male A. gazella.Full size imageAs expected, when applying both targeted and untargeted sequencing methodologies, TES approach resulted in a recovery of many vertebrate viral assignations (Table 1) whereas untargeted sequencing enabled a better detection of viruses known to infect invertebrates (Table 2). To describe the complete A. gazella fecal virome, sequences obtained by both sequencing methodologies were considered all together, representing a total of 2.62 million reads.Table 1 Vertebrate viral assignations obtained from fecal samples sequencing from male A. gazella. Ranges of Genome coverage, nucleotide identity and aminoacidic identity are expressed in percentages.Full size tableTable 2 Invertebrate viral assignations obtained from fecal samples sequencing from male A. gazella. Colours represent the presence of each assignation in the processed pools. Ranges of Coverage, NT ID and AA ID are represented in percentages.Full size table
A. gazella virusesFur seal picorna-like virusFur seal picorna-like virus was firstly described in a fecal sample obtained from A. gazella in King George Island in the South Shetland Islands, Antarctica by Krumbholz and co-workers16.In this study, we report a total of 19 contigs resulting after assembling 2671 reads obtained from 4/4 fecal pools analysed being the most prevalent virus described in this study. One of the contigs covered 96.91% of the fur seal picorna-like virus genome and presented a nucleotide homology of 99.38% with the reference strain described in 2017. The other contigs coverage ranged from 19.75 to 21.22% with a 45.92 to 90.5% nucleotide identity with reference strain NC_035110. Four contigs matching the ORF2 polyprotein are represented in Fig. 2 where differences among them and with the reference strain are showed.Figure 2Nucleotide alignment of ORF2 sequences from the A. gazella picorna-like contigs compared to the ORF2 from RefSeq NC_0351110. In consensus strain, position 1 represents position 6523 from RefSeqs genome.Full size imagePicornaviruses are known to cause a wide variety of diseases in vertebrate hosts, especially mammals29, but the role of Fur seal picorna-like virus in pathogenesis development is still unknown30. Many picornaviruses are transmitted horizontally via fecal–oral or airborne routes29. The fact that these sequences were detected in all the fecal pools obtained from animals with no evidence of disease may that suggest the virus may have a stable endemic relationship within that seal population.Torque teno pinniped virusLambdatorquevirus is a genus within the Anelloviridae family. The genus comprises 8 species named Torque teno pinniped virus 2 to 9 isolated from different pinniped species: A. gazella (Torque teno pinniped virus 6 and 7)17, Phoca vitulina (Torque teno pinniped virus 2, 3, 4)31, Zalophus californianus (Torque teno pinniped virus 5)32 and Leptonychotes weddellii (Torque teno pinniped 8 and 9)33.One contig with a nuleotide similarity of 95.12% against Torque teno pinniped virus 7 was obtained from one of the fecal pools. This virus had been described in these animals inhabiting Livingston Island in 2016, using rolling circle amplification and subsequent Sanger sequencing from buccal swabs17. However, sequences obtained in this study belong to partial ORF2 which is not the optimal genome region for typing purposes or phylogenetic analysis.These members of the Anelloviridae represent the more abundant viruses found in human, animals and environmental samples although their etiological role in any disease has not been clearly identified being considered a persistent virus ubiquitous to several different tissues34,35No Torque teno virus sequences were detected in serum samples which agree with what was observed for Zalophus californianus anellovirus prevalently detected in different tissues, like lung and liver, but not in blood samples. Interestingly, other known anelloviruses are typically found in blood or plasma samples32.MamastrovirusTwo of the fecal pools analyzed presented Mamastrovirus sequences. The presence of these viruses in humans and other mammals is widely known, as well as their involvement in gastroenteritis development36. The four contigs obtained (comprising 1008 sequences) showed homologies against reference genomes, ranging from 45.70% to 59.37% when compared at nucleotide level and 36.69% to 46.69% when compared at aminoacidic level. Phylogenetic analysis of partial OFR2 regions of these contigs indicate its closer similarity with sequences from California Sea Lion astroviruses, a virus that was determined as to be the most prevalent in fecal samples from these animals (Z. californianus)37. This finding suggests that these sequences may belong to a yet unknown virus like Z. californianus astrovirus and may indicate that such virus is prevalent in the sampled area (detected in 2/4 fecal pools studied) and the second more abundant virus (1008 reads) in the studied fecal samples (Fig. 3).Figure 3Phylogenetic consensus tree based on partial ORF2 sequences from the Mamastrovirus contigs sequenced from A. gazella scats (in bold). Bootstrap resampling with 1000 replicates.Full size imageAdeno associated virus 2Two of the studied fecal pools presented 138 sequences, forming 3 contigs with nucleotide identities ranging from 46.91 to 48.04% (Table 1), that matched adeno associated viruses previously described in Z. californianus, humans and other mammals with and unknow etiologic role (Fig. 4). The detected sequences probably correspond to fur seal adeno associated viruses never described before. The detection of these viruses is quite common in other mammals suggesting they could cause persistent infections in their hosts, but no etiological role has been attributed to them38.Figure 4Phylogenetic consensus tree of the Adeno-associated virus contigs sequenced from A. gazella scats (in bold). Bootstrap resampling with 1000 replicates.Full size imageNorovirusA norovirus contig was obtained in one of the four pools analyzed. Noroviruses are the most relevant non-bacterial gastroenteritis etiological agents in humans39, with its presence widely described in other mammals40. The contig detected in the fecal samples, represented the 4.43% of the viral genome, was in the VP1 region and comprised 56 reads with an identity  > 99% to California sea lion norovirus described by Teng and collaborators in 201841 (Fig. 5). Results obtained suggest these sequences belong to a putative new norovirus specie.Figure 5Phylogenetic consensus tree of the Norovirus contig sequenced from A. gazella scats (in bold). Bootstrap resampling with 1000 replicates.Full size imageViruses in serum samplesAll the viral sequences obtained from serum samples (970 reads) matched to CRESS-DNA viral sequences from unknown hosts.The fact that no other viruses were identified in serum samples suggests the animals tested were not under active viremia at the time of sample collection or it was not detectable by the applied methodology.Diet related virusesSeveral virus sequences similar to viruses known to have invertebrate animals as hosts were detected in fecal pools, mainly by UVM although some also by TES. These viruses are probably present in fur seal feces because of dietary habits although, since scats were collected from the ground nearby the animals, environmental cross-contamination should not be ruled out.Sequences with high coverage or similarities to any described virus are showed in Table 2.The high prevalence of virus sequences from crustaceans in the feces analyzed is hardly surprising because A. gazella inhabiting the Antarctic peninsula and the Atlantic sector of the Southern Ocean feed mostly on Antarctic krill Euphasia superba during the summer months42,43,44,45,46,47,48. Sequences from cephalopod viruses were also detected, although were much scarcer than those from crustaceans. This also agrees with current knowledge about the diet of A. gazella in the Atlantic sector of the Southern Ocean, where octopuses and squids are regularly consumed, although in low numbers44,45,46. It is worth noting than not cephalopod beak was recovered from the scats analyzed here48. Among all invertebrate viruses identified, some sequences present low identities with genomes from available databases, probably because Antarctica wildlife has been scarcely explored, forcing bioinformatic analysis to match them with the most similar viruses from these databases.No fish viruses were found in this study. Hard skeletal remains of fishes are often recovered from the scats of A. gazella from the Atlantic sector of the Southern Ocean42,43,44,45,46,47 and occurred indeed in the samples analysed here48, but stable isotope analysis of blood and whiskers revealed a negligible contribution of fish to the assimilate diet of juvenile and subadult male A. gazella49, which likely explain the absence of fish viruses in the samples analized here. Additionaly, no data on the virome present in the fish species regularly consumed by A. gazella has been published to our knowledge, with information limited to the bacteriome32, so even in case fish viruses were sequenced, it might not be correctly assigned to a fish host. Nevertheless, the methodology applied in this study had been successfully applied to the identification of the virome of Atlantic fishes50. Furthermore, Li and coworkers.37 and Wille and coworkers.22 also observed viral sequences probably corresponding to fish when analyzing the fecal virome of the California sea lions and Antarctic penguins.On the other hand, sequences highly similar to Coelho and Khabarov viral polymerases (greater than 98% of aminoacid identity), previously described in chinstrap penguins (Pygoscelis antarcticus) by Wille and coworkers22, were found in this study. The consumption of penguins by A. gazella during the summer months has been reported widely51,52,53,54,55, penguins feathers were reported from the scats analyzed in this study48 and stable isotope analysis of blood and whiskers revealed penguins as the second most relevant prey from juvenile and subadult male A. gazella in the population studied here49. This evidence is consistent with the presence of virus from chinstrap penguins in the samples analysed here. All in all, the study of fecal virome constitutes a very promising tool to explore the consumers’ diet. More