Philippa Kaur

Martin, G., Martin-Clouaire, R. & Duru, M. Farming system design to feed the changing world. A review. Agron. Sustain. Dev. 33, 131–149 (2013).
Google Scholar 
McElwee, G. & Bosworth, G. Exploring the strategic skills of farmers across a typology of farm diversification approaches. J. Farm Manag. 13, 819–838 (2010).
Google Scholar 
Maghrebi, M. et al. Iran’s agriculture in the anthropocene. Earth’s Future. https://doi.org/10.1029/2020EF001547 (2020).Article 

Google Scholar 
Raorane, A. A. & Kulkarni, R. V. Data mining: An effective tool for yield estimation in the agricultural sector. Int. J. Emerg. Trends Technol. Comput. Sci. 1, 1–4 (2012).
Google Scholar 
Gonzalez-Sanchez, A., Frausto-Solis, J. & Ojeda-Bustamante, W. Attribute selection impact on linear and nonlinear regression models for crop yield prediction. Sci. World J. 2014, 509429 (2014).
Google Scholar 
Salman, S. A. et al. Changes in climatic water availability and crop water demand for Iraq region. Sustainability 12, 3437 (2020).
Google Scholar 
Mahmood, N., Arshad, M., Kächele, H., Ullah, A. & Müller, K. Economic efficiency of rainfed wheat farmers under changing climate: Evidence from Pakistan. Environ. Sci. Pollut. Res. 27, 34453–34467 (2020).
Google Scholar 
Pracha, A. S. & Volk, T. A. An edible energy return on investment (EEROI) analysis of wheat and rice in Pakistan. Sustainability 3, 2358–2391 (2011).
Google Scholar 
Canadell, J. et al. Abberton, M., Conant, R., & Batello, C. (Eds.). (2010). Grassland carbon sequestration: Management, policy and economics. Food and Agriculture Organization of the United Nations, Integrated Crop Management, Vol. 11–2010. Ahlstrom, A., Raupach, M., Schurgers. Sensit. A Semi-Arid Grassl. To Extrem. Precip. Events 127, 6 (2021).
Google Scholar 
Canton, H. Food and Agriculture Organization of the United Nations—FAO. In The Europa Directory of International Organizations 2021 (ed. Canton, H.) 297–305 (Routledge, 2021).
Google Scholar 
Abdullah, A. et al. Potential for sustainable utilisation of agricultural residues for bioenergy production in Pakistan: An overview. J. Clean. Prod. 287, 125047 (2020).
Google Scholar 
Mughal, I. et al. Protein quantification and enzyme activity estimation of Pakistani wheat landraces. PLoS ONE 15, e0239375 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Dorosh, P. & Salam, A. Wheat markets and price stabilisation in Pakistan: An analysis of policy options. Pak. Dev. Rev. 47, 71–87 (2008).
Google Scholar 
Fowke, V. The National Policy and the Wheat Economy (University of Toronto Press, 2019).
Google Scholar 
Hussain, S. et al. Study the effects of COVID-19 in Punjab, Pakistan using space-time scan statistic for policy measures in regional agriculture and food supply chain. Environ. Sci. Pollut. Res. Int. 20, 1–14 (2021).
Google Scholar 
Sajjad, S. A. Story of Pakistan’s Elite Wheat (The Express Tribune, 2017).
Google Scholar 
Durgun, Y. Ö., Gobin, A., Duveiller, G. & Tychon, B. A study on trade-offs between spatial resolution and temporal sampling density for wheat yield estimation using both thermal and calendar time. Int. J. Appl. Earth Obs. Geoinf. 86, 101988 (2020).
Google Scholar 
Vannoppen, A. et al. Wheat yield estimation from NDVI and regional climate models in Latvia. Remote Sens. 12, 2206 (2020).ADS 

Google Scholar 
Irmak, A. et al. Artificial neural network model as a data analysis tool in precision farming. Trans. ASABE 49, 2027–2037 (2006).
Google Scholar 
Bannerjee, G., Sarkar, U., Das, S. & Ghosh, I. Artificial intelligence in agriculture: A literature survey. Int. J. Sci. Res. Comput. Sci. Appl. Manag. Stud. 7, 1–6 (2018).
Google Scholar 
Patrício, D. I. & Rieder, R. Computer vision and artificial intelligence in precision agriculture for grain crops: A systematic review. Comput. Electron. Agric. 153, 69–81 (2018).
Google Scholar 
Yaseen, Z. M. et al. Prediction of evaporation in arid and semi-arid regions: A comparative study using different machine learning models. Eng. Appl. Comput. Fluid Mech. 14, 70–89 (2019).
Google Scholar 
Bauer, M. E. The role of remote sensing in determining the distribution and yield of crops. In Advances in Agronomy (ed. Sparks, D. L.) 271–304 (Elsevier, 1975). https://doi.org/10.1016/s0065-2113(08)70012-9.Chapter 

Google Scholar 
Dempewolf, J. et al. Wheat yield forecasting for Punjab Province from vegetation index time series and historic crop statistics. Remote Sens. 6, 9653–9675 (2014).ADS 

Google Scholar 
Hamid, N., Pinckney, T. C., Gnaegy, S. & Valdes, A. The Wheat Economy of Pakistan: Setting and Prospects (IFPRI, 2015).
Google Scholar 
Muhammad, K. Description of the Historical Background of Wheat Improvement in Baluchistan, Pakistan (Agriculture Research Institute (Sariab, Quetta, Baluchistan, Pakistan), 1989).
Google Scholar 
Iqbal, N., Bakhsh, K., Maqbool, A. & Abid Shohab, A. Use of the ARIMA model for forecasting wheat area and production in Pakistan. J. Agric. Soc. Sci. 1, 120–122 (2005).
Google Scholar 
Sher, F. & Ahmad, E. Forecasting wheat production in Pakistan. LAHORE J. Econ. 13, 57–85 (2008).
Google Scholar 
Khan, N. et al. Determination of cotton and wheat yield using the standard precipitation evaporation index in Pakistan. Arab. J. Geosci. 14, 1–16 (2021).
Google Scholar 
Rahman, M. M., Haq, N. & Rahman, R. M. Machine learning facilitated rice prediction in Bangladesh. In 2014 Annual Global Online Conference on Information and Computer Technology. https://doi.org/10.1109/gocict.2014.9 (2014).Chen, C. & Mcnairn, H. A neural network integrated approach for rice crop monitoring. Int. J. Remote Sens. 27, 1367–1393 (2006).
Google Scholar 
Kaul, M., Hill, R. L. & Walthall, C. Artificial neural networks for corn and soybean yield prediction. Agric. Syst. 85, 1–18 (2005).
Google Scholar 
Deo, R. C., Samui, P., Kisi, O. & Yaseen, Z. M. Intelligent Data Analytics for Decision-Support Systems in Hazard Mitigation: Theory and Practice of Hazard Mitigation (Springer Nature, 2020).
Google Scholar 
Sanikhani, H. et al. Survey of different data-intelligent modeling strategies for forecasting air temperature using geographic information as model predictors. Comput. Electron. Agric. 152, 242–260 (2018).
Google Scholar 
Hai, T. et al. Global solar radiation estimation and climatic variability analysis using extreme learning machine based predictive model. IEEE Access 8, 12026–12042 (2020).
Google Scholar 
Ramos, A. P. M. et al. A random forest ranking approach to predict yield in maize with UAV-based vegetation spectral indices. Comput. Electron. Agric. 178, 105791 (2020).
Google Scholar 
Suchithra, M. S. & Pai, M. L. Improving the prediction accuracy of soil nutrient classification by optimizing extreme learning machine parameters. Inf. Process. Agric. 7, 72–82 (2020).
Google Scholar 
Feng, Z., Huang, G. & Chi, D. Classification of the complex agricultural planting structure with a semi-supervised extreme learning machine framework. Remote Sens. 12, 3708 (2020).ADS 

Google Scholar 
Tur, R. & Yontem, S. A comparison of soft computing methods for the prediction of wave height parameters. Knowl. Based Eng. Sci. 2, 31–46 (2021).
Google Scholar 
Yaseen, Z. M., Ali, M., Sharafati, A., Al-Ansari, N. & Shahid, S. Forecasting standardized precipitation index using data intelligence models: Regional investigation of Bangladesh. Sci. Rep. 11, 1–25 (2021).
Google Scholar 
Sharafati, A., Asadollah, S. B. H. S. & Neshat, A. A new artificial intelligence strategy for predicting the groundwater level over the Rafsanjan aquifer in Iran. J. Hydrol. https://doi.org/10.1016/j.jhydrol.2020.125468 (2020).Article 

Google Scholar 
Huang, G.-B., Zhu, Q.-Y. & Siew, C.-K. Extreme learning machine: Theory and applications. Neurocomputing 70, 489–501 (2006).
Google Scholar 
Adnan, R. M. et al. Improving streamflow prediction using a new hybrid ELM model combined with hybrid particle swarm optimization and grey wolf optimization. Knowl. Based Syst. 230, 107379 (2021).
Google Scholar 
Yaseen, Z. M. et al. Stream-flow forecasting using extreme learning machines: A case study in a semi-arid region in Iraq. J. Hydrol. 542, 603–614 (2016).ADS 

Google Scholar 
Prasad, R., Deo, R. C., Li, Y. & Maraseni, T. Ensemble committee-based data intelligent approach for generating soil moisture forecasts with multivariate hydro-meteorological predictors. Soil Tillage Res. https://doi.org/10.1016/j.still.2018.03.021 (2018).Article 

Google Scholar 
Tiyasha, T. et al. Functionalization of remote sensing and on-site data for simulating surface water dissolved oxygen: Development of hybrid tree-based artificial intelligence models. Mar. Pollut. Bull. 170, 112639 (2021).CAS 
PubMed 

Google Scholar 
Ali, M. et al. Variational mode decomposition based random forest model for solar radiation forecasting: New emerging machine learning technology. Energy Rep. 7, 6700–6717 (2021).
Google Scholar 
Khozani, Z. S. et al. Determination of compound channel apparent shear stress: Application of novel data mining models. J. Hydroinform. 21, 798–811 (2019).MathSciNet 

Google Scholar 
Dorigo, M. & Di Caro, G. Ant colony optimization: A new meta-heuristic. In Proceedings of the 1999 Congress on Evolutionary Computation, CEC 1999. https://doi.org/10.1109/CEC.1999.782657 (1999).Mullen, R. J., Monekosso, D., Barman, S. & Remagnino, P. A review of ant algorithms. Expert Syst. Appl. https://doi.org/10.1016/j.eswa.2009.01.020 (2009).Article 

Google Scholar 
Sweetlin, J. D., Nehemiah, H. K. & Kannan, A. Feature selection using ant colony optimization with tandem-run recruitment to diagnose bronchitis from CT scan images. Comput. Methods Prog. Biomed. https://doi.org/10.1016/j.cmpb.2017.04.009 (2017).Article 

Google Scholar 
Cordon, O., Herrera, F. & Stützle, T. A review on the ant colony optimization metaheuristic: Basis, models and new trends. Mathw. Comput. 9, 2–3 (2002).MathSciNet 
MATH 

Google Scholar 
Singh, G., Kumar, N. & Kumar Verma, A. Ant colony algorithms in MANETs: A review. J. Netw. Comput. Appl. https://doi.org/10.1016/j.jnca.2012.07.018 (2012).Article 

