Ecology

Spehar, S. N. et al. Orangutans venture out of the rainforest and into the anthropocene. Sci. Adv. 4, e1701422. https://doi.org/10.1126/sciadv.1701422 (2018).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Suganuma, Y. et al. Magnetostratigraphy of the Miocene Chiang Muan Formation, northern Thailand. Implications for revised chronology of the earliest Miocene hominoid in Southeast Asia. Palaeogeogr. Palaeoclimatol. Plaeoecol. 239, 75–86 (2006).
Google Scholar 
Coster, P. et al. A complete magnetic-polarity stratigraphy of the Miocene continental deposits of Mae Moh Basin, northern Thailand, and a reassessment of the age of hominoid-bearing localities in northern Thailand. Geol. Soc. Am. Bull. 122, 1180–1191 (2010).ADS 

Google Scholar 
Begun, D. R. The Miocene hominoid radiations. In A Companion to Paleoanthropology (ed. Begun, D. R.) 398–416 (Blackwell Publishing, 2013).
Google Scholar 
Pugh, K. D. Phylogenetic analysis of Middle-Late Miocene apes. J. Hum. Evol. 165, 1–33 (2022).
Google Scholar 
Chaimanee, Y. et al. Khoratpithecus piriyai, a Late Miocene Hominoid of Thailand. Am. J. Phys. Anthropol. 131, 311–323 (2006).PubMed 

Google Scholar 
Chavasseau, O. et al. Advances in the biochronology and biostratigraphy of the continental Neogene of Myanmar. In Fossil Mammals in Asia. Neogene Biostratigraphy and Chronology (eds Wang, X. et al.) 461–474 (Columbia University Press, 2013).
Google Scholar 
Patnaik, R. Indian Neogene Siwalik Mammalian biostratigraphy. An overview. In Fossil Mammals in Asia Neogene Biostratigraphy and Chronology (eds Wang, X. et al.) 423–444 (Columbia University Press, 2013).
Google Scholar 
Chaimanee, Y. et al. A middle Miocene hominoid from Thailand and orangutan origins. Nature 422, 61–65 (2003).ADS 
CAS 
PubMed 

Google Scholar 
Chaimanee, Y. et al. A new orang-utan relative from the Late Miocene of Thailand. Nature 427, 439–441 (2004).ADS 
CAS 
PubMed 

Google Scholar 
Chaimanee, Y., Lazzari, V., Chaivanich, K. & Jaeger, J.-J. First maxilla of a late Miocene hominid from Thailand and the evolution of pongine derived characters. J. Hum. Evol. 134, 102636. https://doi.org/10.1016/j.jhevol.2019.06.007 (2019).Article 
PubMed 

Google Scholar 
Jaeger, J.-J. et al. First Hominoid from the Late Miocene of the Irrawaddy formation (Myanmar). PLoS ONE 6, 1–14 (2011).
Google Scholar 
Begun, D. R. European hominoids. In The Primate Fossil Record (ed. Hartwig, W. C.) 339–368 (Cambridge University Press, 2002).
Google Scholar 
Kelley, J. & Gao, F. Juvenile hominoid cranium from the late Miocene of southern China and hominoid diversity in Asia. Proc. Natl. Acad. Sci. U.S.A. 109, 6882–6885 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Kettle, C. J., Maycock, C. R. & Burslem, D. New directions in dipterocarp biology and conservation: A synthesis. Biotropica 44, 658–660. https://doi.org/10.1111/j.1744-7429.2012.00912.x (2012).Article 

Google Scholar 
Cannon, C. H., Morley, R. J. & Bush, A. B. G. The current refugial rainforests of Sundaland are unrepresentative of their biogeographic past and highly vulnerable to disturbance. Proc. Natl. Acad. Sci. U.S.A. 106, 11188–11193. https://doi.org/10.1073/pnas.0809865106 (2009).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Nelson, S. V. Isotopic reconstruction of habitat change surrounding the extinction of Sivapithecus, a Miocene hominoid, in the Siwalik Group of Pakistan. Palaeogeogr. Palaeoclimatol. Palaeoecol. 243, 204–222 (2007).
Google Scholar 
Bender, M. M. Variations in the 13C/12C ratios of plants in relation to the pathway of photosynthetic carbon dioxide fixation. Phytochemistry 10, 1239–1244 (1971).CAS 

Google Scholar 
Kohn, M. J. Carbon isotope compositions of terrestrial C3 plants as indicators of (paleo)ecology and (paleo)climate. Proc. Natl. Acad. Sci. 107, 19691–19695 (2010).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Bonafini, M., Pellegrini, M., Ditchfield, P. & Pollard, A. M. Investigation of the ‘canopy effect’ in the isotope ecology of temperate woodlands. J. Archaeol. Sci. 40, 3926–3935. https://doi.org/10.1016/j.jas.2013.03.028 (2013).Article 

Google Scholar 
Krigbaum, J., Berger, M. H., Daegling, D. J. & McGraw, W. S. Stable isotope canopy effects for sympatric monkeys at Tai Forest, Cote d’Ivoire. Biol. Lett. 9, 20130466. https://doi.org/10.1098/rsbl.2013.0466 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Dansgaard, W. Stable isotopes in precipitation. Tellus 16, 436–468 (1964).ADS 

Google Scholar 
Fannin, L. D. & McGraw, W. S. Does oxygen stable isotope composition in primates vary as a function of vertical stratification or folivorous behaviour?. Folia Primatol. Int. J. Primatol. 91, 219–227. https://doi.org/10.1159/000502417 (2020).Article 

Google Scholar 
Louys, J. & Roberts, P. Environmental drivers of megafauna and hominin extinction in Southeast Asia. Nature 586, 402–406. https://doi.org/10.1038/s41586-020-2810-y (2020).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Zin-Maung-Maung-Thein, et al. Stable isotope analysis of the tooth enamel of Chaingzauk mammalian fauna (late Neogene, Myanmar) and its implication to paleoenvironment and paleogeography. Palaeogeogr. Palaeoclimatol. Palaeoecol. 300, 11–22. https://doi.org/10.1016/j.palaeo.2010.11.016 (2011).Article 

Google Scholar 
Patnaik, R., Cerling, T. E., Uno, K. T. & Fleagle, J. G. Diet and habitat of Siwalik primates Indopithecus, Sivaladapis and Theropithecus. Ann. Zool. Fenn. 51, 123–142. https://doi.org/10.5735/086.051.0214 (2014).Article 

Google Scholar 
Pushkina, D., Bocherens, H., Chaimanee, Y. & Jaeger, J.-J. Stable carbon isotope reconstructions of diet and paleoenvironment from the late Middle Pleistocene Snake Cave in Northeastern Thailand. Naturwissenschaften 97, 299–309 (2010).ADS 
CAS 
PubMed 

Google Scholar 
Nelson, S. V. The paleoecology of early Pleistocene Gigantopithecus blacki inferred from isotopic analyses. Am. J. Phys. Anthropol. 155, 571–578. https://doi.org/10.1002/ajpa.22609 (2014).Article 
PubMed 

Google Scholar 
Qu, Y. et al. Preservation assessments and carbon and oxygen isotopes analysis of tooth enamel of Gigantopithecus blacki and contemporary animals from Sanhe Cave, Chongzuo, South China during the Early Pleistocene. Quat. Int. 354, 52–58. https://doi.org/10.1016/j.quaint.2013.10.053 (2014).Article 

Google Scholar 
Bocherens, H. et al. Flexibility of diet and habitat in Pleistocene South Asian mammals. Implications for the fate of the giant fossil ape Gigantopithecus. Quat. Int. 434, 148–155 (2017).
Google Scholar 
Bacon, A.-M. et al. Nam Lot (MIS 5) and Duoi U’Oi (MIS 4) Southeast Asian sites revisited. Zooarchaeological and isotopic evidences. Palaeogeogr. Palaeoclimatol. Palaeoecol. 512, 132–144. https://doi.org/10.1016/j.palaeo.2018.03.034 (2018).Article 

Google Scholar 
Jiang, Q.-Y., Zhao, L., Guo, L. & Hu, Y.-W. First direct evidence of conservative foraging ecology of early Gigantopithecus blacki (~2 Ma) in Guangxi, southern China. Am. J. Phys. Anthropol. https://doi.org/10.1002/ajpa.24300 (2021).Article 
PubMed 

Google Scholar 
Ma, J. et al. Isotopic evidence of foraging ecology of Asian elephant (Elephas maximus) in South China during the Late Pleistocene. Quat. Int. 443, 160–167. https://doi.org/10.1016/j.quaint.2016.09.043 (2017).Article 

Google Scholar 
Ma, J., Wang, Y., Jin, C., Hu, Y. & Bocherens, H. Ecological flexibility and differential survival of Pleistocene Stegodon orientalis and Elephas maximus in mainland southeast Asia revealed by stable isotope (C, O) analysis. Quat. Sci. Rev. 212, 33–44. https://doi.org/10.1016/j.quascirev.2019.03.021 (2019).ADS 
Article 

Google Scholar 
Janssen, R. et al. Tooth enamel stable isotopes of Holocene and Pleistocene fossil fauna reveal glacial and interglacial paleoenvironments of hominins in Indonesia. Quat. Sci. Rev. 144, 145–154. https://doi.org/10.1016/j.quascirev.2016.02.028 (2016).ADS 
Article 

Google Scholar 
Wang, W. et al. Sequence of mammalian fossils, including hominoid teeth, from the Bubing Basin caves, South China. J. Hum. Evol. 52, 370–379. https://doi.org/10.1016/j.jhevol.2006.10.003 (2007).Article 
PubMed 

Google Scholar 
Suraprasit, K., Bocherens, H., Chaimanee, Y., Panha, S. & Jaeger, J.-J. Late Middle Pleistocene ecology and climate in Northeastern Thailand inferred from the stable isotope analysis of Khok Sung herbivore tooth enamel and the land mammal cenogram. Quat. Sci. Rev. 193, 24–42. https://doi.org/10.1016/j.quascirev.2018.06.004 (2018).ADS 
Article 

Google Scholar 
Bocherens, H., Fizet, M. & Mariotti, A. Diet, physiology and ecology of fossil mammals as inferred from stable carbon and nitrogen biogeochemistry. Implications for Pleistocene bears. Palaeogeogr. Palaeoclimatol. Palaeoecol. 107, 213–225 (1994).
Google Scholar 
Koch, P. L., Tuross, N. & Fogel, M. L. The effects of sample treatment and diagenesis on the isotopic integrity of carbonate in biogenic hydroxylapatite. J. Archaeol. Sci. 24, 417–429 (1997).
Google Scholar 
Wright, L. E. & Schwarcz, H. P. Correspondence between stable carbon, oxygen and nitrogen isotopes in human tooth enamel and dentine. Infant diets at Kaminaljuyú. J. Archaeol. Sci. 26, 1159–1170 (1999).
Google Scholar 
Szpak, P., Metcalfe, J. Z. & Macdonald, R. A. Best practices for calibrating and reporting stable isotope measurments in archaeology. J. Archaeol. Sci. Rep. 13, 609–616 (2017).
Google Scholar 
Coplen, T. B. Guidelines and recommended terms for expression of stable-isotope-ratio and gas-ratio measurement results. Rapid Commun. Mass Spectrom. 25, 2538–2560 (2011).ADS 
CAS 
PubMed 

Google Scholar 
Bond, A. L. & Hobson, K. A. Reporting stable-isotope ratios in ecology. Recommended terminology, guidelines and best practices. Waterbirds 35, 324–331 (2012).
Google Scholar 
Craig, H. Carbon 13 in plants and the relationships between carbon 13 and carbon 14 variations in nature. J. Geol. 62, 115–149. https://doi.org/10.1086/626141 (1954).ADS 
CAS 
Article 

Google Scholar 
Cerling, T. E. & Harris, J. M. Carbon isotope fractionation between diet and bioapatite in ungulate mammals and implications for ecological and paleoecological studies. Oecologia 120, 347–363 (1999).ADS 
PubMed 

Google Scholar 
Passey, B. H. et al. Carbon isotope fractionation between diet, breath CO2, and bioapatite in different mammals. J. Archaeol. Sci. 32, 1459–1470. https://doi.org/10.1016/j.jas.2005.03.015 (2005).Article 

Google Scholar 
Howland, M. R. et al. Expression of the dietary isotope signal in the compound-specific δ13C values of pig bone lipids and amino acids. Int. J. Osteoarchaeol. 13, 54–65. https://doi.org/10.1002/oa.658 (2003).Article 

Google Scholar 
Crowley, B. E. et al. Stable carbon and nitrogen isotope enrichment in primate tissues. Oecologia 164, 611–626. https://doi.org/10.1007/s00442-010-1701-6 (2010).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Keeling, C. D. The Suess effect: 13Carbon–14Carbon interrelations. Environ. Int. 2, 229–300. https://doi.org/10.1016/0160-4120(79)90005-9 (1979).CAS 
Article 

Google Scholar 
Marino, B. D., McElroy, M. B., Salawitch, R. J. & Spaulding, W. G. Glacial-to-interglacial variations in the carbon isotopic composition of atmospheric CO2. Nature 357, 461–466. https://doi.org/10.1038/357461a0 (1992).ADS 
CAS 
Article 

Google Scholar 
Tipple, B. J., Meyers, S. R. & Pagani, M. Carbon isotope ratio of Cenozoic CO2 A comparative evaluation of available geochemical proxies. Paleoceanography https://doi.org/10.1029/2009PA001851 (2010).Article 

Google Scholar 
Zachos, J., Pagani, M., Sloan, L., Thomas, E. & Billups, K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science 292, 686–693 (2001).ADS 
CAS 
PubMed 

Google Scholar 
Cerling, T. E., Harris, J. M., Leakey, M. G., Passey, B. H. & Levin, N. E. Stable carbon and oxygen isotopes in East African Mammals. Modern and fossil. In Cenozoic Mammals of Africa (ed. Werdelin, L.) 941–952 (University of California Press, 2010).
Google Scholar 
Friedli, H., Lötscher, H., Oeschger, H., Siegenthaler, U. & Stauffer, B. Ice core record of the 13C/12C ratio of atmospheric CO2 in the past two centuries. Nature 324, 237–238. https://doi.org/10.1038/324237a0 (1986).ADS 
CAS 
Article 

Google Scholar 
Nelson, S. V. Paleoseasonality inferred from equid teeth and intra-tooth isotopic variability. Palaeogeogr. Palaeoclimatol. Palaeoecol. 222, 122–144 (2005).
Google Scholar 
Komsta, L. Processing data for outliers. R News 6, 10–13 (2006).
Google Scholar 
Hutchinson, G. E. Concluding remarks. In Cold spring Harbor Symposium on Quantitative Biology, edited by Q. Biology (1957).Hutchinson, G. E. An Introduction to Population Ecology (Yale University Press, 1978).MATH 

Google Scholar 
Baumann, C., Bocherens, H., Drucker, D. G. & Conard, N. J. Fox dietary ecology as a tracer of human impact on Pleistocene ecosystems. PLoS ONE 15, e0235692. https://doi.org/10.1371/journal.pone.0235692 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Jackson, A. L., Inger, R., Parnell, A. C. & Bearhop, S. Comparing isotopic niche widths among and within communities: SIBER—Stable Isotope Bayesian Ellipses in R. J. Anim. Ecol. 80, 595–602. https://doi.org/10.1111/j.1365-2656.2011.01806.x (2011).Article 
PubMed 

Google Scholar 
Nelson, S. V. & Hamilton, M. I. Evolution of the human dietary niche. Initial transitions. In Chimpanzees and Human Evolution (eds Muller, M. N. et al.) 286–310 (Harvard University Press, 2017).
Google Scholar 
Sun, F. et al. Paleoenvironment of the late Miocene Shuitangba hominoids from Yunnan, Southwest China: Insights from stable isotopes. Chem. Geol. 569, 120123. https://doi.org/10.1016/j.chemgeo.2021.120123 (2021).ADS 
CAS 
Article 

Google Scholar 
Nelson, S. V. Chimpanzee fauna isotopes provide new interpretations of fossil ape and hominin ecologies. Proc. R. Soc. B: Biol. Sci. 280, 20132324. https://doi.org/10.1098/rspb.2013.2324 (2013).CAS 
Article 

Google Scholar 
Merceron, G., Taylor, S., Scott, R., Chaimanee, Y. & Jaeger, J.-J. Dietary characterization of the hominoid Khoratpithecus (Miocene of Thailand). Evidence from dental topographic and microwear texture analyses. Naturwissenschaften 93, 329–333. https://doi.org/10.1007/s00114-006-0107-0 (2006).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Kay, R. F. The nut-crackers—A new theory of the adaptations of the ramapithecinae. Am. J. Phys. Anthropol. 55, 141–151 (1981).
Google Scholar 
Nelson, S. V. The Extinction of Sivapithecus. Faunal and Environmental Changes Surrounding the Disappearance of a Miocene Hominoid in the Siwaliks of Pakistan (Brill Academic Publishers, 2003).
Google Scholar 
Kanamori, T., Kuze, N., Bernard, H., Malim, T. P. & Kohshima, S. Feeding ecology of Bornean orangutans (Pongo pygmaeus morio) in Danum Valley, Sabah, Malaysia: A 3-year record including two mast fruitings. Am. J. Primatol. 72, 820–840. https://doi.org/10.1002/ajp.20848 (2010).Article 
PubMed 

Google Scholar 
Vogel, E. R. et al. Nutritional ecology of wild Bornean orangutans (Pongo pygmaeus wurmbii) in a peat swamp habitat. Effects of age, sex, and season. Am. J. Primatol. 79, 1–20. https://doi.org/10.1002/ajp.22618 (2017).Article 
PubMed 

Google Scholar 
Louys, J. et al. Sumatran orangutan diets in the Late Pleistocene as inferred from dental microwear texture analysis. Quat. Int. 603, 74–81. https://doi.org/10.1016/j.quaint.2020.08.040 (2021).Article 

Google Scholar 
Quade, J., Cerling, T. E. & Bowman, J. R. Development of Asian monsoon revealed by marked ecological shift during the latest Miocene in northern Pakistan. Nature 342, 163–166 (1989).ADS 

Google Scholar 
Hoorn, C., Ohja, T. & Quade, J. Palynological evidence for vegetation development and climatic change in the sub-Himalayan Zone (Neogene, Central Nepal). Palaeogeogr. Palaeoclimatol. Palaeoecol. 163, 133–161 (2000).
Google Scholar 
Morley, R. J. A review of the Cenozoic palaeoclimate history of Southeast Asia. In Biotic Evolution and Environmental Change in Southeast Asia (eds Gower, D. et al.) 79–114 (Cambridge University Press, 2012).
Google Scholar 
Morley, R. J. Assembly and division of the South and South-East Asian flora in relation to tectonics and climate change. J. Trop. Ecol. 34, 209–234. https://doi.org/10.1017/S0266467418000202 (2018).Article 

Google Scholar 
Sepulchre, P. et al. Mid-tertiary paleoenvironments in Thailand. Pollen evidence. Clim. Past 6, 461–473 (2010).
Google Scholar 
Sepulchre, P., Jolly, D., Ducrocq, S., Chaimanee, Y. & Jaeger, J.-J. Mid-tertiary palaeoenvironments in Thailand. Pollen evidence. Clim. Past Discuss. 5, 709–734 (2009).ADS 

Google Scholar 
Fleagle, J. G., Janson, C. H. & Reed, K. E. Primate Communities (Cambridge University Press, 1999).
Google Scholar 
Fleagle, J. G. Primate Adaptation and Evolution 3rd edn. (Elsevier, 2013).
Google Scholar 
Pilbeam, D. Gigantopithecus and the origins of Hominidae. Nature 225, 516–519. https://doi.org/10.1038/225516a0 (1970).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Jiang, Q.-Y., Zhao, L.-X. & Hu, Y.-W. Isotopic (C, O) variations of fossil enamel bioapatite caused by different preparation and measurement protocols: A case study of Gigantopithecus fauna. Vertebr. PalAsiat. 58, 159–168 (2020).
Google Scholar 
Hunt, K. D. Why are there apes? Evidence for the co-evolution of ape and monkey ecomorphology. J. Anat. 228, 630–685. https://doi.org/10.1111/joa.12454 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Zihlman, A. L., Mcfarland, R. K. & Underwood, C. E. Functional anatomy and adaptation of male gorillas (Gorilla gorilla gorilla) with comparison to male orangutans (Pongo pygmaeus). Anat. Rec. Adv. Integr. Anat. Evol. Biol. 294, 1842–1855. https://doi.org/10.1002/ar.21449 (2011).Article 

Google Scholar 
Thorpe, S. K. & Crompton, R. H. Orangutan positional behavior and the nature of arboreal locomotion in Hominoidea. Am. J. Phys. Anthropol. 131, 384–401. https://doi.org/10.1002/ajpa.20422 (2006).Article 
PubMed 

Google Scholar 
Barry, J. C. The history and chronology of Siwalik cercopithecids. J. Hum. Evol. 2, 47–58 (1987).
Google Scholar 
Jablonski, N. G., Whitfort, M. J., Roberts-Smith, N. & Qinqi, X. The influence of life history and diet on the distribution of catarrhine primates during the Pleistocene in eastern Asia. J. Hum. Evol. 39, 131–157 (2000).CAS 
PubMed 

Google Scholar 
Takai, M., Saegusa, H., Thaung-Htike, & Zin-Maung-Maung-Thein,. Neogene mammalian fauna in Myanmar. Asian Paleoprimatol. 4, 143–172 (2006).
Google Scholar 
Houle, A., Chapman, C. A. & Vickery, W. L. Intratree vertical variation of fruit density and the nature of contest competition in frugivores. Behav. Ecol. Sociobiol. 64, 429–441. https://doi.org/10.1007/s00265-009-0859-6 (2010).Article 

Google Scholar 
Vuille, M., Werner, M., Bradley, R. S. & Keimig, F. Stable isotopes in precipitation in the Asian monsoon region. J. Geophys. Res. 110, D23108 (2005).ADS 

Google Scholar  More

Rafael Almeida and his colleagues estimate that floating solar panels on 5–10% of the area of large reservoirs could help the world to reach electricity decarbonization targets by 2050 (R. M. Almeida et al. Nature 606, 246–249; 2022). On small lakes in Europe and Asia, however, the existing coverage is significantly higher (averaging 50%, according to our unpublished data), with potentially greater ecological impact (G. Exley et al. Solar Energy 219, 24–33; 2021).
Competing Interests
The authors declare no competing interests. More