Google Scholar 
Kumar, S., Solanki, V. K., Choudhary, S. K., Selamat, A. & González Crespo, R. Comparative study on ant colony optimization (ACO) and K-means clustering approaches for jobs scheduling and energy optimization model in internet of things (IoT). Int. J. Interact. Multimed. Artif. Intell. 6, 107 (2020).
Google Scholar 
Paniri, M., Dowlatshahi, M. B. & Nezamabadi-pour, H. MLACO: A multi-label feature selection algorithm based on ant colony optimization. Knowl. Based Syst. 192, 105285 (2020).
Google Scholar 
Yaseen, Z. M., Sulaiman, S. O., Deo, R. C. & Chau, K.-W. An enhanced extreme learning machine model for river flow forecasting: State-of-the-art, practical applications in water resource engineering area and future research direction. J. Hydrol. 569, 387–408 (2019).ADS 

Google Scholar 
Manju Parkavi, R., Shanthi, M. & Bhuvaneshwari, M. C. Recent trends in ELM and MLELM: A review. Adv. Sci. Technol. Eng. Syst. https://doi.org/10.25046/aj020108 (2017).Article 

Google Scholar 
Araba, A. M., Memon, Z. A., Alhawat, M., Ali, M. & Milad, A. Estimation at completion in Civil engineering projects: Review of regression and soft computing models. Knowl. Based Eng. Sci. 2, 1–12 (2021).
Google Scholar 
Tamura, S. & Tateishi, M. Capabilities of a four-layered feedforward neural network: Four layers versus three. IEEE Trans. Neural Netw. 8, 251–255 (1997).CAS 
PubMed 

Google Scholar 
Huang, G.-B. Learning capability and storage capacity of two-hidden-layer feedforward networks. IEEE Trans. Neural Netw. 14, 274–281 (2003).PubMed 

Google Scholar 
Ali, M., Deo, R. C., Downs, N. J. & Maraseni, T. Multi-stage hybridized online sequential extreme learning machine integrated with Markov Chain Monte Carlo copula-Bat algorithm for rainfall forecasting. Atmos. Res. 213, 450–464 (2018).
Google Scholar 
Liang, N.-Y., Huang, G.-B., Saratchandran, P. & Sundararajan, N. A fast and accurate online sequential learning algorithm for feedforward networks. IEEE Trans. Neural Netw. 17, 1411–1423 (2006).PubMed 

Google Scholar 
Lan, Y., Soh, Y. C. & Huang, G.-B. Ensemble of online sequential extreme learning machine. Neurocomputing 72, 3391–3395 (2009).
Google Scholar 
Yadav, B., Ch, S., Mathur, S. & Adamowski, J. Discharge forecasting using an online sequential extreme learning machine (OS-ELM) model: A case study in Neckar River, Germany. Measurement 92, 433–445 (2016).ADS 

Google Scholar 
Breiman, L. Bagging predictors. Mach. Learn. 24, 123–140 (1996).MATH 

Google Scholar 
Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).MATH 

Google Scholar 
Al-Sulttani, A. O. et al. Proposition of new ensemble data-intelligence models for surface water quality prediction. IEEE Access 9, 108527–108541 (2021).
Google Scholar 
Carranza, C., Nolet, C., Pezij, M. & Van Der Ploeg, M. Root zone soil moisture estimation with random forest. J. Hydrol. 593, 125840 (2021).
Google Scholar 
Evans, J. S., Murphy, M. A., Holden, Z. A. & Cushman, S. A. Modeling species distribution and change using random forest. In Predictive Species and Habitat Modeling in Landscape Ecology (eds Ashton Drew, C. et al.) 139–159 (Springer, 2011).
Google Scholar 
Rahmati, O., Pourghasemi, H. R. & Melesse, A. M. Application of GIS-based data driven random forest and maximum entropy models for groundwater potential mapping: A case study at Mehran Region, Iran. CATENA 137, 360–372 (2016).
Google Scholar 
Prasad, R., Ali, M., Kwan, P. & Khan, H. Designing a multi-stage multivariate empirical mode decomposition coupled with ant colony optimization and random forest model to forecast monthly solar radiation. Appl. Energy 236, 778–792 (2019).
Google Scholar 
Sharafati, A. et al. The potential of novel data mining models for global solar radiation prediction. Int. J. Environ. Sci. Technol. https://doi.org/10.1007/s13762-019-02344-0 (2019).Article 

Google Scholar 
Service, A. M. I. District-Wise Area of Wheat Crop. Available at: http://www.amis.pk/Agristatistics/DistrictWise/2010-2012/Wheat.html (2012).Service, A. M. I. District-Wise Area of Wheat Crop. Available at: http://www.amis.pk/Agristatistics/DistrictWise/2012-2014/Wheat.html (2014).Punjab, P. Population. Available at: https://en.wikipedia.org/wiki/Punjab_Pakistan (2015).Steiniger, S. & Hunter, A. J. S. The 2012 free and open source GIS software map—A guide to facilitate research, development, and adoption. Comput. Environ. Urban Syst. 39, 136–150 (2013).
Google Scholar 
Hsu, C.-W. et al. A practical guide to support vector classification. BJU Int. https://doi.org/10.1177/02632760022050997 (2008).Article 
PubMed 

Google Scholar 
Bergmeir, C. & Benítez, J. M. On the use of cross-validation for time series predictor evaluation. Inf. Sci. (NY) 191, 192–213 (2012).
Google Scholar 
Xia, Y., Liu, C., Li, Y. Y. & Liu, N. A boosted decision tree approach using Bayesian hyper-parameter optimization for credit scoring. Expert Syst. Appl. https://doi.org/10.1016/j.eswa.2017.02.017 (2017).Article 

Google Scholar 
Yen, B. C., ASCE Task Committee on Definition of Criteria for Evaluation of Watershed Models of the Watershed Management Committee Irrigation and Drainage Division. Discussion and closure: Criteria for evaluation of watershed models. J. Irrig. Drain. Eng. 121, 130–132 (1995).
Google Scholar 
Yaseen, Z. M. An insight into machine learning models era in simulating soil, water bodies and adsorption heavy metals: Review, challenges and solutions. Chemosphere 277, 130126 (2021).ADS 
CAS 
PubMed 

Google Scholar 
Dawson, C. W., Abrahart, R. J. & See, L. M. HydroTest: A web-based toolbox of evaluation metrics for the standardised assessment of hydrological forecasts. Environ. Model. Softw. 22, 1034–1052 (2007).
Google Scholar 
Legates, D. R. & Mccabe, G. J. Evaluating the use of ‘goodness-of-fit’ measures in hydrologic and hydroclimatic model validation. Water Resour. Res. 35, 233–241 (1999).ADS 

Google Scholar 
Willmott, C. J. & Willmott, C. J. Some comments on the evaluation of model performance. Bull. Am. Meteorol. Soc. https://doi.org/10.1175/1520-0477(1982)063%3c1309:SCOTEO%3e2.0.CO;2 (1982).Article 
MATH 

Google Scholar 
Willmott, C. J. On the validation of models. Phys. Geogr. https://doi.org/10.1080/02723646.1981.10642213 (1981).Article 
MATH 

Google Scholar 
Sharafati, A., Yasa, R. & Azamathulla, H. M. Assessment of stochastic approaches in prediction of wave-induced pipeline scour depth. J. Pipeline Syst. Eng. Pract. 9, 04018024 (2018).
Google Scholar 
Mohammadi, K. et al. A new hybrid support vector machine-wavelet transform approach for estimation of horizontal global solar radiation. Energy Convers. Manag. 92, 162–171 (2015).
Google Scholar 
Willmott, C. J., Robeson, S. M. & Matsuura, K. A refined index of model performance. Int. J. Climatol. 32, 2088–2094 (2012).
Google Scholar 
Nash, J. E. & Sutcliffe, J. V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 10, 282–290 (1970).ADS 

Google Scholar 
Yaseen, Z. M. et al. Hourly river flow forecasting: Application of emotional neural network versus multiple machine learning paradigms. Water Resour. Manag. 34, 1075–1091 (2020).
Google Scholar 
Bhagat, S. K., Tung, T. M. & Yaseen, Z. M. Heavy metal contamination prediction using ensemble model: Case study of Bay sedimentation, Australia. J. Hazard. Mater. 403, 123492 (2021).CAS 
PubMed 

Google Scholar 
Hora, J. & Campos, P. A review of performance criteria to validate simulation models. Expert Syst. 32, 578–595 (2015).
Google Scholar 
Nourani, V., Kisi, Ö. & Komasi, M. Two hybrid Artificial Intelligence approaches for modeling rainfall-runoff process. J. Hydrol. https://doi.org/10.1016/j.jhydrol.2011.03.002 (2011).Article 

Google Scholar 
Ertekin, C. & Yaldiz, O. Comparison of some existing models for estimating global solar radiation for Antalya (Turkey). Energy Convers. Manag. 41, 311–330 (2000).
Google Scholar 
Li, M. F., Tang, X. P., Wu, W. & Liu, H. B. General models for estimating daily global solar radiation for different solar radiation zones in mainland China. Energy Convers. Manag. 70, 139–148. https://doi.org/10.1016/j.enconman.2013.03.004 (2013).Article 

Google Scholar 
Xu, Z., Hou, Z., Han, Y. & Guo, W. A diagram for evaluating multiple aspects of model performance in simulating vector fields. Geosci. Model Dev. 9, 4365–4380 (2016).ADS 

Google Scholar 
Dan Foresee, F. & Hagan, M. T. Gauss–Newton approximation to bayesian learning. In IEEE International Conference on Neural Networks—Conference Proceedings. https://doi.org/10.1109/ICNN.1997.614194 (1997).Akhtar, I. U. H. Pakistan needs a new crop forecasting system (2012).Stathakis, D., Savina, I. & Nègrea, T. Neuro-fuzzy modeling for crop yield prediction. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 34, 1–4 (2006).
Google Scholar 
Kumar, P., Gupta, D. K., Mishra, V. N. & Prasad, R. Comparison of support vector machine, artificial neural network, and spectral angle mapper algorithms for crop classification using LISS IV data. Int. J. Remote Sens. 36, 1604–1617 (2015).
Google Scholar 
Sun, J., Xu, W. & Feng, B. A global search strategy of quantum-behaved particle swarm optimization. In 2004 IEEE Conference on Cybernetics and Intelligent Systems. https://doi.org/10.1109/iccis.2004.1460396 (2004)Naganna, S. et al. Dew point temperature estimation: Application of artificial intelligence model integrated with nature-inspired optimization algorithms. Water. https://doi.org/10.3390/w11040742 (2019).Article 

Google Scholar 
Gilles, J. Empirical wavelet transform. IEEE Trans. Signal Process. 61, 3999–4010 (2013).ADS 
MathSciNet 
MATH 

Google Scholar 
Bokde, N., Feijóo, A., Al-Ansari, N., Tao, S. & Yaseen, Z. M. The hybridization of ensemble empirical mode decomposition with forecasting models: Application of short-term wind speed and power modeling. Energies 13, 1666 (2020).
Google Scholar 
Chau, K. W. & Wu, C. L. A hybrid model coupled with singular spectrum analysis for daily rainfall prediction. J. Hydroinform. 12, 458–473 (2010).
Google Scholar  More

Martin, P. S. & Klein, R. G. Quaternary extinctions: a prehistoric revolution. (University of Arizona Press, 1984).Waguespack, N. M. & Surovell, T. A. Clovis hunting strategies, or how to make out on plentiful resources. Am. Antiq. 68, 333–352 (2003).
Google Scholar 
Surovell, T. A., Pelton, S. R., Anderson-Sprecher, R. & Myers, A. D. Test of Martin’s overkill hypothesis using radiocarbon dates on extinct megafauna. Proc. Natl. Acad. Sci. 113, 886–891 (2016).CAS 
PubMed 
ADS 

Google Scholar 
Martin, P. S. Prehistoric overkill: the global model. In Quaternary extinctions: a prehistoric revolution (eds. Martin, P. S. & Klein, R. G.) 355–403 (University of Arizona Press, 1984).Barnosky, A. D. & Lindsey, E. L. Timing of Quaternary megafaunal extinction in South America in relation to human arrival and climate change. Quatern. Int. 217, 10–29 (2010).
Google Scholar 
Prescott, G. W., Williams, D. R., Balmford, A., Green, R. E. & Manica, A. Quantitative global analysis of the role of climate and people in explaining late Quaternary megafaunal extinctions. Proc. Natl. Acad. Sci. 109, 4527–4531 (2012).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Sandom, C., Faurby, S., Sandel, B. & Svenning, J.-C. Global late Quaternary megafauna extinctions linked to humans, not climate change. Proc. R. Soc. B Biol. Sci. 281, 20133254 (2014).
Google Scholar 
Wolfe, A. L. & Broughton, J. M. A foraging theory perspective on the associational critique of North American Pleistocene overkill. J. Archaeol. Sci. 119, 105162 (2020).
Google Scholar 
Berger, J., Swenson, J. E. & Persson, I. L. Recolonizing carnivores and naïve prey: Conservation lessons from pleistocene extinctions. Science 291, 1036–1039 (2001).CAS 
PubMed 
ADS 

Google Scholar 
Brook, B. W. & Bowman, D. M. J. S. The uncertain blitzkrieg of Pleistocene megafauna. J. Biogeogr. 31, 517–523 (2004).
Google Scholar 
Johnson, C. N. Determinants of loss of mammal species during the Late Quaternary ‘megafauna’ extinctions: life history and ecology, but not body size. Proc. R. Soc. London. Ser. B Biol. Sci. 269, 2221–2227 (2002).CAS 

Google Scholar 
Bourgon, N. et al. Trophic ecology of a Late Pleistocene early modern human from tropical Southeast Asia inferred from zinc isotopes. J. Hum. Evol. 161, 103075 (2021).PubMed 

Google Scholar 
Meltzer, D. J. Overkill, glacial history, and the extinction of North America’s Ice Age megafauna. Proc. Natl. Acad. Sci. 117, 28555–28563 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Stewart, M., Carleton, W. C. & Groucutt, H. S. Climate change, not human population growth, correlates with Late Quaternary megafauna declines in North America. Nat. Commun. 12, 965 (2021).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Nogués-Bravo, D., Rodríguez, J., Hortal, J., Batra, P. & Araújo, M. B. Climate change, humans, and the extinction of the woolly mammoth. PLoS Biol. 6, e79 (2008).PubMed 
PubMed Central 

Google Scholar 
Koch, P. L. & Barnosky, A. D. Late quaternary extinctions: State of the debate. Annu. Rev. Ecol. Evol. Syst. 37, 215–250 (2006).
Google Scholar 
Cardillo, M. Multiple causes of high extinction risk in large mammal species. Science 309, 1239–1241 (2005).CAS 
PubMed 
ADS 

Google Scholar 
Meiri, S. & Liang, T. Rensch’s rule—Definitions and statistics. Glob. Ecol. Biogeogr. 30, 573–577 (2021).
Google Scholar 
Lyons, S. K. et al. The changing role of mammal life histories in Late Quaternary extinction vulnerability on continents and islands. Biol. Lett. 12, 20160342 (2016).PubMed 
PubMed Central 

Google Scholar 
Alroy, J. A multispecies overkill simulation of the end-pleistocene megafaunal mass extinction. Science 292, 1893–1896 (2001).CAS 
PubMed 
ADS 

Google Scholar 
Smaers, J. B. et al. The evolution of mammalian brain size. Sci. Adv. 7, 1–12 (2021).
Google Scholar 
Jerison, H. J. Evolution of the Brain and Intelligence (Academic Press, 1973). https://doi.org/10.2307/4512058.Book 

Google Scholar 
Sol, D., Bacher, S., Reader, S. M. & Lefebvre, L. Brain size predicts the success of mammal species introduced into novel environments. Am. Nat. 172, S63–S71 (2008).PubMed 

Google Scholar 
Møller, A. P. & Erritzøe, J. Brain size in birds is related to traffic accidents. R. Soc. Open Sci. 4, 161040 (2017).PubMed 
PubMed Central 
ADS 

Google Scholar 
Sayol, F., Sol, D. & Pigot, A. L. Brain size and life history interact to predict urban tolerance in birds. Front. Ecol. Evol. 8, 58 (2020).
Google Scholar 
Budd, G. E. & Jensen, S. The origin of the animals and a ‘Savannah’ hypothesis for early bilaterian evolution. Biol. Rev. 92, 446–473 (2017).PubMed 

Google Scholar 
Benoit, J. et al. Brain evolution in Proboscidea (Mammalia, Afrotheria) across the Cenozoic. Sci. Rep. 9, 9323 (2019).PubMed 
PubMed Central 
ADS 

Google Scholar 
Møller, A. P. & Erritzøe, J. Brain size and the risk of getting shot. Biol. Lett. 12, 20160647 (2016).PubMed 
PubMed Central 

Google Scholar 
Di Febbraro, M. et al. Does the jack of all trades fare best? Survival and niche width in Late Pleistocene megafauna. J. Biogeogr. 44, 2828–2838 (2017).
Google Scholar 
Morris, S. D., Kearney, M. R., Johnson, C. N. & Brook, B. W. Too hot for the devil? Did climate change cause the mid-Holocene extinction of the Tasmanian devil Sacrophilus harrisii from mainland Australia? Ecography 2022, (2022).Fillios, M., Crowther, M. S. & Letnic, M. The impact of the dingo on the thylacine in Holocene Australia. World Archaeol. 44, 118–134 (2012).
Google Scholar 
González-Lagos, C., Sol, D. & Reader, S. M. Large-brained mammals live longer. J. Evol. Biol. 23, 1064–1074 (2010).PubMed 

Google Scholar 
Barton, R. A. & Capellini, I. Maternal investment, life histories, and the costs of brain growth in mammals. Proc. Natl. Acad. Sci. U.S.A. 108, 6169–6174 (2011).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Abelson, E. S. Brain size is correlated with endangerment status in mammals. Proc. R. Soc. B Biol. Sci. 283, 20152772 (2016).
Google Scholar 
Gonzalez-Voyer, A., González-Suárez, M., Vilà, C. & Revilla, E. Larger brain size indirectly increases vulnerability to extinction in mammals. Evolution (N.Y.) 70, 1364–1375 (2016).
Google Scholar 
Ives, A. R. & Helmus, M. R. Generalized linear mixed models for phylogenetic analyses of community structure. Ecol. Monogr. 81, 511–525 (2011).
Google Scholar 
Castiglione, S. et al. A new method for testing evolutionary rate variation and shifts in phenotypic evolution. Methods Ecol. Evol. 9, 974–983 (2018).
Google Scholar 
Billet, G. Phylogeny of the Notoungulata (Mammalia) based on cranial and dental characters. J. Syst. Palaeontol. 9, 481–497 (2011).
Google Scholar 
Shultz, S., Bradbury, R. B., Evans, K. L., Gregory, R. D. & Blackburn, T. M. Brain size and resource specialization predict long-term population trends in British birds. Proc. R. Soc. B Biol. Sci. 272, 2305–2311 (2005).
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 

Google Scholar 
Abelson, E. S. Big brains reduce extinction risk in Carnivora. Oecologia 191, 721–729 (2019).PubMed 
ADS 

Google Scholar 
Lundgren, E. J. et al. Introduced herbivores restore Late Pleistocene ecological functions. Proceedings of the National Academy of Sciences 117, 7871–7878 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Shultz, S. & Dunbar, R. Encephalization is not a universal macroevolutionary phenomenon in mammals but is associated with sociality. Proceedings of the National Academy of Sciences 107, 21582–21586 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Gould, S. J. & Vrba, E. S. Exaptation—A missing term in the science of form. Paleobiology 8, 4–15 (1982).
Google Scholar 
Wroe, S. et al. Climate change frames debate over the extinction of megafauna in Sahul (Pleistocene Australia-New Guinea). Proc. Natl. Acad. Sci. U.S.A. 110, 8777–8781 (2013).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
Barnosky, A. D., Koch, P. L., Feranec, R. S., Wing, S. L. & Shabel, A. B. Assessing the Causes of Late Pleistocene Extinctions on the Continents. Science 306, 70–75 (2004).Article 
PubMed 

Google Scholar 
Profico, A., Buzi, C., Melchionna, M., Veneziano, A. & Raia, P. Endomaker, a new algorithm for fully automatic extraction of cranial endocasts and the calculation of their volumes. Am. J. Phys. Anthropol. 172, 511–515 (2020).PubMed 

Google Scholar 
Damuth, J. & Macfadden, B. J. Body Size in Mammalian Paleobiology: Estimation and Biological Implications (Cambridge University Press, 1990).
Google Scholar 
Zagwijn, W. H. The beginning of the Ice Age in Europe and its major subdivisions. Quatern. Sci. Rev. 11, 583–591 (1992).ADS 

Google Scholar 
Hearty, P. J., Hollin, J. T., Neumann, A. C., O’Leary, M. J. & McCulloch, M. Global sea-level fluctuations during the Last Interglaciation (MIS 5e). Quatern. Sci. Rev. 26, 2090–2112 (2007).ADS 

Google Scholar 
Ashwell, K. W. S., Hardman, C. D. & Musser, A. M. Brain and behaviour of living and extinct echidnas. Zoology 117, 349–361 (2014).PubMed 

Google Scholar 
Castiglione, S. et al. The influence of domestication, insularity and sociality on the tempo and mode of brain size evolution in mammals. Biol. J. Linn. Soc. 132, 221–231 (2021).
Google Scholar 
Wilkins, A. S., Wrangham, R. W. & Tecumseh Fitch, W. The ‘domestication syndrome’ in mammals: A unified explanation based on neural crest cell behavior and genetics. Genetics 197, 795–808 (2014).PubMed 
PubMed Central 

Google Scholar 
Sayol, F., Steinbauer, M. J., Blackburn, T. M., Antonelli, A. & Faurby, S. Anthropogenic extinctions conceal widespread evolution of flightlessness in birds. Sci. Adv. 6, eabb6095 (2020).PubMed 
PubMed Central 
ADS 

Google Scholar 
Fromm, A., Meiri, S. & McGuire, J. Big, flightless, insular and dead: Characterising the extinct birds of the Quaternary. J. Biogeogr. 48(9), 2350–2359. https://doi.org/10.1111/jbi.14206 (2021).Article 

Google Scholar 
Meiri, S., Dayan, T. & Simberloff, D. The generality of the island rule reexamined. J. Biogeogr. 33, 1571–1577 (2006).
Google Scholar 
Larramendi, A. & Palombo, M. R. Body Size, Structure, Biology and Encephalization Quotient of Palaeoloxodon ex gr. P. falconeri from Spinagallo Cave (Hyblean plateau, Sicily). Hystrix, the Italian Journal of Mammalogy 26, 102–109 (2015).Article 