Widespread occurrence and environmental stimuli responsiveness of trypsin in marine phytoplanktonTo assess whether trypsin occurs broadly in marine phytoplankton and what ecological functions phytoplankton trypsin genes may play, we investigated the occurrence of trypsin genes and environmental stimuli regulating their expression based on PhyloDB, Tara Oceans unigenes and metatranscriptomes datasets. From Tara Oceans unigenes and metatranscriptomes, trypsin homologs were found at all the sampling stations worldwide and in all major phytoplankton phyla (Fig. 1a and Supplementary Fig. 1). The broad phylogenetic representation is corroborated by the prevalence of trypsin in the individual species’ transcriptomes in the PhyloDB database (Fig. 1b), most notably in Bacillariophyta, Dinophyta, Chlorophyta, Cryptophyta and Haptophyta, the major eukaryotic groups of phytoplankton in the ocean. These indicate that trypsin is widely distributed in phytoplankton both taxonomically and geographically, a finding that advances our knowledge on the distribution of this ancient enzyme. Moreover, phylogenetic and structure alignment analysis showed that phytoplankton trypsins are more closely related with bacterial trypsins than metazoan and fungal counterparts, but contain the conserved important residues and structure typical of animal trypsins (Supplementary Figs. 2–4). These observations suggest some complex evolutionary trajectory that might result in functional innovation of phytoplankton trypsin.Fig. 1: Widespread occurrence and environmental nutrient responsiveness of trypsin in global marine phytoplankton.a Wide geographic distribution of trypsin in phytoplankton found in Tara Oceans. Color scale depicts trypsin mRNA abundance. b Wide taxonomic distribution of trypsin in algae found in PhyloDB. c Environmental nutrient drivers of phytoplankton trypsin abundance. Pairwise comparisons of environmental nutrient concentrations are shown with a color gradient denoting Pearson’s correlation coefficient. The trypsin abundance and taxonomic distribution based on the 5–180 µm size fraction from SRF layer from Tara Ocean datasets. Taxonomic trypsin abundance was related to each nutrient factor by partial (geographic distance-corrected) Mantel tests. Edge width corresponds to the Mantel’s r statistic for the corresponding distance correlations, and edge color denotes the statistical significance based on 9999 permutations. Baci Bacillariophyta, Dino Dinophyta, Chlo Chlorophyta, Cryp Cryptophyta, Hapt Haptophyta. Source data are provided as a Source Data file.Full size imageWe found a large amount of trypsin gene duplication, 5 copies to 65 copies in each algal genome we examined6. The evolution of the gene family, in gene sequence and organization relative to other functional domain, need to be treated in a separate paper6, but the rampant gene duplication suggests that trypsin may have important roles in phytoplankton. Moreover, our correlation analysis for trypsin gene expression with environmental parameters in the Tara Oceans metatranscriptomic data showed that the phytoplankton trypsin transcript abundance was correlated with environmental conditions in some taxa, size fractions, and water depths, evidence that trypsin may be important in phytoplankton to adapt to dynamical environmental conditions6. To further explore specific environmental drivers modulating the expression of trypsin, we analyzed distance-corrected dissimilarities of phytoplankton trypsin transcript abundance with environmental nutrient factors using the partial Mantel test. Analyses were restricted to the 5–20 and 20–180 µm size fractions from surface layer as their trypsin appeared to be more responsive to environmental stimuli. As shown in Fig. 1c, trypsin expression in Bacillariophyta, Dinophyta, Chlorophyta, Cryptophyta and Haptophyta was differentially correlated with nutrient availability, most notably in Bacillariophyta and Chlorophyta. Moreover, nitrate and nitrite (NO3, NO3_5m*, and NO3_NO2) and phosphate (PO4) were the strongest correlates of both Bacillariophyta and Chlorophyta trypsin transcript abundances (Fig. 1c). Hence, we posit that trypsin have important functions in the response of phytoplankton to N and P nutrient conditions.Involvement of trypsin in nitrogen and phosphorus nutrient responsesTo gain mechanistic insights into the function of trypsin in phytoplankton, we conducted experiments on the model diatom Phaeodactylum tricornutum. We identified ten trypsin genes from its genome (Supplementary Table 1), and based on qRT-PCR, we observed their growth stage- and condition-specific expression variations (Fig. 2a and Supplementary Fig. 5). Interestingly, one of these genes (PtTryp2) exhibited opposite directions of expression dynamic under N- and P-depleted conditions: downregulated under N-depleted but upregulated under P-depleted condition (Fig. 2a). Furthermore, PtTryp2 transcript increased with increasing cellular N content but decreased with increasing cellular P content (Fig. 2b, c). These results suggest that PtTryp2 is involved in an opposite-direction regulation of responses to nitrogen and phosphorus nutrient status.Fig. 2: Involvement of PtTryp2 in nitrogen and phosphorus nutrient responses.a PtTryp2 expression in P. tricornutum under different growth stages and conditions based on qRT-PCR. Nutrient-replete, HNHP; N-depletion, LNHP; P-depletion, HNLP. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. b Time-course expression patterns of PtTryp2 when P. tricornutum was grown with different forms of nitrogen nutrients. Data are presented as mean values ± SD (n = 2 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. c PtTryp2 expression pattern after phosphorus supplement. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. Source data are provided as a Source Data file.Full size imageTo interrogate the function of PtTryp2 in N and P nutrient responses, we analyzed the physiology of homologous overexpression and CRISPR/Cas9 knockout lines we generated. A PtTryp2-overexpression cell line with C-terminal eGFP fusion (named PtTryp2-OE) was generated, and the expression of OE cell line was confirmed at a protein level through Western blot (Fig. 3a). Because the function of a protein corresponds with its subcellular location, we first examined where PtTryp2 is located inside P. tricornutum cells. By computational simulation, we find PtTryp2 is potentially localized in the chloroplast via the secretory pathway (Supplementary Table 2), in accordance with the fact that chloroplasts contain a rather high number of proteases and are the main location of nutrients assimilation and remobilization7. To obtain experimental verification of the chloroplast localization, we carried out subcellular localization analysis in the OE and OEC cell lines using confocal fluorescence microscopy. Interestingly, results show PtTryp2-eGFP are localized in both the chloroplast and cytoplasmic endoplasmic reticulum (ER), to the exclusion of the nucleus and Golgi apparatus, whereas the fluorescence from the eGFP blank vector control is outspread in the cell instead of being co-localized with chloroplast and ER (Fig. 3b and Supplementary Figs. 6–8). Further analyses show that PtTryp2 lacks the C-terminal -(K/H) DEL sequences, a typical ER-retention signal that prevents ER-resident proteins from being transported to downstream locations of the secretory system8,9. Hence, PtTtryp2 is evidently transported via the ER to the chloroplast, as in the case of the previously documented light-harvesting chlorophyll a/b-binding protein in Euglena10.Fig. 3: Subcellular localization of PtTryp2.a Detection of the expression of GFP-PtTryp2 by Western blot using anti-GFP primary antibody. Left panel, GFP-PtTryp2 fusion protein. Middle panel, GFP protein. GAPDH (on the right) was detected using anti-GAPDH as the control to indicate equal protein quantities loaded to each lane. The GFP-PtTryp2 was confirmed expressed successfully at protein level in OE cell line. All experiments were repeated independently three times, and similar results were obtained. b Confocal micrographs showing subcellular localization of GFP-PtTryp2 in chloroplast (PAF, showing red autofluorescence) and endoplasmic reticulum (ER, showing blue fluorescent stain by ER-Tracker) but not in nucleus (Hoechst 33342, showing blue fluorescent stain). TL merge, merger of the fluorescence images with transmission light image. Scale bar, 10 µm, applies to all images. All experiments were repeated independently three times, and similar results were obtained. Source data are provided as a Source Data file.Full size imagePtTryp2 contains one trypsin domain and two internal repeats 1 (RPT) (Fig. 4a), offering one single target for trypsin mutagenesis. Using an optimized efficient CRISPR/Cas9 gene editing system11, we obtained three PtTryp2 mutants with different mutation characteristics in the trypsin domain (named KO1, KO2, and KO3, respectively; Fig. 4b). As shown in Fig. 4c, compared with the knockout control cell line (KOC), all three PtTryp2-KO lines exhibited a significantly diminished PtTryp2 expression under both nutrient depletion and repletion; conversely, the OE cell line displayed markedly elevated PtTryp2 expression in comparison to the overexpression control cell line (OEC). Moreover, the PtTryp2 expression level in KOC cell lines strongly responded to the ambient N and P level, but consistently showed a constant and low expression pattern in KO lines (Fig. 4d). These results verified that KO cell lines with the loss of PtTryp2 function, and OE with enhanced function of PtTryp2, can be used for subsequent functional analyses of PtTryp2.Fig. 4: Mutation generations of PtTryp2 and characters of mutants.a Schematic presentation of PtTryp2 protein. The target site (vertical arrow) for CRISPR/Cas9-based knockout is located within the conserved functional domain (green pentagon), with PAM motif shown in orange font. Red rectangle on the left depicts signal peptide; RPT: internal repeat 1; b Alignment of partial PtTryp2 sequences of the CRISPR/Cas9-generated mutants showing frameshift indels compared to wild type. The frequency by which the sequence was detected within the same colony is indicated in parenthesis. Font color coding: Black, WT sequence; Orange, functional domain containing target for CRISPR/Cas9; Purple, PAM sequence; Blue, Inserted bases; Red dashes, deleted bases. c PtTryp2 expression patterns of knockout and overexpression mutants under different conditions. FC fold change. Data are presented as mean values ± SD (n = 3 biologically independent samples). d PtTryp2 expression of knockout mutants exhibited no response to ambient N and P fluctuation. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. e Growth curves of different PtTryp2 mutants under different N and P conditions. Nutrient conditions in c–e are indicated by HNHP (Nutrient-replete), LNHP (N-depleted, P-replete), HNLP (N-replete, P-depleted), and LNLP (Nutrient-depleted). Data are presented as mean values ± SD (n = 3 biologically independent samples). Source data are provided as a Source Data file.Full size imageMoreover, we observed the growth physiology of different PtTryp2 mutants across different nutrient conditions. As shown in Fig. 4e and Supplementary Fig. 9, both of the knockout and overexpression of PtTryp2 resulted in decreases in the exponential growth rates (days 1–4) and maximum cell density across different N and P culture conditions. Taken together, these results demonstrate that both elevation and reduction of PtTryp2 expression result in cell growth repression, evidence that PtTryp2 has a crucial role in modulating cell growth in response to different N and P conditions.
PtTryp2 represses nitrogen assimilation and metabolismTranscriptomic data show that PtTryp2 knockout led to the upregulation of most of the nitrogen assimilation and metabolism genes under both N-depleted and replete conditions (Fig. 5a). The transcriptomic data are confirmed to be reproducible based on the correlation analysis of housekeeping genes (Supplementary Fig. 10 and Supplementary Table 3). Notably, the expression fold change of most N assimilation and metabolism genes under N-depleted, P-replete (LNHP) versus nutrient repete (HNHP) conditions were moderated in the PtTryp2 knockout mutant compared to that in its control (KOC), with the exception of GOGAT, which exhibited larger response to the nutrient changes in KOC (Fig. 5a). All these indicate that the inactivation of PtTryp2 enhanced N assimilation and metabolism to mitigate cell stress and reduce overall transcriptomic swing from N-depletion. Under replete conditions (HNHP), substantial transcriptional reprogramming and a significant increase in nitrate uptake rate and cellular N content was observed in the knockout mutants (KO1, KO2 and KO3) (Fig. 5b). The physiological changes were reversed in the overexpression cell lines: a decline in nitrate uptake rate and cellular N content was noted in PtTryp2-OE (Fig. 5c). All the results demonstrate that PtTryp2 functions as a repressor of nitrogen assimilation and metabolism.Fig. 5: Transcriptomic and physiological evidence that PtTryp2 directly represses nitrogen assimilation and metabolism.a PtTryp2 knockout resulted in upregulation of major nitrate-uptake and N-metabolism genes in PtTryp2 knockout (KO1) and control (KOC) under N-depleted (LNHP), P-depleted (HNLP), and nutrient-replete conditions (HNHP). NRT nitrate transporter, NR nitrate reductase, NiR nitrite reductase, GS glutamine synthetase, GOGAT glutamate synthase, GDH glutamate dehydrogenase, 2OG 2-Oxoglutarate; b NO3− uptake rate and cellular N content, increasing dramatically in PtTryp2-KO under HNHP, but decreasing remarkably under HNLP. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. c NO3− uptake rate and cellular N content, decreasing remarkably in PtTryp2-overexpressing P. tricornutum under HNHP, but increasing under HNLP. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. d Venn diagram showing the number of N-depletion induced DEGs in PtTryp2-KO1 and KOC. In parentheses, total number of DEGs; red font, upregulated; green font, downregulated. e Log2 fold changes (FC) of N-depletion induced differential gene expression in PtTryp2-KO1 against that in KOC. Most data points (93.37%) are distributed in 1,3 quadrants, indicating the same direction of change. Source data are provided as a Source Data file.Full size imageIn addition, when comparing N-depleted with N-replete conditions, 646 differentially expressed genes (DEGs) were identified in the blank vector control (KOC) but only 187 in PtTryp2-KO1, considerably fewer in the knockout mutant (Fig. 5d). Besides, the magnitude of change was smaller in PtTryp2-KO1 than in KOC for the majority (73%) of the DEGs (Fig. 5e). It is thus evident PtTryp2 in the wild type functions as an amplifier of general metabolic response to N-starvation by repressing nitrogen assimilation and metabolism. Notably, the PtTryp2-KO-promoted and PtTryp2-OE-repressed NO3− uptake patterns observed under nutrient repletion were reversed under P-depletion, indicating that PtTryp2’s roles in N and P signaling are not separated, but rather the protein might mediate the cross-talk between N and P signaling.Besides the direction of action (repression or promotion) shown above, the function of PtTryp2 involves another layer of regulation: the direction of its own expression changes. We find that PtTryp2 expression decreased under N-depletion and increased after N-supplement. Under this two-level regulatory scheme, PtTryp2 is a repressor of N uptake and assimilation genes and a promoter of N starvation-responsiveness in general metabolic pathways per se; yet its own expression decreases under N-limitation to upgrade N-uptake and assimilation under N depletion, and increases under N richness to prevent excessive N-uptake and assimilation; meanwhile, the decreased expression of PtTryp2 actually dampens the dynamic swing in the metabolic landscape in response to N-starvation. This PtTryp2-based regulatory mechanism might enable cells to swiftly respond to fluctuating N availability and cellular demand in order to finetune N responses so that N acquisition is optimized.
PtTryp2 promotes P starvation-induced genes and Pi uptakeAs shown above, PtTryp2 expression is downregulated under N-deficiency to release PtTryp2’s repressing effects on N-starvation response and to promote N uptake, thereby the cells achieve N homeostasis, and an opposite expression pattern of PtTryp2 was observed under P-deficiency, suggesting a N-P coregulation. However, the role of PtTryp2 in P-starvation responses and P homeostasis still needs to be unraveled. Toward that goal, we examined the effects of PtTryp2 inactivation on the expression changes of P starvation-induced genes and the inhibitory regulator of P signaling (SPX), which in plants is a typical P starvation response mechanism12. Consistently, most of Pi transporters (PTs) and alkaline phosphatase (APs) exhibited upregulation to P starvation response in KOC, but most of SPX genes showed downregulation (Fig. 6a).Fig. 6: Transcriptomic and physiological evidence that PtTryp2 positively modulates P starvation-induced genes during Pi starvation.a PtTryp2 knockout resulted in a reverse regulation of most P starvation-induced genes relative to that in control (KOC) under N-depleted (LNHP), P-depleted (HNLP), and nutrient-replete (HNHP) conditions. b PtTryp2 knockout caused decreases in Pi uptake and cellular P content under nutrient-replete condition (HNHP) but caused increases under N-depleted condition (LNHP). Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. c PtTryp2 knockout caused increases in Pi uptake rate and cellular P content under HNHP and LNHP. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. d Venn diagram showing the number of P-depletion induced DEGs in PtTryp2-KO1 and KOC. In parentheses, total number of DEGs; red font, upregulated; green font, downregulated. e Log2 fold changes (FC) of P-depletion induced differential gene expression in PtTryp2-KO1 against that in KOC. Most data points (95.69%) are distributed in 1,3 quadrants, indicating the same direction of change. Source data are provided as a Source Data file.Full size imageInterestingly, under P-depletion, PtTryp2 knockout downregulated the expression of most of PTs and APs, but upregulated most of the SPX genes (Fig. 6a), revealing PtTryp2’s role in WT to promote P-starvation responses. Consistent with gene transcription, PtTryp2 knockout lowered Pi uptake rate and cellular P content under the nutrient-replete condition (Fig. 6b), whereas an increase was noted in the overexpression cell line PtTryp2-OE (Fig. 6c). Based on RNA-seq, remarkably more DEGs were found for the P-depleted versus nutrient-replete comparison in PtTryp2-KO1 (1501) than that in KOC (277) (Fig. 6d). Besides, in PtTryp2-KO1, the majority of these DEGs (77.25%) exhibited greater fold changes than that in KOC (Fig. 6e). These results indicate that PtTryp2 upregulation in the wild type would dampen metabolic reprogramming in responses to P-limitation, and PtTryp2 downregulation would prevent cells from over P accumulation after P supplement, as opposed to the response to N-depletion. All these findings are indicative that PtTryp2 in the WT functions to upregulate the P starvation-induced genes and restrict general metabolic reconfiguration in response to P-limitation, a mechanism to maintain P homeostasis. Similar to that the PtTryp2-KO-promoted and PtTryp2-OE-repressed NO3− uptake patterns were reversed under P-depletion, the PtTryp2-KO-repressed Pi uptake pattern was reversed under N-depletion (Fig. 6b), implying that PtTryp2 might mediate the cross-talk between N and P signaling. The PtTryp2-OE-promoted Pi uptake pattern was not reversed under N-depletion, however, because N-depletion downregulated the expression of PtTryp2, resulting in the PtTryp2 expression pattern between OEC and OE similar to that under nutrient repletion.
PtTryp2 coordinately regulate N and P uptake and mediates N-P cross-talkGiven the PtTryp2-mediated cross-talk between N and P signaling in P. tricornutum implied in the results presented above, we were tempted to investigate the nature and the mechanism the cross-talk. Here, we uncover Pi and NO3− antagonistic interactions in P. tricornutum, which resemble that in land plants to achieve an optimal N-P nutrient balance13,14. In wild-type (WT) P. tricornutum, we observed a significant repression of NO3− uptake under P starvation and a significant repression of Pi uptake rate under N starvation. Consequently, cellular N content decreased under the P-depleted condition, and cellular P content decreased under the N-depleted condition, relative to nutrient-replete conditions (Fig. 7a, b). In accordance, the transcription of N assimilation and metabolism genes was repressed by P deficiency, and that of P starvation-induced genes was repressed by N limitation (Supplementary Fig. 11). Moreover, transcriptomic results demonstrated that PtTryp2 knockout led to the magnification of Pi and NO3− antagonistic interaction (Supplementary Fig. 11), linking PtTryp2 inactivation to the disruption of the N-P homeostasis. Taken together, our data reveal that PtTryp2’s function operates in opposite directions for N and P responses, but in a coordinated manner, consistent with a role to coregulate N and P signaling.Fig. 7: Illustration that PtTryp2 coordinately regulates N and P acquisition under fluctuating nutritional conditions.a NO3− uptake and cellular N content repressed under HNLP in wild-type cells (WT). Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. b Pi uptake and cellular P content repressed under LNHP in wild-type cells (WT). Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. c Time-course expression of PtTryp2 showed co-varied with the N/P nutrient ratio. Moreover, PtTryp2 expression fluctuated less at the N/P ratio of 16:1 compared to other N/P ratios. The 4 h after nutrient addition represents nutrient-repletion and 72 h nutrient-depletion. Data are presented as mean values ± SD (n = 3 biologically independent samples). d The cellular N/P ratio was significantly elevated by the inactivation of PtTryp2. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. e The cellular N/P ratio was significantly decreased by the overexpression of PtTryp2. Data are presented as mean values ± SD (n = 3 biologically independent samples). The comparisons between the averages of the two groups were evaluated using the one-tailed Student’s t test. The p values with significance (p ≤ 0.05) are shown. f Hypothetical model depicting the role of PtTryp2 in balancing N and P acquisition. Under N-depletion, PtTryp2 expression is downregulated to promote N-starvation responses and repress P-starvation responses. In contrast, under P-depletion, PtTryp2 expression is upregulated to reinforce P-starvation responses and lessen N-starvation response. By this feedback loop, optimal N-P uptake is achieved to maintain stoichiometric homeostasis. Upregulated genes and enhancement processes are shown in red, downregulated genes and weakened processes colored green. The black arrows depict transcriptional activation. Black bar at line’s end depicts inhibitory regulation. The gray arrows depict possible but unverified interaction between PtTryp2 and the existing P regulating cascade SPX-PHR or an equivalent of the N regulating cascade known in plants (SPX-NLP where NLP stands for NIN-like protein, a transcription factor). Source data are provided as a Source Data file.Full size imageTo further illustrate this, we have carried out PtTryp2 expression pattern analysis across different N/P nutrient stoichiometric ratio conditions, and found that PtTryp2 expression co-varied with the N/P nutrient ratio (Fig. 7c). The time-course analysis showed that PtTryp2 expression fluctuated less under different N or P conditions at the N/P ratio of 16:1 compared to other N/P ratios. The N/P nutrient ratio of 16:1 is considered balanced stoichiometry (Redfield ratio) and appears to be optimal for P. tricornutum growth (Supplementary Fig. 12), as previously documented15, suggesting that at this nutrient stoichiometry there is no need for a significant change in PtTryp2 expression to maintain N/P balance, but other N:P nutrient ratios deviating from 16:1 caused changes in PtTryp2 expression to maintain N/P balance. Moreover, the extent of change in PtTryp2 expression varied between cultures with different levels of N:P nutrient ratios, and between 4 and 72 h after culture inoculation from N- and P-depletion-acclimated parent culture into the experimental nutrient conditions. At 72 h PtTryp2 expression level increased with the degree of P stress (the higher the N:P ratio, the more P stressed the cultures were), except for the N:P = 1:1 condition, an extreme N-limited condition that seemed to cause PtTryp2 expression not to respond according to the general trend. Overall, all these data indicate that PtTryp2 responds strongly to the variability of the N:P ratio. Correspondingly, the cellular N/P ratio under nutrient-repletion also seems to be influenced by PtTryp2 expression level: the cellular N/P ratio was significantly elevated by PtTryp2 knockout, but conversely, was significantly decreased by the overexpression of PtTryp2 (Fig. 7d). Evidently, PtTryp2 serves to coordinate N and P uptake and metabolism to dampen the amplitude of N:P ratio changes that occur when the P. tricornutum cells experience fluctuations in nutrient conditions16,17. That is, PtTryp2 in P. tricornutum acts like an amplitude reducer of the N-P seesaw to achieve the N and P stoichiometric homeostasis (Fig. 7f).As critical nutrients for phytoplankton and plants, the balance and homeostasis of N and P are crucial to the growth of the organisms. For plants, nutrient supply in the soil is highly variable; therefore, to achieve optimal and coordinated utilization of N and P, integration of N and P signaling into an integrated network is required18. Recent studies have revealed the critical components of the network in the model plants Arabidopsis thaliana and Oryza sativa12,19,20,21. Similarly, phytoplankton in the ocean face remarkable environmental nutrient variations, and N and P nutrients are often limited22,23. Although the respective responses to N and P deficiencies have been extensively studied in phytoplankton24,25, an integrative signaling pathway of N-P nutrition cross-talk has remained unknown until now. It is striking to find that trypsin, rather than homologs of plant NRT1.1 and NIGT114,19, mediates and regulates the nitrate-phosphate signaling cross-talk.The two-level model of PtTryp2 function (Fig. 7f), including the direction of PtTryp2 action and the direction of PtTryp2 expression changes, demonstrate that PtTryp2 functions by shifting the setpoints, by tuning its own expression level, at which N signaling or P signaling is triggered in response to environmental nutrient fluctuations so that cells commit to appropriate responses. However, much of the mechanics in the regulatory cascade, from environmental nutrient sensing, PtTryp2-mediated signaling, to the regulation of the effectors such as N- and Pi-transporters and assimilatory genes, remains to be elucidated. Although the interplay between N and P nutrition based on SPX-NLP-NIGT1 and SPX-PHR-NIGT1 cascades, respectively have been uncovered in plants12,19, how PtTryp2 interacts with the SPX-PHR cascade26 and whether a SPX-NLP cascade or other regulatory cascades exist and interact with PtTryp2 for P and N nutrient regulation in phytoplankton remain to be addressed.As an initial attempt, we have performed transcriptional regulatory interaction analysis based on the Inferelator algorithm27 to predict the potential co-regulated genes in the PtTryp2-dependent regulatory cascade. Consequently, a set of 1034 genes co-regulated with PtTryp2 were identified, including 10 transcription factors (Supplementary Table 4), 10 N metabolism and assimilation genes, and a P responsive gene (Supplementary Fig. 13). Moreover, the functional enrichment of the gene set showed that PtTryp2 is possibly involved in post-transcriptional regulation, intracellular signal transduction pathway and kinase-based phosphorus metabolism and recycle pathway (Supplementary Fig. 14). The results hint on a potentially complex regulatory network that requires much more transcriptomes derived from more growth conditions than just the N and P conditions used in this study and other experimental approaches to unravel.We used the potential co-regulated gene list identified in this study in a comparative analysis with the published co-regulatory analysis datasets that contained hundreds of public RNA-seq datasets: DiatomPortal28 and PhaeoNet29. Interestingly, based on the DiatomPortal dataset, the PtTryp2 was found in the Phatr_hclust_0381 hierarchical cluster that consists of 10 genes, which has been identified as the GO term of ubiquitin-dependent protein catabolism. In terrestrial plants, the ubiquitination and degradation of SPX4 was found to mediate the nitrate-phosphate interaction signaling pathway by enabling the release of PHR2 and NLP3 into the nucleus to activate the expression of both phosphate- and nitrate-responsive genes12,19. In addition, we found 120 genes that were common in our gene list and PhaeoNet, some of which are transcription factors.Taken together, our analyses showed that the deletion and overexpression of PtTryp2 simultaneously impacted nitrogen and phosphorus uptake, nitrogen and phosphorus contents of the cell, and the N:P ratio. The simultaneous impact on N and P in opposite directions suggests that this protein either directly regulates the N and P uptake machinery or is close to the direct regulator, e.g., functioning through the ubiquitination and degradation of the direct regulators as in terrestrial plants. Furthermore, it is conceivable that one or more intermediate relays between PtTryp2 and the direct regulator would make it extremely challenging, if not impossible, to exert such precise and coordinated bidirectional regulation on N and P. To understand the mechanics of the regulatory mechanism, co-immunoprecipitation and Chromation immunoprecipitation sequencing are underway in our laboratory to experimentally identify the potential proteins and DNAs interacting with PtTryp2. Further studies on multiple fronts surrounding trypsin and its regulatory pathway are required for gaining an in-depth understanding of the interplay between N and P nutrition in phytoplankton and how phytoplankton will adapt to the potentially more variable and skewed N-P environment in the Anthropocene oceans. More

Brazil’s government risks fuelling the climate and biodiversity crisis by offsetting the fertilizer shortage resulting from Russia’s invasion of Ukraine this year (J. Liu et al. Nature 604, 425 (2022); S. Osendarp et al. Nature 604, 620–624; 2022). To produce an alternative fertilizer, it plans to mine up to 12 million tonnes annually of rhodoliths taken from an area in the South Atlantic that is roughly the size of the United Kingdom (see go.nature.com/3yhiyio).A full list of co-signatories to this letter appears in Supplementary Information.
Competing Interests
The author declares no competing interests. More

Faegri, K. & van der Pijl, L. The Principles of Pollination Ecology 2nd edn. (Pergamon Press, 1971).
Google Scholar 
Crepet, W. L. The role of insect pollination in the evolution of the angiosperms. In Pollination Biology (ed. Real, L.) 29–50 (Academic Press, 1983).Chapter 

Google Scholar 
Dressler, R. L. Biology of the orchid bees (Euglossini). Annu. Rev. Ecol. Syst. 13, 373–394 (1982).Article 

Google Scholar 
Ollerton, J. et al. A global test of the pollination syndrome hypothesis. Ann. Bot. 103, 1471–1480 (2009).PubMed 
PubMed Central 
Article 

Google Scholar 
IPBES. The assessment report of the Intergovernmental Science-Policy Platform on Biodiversity and Ecosystem Services on pollinators, pollination and food production. (IPBES Secretariat, 2016).Kearns, C. A., Inouye, D. W. & Waser, N. M. Endangered mutualisms : The conservation of plant-pollinator interactions. Annu. Rev. Ecol. Syst. 29, 83–112 (1998).Article 

Google Scholar 
Ollerton, J. Pollinator diversity: Distribution, ecological function, and conservation. Annu. Rev. Ecol. Evol. Syst. 48, 353–376 (2017).Article 

Google Scholar 
Irwin, R. E., Bronstein, J. L., Manson, J. S. & Richardson, L. Nectar robbing: Ecological and evolutionary perspectives. Annu. Rev. Ecol. Evol. Syst. 41, 271–292 (2010).Article 

Google Scholar 
Hargreaves, A. L., Harder, L. D. & Johnson, S. D. Consumptive emasculation: The ecological and evolutionary consequences of pollen theft. Biol. Rev. 84, 259–276 (2009).PubMed 
Article 

Google Scholar 
McCall, A. C. & Irwin, R. E. Florivory: The intersection of pollination and herbivory. Ecol. Lett. 9, 1351–1365 (2006).PubMed 
Article 

Google Scholar 
King, C., Ballantyne, G. & Willmer, P. G. Why flower visitation is a poor proxy for pollination: Measuring single-visit pollen deposition, with implications for pollination networks and conservation. Methods Ecol. Evol. 4, 811–818 (2013).Article 

Google Scholar 
Komamura, R., Koyama, K., Yamauchi, T., Konno, Y. & Gu, L. Pollination contribution differs among insects visiting cardiocrinum cordatum flowers. Forests 12, 452 (2021).Article 

Google Scholar 
Aizen, M. A. Flower sex ratio, pollinator abundance, and the seasonal pollination dynamics of a protandrous plant. Ecology 82, 127–144 (2001).Article 

Google Scholar 
Forrest, J. R. K., Ogilvie, J. E., Gorischek, A. M. & Thomson, J. D. Seasonal change in a pollinator community and the maintenance of style length variation in Mertensia fusiformis (Boraginaceae). Ann. Bot. 108, 1–12 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
Burrill, R. M. & Dietz, A. The response of honey bees to variations in solar radiation and temperature. Apidologie 12, 319–328 (1981).Article 

Google Scholar 
Corbet, S. A. et al. Temperature and the pollinating activity of social bees. Ecol. Entomol. 18, 17–30 (1993).Article 

Google Scholar 
Herrera, C. M. Daily patterns of pollinator activity, differential pollinating effectiveness, and floral resource availability, in a summer-flowering Mediterranean shrub. Oikos 58, 277–288 (1990).Article 

Google Scholar 
Zoller, L., Bennett, J. M. & Knight, T. M. Diel-scale temporal dynamics in the abundance and composition of pollinators in the Arctic summer. Sci. Rep. 10, 21187 (2020).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kirk, W. D. J., Ali, M. & Breadmore, K. N. The effects of pollen beetles on the foraging behaviour of honey bees. J. Apic. Res. 34, 15–22 (1995).Article 

Google Scholar 
Tan, K. et al. Fearful foragers: Honey bees tune colony and individual foraging to multi-predator presence and food quality. PLoS ONE 8, e75841 (2013).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Baldock, K. C. R., Memmott, J., Carlos Ruiz-Guajardo, J., Roze, D. & Stone, G. N. Daily temporal structure in African savanna flower visitation networks and consequences for network sampling. Ecology 92, 687–698 (2011).PubMed 
Article 

Google Scholar 
Edwards, J., Smith, G. P. & McEntee, M. H. F. Long-term time-lapse video provides near complete records of floral visitation. J. Pollinat. Ecol. 16, 91–100 (2015).CAS 
Article 

Google Scholar 
Knop, E. et al. Rush hours in flower visitors over a day–night cycle. Insect Conserv. Divers. 11, 267–275 (2018).Article 

Google Scholar 
Herrera, J. Pollination relationships in southern Spanish Mediterranean shrublands. J. Ecol. 76, 274–287 (1988).Article 

Google Scholar 
Eckhart, V. M. Spatio-temporal variation in abundance and variation in foraging behavior of the pollinators of gynodioecious Phacelia linearis (Hydrophyllaceae). Oikos 64, 573–586 (1992).Article 

Google Scholar 
Nakano, C. & Washitani, I. Variability and specialization of plant-pollinator systems in a northern maritime grassland. Ecol. Res. 18, 221–246 (2003).Article 

Google Scholar 
Moeller, D. A. Pollinator community structure and sources of spatial variation in plant-pollinator interactions in Clarkia xantiana ssp. xantiana. Oecologia 142, 28–37 (2005).ADS 
PubMed 
Article 

Google Scholar 
Lortie, C. J., Budden, A. E. & Reid, A. M. From birds to bees: Applying video observation techniques to invertebrate pollinators. J. Pollinat. Ecol. 6, 125–128 (2012).
Google Scholar 
Pegoraro, L., Hidalgo, O., Leitch, I. J., Pellicer, J. & Barlow, S. E. Automated video monitoring of insect pollinators in the field. Emerg. Top. Life Sci. 4, 87–97 (2020).PubMed 
Article 

Google Scholar 
Steen, R. & Aase, A. L. T. O. Portobale digital video surveillance system for monitoring flower-visiting bumblebees. J. Pollinat. Ecol. 5, 90–94 (2011).Article 

Google Scholar 
Sakamoto, R. L., Morinaga, S. I., Ito, M. & Kawakubo, N. Fine-scale flower-visiting behavior revealed by using a high-speed camera. Behav. Ecol. Sociobiol. 66, 669–674 (2012).Article 

Google Scholar 
Georgian, E., Fang, Z., Emshwiller, E. & Pidgeon, A. The pollination ecology of Rhododendron floccigerum Franchet (Ericaceae) in Weixi, Yunnan Province, China. J. Pollinat. Ecol. 16, 72–81 (2015).Article 

Google Scholar 
Suetsugu, K., Nakahama, N., Ito, A. & Isagi, Y. Time-lapse photography reveals the occurrence of unexpected bee-pollination in Calanthe izuinsularis, an endangered orchid endemic to the Izu archipelago. J. Nat. Hist. 51, 783–792 (2017).Article 

Google Scholar 
Steen, R. Diel activity, frequency and visit duration of pollinators in focal plants: in situ automatic camera monitoring and data processing. Methods Ecol. Evol. 8, 203–213 (2017).Article 

Google Scholar 
Burton, A. C. et al. Wildlife camera trapping: A review and recommendations for linking surveys to ecological processes. J. Appl. Ecol. 52, 675–685 (2015).Article 

Google Scholar 
Edwards, J., Griffin, A. J. & Knoedler, M. R. Simultaneous recordings of insect visitors to flowers show spatial and temporal heterogeneity. Ann. Entomol. Soc. Am. 112, 93–98 (2019).Article 

Google Scholar 
Droissart, V. et al. PICT: A low-cost, modular, open-source camera trap system to study plant–insect interactions. Methods Ecol. Evol. 12, 1389–1396 (2021).Article 

Google Scholar 
Li, Y. et al. Paleobiogeography of the lotus plant (Nelumbonaceae: Nelumbo) and its bearing on the paleoclimatic changes. Palaeogeogr. Palaeoclimatol. Palaeoecol. 399, 284–293 (2014).Article 

Google Scholar 
Li, J. & Huang, S. Effective pollinators of Asian sacred lotus (Nelumbo nucifera): Contemporary pollinators may not reflect the historical pollination syndrome. Ann. Bot. 104, 845–851 (2009).PubMed 
PubMed Central 
Article 

Google Scholar 
Lee, T. D. Patterns of fruit and seed production. In Plant Reproductive Ecology: Patterns and Strategies (eds Doust, J. L. & Doust, L. L.) 179–202 (Oxford University Press, 1988).
Google Scholar 
Moro, C. F., Yonekura, M., Kouzuma, Y., Agrawal, G. K. & Rakwal, R. Lotus—a source of food and medicine: Current status and future perspectives in context of the seed proteomics. Int. J. Life Sci. 7, 1–5 (2013).CAS 
Article 

Google Scholar 
Zhu, M., Liu, T. & Guo, M. Current advances in the metabolomics study on lotus seeds. Front. Plant Sci. 7, 891 (2017).
Google Scholar 
Guo, H. B. Cultivation of lotus (Nelumbo nucifera Gaertn. Ssp. nucifera) and its utilization in China. Genet. Resour. Crop Evol. 56, 323–330 (2009).Article 

Google Scholar 
Vogel, S. & Hadacek, F. Contributions to the functional anatomy and biology of Nelumbo nucifera (Nelumbonaceae) III. An ecological reappraisal of floral organs. Plant Syst. Evol. 249, 173–189 (2004).Article 

Google Scholar 
Seymour, R. S. & Schultze-Motel, P. Thermoregulating lotus flowers. Nature 383, 305 (1996).ADS 
CAS 
Article 

Google Scholar 
Dieringer, G., Leticia Cabrera, R. & Mottaleb, M. Ecological relationship between floral thermogenesis and pollination in Nelumbo lutea (Nelumbonaceae). Am. J. Bot. 101, 357–364 (2014).PubMed 
Article 

Google Scholar 
Johnson, J. B. & Omland, K. S. Model selection in ecology and evolution. Trends Ecol. Evol. 19, 101–108 (2004).PubMed 
Article 

Google Scholar 
R Development Core Team. R: A language and environment for statistical computing. (R Foundation for Statistical Computing, 2020).Inouye, D. W., Gill, D. E., Dudash, M. R. & Fenster, C. B. A model and lexicon for pollen fate. Am. J. Bot. 81, 1517–1530 (1994).Article 

Google Scholar 
Ne’Eman, G., Jürgens, A., Newstrom-Lloyd, L., Potts, S. G. & Dafni, A. A framework for comparing pollinator performance: Effectiveness and efficiency. Biol. Rev. 85, 435–451 (2010).PubMed 

Google Scholar 
Kearns, C. A. & Inouye, D. W. Techniques for Pollination Biologists (University Press of Colorado, 1993).
Google Scholar 
Delaplane, K. S., Mayer, D. R. & Mayer, D. F. Crop Pollination by Bees (CABI Publishing, 2000).
Google Scholar 
Oldroyd, B. P. & Nanork, P. Conservation of Asian honey bees. Apidologie 40, 296–312 (2009).Article 

Google Scholar 
Abou-Shaara, H. F. The foraging behaviour of honey bees Apis mellifera: A review. Vet. Med. (Praha) 59, 1–10 (2014).Article 

Google Scholar 
Reyes-Carrillo, J. L., Eischen, F. A., Cano-Rios, P., Rodríguez Martínez, R. & Nava Camberos, U. Pollen collection and honey bee forager distribution in cantaloupe. Acta Zool. Mex. 23, 29–36 (2007).
Google Scholar 
Bloch, G., Bar-Shai, N., Cytter, Y. & Green, R. Time is honey: Circadian clocks of bees and flowers and how their interactions may influence ecological communities. Philos. Trans. R. Soc. B Biol. Sci. 372, 20160256 (2017).Article 
CAS 

Google Scholar 
Lawson, D. A. & Rands, S. A. The effects of rainfall on plant–pollinator interactions. Arthropod. Plant. Interact. 13, 561–569 (2019).Article 

Google Scholar 
Antiqueira, P. A. P. et al. Precipitation and predation risk alter the diversity and behavior of pollinators and reduce plant fitness. Oecologia 192, 745–753 (2020).ADS 
PubMed 
Article 

Google Scholar 
Matsuura, M. Comparative biology of the five Japanese species of the genus Vespa (Hymenoptera, Vespidae). Bull. Fac. Agric. Mie Univ. 69, 1–131 (1984).
Google Scholar 
Dukas, R. Effects of perceived danger on flower choice by bees. Ecol. Lett. 4, 327–333 (2001).Article 

Google Scholar 
Romero, G. Q., Antiqueira, P. A. P. & Koricheva, J. A meta-analysis of predation risk effects on pollinator behaviour. PLoS ONE 6, e20689 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kluser, S. & Peduzzi, P. Global pollinator decline : A literature review. UNEP/GRID-Europe (2007).Potts, S. G. et al. Declines of managed honey bees and beekeepers in Europe. J. Apic. Res. 49, 15–22 (2010).Article 

Google Scholar 
Paudel, Y. P., Mackereth, R., Hanley, R. & Qin, W. Honey bees (Apis mellifera L.) and pollination issues: current status, impacts and potential drivers of decline. J. Agric. Sci. 7, 93–109 (2015).
Google Scholar 
Theisen-Jones, H. & Bienefeld, K. The Asian honey bee (Apis cerana) is significantly in decline. Bee World 93, 90–97 (2016).Article 

Google Scholar 
Cameron, S. A. et al. Patterns of widespread decline in North American bumble bees. Proc. Natl. Acad. Sci. U.S.A. 108, 662–667 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Wood, T. J. et al. Managed honey bees as a radar for wild bee decline?. Apidologie 51, 1100–1116 (2020).Article 

Google Scholar 
Thapa, R. Honeybees and other insect pollinators of cultivated plants: A review. J. Inst. Agric. Anim. Sci. 27, 1–23 (2006).Article 

Google Scholar 
Nicholls, C. I. & Altieri, M. A. Plant biodiversity enhances bees and other insect pollinators in agroecosystems. A review. Agron. Sustain. Dev. 33, 257–274 (2013).Article 

Google Scholar 
Rader, R., Howlett, B. G., Cunningham, S. A., Westcott, D. A. & Edwards, W. Spatial and temporal variation in pollinator effectiveness: Do unmanaged insects provide consistent pollination services to mass flowering crops?. J. Appl. Ecol. 49, 126–134 (2012).Article 

Google Scholar 
Nikkeshi, A., Inoue, H., Arai, T., Kishi, S. & Kamo, T. The bumblebee Bombus ardens ardens (Hymenoptera: Apidae) is the most important pollinator of Oriental persimmon, Diospyros kaki (Ericales: Ebenaceae), in Hiroshima, Japan. Appl. Entomol. Zool. 54, 409–419 (2019).Article 

Google Scholar 
Macgregor, C. J. & Scott-Brown, A. S. Nocturnal pollination: an overlooked ecosystem service vulnerable to environmental change. Emerg. Top. Life Sci. 4, 19–32 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Knop, E. et al. Artificial light at night as a new threat to pollination. Nature 548, 206–209 (2017).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Gilpin, A. M., Denham, A. J. & Ayre, D. J. The use of digital video recorders in pollination biology. Ecol. Entomol. 42, 383–388 (2017).Article 

Google Scholar 
Barlow, S. E. & O’Neill, M. A. Technological advances in field studies of pollinator ecology and the future of e-ecology. Curr. Opin. Insect Sci. 38, 15–25 (2020).PubMed 
Article 

Google Scholar 
Suzuki-Ohno, Y. et al. Deep learning increases the availability of organism photographs taken by citizens in citizen science programs. Sci. Rep. 12, 1–10 (2022).Article 
CAS 

Google Scholar 
Høye, T. T. et al. Deep learning and computer vision will transform entomology. Proc. Natl. Acad. Sci. U.S.A. 118, 1–10 (2021).
Google Scholar  More

Mal, T. K. & Lovett-Doust, J. Phenotypic plasticity in vegetative and reproductive traits in an invasive weed, Lythrum salicaria (Lythraceae), in response to soil moisture. Am. J. Bot. 92, 819–825 (2005).PubMed 
Article 

Google Scholar 
Wang, S., Li, L. & Zhou, D.-W. Morphological plasticity in response to population density varies with soil conditions and growth stage in Abutilon theophrasti (Malvaceae). Plant Ecol. 218, 785–797 (2017).Article 

Google Scholar 
Eid, E. M., Shaltout, K. H., Al-Sodany, Y. M., Haroun, S. A. & Jensen, K. A comparison of the functional traits of Phragmites australis in Lake Burullus (a Ramsar site in Egypt): Young vs. old populations over the nutrient availability gradient. Ecol. Eng. 166, 106244 (2021).Article 

Google Scholar 
Hassan, M. O. et al. Habitat variations affect morphological, reproductive and some metabolic traits of Mediterranean Centaurea glomerata Vahl populations. Heliyon 6, e04173 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Więcław, H. Within-species variation among populations of the Carex flava complex as a function of habitat conditions. Plant Ecol. Divers. 10, 443–451 (2017).Article 

Google Scholar 
MacLeod, N. Phylogenetic signals in morphometric data. In Morphology, Shape and Phylogeny (eds MacLeod, N. & Forey, P.) 100–138 (Taylor & Francis, Routledge, 2002).Chapter 

Google Scholar 
Gratani, L. Plant phenotypic plasticity in response to environmental factors. Adv. Bot. 2014, 208747 (2014).
Google Scholar 
Koopman, J. et al. Global distribution of Carex buekii (Cyperaceae) reappraised. Phytotaxa 358, 139–161 (2018).Article 

Google Scholar 
Egorova, T. V. The Sedges (Carex L.) of Russia and Adjacent States (Within the Limits of the Former USSR) (St.-Petersburg State Chemical-Pharmaceutical Academy, St.-Petersburg, 1999).
Google Scholar 
Burkart, M. River corridor plants (Stromtalpflanzen) in Central European lowland: A review of a poorly understood plant distribution pattern: River corridor plants. Glob. Ecol. Biogeogr. 10, 449–468 (2001).Article 

Google Scholar 
Więcław, H. et al. Ecology, threats and conservation status of Carex buekii (Cyperaceae) in Central Europe. Sci. Rep. 9, 11162 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Nobis, A. & Skórka, P. River corridor plants revisited: What drives their unique distribution patterns?. Plant Biosyst. 150, 244–253 (2016).Article 

Google Scholar 
Spink, A., Sparks, R. E., Van Oorschot, M. & Verhoeven, J. T. A. Nutrient dynamics of large river floodplains. Regul. Rivers. Res. Manag. 14, 203–216 (1998).Article 

Google Scholar 
Myśliwy, M. Diversity and environmental variability of riparian tall herb fringe communities of the order Convolvuletalia sepium in Polish river valleys. Monographiae Botaniceae 108, 1–129 (2019).Article 

Google Scholar 
Fischer, W. Die Stromtalpflanzen Brandenburgs. Unter Havel. Naturkundliche Berichte 5, 4–13 (1996).
Google Scholar 
Thiers, B. Index Herbariorum: A global directory of public herbaria and associated staff. New York Botanical Garden’s Virtual Herbarium. Available from: http://sweetgum.nybg.org/science/ih/ (assessed: 08 March 2022).Methods of Soil Analysis: Part 3 Chemical Methods. (Soil Science Society of America, American Society of Agronomy, 1996). https://doi.org/10.2136/sssabookser5.3.StatSoft Inc. Electronic Statistics Textbook. (Tulsa, OK: StatSoft, 2013).ter Braak, C. J. F. & Smilauer, P. CANOCO Reference Manual and CanoDraw for Windows User’s Guide: Software for Canonical Community Ordination (version 4.5). (Ithaca NY, 2002).Xu, X. et al. Effects of potassium levels on plant growth, accumulation and distribution of carbon, and nitrate metabolism in apple dwarf rootstock seedlings. Front. Plant Sci. 11, 904 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Sardans, J. & Peñuelas, J. Potassium control of plant functions: Ecological and agricultural implications. Plants 10, 419 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Broadley, M. R. et al. Phylogenetic variation in the shoot mineral concentration of angiosperms. J. Exp. Bot. 55, 321–336 (2004).CAS 
PubMed 
Article 

Google Scholar 
Watanabe, T. et al. Evolutionary control of leaf element composition in plants. New Phytol. 174, 516–523 (2007).CAS 
PubMed 
Article 

Google Scholar 
Konings, H., Koot, E. & Wolf, A. T. Growth characteristics, nutrient allocation and photosynthesis of Carex species from floating fens. Oecologia 80, 111–121 (1989).ADS 
CAS 
PubMed 
Article 

Google Scholar 
Busch, J. Characteristic values of key ecophysiological parameters in the genus Carex. Flora 196, 405–430 (2001).Article 

Google Scholar 
Zhang, D. et al. Effect of hydrological fluctuation on nutrient stoichiometry and trade-offs of Carex schmidtii. Ecol. Ind. 120, 106924 (2021).Article 
CAS 

Google Scholar 
Zhang, D. et al. Growth and physiological responses of Carex schmidtii to water-level fluctuation. Hydrobiologia 847(3), 967–981 (2020).CAS 
Article 

Google Scholar 
Yan, H. et al. Growth and physiological responses to water depths in Carex schmidtii Meinsh. PLoS ONE 10(5), e0128176 (2015).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Luo, W. & Xie, Y. Growth and morphological responses to water level and nutrient supply in three emergent macrophyte species. Hydrobiologia 624(1), 151–160 (2009).CAS 
Article 

Google Scholar 
Lu, Y. Growth and morphological responses to water level variations in two Carex species from Sanjiang Plain, China. Afr. J. Agric. Res. 6, 28–34 (2011).
Google Scholar 
Cao, Y. et al. Flooding influences on the C, N and P stoichiometry in terrestrial ecosystems: A meta-analysis. CATENA 215, 106287 (2022).CAS 
Article 

Google Scholar 
Sardans, T., Peñuelas, T., Prieto, P. & Estiarte, M. Drought and warming induced changes in P and K concentration and accumulation in plant biomass and soil in a Mediterranean shrubland. Plant Soil 306, 261–271 (2007).Article 
CAS 

Google Scholar 
Flórez-Flórez, C. P., León-Peláez, J. D., Osorio-Vega, N. W. & Restrepo-Llano, M. F. Nutrient dynamics in forest plantations of Azadirachta indica (Meliaceae) established for restoration of degraded lands in Colombia. Rev. Biol. Trop. 61, 515–529 (2013).PubMed 
Article 

Google Scholar 
Jordan-Meille, L. & Pellerin, S. Leaf area establishment of a maize (Zea mays L.) field crop under potassium deficiency. Plant Soil 265, 75–92 (2004).CAS 
Article 

Google Scholar 
Gerardeaux, E., Jordan-Meille, L., Constantin, J., Pellerin, S. & Dingkuhn, M. Changes in plant morphology and dry matter partitioning caused by potassium deficiency in Gossypium hirsutum L. Environ. Exp. Bot. 67, 451–459 (2010).CAS 
Article 

Google Scholar 
Bailey, J. S. & Laidlaw, A. S. Growth and development of white clover (Trifolium repens L.) as influenced by P and K nutrition. Ann. Bot. 81, 783–786 (1998).Article 

Google Scholar 
White, P. Relationship between the development and growth of rye (Secale cereale L.) and the potassium concentration in solution. Ann. Bot. 72, 349–358 (1993).CAS 
Article 

Google Scholar 
Pujos, A. & Morard, P. Effects of potassium deficiency on tomato growth and mineral nutrition at the early production stage. Plant Soil 189, 189–196 (1997).CAS 
Article 

Google Scholar 
Osakabe, Y. et al. Osmotic stress responses and plant growth controlled by potassium transporters in Arabidopsis. Plant Cell 25, 609–624 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Lebaudy, A. et al. Plant adaptation to fluctuating environment and biomass production are strongly dependent on guard cell potassium channels. Proc. Natl. Acad. Sci. 105, 5271–5276 (2008).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Tränkner, M., Tavakol, E. & Jákli, B. Functioning of potassium and magnesium in photosynthesis, photosynthate translocation and photoprotection. Physiol. Plantarum 163, 414–431 (2018).Article 
CAS 

Google Scholar 
Du, Q. et al. Effect of potassium deficiency on root growth and nutrient uptake in maize (Zea mays L.). Agric. Sci. 8, 1263–1277 (2017).CAS 

Google Scholar 
Hu, W., Coomer, T. D., Loka, D. A., Oosterhuis, D. M. & Zhou, Z. Potassium deficiency affects the carbon-nitrogen balance in cotton leaves. Plant Physiol. Biochem. 115, 408–417 (2017).CAS 
PubMed 
Article 

Google Scholar 
Reisch, C., Meier, S., Schmid, C. & Bartelheimer, M. Clonal diversity and genetic variation of the sedge Carex nigra in an alpine fen depend on soil nutrients. PeerJ 8, e8887 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Lenssen, J. P. M., Menting, F. B. J. & Van der Putten, W. H. Plant responses to simultaneous stress of waterlogging and shade: Amplified or hierarchical effects?. New Phytol. 157, 281–290 (2003).CAS 
PubMed 
Article 

Google Scholar 
Liu, Z. G. & Li, Z. Q. Effects of different grazing regimes on the morphological traits of Carex duriuscula on the Inner Mongolia steppe, China. N. Z. J. Agric. Res. 53(1), 5–12 (2010).Article 

Google Scholar 
Więcław, H. et al. Morphological variability and genetic diversity in Carex buxbaumii and Carex hartmaniorum (Cyperaceae) populations. PeerJ 9, e11372 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
Więcław, H., Kurnicki, B., Bihun, M., Białecka, B. & Koopman, J. Carex section Racemosae (Cyperaceae) in Europe: Morphological diversity, taxonomy and phylogenetic relationships. Bot. J. Linn. Soc. 183, 124–145 (2017).
Google Scholar 
Jiménez-Mejías, P., Benítez-Benítez, C., Fernández-Mazuecos, M. & Martín-Bravo, S. Cut from the same cloth: The convergent evolution of dwarf morphotypes of the Carex flava group (Cyperaceae) in Circum-Mediterranean mountains. PLoS ONE 12(12), e0189769 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Więcław, H. Carex flava agg. (section Ceratocystis, Cyperaceae) in Poland: taxonomy, morphological variation, and soil conditions. Biodivers. Res. Conserv. 33, 3–51 (2014).Article 

Google Scholar 
Kalela, A. Systematische und Pflanzengeographische Studien an der Carex-Subsektion Alpinae Kalela. Annales Botanici Societatis Zoologicae-Botanicae Fennicae 19, 1–218 (1944).
Google Scholar 
Wallnöfer, B. Uber Carex melanostachya, C. norvegica, C. cespitosa und C. hartmanii in Südtirol. Gredleriana 4, 413–418 (2004).
Google Scholar 
Gebauer, S., Röser, M. & Hoffmann, M. H. Molecular phylogeny of the species-rich Carex sect. Racemosae (Cyperaceae) based on four nuclear and chloroplast markers. Syst. Bot. 40, 433–447 (2015).Article 

Google Scholar 
Molina, A., Acedo, C. & Llamas, F. Taxonomy and new taxa in Eurasian Carex (Section Phaestoglochin, Cyperaceae). Syst. Bot. 33, 237–250 (2008).Article 

Google Scholar 
Molina, A., Acedo, C. & Llamas, F. Taxonomy and new taxa of the Carex divulsa aggregate in Eurasia (section Phaestoglochin, Cyperaceae). Botan. J. Linn. Soc. 156, 385–409 (2008).Article 

Google Scholar 
Jiménez-Mejías, P. & Luceño, M. Cyperaceae. in Euro+Med. etc. Plantbase – the information resource for Euro-Mediterranean plant diversity. Available from: http://www.emplantbase.org/home.html (accessed 07 January 2022). (eds Greuter, W. & Raab-Straube, E. von) (2011).Míguez, M., Martín-Bravo, S. & Jiménez-Mejías, P. Reconciling morphology and phylogeny allows an integrative taxonomic revision of the giant sedges of Carex section Rhynchocystis (Cyperaceae). Botan. J. Linn. Soc. 188, 34–58 (2018).Article 

Google Scholar 
Kaplan, Z. et al. Distributions of vascular plants in the Czech Republic. Preslia 93, 255–304 (2021).Article 

Google Scholar  More

Shultz, S. & Dunbar, R. Encephalization is not a universal macroevolutionary phenomenon in mammals but is associated with sociality. Proc. Natl Acad. Sci. USA 107, 21582–21586 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Jerison, H. J. Animal intelligence as encephalization. Phil. Trans. R. Soc. Lond. B 308, 21–35 (1985).CAS 
Article 

Google Scholar 
Roth, G. & Dicke, U. Evolution of the brain and intelligence. Trends Cogn. Sci. 9, 250–257 (2005).PubMed 
Article 

Google Scholar 
Lefebvre, L., Whitle, P., Lascaris, E. & Finkelstein, A. Feeding innovations and forebrain size in birds. Anim. Behav. 53, 549–560 (1997).Article 

Google Scholar 
Overington, S. E., Morand-Ferron, J., Boogert, N. J. & Lefebvre, L. Technical innovations drive the relationship between innovativeness and residual brain size in birds. Anim. Behav. 78, 1001–1010 (2009).Article 

Google Scholar 
Reader, S. M., Hager, Y. & Laland, K. N. The evolution of primate general and cultural intelligence. Phil. Trans. R. Soc. B 366, 1017–1027 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
Benson-Amram, S., Dantzer, B., Stricker, G., Swanson, E. M. & Holekamp, K. E. Brain size predicts problem-solving ability in mammalian carnivores. Proc Natl Acad. Sci. USA 113, 2532–2537 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Reader, S. M. & Laland, K. N. Social intelligence, innovation, and enhanced brain size in primates. Proc. Natl Acad. Sci. USA 99, 4436–4441 (2002).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Fristoe, T. S., Iwaniuk, A. N. & Botero, C. A. Big brains stabilize populations and facilitate colonization of variable habitats in birds. Nat. Ecol. Evol. 1, 1706–1715 (2017).PubMed 
Article 

Google Scholar 
van Woerden, J. T., van Schaik, C. P. & Isler, K. Effects of seasonality on brain size evolution: evidence from Strepsirrhine primates. Am. Nat. 176, 758–767 (2010).PubMed 
Article 

Google Scholar 
Ducatez, S., Sol, D., Sayol, F. & Lefebvre, L. Behavioural plasticity is associated with reduced extinction risk in birds. Nat. Ecol. Evol. 4, 788–793 (2020).PubMed 
Article 

Google Scholar 
Herculano-Houzel, S. Brains matter, bodies maybe not: the case for examining neuron numbers irrespective of body size. Ann. NY Acad. Sci. 1225, 191–199 (2011).PubMed 
Article 

Google Scholar 
Logan, C. J. et al. Beyond brain size: uncovering the neural correlates of behavioral and cognitive specialization. Comp. Cogn. Behav. Rev. 13, 55–89 (2018).Article 

Google Scholar 
Jerison, H. Evolution of the Brain and Intelligence (Academic Press, 1973).Herculano-Houzel, S. Numbers of neurons as biological correlates of cognitive capability. Curr. Opin. Behav. Sci. 16, 1–7 (2017).Article 

Google Scholar 
Van Schaik, C. P., Triki, Z., Bshary, R. & Heldstab, S. A. A farewell to the encephalization quotient: a new brain size measure for comparative primate cognition. Brain Behav. Evol. 96, 1–12 (2021).PubMed 
Article 

Google Scholar 
Striedter, G. F. Principles of Brain Evolution (Sinauer Associates, 2005).
Google Scholar 
MacLean, E. L. et al. The evolution of self-control. Proc. Natl Acad. Sci. USA 111, E2140–E2148 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Matějů, J. et al. Absolute, not relative brain size correlates with sociality in ground squirrels. Proc. R. Soc. B 283, 20152725 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Deaner, R. O., Isler, K., Burkart, J. & Van Schaik, C. Overall brain size, and not encephalization quotient, best predicts cognitive ability across non-human primates. Brain Behav. Evol. 70, 115–124 (2007).PubMed 
Article 

Google Scholar 
Smaers, J. B., Dechmann, D. K. N., Goswami, A., Soligo, C. & Safi, K. Comparative analyses of evolutionary rates reveal different pathways to encephalization in bats, carnivorans, and primates. Proc. Natl Acad. Sci. USA 109, 18006–18011 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Smaers, J. B. et al. The evolution of mammalian brain size. Sci. Adv. 7, eabe2101 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
Němec, P. & Osten, P. The evolution of brain structure captured in stereotyped cell count and cell type distributions. Curr. Opin. Neurobiol. 60, 176–183 (2020).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Olkowicz, S. et al. Birds have primate-like numbers of neurons in the forebrain. Proc. Natl Acad. Sci. USA 113, 7255–7260 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Kverková, K. et al. The evolution of brain neuron numbers in amniotes. Proc. Natl Acad. Sci. USA 119, e2121624119 (2022).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Iwaniuk, A. N. & Hurd, P. L. The evolution of cerebrotypes in birds. Brain Behav. Evol. 65, 215–230 (2005).PubMed 
Article 

Google Scholar 
Timmermans, S., Lefebvre, L., Boire, D. & Basu, P. Relative size of the hyperstriatum ventrale is the best predictor of feeding innovation rate in birds. Brain Behav. Evol. 56, 196–203 (2000).CAS 
PubMed 
Article 

Google Scholar 
Sayol, F., Lefebvre, L. & Sol, D. Relative brain size and its relation with the associative pallium in birds. Brain Behav. Evol. 87, 69–77 (2016).PubMed 
Article 

Google Scholar 
Healy, K. et al. Ecology and mode-of-life explain lifespan variation in birds and mammals. Proc. R. Soc. B 281, 20140298 (2014).Deaner, R. O., Barton, R. A. & van Schaik, C. P. in Primate Life Histories and Socioecology (eds Kappeler, P. M. & Pereira, M. E.) 233–265 (Univ. of Chicago Press, 2003).Sol, D., Sayol, F., Ducatez, S. & Lefebvre, L. The life-history basis of behavioural innovations. Phil. Trans. R. Soc. B 371, 20150187 (2016).PubMed 
PubMed Central 
Article 

Google Scholar 
Dukas, R. Evolutionary biology of animal cognition. Ann. Rev. Ecol. Evol. Syst. 35, 347–374 (2004).Article 

Google Scholar 
Ricklefs, R. E. The cognitive face of life histories. Wilson Bull. 116, 119–133 (2004).Article 

Google Scholar 
Martin, T. E., Oteyza, J. C., Boyce, A. J., Lloyd, P. & Ton, R. Adult mortality probability and nest predation rates explain parental effort in warming eggs with consequences for embryonic development time. Am. Nat. 186, 223–236 (2015).PubMed 
Article 

Google Scholar 
Unzeta, M., Martin, T. E. & Sol, D. Daily nest predation rates decrease with body size in passerine birds. Am. Nat. 196, 743–754 (2020).PubMed 
Article 

Google Scholar 
Charvet, C. J. & Striedter, G. F. Developmental modes and developmental mechanisms can channel brain evolution. Front. Neuroanat. 5, 4 (2011).Finlay, B. L. & Darlington, R. B. Linked regularities in the development and evolution of mammalian brains. Science 268, 1578–1584 (1995).CAS 
PubMed 
Article 

Google Scholar 
Herculano-Houzel, S. Isotropic fractionator: a simple, rapid method for the quantification of total cell and neuron numbers in the brain. J. Neurosci. 25, 2518–2521 (2005).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Massen, J. J. M. et al. Brain size and neuron numbers drive differences in yawn duration across mammals and birds. Commun. Biol. 4, 1–10 (2021).Article 

Google Scholar 
Ramsey, G., Bastian, M. L. & Schaik, C. Van Animal innovation defined and operationalized. Behav. Brain Sci. 30, 393–437 (2007).PubMed 
Article 

Google Scholar 
Lefebvre, L. A global database of feeding innovations in birds. Wilson J. Ornithol. 132, 803–809 (2021).Article 

Google Scholar 
Barton, R. A. Embodied cognitive evolution and the cerebellum. Phil. Trans. R. Soc. B 367, 2097–2107 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
Gutiérrez-Ibáñez, C., Iwaniuk, A. N. & Wylie, D. R. Parrots have evolved a primate-like telencephalic–midbrain–cerebellar circuit. Sci. Rep. 8, 9960 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Brieuc, M. S. O. O., Waters, C. D., Drinan, D. P. & Naish, K. A. A practical introduction to random forest for genetic association studies in ecology and evolution. Mol. Ecol. Res. 18, 755–766 (2018).Article 

Google Scholar 
Hadfield, J. D. & 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 
Article 

Google Scholar 
Güntürkün, O., Ströckens, F., Scarf, D. & Colombo, M. Apes, feathered apes, and pigeons: differences and similarities. Curr. Opin. Behav. Sci. 16, 35–40 (2017).Article 

Google Scholar 
Ströckens, F. et al. High associative neuron numbers could drive cognitive performance in corvid species. J. Comp. Neurol. 530, 1588–1605 (2022).PubMed 
Article 

Google Scholar 
Shanahan, M., Bingman, V. P., Shimizu, T., Wild, M. & Güntürkün, O. Large-scale network organisation in the avian forebrain: a connectivity matrix and theoretical analysis. Front. Comput. Neurosci. 7, 89 (2013).Emery, N. J. Cognitive ornithology: the evolution of avian intelligence. Phil. Trans. R. Soc. B 361, 23–43 (2006).PubMed 
Article 

Google Scholar 
Lambert, M. L., Jacobs, I., Osvath, M. & von Bayern, A. M. P. Birds of a feather? Parrot and corvid cognition compared. Behaviour 156, 505–594 (2019).Article 

Google Scholar 
Ksepka, D. T. et al. Tempo and pattern of avian brain size evolution. Curr. Biol. 30, 2026–2036 (2020).CAS 
PubMed 
Article 

Google Scholar 
Herculano-Houzel, S., Manger, P. R. & Kaas, J. H. Brain scaling in mammalian evolution as a consequence of concerted and mosaic changes in numbers of neurons and average neuronal cell size. Front. Neuroanat. 8, 77 (2014).Smaers, J. B., Mongle, C. S., Safi, K. & Dechmann, D. K. N. Allometry, evolution and development of neocortex size in mammals. Prog. Brain Res. 250, 83–107 (2019).PubMed 
Article 

Google Scholar 
Cárdenas, A. & Borrell, V. Molecular and cellular evolution of corticogenesis in amniotes. Cell Mol. Life Sci. 77, 435–1460 (2020).Article 
CAS 

Google Scholar 
García-Moreno, F. & Molnár, Z. Variations of telencephalic development that paved the way for neocortical evolution. Prog. Neurobiol. 194, 101865 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
Charvet, C. J. & Striedter, G. F. Developmental basis for telencephalon expansion in waterfowl: enlargement prior to neurogenesis. Proc. R. Soc. B 276, 3421–3427 (2009).PubMed 
PubMed Central 
Article 

Google Scholar 
Striedter, G. F. & Charvet, C. J. Developmental origins of species differences in telencephalon and tectum size: morphometric comparisons between a parakeet (Melopsittacus undulatus) and a quail (Colinus virgianus). J. Comp. Neurol. 507, 1663–1675 (2008).PubMed 
Article 

Google Scholar 
Sibly, R. M. & Brown, J. H. Effects of body size and lifestyle on evolution of mammal life histories. Proc. Natl Acad. Sci. USA 104, 17707–17712 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Uomini, N., Fairlie, J., Gray, R. D. & Griesser, M. Extended parenting and the evolution of cognition. Phil. Trans. R. Soc. Lond. B 375, 20190495 (2020).Article 

Google Scholar 
Reiner, A. et al. Revised nomenclature for avian telencephalon and some related brainstem nuclei. J. Comp. Neurol. 473, 377–414 (2004).PubMed 
PubMed Central 
Article 

Google Scholar 
Mullen, R. J., Buck, C. R. & Smith, A. M. NeuN, a neuronal specific nuclear protein in vertebrates. Development 116, 201–211 (1992).CAS 
PubMed 
Article 

Google Scholar 
Mezey, S. et al. Postnatal changes in the distribution and density of neuronal nuclei and doublecortin antigens in domestic chicks (Gallus domesticus). J. Comp. Neurol. 520, 100–116 (2012).CAS 
PubMed 
Article 

Google Scholar 
Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217–223 (2012).Article 

Google Scholar 
Ducatez, S. & Lefebvre, L. Patterns of research effort in birds. PLoS ONE 9, e89955 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Sheard, C. et al. Ecological drivers of global gradients in avian dispersal inferred from wing morphology. Nat. Commun. 11, 2463 (2020).Cooney, C. R. et al. Ecology and allometry predict the evolution of avian developmental durations. Nat. Commun. 11, 2383 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
Botelho, J. F. & Faunes, M. The evolution of developmental modes in the new avian phylogenetic tree. Evol. Dev. 17, 221–223 (2015).PubMed 
Article 

Google Scholar 
Bürkner, P.-C. Brms: an R package for Bayesian multilevel models using Stan. J. Stat. Softw. 80, 1–28 (2017).Article 

Google Scholar 
Pigot, A. L. et al. Macroevolutionary convergence connects morphological form to ecological function in birds. Nat. Ecol. Evol. 4, 230–239 (2020).PubMed 
Article 

Google Scholar 
Jetz, W., Thomas, G. H., Joy, J. B., Hartmann, K. & Mooers, A. O. The global diversity of birds in space and time. Nature 491, 444–448 (2012).CAS 
PubMed 
Article 

Google Scholar 
Berk, R. A. Statistical Learning from a Regression Perspective (Springer International, 2017).Liaw, A. & Wiener, M. Classification and regression by randomForest. R News 2, 18–22 (2002).
Google Scholar 
Lleonart, J., Salat, J. & Torres, G. J. Removing allometric effects of body size in morphological analysis. J. Theor. Biol. 205, 85–93 (2000).CAS 
PubMed 
Article 

Google Scholar 
Sayol, F., Downing, P. A., Iwaniuk, A. N., Maspons, J. & Sol, D. Predictable evolution towards larger brains in birds colonizing oceanic islands. Nat. Commun. 9, 2820 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
Torres, C. R., Norell, M. A. & Clarke, J. A. Bird neurocranial and body mass evolution across the end-Cretaceous mass extinction: the avian brain shape left other dinosaurs behind. Sci. Adv. 7, eabg7099 (2021). More

Satellite tracking has yielded key information about the movements and behaviour of marine vertebrates in ways that were previously logistically impossible34. In the current study, we tagged the first 18 angler-caught ABFT in Skagerrak, and tracked their movements for up to one year. Despite the majority of tags detaching prematurely, our data provides new insights regarding the migration behaviour and habitat use of this species, both locally within the Nordic region and more widely throughout the northeast Atlantic and western Mediterranean Sea. Most fish (N = 9) left Skagerrak via the Norwegian Trench, heading north before exiting into the Atlantic. In addition, the two tags which remained deployed for approximately one full year showed a return migration into the Skagerrak from the northern North Sea and southern Norwegian Sea regions, re-entering north of the British Isles and through the Norwegian Trench. No fish exited or re-entered through the English Channel or the southern North Sea. These observations of entry/exit from the Skagerrak are similar to migration behaviour inferred from historical commercial fishery data in the region during the 1950s–1960s16,19. These historical records also demonstrated that some individuals migrated from the southern Norwegian Sea into the Skagerrak, Kattegat and Øresund, before leaving the area several weeks later, potentially indicating exploratory feeding on herring and mackerel, abundant in the area during this time of year. Our new tagging results confirm this behaviour among at least some of the ABFT migrating to these areas.The migration patterns revealed by our tagging study exposes tuna entering and exiting the Skagerrak, Kattegat and Øresund to targeted exploitation by regional commercial fishing vessels. Presently, these vessels catch ABFT under a Norwegian quota (315 tonnes in 2021) but additional countries in the region may acquire a quota in the future. Moreover, the relatively narrow size distribution of tunas caught indicates that this migratory behaviour may only be performed by a limited number of year classes35, meaning that the continued long-term migration of ABFT to these waters is highly dependent on recruitment and survival of younger year classes. These younger year classes, perhaps once they reach a certain size, could then also undertake a migration to Skagerrak–Kattegat–Øresund. However, the combination of local exploitation pressures, and the presently limited number of year classes found in Skagerrak could result in ABFT migrating into Skagerrak–Kattegat–Øresund being a short-lived phenomenon if those year classes are subject to a large yearly fishing mortality (both regionally within the Nordic region, and more generally throughout the population range) and no younger year classes appear. Additionally, currently there is no scientifically-derived estimates of ABFT abundance for this region. We suggest to monitor the size distribution and abundance of ABFT in Scandinavian waters in the coming years to (1) confirm that visiting ABFT consist of only a few year classes, and clarify if younger year classes begin to appear, (2) evaluate how the numbers migrating to the region annually may change over time (e.g., under different levels of exploitation, or in relation to environmental factors).While most of our tagged ABFT went north after exiting the Skagerrak, one individual turned south into the south-central North Sea before eventually leaving through the northern part of the North Sea. The region to which it migrated in the North Sea is congruent with earlier commercial catches and sightings in this region, including the Dogger Bank vicinity15,16. Although the exact routes that tagged individuals followed were not identical, no individuals used the shortest route to reach the Atlantic: from the Skagerrak through the North Sea to the English Channel, and further south to the Bay of Biscay and other southern regions. Migration along a northerly route probably reflects a trade-off between the potential for higher energetic gain from more abundant food and higher energy resources, and the longer migration distance. This could suggest that ABFT either follow the food, or simply follow the same route by which they came through learned behaviour.Three tags remained attached long enough to explore long-term migration patterns and showed widely different behaviours. One fish crossed the Atlantic and utilized areas near the Grand Banks, crossing the ICCAT management boundary between the Western and Eastern stocks of ABFT (the 45° meridian), while the other two fish remained in the eastern Atlantic. The area west of Ireland, the Bay of Biscay and the area west of Portugal appear to be important feeding areas when the fish are not in Skagerrak or the Norwegian Sea. These results reflect interconnected seascapes for foraging through the NE Atlantic. Connecting foraging grounds off Ireland and the Bay of Biscay, which was previously suggested by Ref.24 is further corroborated by one of the fish tagged in this study, which passed over the Irish continental shelf when returning to Skagerrak in 2018.Depth and temperature useWithin ICES Area 3a, ABFT were predominantly roaming the upper water column, with most observations in the upper 100 m. However, some ABFT did dive to much deeper depths, with the maximum depth recorded being 520 m, showing that they can use the majority of the depth range available in the area (max. depth in the Norwegian Trench is app. 725 m, but represents a relatively small area). The behaviour likely reflects foraging, as ABFT were also observed by both the scientific tagging crews and the anglers to actively chase prey fish, like garfish and mackerel, at the surface during the tagging operations. The temperature ranges recorded varied between 7 and 17 °C. Both the depths and temperatures recorded are well within the thermal and depth limits reported in the literature for ABFT36.SpawningABFT have been shown to successfully spawn at temperatures above 20 °C at night30,31, and to display a distinct dive pattern thought to represent courtship and spawning behaviour29. When matching this described behaviour with the data from fish 34859 in the Mediterranean Sea, almost identical behavioural patterns were detected on specific days (Fig. 4). In total, seven days aligned with temperatures above 20 °C and oscillatory movement past the thermocline. All detected spawning events occurred west of Sardinia, where fishing for mature ABFT has been conducted for centuries37.In light of the recently proposed third spawning area in the Slope Sea of northeast United States38 and other proposed areas outside the Mediterranean19, it is relevant to look for similar temperature and behavioural patterns for fish 34840, which did not enter the Mediterranean Sea, and instead stayed in the eastern Atlantic. We found that this fish did not display a similar oscillatory behaviour, and the temperature experienced during the alleged spawning period (June–July) was above 20 °C only once (20.4 °C on 11 July). In this period, the fish was on the continental shelf west of Ireland, likely feeding and not spawning. Due to the size of the fish (247 cm CFL), reflecting a likely age of 14–16 years (matching the strong 2003 cohort), and the assumption that all eastern ABFT above five years and western ABFT above eight years are mature, we find it unlikely that this fish was immature. As such, these observations may suggest that this fish skipped spawning in 2018. Fish 34861 surfaced on 25 April and the tag was not recovered. The transmitted data does not allow for a detailed analysis of potential spawning behaviour for this fish. It did however, display 6 days where maximum temperatures from the transmitted dataset reached 20 °C (observations from 15. March to 20 April, with temperatures ranging from 20 to 20.6 °C). Given the lack of detailed behaviour and the fact that this time is well outside the normal spawning time for Mediterranean ABFT, we propose that this ABFT did not spawn in that period. However, the documentation of spawning depends on the general applicability of the temperature limits and nightly spawning behaviour30,31. More studies documenting spawning behaviour will be needed to corroborate if this pattern is consistent among locations and stocks. We also suggest more studies with longer lasting tags to elucidate if skipped spawning is a common behaviour and if fish skip one or more consecutive spawning seasons. Skipped spawning has been demonstrated in many fish species, including both freshwater and marine fish39, and likely reflects physiological condition40. If a considerable proportion of the adult population skips spawning every season, current population models, which assume annual spawning by all adult fish, should be modified to more accurately reflect population egg production and reproductive output. Current population modelling may be even further challenged if the proportion of adults that skip spawning varies over time, perhaps depending on environmental conditions. However, we acknowledge that only one of two fish followed through the spawning season appeared to skip spawning, and therefore caution against broad general interpretations. More studies are needed to verify that skipped spawning is a common behaviour, and if so, to estimate just how common that behaviour is.
Return migrationIn exploited fish populations, large adults are hypothesized to be important components of the spawning population because they contribute more to recruitment than smaller individuals due to a variety of maternal effects including higher fecundity, better quality of eggs and differences in spawning behaviour (e.g. time, location)41. Although such effects remain to be documented for ABFT, it may be prudent to conserve these large individuals as a precautionary measure, to maximize their potential contributions to reproduction and recruitment.In order to protect these fish, new knowledge about their movements and distribution is required. Data from ABFT deployed with long-term electronic tags suggests that after spawning in the Gulf of Mexico, the fish return to the feeding grounds where they were initially tagged, indicating a return feeding migration7. The same has been observed more recently from ABFT tagged in Ireland24, and other large highly migratory fish species (e.g., swordfish, Xiphias gladius42). In the current study, both ABFT that retained the tag for one year also returned to the same area, suggesting a similar seasonal return feeding migration. We also note that ABFT appeared to perform recurrent visits to the Norwegian Sea, Ireland and the Bay of Biscay on their way from Nordic waters and upon their return to the latter. Hence, we hypothesize that large ABFT in Nordic waters generally return to the same feeding area the following year, given suitable habitat features (e.g., food and temperature conditions), and follow a similar migration route as they do so. More studies are nonetheless needed to confirm this hypothesis, given few long-term deployments in the current study. For a deeper understanding of behavioural repeatability, and if/when shifts in the behaviour occur, it will be necessary to follow the same fish over multiple years. Such studies would also act as a highly valuable indicator of survival, independent of stock assessment-derived mortality estimates, and could be used to estimate the local abundance of larger ABFT43. Thus, a promising avenue for future research would be to deploy long-lasting ( > 5 to 10 years) acoustic tags and use existing infrastructure from networks such as the European Tracking Network to track these large fish over the next decade44. Given that ABFT appear to return to the area annually, we suggest that Skagerrak is a promising area for the future deployment and retrieval of PSATs and other long-lasting tags, because of the relatively easy access to locate and recover detached floating tags, given that the area is reachable from land within a few hours by boat. Retrieving PSATs that have detached from animals enables scientists to access full datasets (in the present case with 5 s resolution, rather than the much coarser and variable resolution typically transmitted). This much higher resolution enables much more detailed analysis, as shown in our analysis of spawning behaviour. Additionally, floating Pop-off Data Storage Tags (PDST) tags may also be a prominent and less costly avenue forward as the geographical region is densely populated and contains many sandy beaches and highly visited coastal areas, giving ample opportunity for tag recovery. Previous studies with floating DSTs in this area have shown remarkably high return rates45.The evidence that ABFT have returned to Nordic waters following many years of rarity or absence, and our findings that at least some individuals return to the same site for feeding in consecutive years, raises new questions about the mechanisms that underlie habitat discovery—or the return to previously used habitats—by highly migratory fish species. How individuals or entire schools have discovered this region again as a suitable feeding area after an absence of more than 50 years is unclear. In light of the positive stock development in the last 1–2 decades22 and modelling studies showing suitable habitat in the area46, density-dependent foraging and exploratory behaviour for new feeding areas may be a prominent hypothesis for their return, potentially accompanied by complex social learning interactions among individuals within the population47,48. New tagging data which documents the use of new or formerly occupied habitats will be essential for understanding these processes and how they might be affected by human pressures (e.g., exploitation, climate change). Such data can help to parameterize and validate advanced conceptual models of group movement behaviour, collective memory and habitat use49,50,51, as well as to inform modern stock assessment models used for management.
Tag deploymentFollowing recommendations from experienced taggers previously operating in the Mediterranean, most fish were tagged in the water alongside the boat. All these tags surfaced prematurely, while two (out of three) tags deployed on tunas brought on board the tagging boat surfaced after approximately one year. Depending on the conditions at sea, tagging along the side of the boat may not be as precise as on-board tagging, and the quality of the tag anchoring cannot be properly assessed. We therefore suggest that tagging on-board a boat is superior to tagging in the water alongside the boat for the deployment of long-lasting tags. This was also suggested in Ref.24. Furthermore, on-board tagging makes biological sampling fast and feasible, as opposed to tagging in the water alongside the boat. However, our advice is limited by a small sample size, making it difficult to draw formal conclusions; more studies are necessary to assess the best method to tag large ABFT. More