Google Scholar 
Slavenko, A., Tallowin, O. J. S., Itescu, Y., Raia, P. & Meiri, S. Late Quaternary reptile extinctions: Size matters, insularity dominates. Glob. Ecol. Biogeogr. 25, 1308–1320 (2016).
Google Scholar 
Tracy, C. R. & George, T. L. On the determinants of extinction. Am. Nat. 139, 102–122 (1992).
Google Scholar 
Manne, L. L., Brooks, T. M. & Pimm, S. L. Relative risk of extinction of passerine birds on continents and islands. Nature 399, 258–261 (1999).CAS 
ADS 

Google Scholar 
Turvey, S. T. In the shadow of the megafauna: prehistoric mammal and bird extinctions across the Holocene. in Holocene Extinctions 17–40 (Oxford University Press, 2009). https://doi.org/10.1093/acprof:oso/9780199535095.003.0002Ebinger, P. A cytoarchitectonic volumetric comparison of brains in wild and domestic sheep. Zeitschrift für Anat. und Entwicklungsgeschichte 144, 267–302 (1974).CAS 

Google Scholar 
Röhrs, M. & Ebinger, P. Welche quantitativen beziehungen bestehen bei säugetieren zwischen schädelkapazität und hirnvolumen? Mammalian Biology 66, 102–110 (2001).Köhler, M. & Moyà-Solà, S. Reduction of brain and sense organs in the fossil insular bovid Myotragus. Brain Behav. Evol. 63, 125–140 (2004).PubMed 

Google Scholar 
de Bello, F. et al. On the need for phylogenetic ‘corrections’ in functional trait-based approaches. Folia Geobot. 50, 349–357 (2015).
Google Scholar 
Bates, D., Sarkar, D., Bates, M. D. & Matrix, L. The lme4 Package. October (2007).Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in linear mixed effects models. J. Stat. Softw. 82, 1–26 (2017).
Google Scholar 
Raia, P. & Meiri, S. The tempo and mode of evolution: Body sizes of island mammals. Evolution 65, 1927–1934 (2011).

Google Scholar 
Montgomery, S. H. et al. The evolutionary history of cetacean brain and body size. Evolution 67, 3339–3353 (2013).
PubMed 

Google Scholar 
Li, D., Dinnage, R., Nell, L. A., Helmus, M. R. & Ives, A. R. phyr: An r package for phylogenetic species-distribution modelling in ecological communities. Methods Ecol. Evol. 11, 1455–1463 (2020).
Google Scholar 
Melchionna, M. et al. Macroevolutionary trends of brain mass in Primates. Biological Journal of the Linnean Society 129, 14–25 (2020).Article 

Google Scholar 
Serio, C. et al. Macroevolution of toothed whales exceptional relative brain size. Evol. Biol. 46, 332–342 (2019).
Google Scholar 
Wickham, H. et al. Welcome to the Tidyverse. Journal of Open Source Software 4, 1686 (2019).Barton, K. Package ‘MuMIn’ Title Multi-Model Inference. CRAN-R (2018). More

Moberg, F. & Folke, C. Ecological goods and services of coral reef ecosystems. Ecol. Econ. 29, 215–233 (1999).
Google Scholar 
Poloczanska, E. S. et al. Responses of marine organisms to climate change across oceans. Front. Mar. Sci. 3, 62 (2016).Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. Nature 543, 373–377(2017).Cornwall, C. E. et al. Resistance of corals and coralline algae to ocean acidification: physiological control of calcification under natural pH variability. Proc. R. Soc. B Biol. Sci. 285, 20181 (2018).
Google Scholar 
Schoepf, V. et al. Coral energy reserves and calcification in a high-CO2 world at two temperatures. PLoS One 8, e75049 (2013).CAS 

Google Scholar 
Kroeker, K. J. et al. Impacts of ocean acidification on marine organisms: quantifying sensitivities and interaction with warming. Glob. Chang. Biol. 19, 1884–1896 (2013).Heron, S. F., Maynard, J. A., Van Hooidonk, R. & Eakin, C. M. Warming trends and bleaching stress of the world’s coral reefs 1985-2012. Sci. Rep. 6, 38402, 1–14 (2016).
Google Scholar 
Marshall, A. T. & Clode, P. Calcification rate and the effect of temperature in a zooxanthellate and an azooxanthellate scleractinian reef coral. Coral Reefs. 23, 218–224 (2004).
Google Scholar 
Jokiel, P. L. & Coles, S. L. Effects of temperature on the mortality and growth of Hawaiian reef corals. Mar. Biol. 208, 201–208 (1977).
Google Scholar 
Rodolfo-Metalpa, R., Huot, Y. & Ferrier-Pagès, C. Photosynthetic response of the Mediterranean zooxanthellate coral Cladocora caespitosa to the natural range of light and temperature. J. Exp. Biol. 211, 1579–1586 (2008).CAS 

Google Scholar 
Cohen, A. L. & McConnaughey, T. A. Geochemical perspectives on coral mineralization. Rev. Mineral. Geochemistry 54, 151–187 (2003).CAS 

Google Scholar 
McCulloch, M. T., Falter, J. L., Trotter, J. & Montagna, P. Coral resilience to ocean acidification and global warming through pH up-regulation. Nat. Clim. Chang. 2, 1–5 (2012).
Google Scholar 
Venn, A. A., Tambutté, É., Holcomb, M., Allemand, D. & Tambutté, S. Live tissue imaging shows reef corals elevate pH under their calcifying tissue relative to seawater. PLoS One 6, e20013 (2011).CAS 

Google Scholar 
Cai, W.-J. et al. Microelectrode characterization of coral daytime interior pH and carbonate chemistry. Nat. Commun. 7, 11144 (2016).CAS 

Google Scholar 
Al-Horani, F. A., Al-Moghrabi, S. M. & de Beer, D. The mechanism of calcification and its relation to photosynthesis and respiration in the scleractinian coral Galaxea fascicularis. Mar. Biol. 142, 419–426 (2003).CAS 

Google Scholar 
Holcomb, M. et al. Coral calcifying fluid pH dictates response to ocean acidification. Sci. Rep. 4, 5207 (2014).CAS 

Google Scholar 
Holcomb, M., DeCarlo, T. M., Gaetani, G. A. & McCulloch, M. T. Factors affecting B/Ca ratios in synthetic aragonite. Chem. Geol. 437, 67–76 (2016).CAS 

Google Scholar 
McCulloch, M. T., D’Olivo, J. P., Falter, J., Holcomb, M. & Trotter, J. A. Coral calcification in a changing World: the interactive dynamics of pH and DIC up-regulation. Nat. Commun. 8, 15686 (2017).Tambutté, S. et al. Calcein labelling and electrophysiology: insights on coral tissue permeability and calcification. Proc. R. Soc. B 279, 19–27 (2012).
Google Scholar 
Trotter, J. et al. Quantifying the pH ‘vital effect’ in the temperate zooxanthellate coral Cladocora caespitosa: Validation of the boron seawater pH proxy. Earth Planet. Sci. Lett. 303, 163–173 (2011).CAS 

Google Scholar 
Schoepf, V., Jury, C. P., Toonen, R. & McCulloch, M. Coral calcification mechanisms facilitate adaptive response to ocean acidification. Proc. R. Soc. B 284, 2117 (2017).
Google Scholar 
Comeau, S., Cornwall, C. E. & McCulloch, M. T. Decoupling between the response of coral calcifying fluid pH and calcification to ocean acidification. Sci. Rep. 7, 7573 (2017).CAS 

Google Scholar 
Schoepf, V., D’Olivo, J. P., Rigal, C., Jung, E. M. U. & Mcculloch, M. T. Heat stress differentially impacts key calcification mechanisms in reef-building corals. Coral Reefs. https://doi.org/10.1007/s00338-020-02038-x (2021).Article 

Google Scholar 
Schoepf, V. et al. Short-term coral bleaching is not recorded by skeletal boron isotopes. PLoS One 9, e112011 (2014).
Google Scholar 
D’Olivo, J. P. & McCulloch, M. T. Response of coral calcification and calcifying fluid composition to thermally induced bleaching stress. Sci. Rep. 7, 2207 (2017).
Google Scholar 
Dishon, G. et al. A novel paleo-bleaching proxy using boron isotopes and high-resolution laser ablation to reconstruct coral bleaching events. Biogeosciences 12, 5677–5687 (2015).
Google Scholar 
Guillermic, M. et al. Thermal stress reduces pocilloporid coral resilience to ocean acidification by impairing control over calcifying fluid chemistry. Sci. Adv. 7, 20172117(2021).Ross, C. L., Falter, J. L. & McCulloch, M. T. Active modulation of the calcifying fluid carbonate chemistry (δ11B, B/Ca) and seasonally invariant coral calcification at sub-tropical limits. Sci. Rep. 7, 1–11 (2017). 13830.
Google Scholar 
D’Olivo, J. P., Ellwood, G., Decarlo, T. M. & Mcculloch, M. T. Deconvolving the long-term impacts of ocean acidification and warming on coral biomineralisation. Earth Planet. Sci. Lett. 526, 115785 (2019).
Google Scholar 
Ross, C. L., DeCarlo, T. M. & McCulloch, M. T. Environmental and physiochemical controls on coral calcification along a latitudinal temperature gradient in Western Australia. Glob. Chang. Biol. 25, 431–447 (2019).
Google Scholar 
Guo, W. Seawater temperature and buffering capacity modulate coral calcifying pH. Sci. Rep. 9, 1–13 (2019).
Google Scholar 
Reynaud, S., Ferrier-Pagès, C., Boisson, F., Allemand, D. & Fairbanks, R. G. Effect of light and temperature on calcification and strontium uptake in the scleractinian coral Acropora verweyi. Mar. Ecol. Prog. Ser. 279, 105–112 (2004).CAS 

Google Scholar 
Dissard, D. et al. Light and temperature effects on δ11B and B/Ca ratios of the zooxanthellate coral Acropora sp.: results from culturing experiments. Biogeosciences 9, 4589–4605 (2012).CAS 

Google Scholar 
Hönisch, B. et al. Assessing scleractinian corals as recorders for paleo-pH: Empirical calibration and vital effects. Geochim. Cosmochim. Acta 68, 3675–3685 (2004).
Google Scholar 
Comeau, S. et al. Flow-driven micro-scale pH variability affects the physiology of corals and coralline algae under ocean acidification. Sci. Reports 9, 1–12 (2019). 2019 91.
Google Scholar 
DeCarlo, T. M., Ross, C. L. & McCulloch, M. T. Diurnal cycles of coral calcifying fluid aragonite saturation state. Mar. Biol. 166, 1–6 (2019).CAS 

Google Scholar 
Ross, C. L., Schoepf, V., DeCarlo, T. M. & McCulloch, M. T. Mechanisms and seasonal drivers of calcification in the temperate coral Turbinaria reniformis at its latitudinal limits. Proc. R. Soc. B 285, 20180 (2018).
Google Scholar 
Krief, S. et al. Physiological and isotopic responses of scleractinian corals to ocean acidification. Geochim. Cosmochim. Acta 74, 4988–5001 (2010).CAS 

Google Scholar 
Coles, S. L. & Jokiel, P. L. Effects of temperature on photosynthesis and respiration in hermatypic corals. Mar. Biol. 43, 209–216 (1977).CAS 

Google Scholar 
Gattuso, J.-P., Allemand, D. & Frankignoulle, M. Photosynthesis and calcification at cellular, organismal and community levels in coral reefs: a review on interactions and control by carbonate chemistry. Am. Zool. 39, 160–183 (1999).CAS 

Google Scholar 
Kajiwara, K., Nagai, A., Ueno, S. & Yokochi, H. Examination of the effect of temperature, light intensity and zooxanthellae concentration on calcification and photosynthesis of scleractinian coral Acropora pulchra. J. Sch. Mar. Sci. Technol. Tokai Univ. 40, 95–103 (1995).
Google Scholar 
Reynaud, S. et al. Interacting effects of CO2 partial pressure and temperature on photosynthesis and calcification in a scleractinian coral. Glob Chang. Biol. 9, 1660–1668 (2003).
Google Scholar 
Furla, P., Galgani, I., Durand, I. & Allemand, D. Sources and mechanisms of inorganic carbon transport for coral calcification and photosynthesis. J. Exp. Biol. 203, 3445–3457 (2000).CAS 

Google Scholar 
Zoccola, D. et al. Bicarbonate transporters in corals point towards a key step in the evolution of cnidarian calcification. Sci. Rep. 5, 9983 (2015).CAS 

Google Scholar 
Allison, N. et al. Corals concentrate dissolved inorganic carbon to facilitate calcification. Nat. Commun. 5, 5741 (2014).CAS 

Google Scholar 
Vajed Samiei, J. et al. Variation in calcification rate of Acropora downingi relative to seasonal changes in environmental conditions in the northeastern Persian Gulf. Coral Reefs. https://doi.org/10.1007/s00338-016-1464-6 (2016).Article 

Google Scholar 
Kuffner, I. B., Hickey, T. D. & Morrison, J. M. Calcification rates of the massive coral Siderastrea siderea and crustose coralline algae along the Florida Keys (USA) outer-reef tract. Coral Reefs. 32, 987–997 (2013).
Google Scholar 
Courtney, T. A. et al. Environmental controls on modern scleractinian coral and reef-scale calcification. Sci. Adv. 3, e170135 (2017).
Google Scholar 
Burton, E. A. & Walter, L. M. Relative precipitation rates of aragonite and Mg calcite from seawater: Temperature or carbonate ion control? Geology 15, 111 (1987).CAS 

Google Scholar 
Lough, J. M. & Barnes, D. Environmental controls on growth of the massive coral Porites. J. Exp. Mar. Bio. Ecol. 245, 225–243 (2000).CAS 

Google Scholar 
Fitt, W. K., Brown, B., Warner, M. E. & Dunne, R. Coral bleaching: interpretation of thermal tolerance limits and thermal thresholds in tropical corals. Coral Reefs. 20, 51–65 (2001).
Google Scholar 
Fisher, R., Bessell-Browne, P. & Jones, R. Synergistic and antagonistic impacts of suspended sediments and thermal stress on corals. Nat. Commun. 10, 1–9 (2019).CAS 

Google Scholar 
Teixeira, C. D. et al. Sustained mass coral bleaching (2016–2017) in Brazilian turbid-zone reefs: taxonomic, cross-shelf and habitat-related trends. Coral Reefs. 38, 801–813 (2019).
Google Scholar 
Bonesso, J. L., Leggat, W. & Ainsworth, T. D. Exposure to elevated sea-surface temperatures below the bleaching threshold impairs coral recovery and regeneration following injury. PeerJ 5, e3719 (2017).
Google Scholar 
Ulstrup, K. E., Kühl, M., van Oppen, M. J. H., Cooper, T. F. & Ralph, P. J. Variation in photosynthesis and respiration in geographically distinct populations of two reef-building coral species. Aquat. Biol. 12, 241–248 (2011).
Google Scholar 
Lough, J. M. & Cantin, N. E. Perspectives on massive coral growth rates in a changing ocean. Biol. Bull. 226, 187–202 (2014).
Google Scholar 
Howells, E. J., Berkelmans, R., van Oppen, M. J. H., Willis, B. L. & Bay, L. K. Historical thermal regimes define limits to coral acclimatization. Ecology 94, 1078–1088 (2013).
Google Scholar 
Veron, J. E. N. Corals of the world. Townsville, Australia (Australian Institute of Marine Science, 2000).Foster, T., Short, J., Falter, J. L., Ross, C. & McCulloch, M. T. Reduced calcification in Western Australian corals during anomalously high summer water temperatures. J. Exp. Mar. Bio. Ecol. 461, 133–143 (2014).CAS 

Google Scholar 
Ross, C. L., Falter, J. L., Schoepf, V. & McCulloch, M. T. Perennial growth of hermatypic corals at Rottnest Island, Western Australia (32°S). PeerJ 3, e781 (2015).
Google Scholar 
McCulloch, M. T., Holcomb, M., Rankenburg, K. & Trotter, J. A. Rapid, high-precision measurements of boron isotopic compositions in marine carbonates. Rapid Commun. Mass Spectrom. RCM 28, 2704–2712 (2014).CAS 

Google Scholar 
Okai, T., Suzuki, A., Kawahata, H., Terashima, S. & Imai, N. Preparation of a new Geological Survey of Japan geochemical reference material: Coral JCp-1. Geostand. Newsl 26, 95–99 (2002).CAS 

Google Scholar 
Dickson, A. G. Thermodynamics of the dissociation of boric acid in synthetic seawater from 273.15 to 318.15 K. Deep Sea Res. Part A. Oceanogr. Res. Pap. 37, 755–766 (1990).CAS 

Google Scholar 
McCulloch, M. T. et al. Resilience of cold-water scleractinian corals to ocean acidification: Boron isotopic systematics of pH and saturation state up-regulation. Geochim. Cosmochim. Acta 87, 21–34 (2012).CAS 

Google Scholar 
Cornwall, C. E. & Hurd, C. L. Experimental design in ocean acidification research: problems and solutions. ICES J. Mar. Sci. 73, 572–581 (2015). More

Van Goethem MW, Pierneef R, Bezuidt OKI, Van De Peer Y, Cowan DA, Makhalanyane TP. A reservoir of ‘historical’ antibiotic resistance genes in remote pristine Antarctic soils. Microbiome. 2018;6:40.Article 

Google Scholar 
D’Costa VM, McGrann KM, Hughes DW, Wright GD. Sampling the antibiotic resistome. Science. 2006;311:374–7.Article 

Google Scholar 
Allen HK, Donato J, Wang HH, Cloud-Hansen KA, Davies J, Handelsman J. Call of the wild: antibiotic resistance genes in natural environments. Nat Rev Microbiol. 2010;8:251–9.CAS 
Article 

Google Scholar 
Martinez JL, Coque TM, Baquero F. What is a resistance gene? Ranking risk in resistomes. Nat Rev Microbiol. 2015;13:116–23.CAS 
Article 

Google Scholar 
Genilloud O. Actinomycetes: still a source of novel antibiotics. Nat Prod Rep. 2017;34:1203–32.CAS 
Article 

Google Scholar 
Ochoa-Hueso R, Plaza C, Moreno-Jimenez E, Delgado-Baquerizo M. Soil element coupling is driven by ecological context and atomic mass. Ecol Lett. 2021;24:319–26.Article 

Google Scholar 
Wardle DA, Walker LR, Bardgett RD. Ecosystem properties and forest decline in contrasting long-term chronosequences. Science. 2004;305:509–13.CAS 
Article 

Google Scholar 
Crews TE, Kitayama K, Fownes JH, Riley RH, Herbert DA, Mueller-Dombois D, et al. Changes in soil phosphorus fractions and ecosystem dynamics across a long chronosequence in Hawaii. Ecology. 1995;76:1407–24.Article 

Google Scholar 
Walker LR, Wardle DA, Bardgett RD, Clarkson BD. The use of chronosequences in studies of ecological succession and soil development. J Ecol. 2010;98:725–36.Article 

Google Scholar 
Delgado-Baquerizo M, Reich PB, Bardgett RD, Eldridge DJ, Lambers H, Wardle DA, et al. The influence of soil age on ecosystem structure and function across biomes. Nat Commun. 2020;11:4721.CAS 
Article 

Google Scholar 
Andersson DI, Hughes D. Antibiotic resistance and its cost: is it possible to reverse resistance? Nat Rev Microbiol. 2010;8:260–71.CAS 
Article 

Google Scholar 
Zhu YG, Johnson TA, Su JQ, Qiao M, Guo GX, Stedtfeld RD, et al. Diverse and abundant antibiotic resistance genes in Chinese swine farms. Proc Natl Acad Sci USA. 2013;110:3435–40.CAS 
Article 

Google Scholar 
Zhu YG, Zhao Y, Li B, Huang CL, Zhang SY, Yu S, et al. Continental-scale pollution of estuaries with antibiotic resistance genes. Nat Microbiol. 2017;2:16270.CAS 
Article 

Google Scholar 
Li J, Cao J, Zhu YG, Chen QL, Shen F, Wu Y, et al. Global survey of antibiotic resistance genes in air. Environ Sci Technol. 2018;52:10975–84.CAS 
Article 

Google Scholar 
Delgado-Baquerizo M, Bardgett RD, Vitousek PM, Maestre FT, Williams MA, Eldridge DJ, et al. Changes in belowground biodiversity during ecosystem development. Proc Natl Acad Sci USA. 2019;116:6891–6.CAS 
Article 

Google Scholar 
Ortiz-Álvarez R, Fierer N, de Los Ríos A, Casamayor EO, Barberán A. Consistent changes in the taxonomic structure and functional attributes of bacterial communities during primary succession. ISME J. 2018;12:1658–67.Article 

Google Scholar 
Shen J, Li Z-M, Hu H, Zeng J, Zhang L-M, He J, et al. Distribution and succession feature of antibiotic resistance genes along a soil development chronosequence in Urumqi No. 1 Glacier of China. Front Microbiol. 2019;10:1569.Article 

Google Scholar 
Drenovsky RE, Vo D, Graham KJ, Scow KM. Soil water content and organic carbon availability are major determinants of soil microbial community composition. Micro Ecol. 2004;48:424–30.CAS 
Article 

Google Scholar 
Bastida F, Eldridge DJ, Garcia C, Kenny Png G, Bardgett RD, Delgado-Baquerizo M. Soil microbial diversity-biomass relationships are driven by soil carbon content across global biomes. ISME J. 2021;15:2081–91.CAS 
Article 

Google Scholar  More

Co-culture dynamicsThis study was designed to enhance understanding of metabolite release and utilization across bloom stages in a simple community of phytoplankton and heterotrophic bacteria. The synthetic community was established with the diatom T. pseudonana and the bacterial strains R. pomeroyi DSS-3, Stenotrophomonas sp. SKA14, and P. dokdonensis MED152. These bacterial strains have high genetic similarity to isolates from phytoplankton cultures [14] and represent taxa that are common in phytoplankton blooms. Metabolites derived from the diatom were the sole source of carbon available for the bacteria, since no organic substrates were added. In addition, none of the bacteria can assimilate nitrate, and usable nitrogen was only available as diatom or bacterial extracellular products. The diatom had its highest specific growth rate of 1.65 d−1 during days 0–3, after which the rate declined (Fig. 1A). The total abundance of heterotrophic bacteria increased steadily but there was a succession that favored P. dokdonensis through day 15, and then R. pomeroyi by day 20; Stenotrophomonas disappeared from the model system by day 3 (Fig. 1B). The presence of bacteria did not affect the growth of diatoms based on comparisons of abundance in co-cultures versus axenic cultures at day 15 (Fig. 1A), as has been found previously [14, 26]. Inorganic nutrients were not limiting ( >5 μM at day 15; Table S1).Fig. 1: Time course of microbial abundances.A Cell abundance based on flow cytometric analysis for co-cultures (5 time points) and axenic cultures (day 15 only) (n = 3). The intensive sampling dates for the early and late bloom comparisons are marked with gray boxes. B Mean relative abundance of bacterial species is based on CFUs (n = 3). The day 0 samples were collected 8 h after inoculation.Full size imageDiatom endometabolite shiftsAnalyses focused on the day 3 (early bloom) and day 15 (late bloom) co-culture time points, for which a complete set of metabolomic and transcriptomic data were collected. Twenty-two diatom endometabolites that were annotated with high confidence by NMR analysis (Table S2) and quantified after normalizing to diatom cell number revealed that endometabolome composition differed substantially between bloom stages. Metabolites with significantly different cellular concentrations included nine compounds that were higher in intracellular concentration during the late bloom; these were arginine, valine, lysine, DHPS, glycerol-3-phosphate, phosphorylcholine, DMSP, glycine betaine, and homarine (T-test; P  More

Discovery of closely related phage sequences with the conserved genetic context of bS21Multiple phage-related sequences with a conserved genomic context were detected from several freshwater metagenome-assembled datasets (see Methods). Genes for bS21, TerL, PVP, prohead core scaffolding, and protease protein (hereafter prohead protease for short), and MCP are encoded in the genomic region. BLASTp search of the TerL sequences against the ggKbase sequences (ggkbase.berkeley.edu) obtained a total of 47 unique scaffolds with the conserved genomic region (Supplementary Table 1). Two related phages were included as outgroups for comparative analyses. The corresponding samples were collected from freshwater lakes or reservoirs (one from a wastewater treatment plant), and all but three were from the oxic layer (see Methods for details).General features of manually curated genomesAll the 49 phage sequences were manually curated to fill scaffolding gaps and fix the assembly errors, and nine of them (including one outgroup phage) were curated to completion (circular and no gaps or local assembly errors) (Supplementary Table 1). A total of 14 related phage genomes from IMG/VR were also included for further analyses. The eight bS21-encoding complete genomes had genome lengths of 293–331 kbp, GC contents of 31.0–33.7% and encoded 350–413 protein-coding genes (coding density, 91.1–94.9%), with 5–25 (average 17) tRNA genes. No alternative coding signal (i.e., stop codon reassignment) was detected in any genome. In comparison, the outgroup complete genome has a size of 308 kbp (450 protein-coding genes, 6 tRNAs, 94.7% coding density) and GC content of 27.3%.Genomic context of bS21 in phagesGenomic context analyses for bS21 genes showed a highly conserved gene architecture across phage genomes in proximity to the region encoding bS21 (see Fig. 1a for example). Specifically, we found that bS21 was consistently located in between two hypothetical protein families (positions 1 and –1 in Fig. 1b and Supplementary Table 2), with core structural proteins—including the TerL, PVP, prohead protease, and MCP—generally located within five genes in both the upstream and downstream DNA. Other hypothetical proteins were also consistently found in this region, although their positions were more variable upstream (positions –4 through –10, Fig. 1b). Importantly, the bS21 gene was consistently encoded in the reverse strand relative to the conserved hypothetical and structural protein genes (Fig. 1a and Supplementary Fig. 1).Fig. 1: Genetic context of the genes encoding bS21 in the phage genomes.a Examples of genetic context of phage genomes with and without bS21. The annotation of protein-coding genes is the same as indicated in b by different colors. Those in white are genes not shown in subfigure (b). b Summary of genetic context of all phage genomes encoding bS21. The relative position of genes near the bS21 gene is shown, and the size of circles indicates the number of phages with a gene belonging to a given protein family (annotation shown on right) at that relative position. Only the 12 most frequent families are shown. The details of the genetic context are shown in Supplementary Fig. 1.Full size imagePhylogeny of bS21-encoding phagesPhylogenetic analyses based on TerL suggested the phages belonging to several groups, we thus assigned them to clades a–e (Fig. 2 and Supplementary Table 1). Most of the phages belong to clades c, d, and e, and they have a broader environmental distribution than clades a and b. Interestingly, we found that some phages within a single clade were from distant sampling sites. Closer inspection indicated they also shared large genomic fragments with high similarity (82–98% for nucleotide sequences; Supplementary Fig. 2). Comparative genome-wide analyses of the complete genomes from the same site but sampled at different time points showed sequence variations in some genes (Supplementary Fig. 3).Fig. 2: The phylogeny of bS21 phages based on the large terminal (TerL) protein sequences.Two closely related phages without bS21 encoded were included as outgroups (shown at the top of the tree). The genomes are assigned to five clades (a, b, c, d, and e) based on the topology of the phylogenetic tree. The numbers in the brackets following the scaffold names indicate the total counts of the same scaffold detected from the corresponding sampling sites. The genomes that were manually curated to completion (circular and no gap) are indicated by squares, and the genome sizes are shown in brackets.Full size imageTerL phylogeny, constructed using sequences from this study and NCBI RefSeq sequences, indicated the most closely related classified phages belong to Caudovirales of either the Myoviridae or Ackermannviridae (Supplementary Fig. 4). A phage baseplate assembly protein was encoded in most curated genomes. This is an important building block for members of Siphoviridae and Myoviridae [8], so we concluded that the bS21-encoding phages are myoviruses.Predicted bacterial hosts of bS21-encoding phagesTo predict host-phage relationships we first used CRISPR-Cas spacers targeting. While none of the 16.5k unique spacers from the relevant metagenomes targeted any of the curated phage genomes from the same sampling sites, a single cross-site target was detected. Specifically, MIW1_072018_0_1um_scaffold_78 was targeted by a spacer (24 nt and no mismatch) from a MIW2 Flavobacterium genome (affiliation: Bacteroidetes, Flavobacteria). We then predicted the bacterial hosts based on the bacterial taxonomic affiliations of the phage gene inventories as previously described [2] (Supplementary Table 3). The results indicated that all of the phages infect members of Bacteroidetes, which were detected in 43 out of 45 samples (Fig. 3 and Supplementary Table 4). The two metagenomic samples without Bacteroidetes identified were both collected via filtering through 0.2 μm and onto 0.1 μm pore size filters. Bacteroidetes were detected in both of the corresponding 0.2 μm fraction samples (Fig. 3).Fig. 3: The relative abundance of the Bacteroidetes classes in all the analyzed samples in this study.The microbial communities were profiled based on ribosomal protein S3 (rpS3) assigned to the Bacteroidetes classes. The sampling sites were indicated by colored names, and the filter sizes used during sampling are shown by circles. The three pairs of filter samples are indicated by colored stars.Full size imageWe profiled the co-detection of phage clades and Bacteroidetes classes to test for specific connections (Supplementary Fig. 5). However, this was uninformative because most samples contained more than one class. However, phages from clades a and b are unlikely to infect class Bacteroidia members, as they did not co-occur in any sample.Comparison of bacterial and phage-encoded bS21Phylogenetic analyses revealed that bS21 protein sequences from phages (this study) and the bacterial bS21 sequences (from the corresponding samples and NCBI RefSeq) clustered separately (Supplementary Fig. 6). The bacterial bS21 sequences that are most similar to phage bS21 were from Bacteroidetes, mostly from the Flavobacteriia class (Supplementary Table 5). We aligned and compared the Bacteroidetes and phage bS21 sequences and mapped the divergent and non-divergent residues to the model of the ribosome of Flavobacterium johnsoniae (Fig. 4a). Multiple divergent positions are located at the beginning of the bS21 sequences and four residues (Arg21, Phe23, Asp25, and Thr28) were significantly divergent (Fig. 4b).Fig. 4: Conservation and differences between phage and bacterial bS21.a Location of bS21 (blue) within the 16S rRNA (green) and the ASD (magenta) of the F. johnsoniae ribosome (PDB ID: 7JIL) [9]. bS21 is in the neck region of the 16S rRNA, interacting closely with the 3’ end of the 16S rRNA, where the ASD is located. The 16S rRNA is shown from the subunit interface direction. b Zebra2 divergency results from an alignment of phage and bacterial bS21 sequences mapped on F. johnsoniae bS21. Divergent positions between phage and bacterial bS21 are shown with red. c Zebra2 conservation results from the same alignment as in (b) mapped on F. johnsoniae bS21 with conserved residues shown in yellow. The stacking interaction between Tyr54 and Adenine 1534 is indicated. d The sequence logo and consensus sequences of phage and bacterial bS21 alignments and the corresponding position of Tyr54 in F. johnsoniae bS21 in the alignment are highlighted. The C-terminal parts are highlighted with gray backgrounds.Full size imageBacteroidetes usually lack the SD sequences. It was recently reported that the bS21 Tyr54 (numbering in F. johnsoniae) is an important residue for blocking the ASD in the 16S rRNA within the ribosome [9]. Our analyses predict that all the analyzed bacterial and phage bS21 in this study have an amino acid with an aromatic ring (often Tyr54 but in a few cases His54, and in one case Phe54) at the position of Tyr54 in F. johnsoniae (Fig. 4c, d and Supplementary Fig. 6). This conservation of the aromatic property in phage bS21 should ensure stacking interaction with Adenine 1534 (numbering in F. johnsoniae 16S) from the ASD. In that way, phage bS21 mimics Bacteroidetes bS21 in the region where it binds the ribosome but differs from it in the region where the mRNA would bind.In contrast, the C-terminal regions of both the bacterial and phage bS21 sets were highly divergent (Fig. 4d). However, the phage C-terminal regions are generally conserved within the clades defined based on TerL phylogeny (Fig. 2 and Supplementary Fig. 7).Metabolic potentials of bS21-encoding phagesFunctional annotation of the predicted protein-coding genes revealed that in addition to bS21, these phages carry other genes related to protein production and stability (Supplementary Table 6). Examples include protein folding chaperones and Clp protease, suggesting the importance of controlling the proteostasis network of the cell. Interestingly, we also identified many genes involved in sugar-related chemistry and polysaccharide biosynthesis. Many of these genes were predicted to perform chemical transformations related to the biosynthesis of lipopolysaccharide, a major component of the Gram-negative bacterial outer membrane. We interpret this as a potential mechanism to remodel the cell surface and prevent superinfection by competitor phages, a strategy common to the phage lysogenic cycle. These phages lack detectable integration machinery (no gene for integrase or resolvase was detected), suggesting the possibility of a non-integrative long-term infection state such as pseudolysogeny [10].Clustering analyses of 22 phages with a minimum genome size of 100 kbp (including the two outgroup genomes) based on the presence/absence of protein families indicated they shared a total of 16 protein families (Supplementary Fig. 8 and Supplementary Table 7). Phosphate starvation-inducible protein PhoH (“fam582”) was the only predicted protein detected in all 22 phages (excluding the shared predicted proteins in the conserved rpS21-encoding region described above). Other common protein families include those related to DNA replication (e.g., DNA primase/helicase, DNA polymerase, HNH endonuclease, thymidylate synthase (EC:2.1.1.45), deoxyuridine 5’-triphosphate nucleotidohydrolase (EC:3.6.1.23)), those associated with virion assembly (e.g., a phage tail sheath protein, phage baseplate assembly protein W), and those for other functions (e.g., chaperone ATPase, alpha-amylase, DegT/DnrJ/EryC1/StrS aminotransferase).Temporal and spatial distribution and activity of bS21-encoding phages in Lake RotseeTo reveal the spatial and temporal distribution of the bS21-encoding phages, we focused on the Lake Rotsee data and profiled phage occurrence based on the sequencing coverage in the metagenomic datasets. The Lake Rotsee samples were collected from the oxic (7 samples) and anoxic (3 samples) layers of the water column. The bS21-encoding phages were readily detected in oxic samples, especially in the under-ice samples when the whole water column was oxic (Fig. 5a).Fig. 5: The spatial and temporal distribution and activity of bS21 phages at Lake Rotsee.a The sequencing coverage of each phage genome in each metagenomic dataset is shown in the heatmaps. The phages are phylogenetically clustered based on their TerL protein sequences (bootstraps shown in numbers), the colored backgrounds are the same as shown in Fig. 2 for different clades. The sampling time points and depths are shown on the left, and the oxygen conditions are indicated by colored circles on the right. Two replicates were sequenced from the 15 m sample collected in 2018. b The percentage of mapped RNA reads to the phage genomes in the corresponding samples (rows labeled in (a)). The mapped RNA reads had a minimum similarity of 98% to the phage genomes. No RNA data were generated for the three samples collected on October 10, 2017. See the figure legend for each genome in the upper right, the circular genomes have names in bold font.Full size imageRotsee Lake RNA reads were mapped to the phage genomes curated from this site to reveal the transcriptional activities of bS21-encoding phages (Fig. 5b). In general, the phages were likely to be most transcriptionally active in the oxic water columns. A total of 736 genes were transcribed in at least one sample (Supplementary Table 8), those for MCP, an AAA ATPase, tail sheath protein, bS21, FKBP-type peptidyl-prolyl cis-trans isomerase, and a methyltransferase FkbM domain protein are among the top 100 most highly transcribed. The high transcriptional activities of MCP in five phages indicated they were in the late stage of replication at the time of sampling.The transcriptional behavior of phage bS21 genesTo seek evidence of a transcriptional relationship involving bS21 and other genes we focused on the three phages that were most active based on the transcriptional level of their 19 shared single-copy genes (Fig. 6a). bS21 had very similar (but slightly lower) transcriptional activities as a neighboring gene (hereafter, bS21_CN gene) encoded on the opposite strand. The bS21_CN gene encodes a hypothetical protein (protein family: fam498) and was not detected in the two outgroup phages without bS21 (Supplementary Table 6). Interestingly, a comparison of the phylogenies of bS21 and bS21_CN showed a very similar evolutionary pattern (Supplementary Fig. 9), likely suggesting their potential functional relationship in the bS21-encoding phages.Fig. 6: The transcription levels of bS21 and core structural protein genes.a The normalized transcriptional level (NTL) of shared single-copy protein families of three phages (indicated by arrows in Fig. 5b) with ≥1000 RNA reads mapped. Two families (including MCP) are listed on a different scale due to their much higher transcription levels. Refer to Fig. 5 for shape symbols that designate phage genomes and samples. b Examples of RNA mapping profiles indicating the co-transcription of some genes neighboring bS21. Hypothetical protein genes are shown in white.Full size imageInspection of the RNA reads mapping profiles indicated that the conserved region encoding bS21 and core structural proteins was not transcribed as an operon, whereas bS21 and bS21_CN, MCP and its upstream hypothetical protein gene, and prohead protease and its downstream hypothetical protein gene may each be transcribed together (Fig. 6b). Given the observed RNA expression patterns, we conclude that the phage-encoded bS21 genes were actively transcribed during late-stage replication, along with other core structural proteins.Genomic context of bS21 genes in published phage genomesTo determine whether the phage bS21 genes are generally co-located with those for core structural proteins in diverse phages, we profiled the genomic context of bS21 in 900 published bS21-encoding phages [2, 11] (Supplementary Table 9). Functional annotations were performed for the upstream and downstream ten genes of the bS21 genes using pVOG (Supplementary Table 10). Of the 20 most abundant pVOGs, 6 were related to core structural assembly (Fig. 7a), i.e., prohead protease (n = 310), MCP (n = 154), PVP (n = 120), TerL (n = 78), neck protein (n = 70), and a tail sheath protein (n = 29). A total of 388 genomes contained at least one of these genes within ten genes of bS21, and eight had all of these six core structural proteins in close proximity. Three pVOGs were related to DNA processing, i.e., an exonuclease (n = 37), an endonuclease (n = 32), DNA helicase (n = 30). Other pVOGs included Hsp20 heat shock protein (n = 127), two ATP-dependent CLP proteases (n = 50 and 47, respectively), and lysozyme (for lysis; n = 29). Interestingly, the prohead protease and the MCP pVOG genes are very close to the bS21 gene (generally 2–4 genes; Fig. 7b), as in the bS21-encoding phage genomes analyzed in this study (2–6 genes away; Fig. 1 and Supplementary Fig. 1).Fig. 7: Neighboring genes within 10 genes of bS21 in published bS21-encoding phage genomes.a The annotation and corresponding functional category (if assigned) of the 20 most commonly detected pVOG genes and their predicted functions are shown on the left, the total number of genomes with the gene are shown on the right. b The distribution of the distance of each gene to bS21 in the genomes. The position of genes next to bS21 (thus distance = 1) is highlighted using a red dashed line. The average distance of each gene to bS21 is shown on the left. c The predicted hosts of bS21-encoding phages with the top 4 most abundant genes detected within 10 genes of bS21. The total count of hosts is shown on the right.Full size imageWe respectively predicted the hosts of the bS21-encoding phages with the four most dominant pVOGs within ten genes of bS21 (Fig. 7c and Supplementary Table 11). The bacterial hosts are diverse and include Proteobacteria, Bacteroidetes, and Firmicutes. More

Storage and displayAll collected datasets (directly downloadable as tabular files), the bibtex file with all references, all reports and all kinetic bioaccumulation metric estimates are publicly available on Zenodo17. An rmarkdown file18,19 was created to build the overview table with information collected from the name of the dataset and from the dataset itself (e.g., column headers, number of data, number of replicates), as well as from the bibtex file. The R package DT was additionally used20 to combine all collected information in a user-friendly manner including a convenient search tool, and the rmarkdown file was finally compiled19 in HTML format for display to the user in packs of 10 lines by default. In such a way, each new dataset added into the repository will compile the rmarkdown file automatically for update.In parallel, the database can also be accessed directly via http://lbbe-shiny.univ-lyon1.fr/mosaic-bioacc/data/database/TK_database.html, or from MOSAICbioacc clicking on the “More scientific TK data” button. An example of the output of the overview table is shown in Fig. 2, while the full table is provided in the supplementary information (Table S2). The collected raw TK data of the database consist in the time-course of several types of chemical substances bioaccumulated in various species via different exposure routes.Fig. 2Screenshot of the first page of the overview table of the database available from MOSAICbioacc.Full size imageDatasets overviewEach dataset is summarized by:

the file name (raw data directly downloadable by clicking on the file name, in text or CSV format),

the genus of the tested organism,

the category of the organism (e.g., aquatic, terrestrial, etc.),

the tested chemical substance,

the duration of the accumulation phase,

the tested exposure routes (e.g., water, sediment, food, pore water),

the total number of observations in the dataset (plus the number of replicate(s) in brackets),

the kinetic bioaccumulation metric median value with its 95% uncertainty interval,

the report which contains all the outputs from MOSAICbioacc (in PDF format),

the link to the reference or the source of the data,

some additional comments (e.g., lipid fraction, growth, biotransformation, if exposure was done for chemical mixtures or not, if total radioactivity was used or not, etc.).

A summary of all datasets is presented in Table 1. Genus were separated in 12 categories: aquatic invertebrates (n = 105), fish (n = 42), insects (n = 17), aquatic worms (n = 10), terrestrial worms (n = 16), seawater sponges (n = 2), seawater plants (n = 1), aquatic algae (n = 1), terrestrial invertebrates (n = 1), vertebrates other than fish (n = 4), marine invertebrates (n = 8), and heterotrichea (n = 4). The most represented genus in the database are Gammarus (aquatic invertebrate, n = 43) and Daphnia (aquatic invertebrate, n = 27), followed by Oncorhynchus (fish, n = 15), genus that are classically used in ecotoxicological experiments. Recommended genus by OECD guidelines for bioaccumulation tests are Eisenia and Enchytraeus for terrestrial organisms (OECD 317)21, and Tubifex or Lumbriculus for aquatic invertebrates exposed to sediment (OECD 315)22; some datasets for these specific species are available in the database (n = 24).Table 1 Summary of the collected TK datasets.Full size tableChemical substances were divided in 10 classes following at the best the nomenclature used in Standartox23: pesticides (n = 71), hydrocarbons (n = 37), metals (n = 20), nanoparticules (n = 23), polychlorobiphenyls (PCB, n = 22), flame retardants (brominated or chlorinated, n = 8), pharmaceutical products (n = 14), PFAS (n = 7), octyphenol (n = 2) and other (n = 7). Among all datasets, the majority of bioaccumulation tests were performed via spiked water (n = 137). Besides, 34 datasets account for biotransformation processes, considering from 1 to 8 metabolites.According to ECHA (2017)2, BCF below 1,000 means that the chemical substance is not bioaccumulative, whereas one ranging between 1,000 and 5,000 corresponds to a bioaccumulative chemical substance: low bioaccumulative if BCF ∈]1,000; 2,000]; mid-bioaccumulative if BCF ∈]2,000; 5,000]. If BCF is >5000, the chemical substance is classified as very bioaccumulative. These ranges are reported in Table 1, where BCF median estimates are >5000 for 25 datasets, indicating a very bioaccumulative capacity of the corresponding chemical substances for the corresponding genus. Concerning BSAF and BMF estimates, their value must be compared to threshold 1. A median BSAF estimate >1 indicates that the corresponding chemical substance can bioaccumulate from soil or sediment into organisms at the base of the non-aquatic food chain24,25; a median BMF estimate >1 indicates that the corresponding chemical substance can biomagnify in the trophic relationship under consideration26. In the database, 16 datasets in 36 led to BSAF >1, for genus Eisenia (n = 2), Enchytraeus (n = 6), Gallus (n = 1), Lumbriculus (n = 2), Metaphire (n = 2), Physa (n = 1), Radix (n = 2)), while 8 datasets in 38 led to BMF >1, for genus Gallus (n = 1), Oncorhynchus (n = 5) and Perca (n = 2). On an ecotoxicological point of view, the highest BCF estimates were obtained for genus Culex and Sialis exposed to chlorpyrifos due to a very low estimate of the elimination rate, for genus Gammarus and Calanus exposed to hydrocarbons, and several aquatic invertebrates exposed to pesticides, especially chlorpyrifos (n = 4), attesting to the potential high bioaccumulation capacity and high risk of toxicity associated with this chemical substance for aquatic organisms. Overall, aquatic invertebrates seem to be the most sensitive category of organisms in terms of bioaccumulation of chemical substances representing 20 in the 25 datasets with a BCF estimates >5000. More

Kemper N. Veterinary antibiotics in the aquatic and terrestrial environment. Ecol Indic. 2008;8:1–13.CAS 
Article 

Google Scholar 
Jechalke S, Heuer H, Siemens J, Amelung W, Smalla K. Fate and effects of veterinary antibiotics in soil. Trends Microbiol. 2014;22:536–45. Available from: https://doi.org/10.1016/j.tim.2014.05.005.CAS 
Article 
PubMed 

Google Scholar 
Kalasseril S, Paul R, J RK V, Pillai D. Investigating the impact of hospital antibiotic usage on aquatic environment and aquaculture systems: A molecular study of quinolone resistance in Escherichia coli. Sci Total Environ. 2020;748:141538. Available from: https://doi.org/10.1016/j.scitotenv.2020.141538.CAS 
Article 

Google Scholar 
Ashbolt NJ. Human Health Risk Assessment (HHRA) for Environmental Development and Transfer of Antibiotic Resistance. Environ Health Perspect. 2013;121:993–1002.Article 

Google Scholar 
Bengtsson-Palme J, Kristiansson E, Larsson DGJ Environmental factors influencing the development and spread of antibiotic resistance. FEMS Microbiol Rev. 2017;(October 2017):68–80. Available from: http://academic.oup.com/femsre/advance-article/doi/10.1093/femsre/fux053/4563583Manaia CM Assessing the Risk of Antibiotic Resistance Transmission from the Environment to Humans: Non-Direct Proportionality between Abundance and Risk. Vol. 25, Trends in Microbiology. 2017.Manaia CM, Macedo G, Fatta-Kassinos D, Nunes OC. Antibiotic resistance in urban aquatic environments: can it be controlled? Appl Microbiol Biotechnol. 2016;100:1543–57.CAS 
Article 

Google Scholar 
Durso LM, Cook KL. Impacts of antibiotic use in agriculture: what are the benefits and risks? Curr Opin Microbiol. 2014;19:37–44. https://doi.org/10.1016/j.mib.2014.05.019. Available fromArticle 
PubMed 

Google Scholar 
Almakki A, Jumas-Bilak E, Marchandin H, Licznar-Fajardo P. Antibiotic resistance in urban runoff. Sci Total Environ. 2019;667:64–76. https://linkinghub.elsevier.com/retrieve/pii/S0048969719306710.CAS 
Article 

Google Scholar 
Andersson DI, Hughes D. Microbiological effects of sublethal levels of antibiotics. Nat Rev Microbiol. 2014;12:465–78. Available from: https://doi.org/10.1038/nrmicro3270.CAS 
Article 
PubMed 

Google Scholar 
Gullberg E, Cao S, Berg OG, Ilbäck C, Sandegren L, Hughes D, et al. Selection of resistant bacteria at very low antibiotic concentrations. PLoS Pathog. 2011;7:1–9.Article 

Google Scholar 
Murray AK, Zhang L, Yin X, Zhang T, Buckling A, Snape J, et al. Novel insights into selection for antibiotic resistance in complex microbial communities. MBio. 2018;9:1–12. http://mbio.asm.org/lookup/doi/10.1128/mBio.00969-18.CAS 
Article 

Google Scholar 
Chow L, Waldron L, Gillings MR. Potential impacts of aquatic pollutants: sub-clinical antibiotic concentrations induce genome changes and promote antibiotic resistance. Front Microbiol. 2015;6:1–10.
Google Scholar 
Bruchmann J, Kirchen S, Schwartz T. Sub-inhibitory concentrations of antibiotics and wastewater influencing biofilm formation and gene expression of multi-resistant Pseudomonas aeruginosa wastewater isolates. Environ Sci Pollut Res. 2013;20:3539–49.CAS 
Article 

Google Scholar 
Gullberg E, Albrecht LM, Karlsson C, Sandegren L, Andersson DI. Selection of a Multidrug Resistance Plasmid by Sublethal Levels of Antibiotics and Heavy Metals. mBio. 2014;5:19–23.Article 

Google Scholar 
Choung S, Yun Z, Kwon EE, Cho Y, Ha U-H, Oh J, et al. Transfer of antibiotic resistance plasmids in pure and activated sludge cultures in the presence of environmentally representative micro-contaminant concentrations. Sci Total Environ. 2014;468–469:813–20. https://doi.org/10.1016/j.scitotenv.2013.08.100.CAS 
Article 
PubMed 

Google Scholar 
Shun-Mei E, Zeng JM, Yuan H, Lu Y, Cai RX, Chen C. Sub-inhibitory concentrations of fluoroquinolones increase conjugation frequency. Microb Pathog. 2018;114:57–62.CAS 
Article 

Google Scholar 
Jutkina J, Rutgersson C, Flach CF, Joakim Larsson DG. An assay for determining minimal concentrations of antibiotics that drive horizontal transfer of resistance. Sci Total Environ. 2016;548–549:131–8. https://doi.org/10.1016/j.scitotenv.2016.01.044.CAS 
Article 
PubMed 

Google Scholar 
Jutkina J, Marathe NP, Flach CF, Larsson DGJ. Antibiotics and common antibacterial biocides stimulate horizontal transfer of resistance at low concentrations. Sci Total Environ. 2018;616–617:172–8. https://doi.org/10.1016/j.scitotenv.2017.10.312.CAS 
Article 
PubMed 

Google Scholar 
Murray AK, Zhang L, Yin X, Zhang T, Buckling A, Snape J, et al. Novel insights into selection for antibiotic resistance in complex microbial communities. MBio. 2018;9:1–12.CAS 
Article 

Google Scholar 
Le-minh N, Khan SJ, Drewes JE, Stuetz RM. Fate of antibiotics during municipal water recycling treatment processes. Water Res. 2010;44:4295–323. https://doi.org/10.1016/j.watres.2010.06.020.CAS 
Article 
PubMed 

Google Scholar 
George J, Halami PM. Sub-inhibitory concentrations of gentamicin triggers the expression of aac(6′)Ie-aph(2″)Ia, chaperons and biofilm related genes in Lactobacillus plantarum MCC 3011. Res Microbiol. 2017;168:722–31. https://doi.org/10.1016/j.resmic.2017.06.002.CAS 
Article 
PubMed 

Google Scholar 
Zhang AN, Li LG, Ma L, Gillings MR, Tiedje JM, Zhang T. Conserved phylogenetic distribution and limited antibiotic resistance of class 1 integrons revealed by assessing the bacterial genome and plasmid collection. Microbiome. 2018;6:1–14.Article 

Google Scholar 
Gillings MR. Integrons: Past, Present, and Future. Microbiol Mol Biol Rev. 2014;78:257–77.Article 

Google Scholar 
Guironnet A, Sanchez-Cid C, Vogel TM, Wiest L, Vulliet E Aminoglycosides analysis optimization using Ion pairing Liquid Chromatography coupled to tandem Mass Spectrometry and application on wastewater samples. J Chromatogr. 2021;1651.Muyzer G, Hottentrager S, Teske A, Wawer C Denaturing gradient gel electrophoresis of PCR-amplified 16S rDNA—a new molecular approach to analyse the genetic diversity of mixed microbial communities. In: Akkermans A, van Elsas J, de Bruijn F, editors. Molecular microbial ecology manual. Dordrecht, The Netherlands: Kluwer Academic Publishers; 1995. p. 1–23.Watanabe K, Kodama Y, Harayama S. Design and evaluation of PCR primers to amplify bacterial 16S ribosomal DNA fragments used for community fingerprinting. J Microbiol Methods. 2001;44:253–62.CAS 
Article 

Google Scholar 
Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, et al. Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013;41:1–11.Article 

Google Scholar 
Masella AP, Bartram AK, Truszkowski JM, Brown DG, Neufeld JD. PANDAseq: paired-end assembler for illumina sequences. BMC Bioinformatics. 2012;13:31 http://www.biomedcentral.com/1471-2105/13/31.CAS 
Article 

Google Scholar 
Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.CAS 
Article 

Google Scholar 
Holmes AJ, Gillings MR, Nield BS, Mabbutt BC, Nevalainen KMH, Stokes HW. The gene cassette metagenome is a basic resource for bacterial genome evolution. Environ Microbiol. 2003;5:383–94.CAS 
Article 

Google Scholar 
Gillings MR, Xuejun D, Hardwick SA, Holley MP, Stokes HW. Gene cassettes encoding resistance to quaternary ammonium compounds: a role in the origin of clinical class 1 integrons? ISME J. 2009;3:209–15.CAS 
Article 

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

Google Scholar 
Minoche AE, Dohm JC, Himmelbauer H Evaluation of genomic high-throughput sequencing data generated on Illumina HiSeq and Genome Analyzer systems. Genome Biol. 2011;12.Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13:1–22.Article 

Google Scholar 
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. Available from: https://doi.org/10.1038/nmeth.1923.CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Eren AM, Esen OC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi’o: An advanced analysis and visualization platformfor’omics data. PeerJ. 2015;2015:1–29.
Google Scholar 
Alcock BP, Raphenya AR, Lau TTY, Tsang KK, Bouchard M, Edalatmand A, et al. CARD 2020: Antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res. 2020;48:D517–25.CAS 
Article 

Google Scholar 
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotechnol. 2017;35:725–31.CAS 
Article 

Google Scholar 
Menzel P, Ng KL, Krogh A Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat Commun. 2016;7.Ramirez SM, Tolmasky EM. Aminoglycoside modifing enzymes. Drug Resist Updat. 2011;13:151–71. Available from: https://doi.org/10.1016/j.drup.2010.08.003.CAS 
Article 

Google Scholar 
Ben Y, Fu C, Hu M, Liu L, Wong MH, Zheng C. Human Health Risk Assessment of Antibiotic Resistance Associated with Antibiotic Residues in the Environment: A Review. Environ Res. 2018;169:483–93. https://www.sciencedirect.com/science/article/pii/S0013935118304298.Article 

Google Scholar 
Bengtsson-Palme J, Larsson DGJ. Concentrations of antibiotics predicted to select for resistant bacteria: Proposed limits for environmental regulation. Environ Int. 2016;86:140–9. https://doi.org/10.1016/j.envint.2015.10.015.CAS 
Article 
PubMed 

Google Scholar 
Sultan I, Rahman S, Jan AT, Siddiqui MT, Mondal AH, Haq QMR Antibiotics, Resistome and Resistance Mechanisms: A Bacterial Perspective. Front Microbiol. 2018;9(September). Available from: https://www.frontiersin.org/article/10.3389/fmicb.2018.02066/fullCasin I, Bordon F, Bertin P, Coutrot A, Podglajen I, Brasseur R, et al. Aminoglycoside 6’-N-acetyltransferase variants of the Ib type with altered substrate profile in clinical isolates of Enterobacter cloacae and Citrobacter freundii. Antimicrob Agents Chemother. 1998;42:209–15.CAS 
Article 

Google Scholar 
Berendonk TU, Manaia CM, Merlin C, Fatta-Kassinos D, Cytryn E, Walsh F, et al. Tackling antibiotic resistance: the environmental framework. Nat Rev Microbiol. 2015;13:310–7.CAS 
Article 

Google Scholar 
Chow LKM, Ghaly TM, Gillings MR. A survey of sub-inhibitory concentrations of antibiotics in the environment. J Environ Sci (China). 2021;99:21–7. https://doi.org/10.1016/j.jes.2020.05.030.Article 

Google Scholar 
Gillings MR. Class 1 integrons as invasive species. Curr Opin Microbiol. 2017;38:10–5. https://doi.org/10.1016/j.mib.2017.03.002.CAS 
Article 
PubMed 

Google Scholar 
Ma L, Li AD, Yin XL, Zhang T. The prevalence of integrons as the carrier of antibiotic resistance genes in natural and man-made environments. Environ Sci Technol. 2017;51:5721–8.CAS 
Article 

Google Scholar 
Gillings M, Boucher Y, Labbate M, Holmes A, Krishnan S, Holley M, et al. The evolution of class 1 integrons and the rise of antibiotic resistance. J Bacteriol. 2008;190:5095–100.CAS 
Article 

Google Scholar 
Bürgmann H, Frigon D, Gaze WH, Manaia CM, Pruden A, Singer AC, et al. Water and sanitation: An essential battlefront in the war on antimicrobial resistance. FEMS Microbiol Ecol. 2018;94.Pena-Miller R, Laehnemann D, Jansen G, Fuentes-Hernandez A, Rosenstiel P, Schulenburg H, et al. When the most potent combination of antibiotics selects for the greatest bacterial load: the smile-frown transition. PLoS Biol. 2013;11:14–6.Article 

Google Scholar  More