Philippa Kaur

Chromatographic separation and mass spectrometric optimizationTo obtain the best monitoring conditions for each compound, a 0.5 mg/L mixed standard solution of 11 pesticides was mixed with the mobile phase through a syringe pump and then injected into the mass spectrometer for tuning. The precursor ion of the compound to be tested was determined by the primary mass spectrometry scan under ESI+ and ESI- modes, and then the product ion was scanned by the secondary mass spectrometry. Two groups of ion pairs with the best sensitivity were selected for detection; one group was used for quantification, and another, for qualitative analysis. The optimization results showed high sensitivity of all the 11 pesticides under the ESI+ mode. Among them, abamectin (B1a), β-cypermethrin, deltamethrin, fenpropathrin, and bifenthrin were [M + NH4]+, and other compounds were [M + H]+. MS parameters of 11 pesticides are mentioned in Table S2.Formic acid and ammonium acetate are commonly used reagents to enhance the ionization of target compounds [M + H]+ and [M + NH4]+ under the ESI+ mode, and they can effectively improve the peak pattern, making the peak sharper and more symmetrical; therefore, they need to be added during gradient elution38. To improve work efficiency, it is necessary to separate and complete the monitoring of 11 pesticides in the shortest possible time; therefore, we selected two different types of chromatographic columns (ACQUITY UPLC HSS C18 and ACQUITY UPLC HSS T3) and three different mobile phases (Ι: 0.1% formic acid aqueous solution—ACN, II: 0.05% formic acid aqueous solution—ACN, and III: 0.1% formic acid/5 mmol/L ammonium acetate aqueous solution—ACN) for optimization experiments. We observed that when using the HSS T3 chromatographic column, β-cypermethrin, deltamethrin, fenpropathrin, and bifenthrin did not show a good retention effect under the three mobile phase systems, and there was substantial tailing of the chromatographic peak. The shape of the chromatographic peak and sensitivity of the target compound were used as evaluation indicators. Compared with Ι and II, mobile phase III produced better sensitivity for all target compounds (Fig. 1), with sharper and more symmetrical peaks of β-cypermethrin, deltamethrin, fenpropathrin, and bifenthrin. This may be because the addition of 5 mmol/L ammonium acetate improved the retention performance of the HSS C18 chromatography columns without affecting the ionization efficiency of all target compounds. In summary, we selected the HSS C18 column for chromatographic separation and used 0.1% formic acid/5 mmol/L ammonium acetate aqueous solution—ACN as the mobile phase to further optimize the gradient elution procedure and effectively separate and detect all the target compounds within 8 min.Figure 1When using HSS C18, the peak areas of 11 pesticides in three different mobile phases.Full size imageOptimization of purification materialsThe flesh of apricot contains sugar, protein, calcium, phosphorus, carotene, thiamine, riboflavin, niacin, and vitamin C. Due to these diverse impurities, the analysis of the sample matrix becomes highly complex. Therefore, these impurities need to be removed from the matrix samples before analysis. Currently, PSA, C18, and MWCNTs are widely used to adsorb to the fruit substrate39. PSA has a strong adsorption capacity for metal ions, fatty acids, sugars, and fat-soluble pigments, C18 has a strong adsorption capacity for non-polar impurities (such as fat, sterol, and volatile oil), while MWCNTs have a strong adsorption capacity for pigments, which can effectively remove chlorophyll, lutein, and carotene. However, C18 and MWCNTs can also simultaneously adsorb pesticides, resulting in poor recovery. Nano-ZrO2 has a large specific surface area and good adsorption stability and has recently been used to purify substrates. It can selectively remove fats and pigments from samples compared to conventional C18 fillers.In the current study, different purification materials were combined for the analysis of 11 pesticide residues and to propose the best purification strategy in the pretreatment of apricot samples. As displayed in Fig. 2, the average recovery of 11 pesticides in the apricot was higher using the C18/nano-ZrO2/MWCNTs than other combinations. Nano-ZrO2 showed better adsorption than PSA in purifying fatty acids, organic acids, polar pigments, and sugars in apricot, owing to its larger specific surface area, better adsorption capacity, and stability. To conclude, the combination of 10 mg C18, 30 mg nano-ZrO2, and 5 mg MWCNTs demonstrated the best recovery for 11 pesticides, with recovery in the range of 72% to 114%, at a pesticide spiking level of 0.01 mg/kg. In summary, we finally determined that among the tested combinations, C18/nano-ZrO2/MWCNTs (10 mg/ 30 mg/5 mg) is the best purification combination for the pre-treatment of apricot samples.Figure 2The recoveries of 11 pesticides in apricot matrix under different scavenger combinations (2–1 C18/nano-ZrO2/MWCNTs, 2–2 PSA/C18/MWCNTs, 2–3 nano-ZrO2/PSA/MWCNTs; 0.01 mg/kg, n = 3).Full size imageLinearity, matrix effects, limit of detection and limit of quantificationThe standard curve obtained from the standard working solutions of 11 pesticides and the calibration curve from blank apricot matrix spiked with 11 pesticides showed good linearity (0.001, 0.005, 0.01, 0.05, 0.1, and 0.5 mg/L), with R2 ≥ 0.9959 for all tested samples (Table 1).Table 1 The standard curves, R2 and MEs of 11 pesticides in apricot.Full size tableTo evaluate MEs, the slopes of matching 11 pesticide standards with solvent and apricot matrix were calculated at the same concentration. According to the derived slope of the matrix-matched calibration curve, MEs of 11 pesticides in apricot were between 89 and 113% (Table 1), well within the range of 80% to 120%, indicating that the MEs could be ignored. It also suggests that the current pre-treatment method has a good purification effect and eliminates the matrix effect very well, laying a robust foundation for the subsequent step of quantitative analysis of samples. We next used the standard solution curve to quantify the 11 pesticide residues in apricot.The LOD refers to the minimum concentration or minimum amount of a component to be tested that can be detected from a test sample under a given confidence level by an analytical method. Its physical meaning is the amount of the measured component when the signal is 3 times the standard deviation (S = 3σ) of the reagent blank signal (background signal). Sometimes it also refers to the amount of the measured component corresponding to when the signal is three times the background signal generated by the reagent blank (S = 3 N). The LOQ refers to the minimum amount of the analyte in the sample that can be quantitatively determined, and the determination result should have a certain accuracy40. The LOQ reflects whether the analytical method has the sensitive quantitative detection ability. The LOQ is the lowest validated level with sufficient recovery and precision, which was estimated to be 0.001 mg/L, while the LOD is the lowest calibration level, which was 2 µg/kg, according to SANTE/12,682/2020.Accuracy and precisionIn the matrix, 11 pesticides were spiked at four levels (0.002, 0.02, 0.1, and 1 mg/kg), and for each spiked sample, there were six replicates. The recoveries of 11 pesticides in apricot at all levels ranged between 72 and 119%. The inter- and intra-level relative standard deviations (RSDs, %) of 11 pesticides in apricot were  More

Lim, C. & Putnam, R. D. Religion, social networks, and life satisfaction. Am. Sociol. Rev. 75, 914–933 (2010).
Google Scholar 
Fox, R. Kinship and marriage: an anthropological perspective/by Robin Fox. (1967).Lévi-Strauss, C. The elementary structures of kinship. (Beacon Press, 1969).Murdock, G. P. Social structure. Macmillan 387 (1949).Chapais, B. Primeval kinship: how pair-bonding gave birth to human society. (Harvard University Press, 2009).Walker, R. S. & Hill, K. R. Causes, consequences, and kin bias of human group fissions. Hum. Nat. 25, 465–475 (2014).PubMed 

Google Scholar 
Shenk, M. K., Towner, M. C., Voss, E. A. & Alam, N. Consanguineous marriage, kinship ecology, and market transition. Curr. Anthropol. 57, S167–S180 (2016).
Google Scholar 
Swann, W. B. Jr., Gómez, A., Seyle, D. C., Morales, J. F. & Huici, C. Identity fusion: The interplay of personal and social identities in extreme group behavior. J. Pers. Soc. Psychol. 96, 995–1011 (2009).PubMed 

Google Scholar 
Richerson, P. J. & Boyd, R. Complex societies. Hum. Nat. 10, 253–289 (1999).CAS 
PubMed 

Google Scholar 
Zelinsky, W. The hypothesis of the mobility transition. Geogr. Rev. 61, 219–249 (1971).
Google Scholar 
Gurven, M., Jaeggi, A. V., von Rueden, C., Hooper, P. L. & Kaplan, H. Does market integration buffer risk, erode traditional sharing practices and increase inequality? A test among Bolivian forager-farmers. Hum. Ecol. Interdiscip. J. 43, 515–530 (2015).PubMed 
PubMed Central 

Google Scholar 
Godoy, R. A. et al. Do markets worsen economic inequalities? Kuznets in the Bush. Hum. Ecol. 32, 339–364 (2004).
Google Scholar 
Kaplan, H. A theory of fertility and parental investment in traditional and modern human societies. Am. J. Phys. Anthropol. 101, 91–135 (1996).
Google Scholar 
Duernecker, G. & Vega-Redondo, F. Social Networks and the Process of Globalization. Rev. Econ. Stud. 85, 1716–1751 (2017).MathSciNet 
MATH 

Google Scholar 
Colleran, H. Market integration reduces kin density in women’s ego-networks in rural Poland. Nat. Commun. 11, 266 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Wilding, R. Families, intimacy and globalization. (Macmillan International Higher Education, 2018).Hackman, J. V. & Kramer, K. L. Kin Ties and market integration in a Yucatec Mayan Village. Soc. Sci. 10, 216 (2021).
Google Scholar 
Norenzayan, A. Big gods: How religion transformed cooperation and conflict. (Princeton University Press, 2013).Lauder, W., Mummery, K. & Sharkey, S. Social capital, age and religiosity in people who are lonely. J. Clin. Nurs. 15, 334–340 (2006).PubMed 

Google Scholar 
Agate, S. T., Zabriskie, R. B. & Eggett, D. L. Praying, playing, and successful families. Marriage Fam. Rev. 42, 51–75 (2007).
Google Scholar 
Day, R. D. et al. Family processes and adolescent religiosity and religious practice: View from the NLSY97. Marriage Fam. Rev. 45, 289–309 (2009).
Google Scholar 
Fagan, P. F. Why religion matters even more: The impact of religious practice on social stability. Backgrounder 1992, 1–19 (2006).
Google Scholar 
Ellison, C. G. & George, L. K. Religious involvement, social ties, and social support in a Southeastern Community. J. Sci. Study Relig. 33, 46–61 (1994).
Google Scholar 
Ellison, C. G. & Xu, X. Religion and families. The Wiley Blackwell companion to the sociology of families 277–299 (2014).Ginges, J., Hansen, I. & Norenzayan, A. Religion and support for suicide attacks. Psychol. Sci. 20, 224–230 (2009).PubMed 

Google Scholar 
Lynch, R., Palestis, B. G. & Trivers, R. Religious devotion and extrinsic religiosity affect in-group altruism and out-group hostility oppositely in rural Jamaica. Evol. Psychol. Sci. 3, 335 (2017).
Google Scholar 
Walker, R. S. & Bailey, D. H. Marrying kin in small-scale societies. Am. J. Hum. Biol. 26, 384–388 (2014).PubMed 

Google Scholar 
Putnam, R. D., Leonardi, R. & Nanetti, R. Y. Making Democracy Work: Civic Traditions in Modern Italy. (Princeton University Press, 1994).Coleman, J. Foundations of Social Theory. (Belknap Press of Harvard University Press, Cambridge, Mass, 1990).Wuthnow, R. The Left Behind: Decline and Rage in Rural America. (Princeton University Press, 2018).Sunstein, C. R. # Republic: Divided democracy in the age of social media. (Princeton University Press, 2018).Putnam, R. D. E Pluribus Unum: Diversity and Community in the Twenty-first Century The 2006 Johan Skytte Prize Lecture. Scan. Polit. Stud. 30, (2007).Putnam, R. Bowling alone: The collapse and revival of American community. (Simon and Schuster, 2000).Olson, M. The Logic of Collective Action: Public Goods and the Theory of Groups, Second printing with new preface and appendix (Harvard Economic Studies). Harvard economic studies, v. 124 (Harvard University Press, 1971).Granovetter, M. S. The strength of weak ties. Am. J. Sociol. (1973).Lynch, R., Lummaa, V. & Panchanathan, K. Integration involves a trade-off between fertility and status for World War II evacuees. Nature Human Behaviour (2019).Beyerlein, K. & Hipp, J. R. Social capital, too much of a good thing? American Religious Traditions and Community Crime. Soc. Forces 84, 995–1013 (2005).
Google Scholar 
Lewis, V. A., Macgregor, C. A. & Putnam, R. D. Religion, networks, and neighborliness: The impact of religious social networks on civic engagement. Soc. Sci. Res. 42, 331–346 (2013).PubMed 

Google Scholar 
Yu, M. & Stiffman, A. R. Positive family relationships and religious affiliation as mediators between negative environment and illicit drug symptoms in American Indian adolescents. Addict. Behav. 35, 694–699 (2010).PubMed 

Google Scholar 
Regnerus, M. D. & Burdette, A. Religious change and adolescent family dynamics. Sociol. Q. 47, 175–194 (2006).
Google Scholar 
Marks, L. Religion and family relational health: An overview and conceptual model. J. Relig. Health (2006).Thornton, A. Reciprocal Influences of Family and Religion in a Changing World. J. Marriage Fam. Couns. 47, 381–394 (1985).
Google Scholar 
Mahoney, A., Pargament, K. I., Murray-Swank, A. & Murray-Swank, N. Religion and the Sanctification of Family Relationships. Rev. Relig. Res. 44, 220–236 (2003).
Google Scholar 
Mahoney, A. Religion in families 1999 to 2009: A relational spirituality framework. J. Marriage Fam. 72, 805–827 (2010).PubMed 
PubMed Central 

Google Scholar 
Ebstyne King, P. & Furrow, J. L. Religion as a resource for positive youth development: religion, social capital, and moral outcomes. Dev. Psychol. 40, 703–713 (2004).PubMed 

Google Scholar 
Dudley, M. G. & Kosinski, F. A. Religiosity and marital satisfaction: A research note. Rev. Relig. Res. 32, 78–86 (1990).
Google Scholar 
Milevsky, A., Smoot, K., Leh, M. & Ruppe, A. Familial and contextual variables and the nature of sibling relationships in emerging adulthood. Marriage Fam. Rev. 37, 123–141 (2005).
Google Scholar 
Galbraith, D. & Shaver, J. H. Religion and Fertility Bibliography. evolutionarydemographyofreligion.Shaver, J. H., Sibley, C. G., Sosis, R., Galbraith, D. & Bulbulia, J. Alloparenting and religious fertility: A test of the religious alloparenting hypothesis. Evol. Hum. Behav. 40, 315–324 (2019).
Google Scholar 
Kaufmann, E. Shall the Religious Inherit the Earth?: Demography and Politics in the Twenty-First Century. (Profile Books, 2010).Ebaugh, H. R. & Curry, M. Fictive Kin as social capital in new immigrant communities. Sociol. Perspect. 43, 189–209 (2000).
Google Scholar 
Taylor, R. J., Chatters, L. M., Woodward, A. T. & Brown, E. Racial and ethnic differences in extended family, friendship, fictive kin and congregational informal support networks. Fam. Relat. 62, 609–624 (2013).PubMed 
PubMed Central 

Google Scholar 
Durkheim, E. The elementary forms of the religious life. Preprint at (1915).Rappaport, R. A. Ritual and Religion in the Making of Humanity. vol. 110 (Cambridge University Press, 1999).Hastings, O. P. Not a lonely crowd? Social connectedness, religious service attendance, and the spiritual but not religious. Soc. Sci. Res. 57, 63–79 (2016).PubMed 

Google Scholar 
Putnam, R. & Campbell, D. E. American grace: How religion is reshaping our civic and political lives. Preprint at (2010).Turke, P. W. Evolution and the demand for children. Popul. Dev. Rev. 15, 61–90 (1989).
Google Scholar 
Sear, R. & Coall, D. How much does family matter? Cooperative breeding and the demographic transition. Popul. Dev. Rev. 37, 81–112 (2011).PubMed 

Google Scholar 
Jenkins, P. Fertility and Faith: The Demographic Revolution and the Transformation of World Religions. (Baylor University Press, 2020).Rothstein, B. Corruption and social trust: Why the fish rots from the head down. Soc. Res. 80, 1009–1032 (2013).
Google Scholar 
Lynch, R, Schaffnit, S. and Shenk, M. OSF preregistration – Does religion help to preserve the density of kin networks often disrupted by globalization? Open Science Framework Registries. https://osf.io/xvyqm/registrations (2020).Alam, N. et al. Health and demographic surveillance system (HDSS) in Matlab, Bangladesh. Int. J. Epidemiol. 46, 809–816 (2017).PubMed 

Google Scholar 
Icddr, B. Health and Demographic Surveillance System-Matlab. 2005 Socioeconomic Census (2007).Imf. International Monetary Fund. World Economic Outlook Database. (2016).Razzaque, A., Streatfield, P. K. & Evans, A. Family size and children’s education in Matlab, Bangladesh. J. Biosoc. Sci. 39, 245–256 (2007).PubMed 

Google Scholar 
Afsar, R. Unravelling the vicious cycle of recruitment: Labour migration from Bangladesh to the gulf states. http://ilo.org/wcmsp5/groups/public/—ed_norm/—declaration/documents/publication/wcms_106536.pdf (2009).Kabeer, N. Ideas, economics and ‘the sociology of supply’: Explanations for fertility decline in Bangladesh. J. Dev. Stud. 38, 29–70 (2001).
Google Scholar 
Novak, J. J. Bangladesh: Reflections on the water. (Indiana University Press, 1993).Shenk, M. K., Towner, M. C., Kress, H. C. & Alam, N. A model comparison approach shows stronger support for economic models of fertility decline. Proc. Natl. Acad. Sci. USA 110, 8045–8050 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Devine, J., Hinks, T. & Naveed, A. Happiness in Bangladesh: The role of religion and connectedness. J. Happiness Stud. 20, 351–371 (2019).
Google Scholar 
Henrich, J. Market incorporation, agricultural change, and sustainability among the Machiguenga Indians of the Peruvian Amazon. Hum. Ecol. 25, 319–351 (1997).
Google Scholar 
Lu, F. Integration into the market among indigenous peoples: A cross-cultural perspective from the Ecuadorian Amazon. Curr. Anthropol. 48, 593–602 (2007).
Google Scholar 
Bürkner, P.-C. Advanced Bayesian Multilevel Modeling with the R Package brms. arXiv [stat.CO] (2017).Team, R. C. & Others. R: A language and environment for statistical computing. (2013).Lynch, R. Kin_density_and-religiosity. (2021).McElreath, R. Statistical rethinking. (2017).Clarke, M. New kinship, Islam, and the liberal tradition: sexual morality and new reproductive technology in Lebanon. J. R. Anthropol. Inst. 14, 153–169 (2008).
Google Scholar 
Swann, W. B. et al. What makes a group worth dying for? Identity fusion fosters perception of familial ties, promoting self-sacrifice. J. Pers. Soc. Psychol. 106, 912–926 (2014).PubMed 

Google Scholar 
Benítez, D. M. Bangladesh: Economy Overview and Structural Changes. (2018).Viry, G. Residential mobility and the spatial dispersion of personal networks: Effects on social support. Soc. Networks 34, 59–72 (2012).
Google Scholar 
Mok, D., Wellman, B. & Carrasco, J. Does distance matter in the age of the internet?. Urban Stud. 47, 2747–2783 (2010).
Google Scholar 
Rivera, M. T., Soderstrom, S. B. & Uzzi, B. Dynamics of dyads in social networks: Assortative, relational, and proximity mechanisms. Annu. Rev. Sociol. 36, 91–115 (2010).
Google Scholar 
Pollet, T. V., Roberts, S. G. B. & Dunbar, R. I. M. Going that extra mile: Individuals travel further to maintain face-to-face contact with highly related kin than with less related kin. PLoS ONE 8, e53929 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Madhavan, S., Clark, S., Araos, M. & Beguy, D. Distance or location? How the geographic distribution of kin networks shapes support given to single mothers in urban Kenya. Geogr. J. 184, 75–88 (2018).
Google Scholar 
Curry, O., Roberts, S. G. B. & Dunbar, R. I. M. Altruism in social networks: evidence for a ‘kinship premium’. Br. J. Psychol. 104, 283–295 (2013).PubMed 

Google Scholar 
Sullivan, K. & Sullivan, A. Adolescent–parent separation. Dev. Psychol. 16, 93 (1980).
Google Scholar 
Roberts, S. G. B. & Dunbar, R. I. M. Communication in social networks: Effects of kinship, network size, and emotional closeness. Pers. Relatsh. 18, 439–452 (2011).
Google Scholar 
Shenk, M. K. et al. Social support, nutrition and health among women in rural Bangladesh: complex tradeoffs in allocare, kin proximity and support network size. Philos. Trans. R. Soc. Lond. B Biol. Sci. 376, 207 (2021).Snopkowski, K. & Sear, R. Grandparental help in Indonesia is directed preferentially towards needier descendants: A potential confounder when exploring grandparental influences on child health. Soc. Sci. Med. 128, 105–114 (2015).PubMed 

Google Scholar 
Schaffnit, S. B. & Sear, R. Support for new mothers and fertility in the United Kingdom: Not all support is equal in the decision to have a second child. Popul. Stud. 71, 345–361 (2017).
Google Scholar 
Boyer, P. The Naturalness of Religious Ideas: A Cognitive Theory of Religion. (University of California Press, 1994).Thomas, M. G. et al. Kinship underlies costly cooperation in Mosuo villages. R Soc Open Sci 5, 171535 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
Maqsood, A. Love as understanding. Am. Ethnol. https://doi.org/10.1111/amet.13000 (2021).Article 

Google Scholar 
Schurmann, A. T. & Mahmud, S. Civil society, health, and social exclusion in Bangladesh. J. Health Popul. Nutr. 27, 536–544 (2009).PubMed 
PubMed Central 

Google Scholar 
Haque, M. R., Hasan, M. S., Alam, N., Barkat, S. & Others. Fertility preferences in Bangladesh. in Family Demography in Asia (Edward Elgar Publishing, 2018).Mattison, S. M. Economic impacts of tourism and erosion of the visiting system among the Mosuo of Lugu Lake. Asia Pac. J. Anthropol. 11, 159–176 (2010).
Google Scholar 
Mattison, S. M. et al. Context specificity of ‘market integration’ among the matrilineal Mosuo of Southwest China. Curr. Anthropol. 63, 118–124 (2022).
Google Scholar 
Uchida, Y., Kitayama, S., Mesquita, B., Reyes, J. A. S. & Morling, B. Is perceived emotional support beneficial? Well-being and health in independent and interdependent cultures. Pers. Soc. Psychol. Bull. 34, 741–754 (2008).PubMed 

Google Scholar 
Reblin, M. & Uchino, B. N. Social and emotional support and its implication for health. Curr. Opin. Psychiatry 21, 201–205 (2008).PubMed 
PubMed Central 

Google Scholar 
Inglehart, R. Faith and freedom: Traditional and modern ways to happiness. Int. Differ. Well-being 351, 397 (2010).
Google Scholar 
Ferriss, A. L. Religion and the Quality of Life. J. Happiness Stud. 3, 199–215 (2002).
Google Scholar 
Greeley, A. & Hout, M. Happiness and lifestyle among conservative Christians. The truth about conservative Christians 1, 150–161 (2006).
Google Scholar 
Pilisuk, M. Kinship, social networks, social support and health. Soc. Sci. Med. 12, 273–280 (1978).CAS 
PubMed 

Google Scholar 
Schaffnit, S. B. & Sear, R. Supportive families versus support from families: The decision to have a child in the Netherlands. Demogr. Res. 37, 417–454 (2017).
Google Scholar 
Hassan, A., Lawson, D., Schaffnit, S. B., Urassa, M. & Sear, R. Childcare in transition: evidence that patterns of childcare differ by degree of market integration in north-western Tanzania. (2021).https://doi.org/10.31219/osf.io/gtc6kPutnam, R. D. Democracies in Flux: The Evolution of Social Capital in Contemporary Society. (Oxford University Press, 2004). More

Gibney, E. R. & Nolan, C. M. Epigenetics and gene expression. Heredity 105, 4–13 (2010).CAS 
PubMed 

Google Scholar 
Vogt, G. Facilitation of environmental adaptation and evolution by epigenetic phenotype variation: Insights from clonal, invasive, polyploid, and domesticated animals. Environ. Epigenet. 3, 1–17 (2017).
Google Scholar 
Beal, A., Rodriguez-Casariego, J., Rivera-Casas, C., Suarez-Ulloa, V. & Eirin-Lopez, J. M. Environmental Epigenomics and Its Applications in Marine Organisms 325–359 (Springer, 2018). https://doi.org/10.1007/13836_2018_28.Book 

Google Scholar 
Hofmann, G. E. Ecological epigenetics in marine metazoans. Front. Mar. Sci. 4, 1–7 (2017).CAS 

Google Scholar 
Richards, C. L. et al. Ecological plant epigenetics: Evidence from model and non-model species, and the way forward. Ecol. Lett. 20, 1576–1590 (2017).PubMed 

Google Scholar 
Ryu, T., Veilleux, H. D., Donelson, J. M., Munday, P. L. & Ravasi, T. The epigenetic landscape of transgenerational acclimation to ocean warming. Nat. Clim. Chang. 8, 504–509 (2018).ADS 

Google Scholar 
Liew, Y. J. et al. Epigenome-associated phenotypic acclimatization to ocean acidification in a reef-building coral. Sci. Adv. 4, 6 (2018).
Google Scholar 
Anastasiadi, D., Díaz, N. & Piferrer, F. Small ocean temperature increases elicit stage-dependent changes in DNA methylation and gene expression in a fish, the European sea bass. Sci. Rep. 7, 1–12 (2017).CAS 

Google Scholar 
Strader, M. E., Wong, J. M., Kozal, L. C., Leach, T. S. & Hofmann, G. E. Parental environments alter DNA methylation in offspring of the purple sea urchin, Strongylocentrotus purpuratus. J. Exp. Mar. Bio. Ecol. 517, 54–64 (2019).
Google Scholar 
Rey, O. et al. Linking epigenetics and biological conservation: Towards a conservation epigenetics perspective. Funct. Ecol. 34, 414–427 (2020).
Google Scholar 
Eirin-Lopez, J. M. & Putnam, H. Editorial: Marine environmental epigenetics. Front. Mar. Sci. 8, 5 (2021).
Google Scholar 
Herrera, C. M. & Bazaga, P. Untangling individual variation in natural populations: Ecological, genetic and epigenetic correlates of longterm inequality in herbivory. Mol. Ecol. 20, 1675–1688 (2011).CAS 
PubMed 

Google Scholar 
Varriale, A. DNA methylation, epigenetics, and evolution in vertebrates: Facts and challenges. Int. J. Evol. Biol. 2014, 1–7 (2014).
Google Scholar 
Liebl, A. L., Wesner, J. S., Russell, A. F. & Schrey, A. W. Methylation patterns at fledging predict delayed dispersal in a cooperatively breeding bird. PLoS ONE 16, e0252227 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Metzger, D. C. H. & Schulte, P. M. Persistent and plastic effects of temperature on DNA methylation across the genome of threespine stickleback (Gasterosteus aculeatus). Proc. R. Soc. B Biol. Sci. 284, 5 (2017).
Google Scholar 
Putnam, H. M., Davidson, J. M. & Gates, R. D. Ocean acidification influences host DNA methylation and phenotypic plasticity in environmentally susceptible corals. Evol. Appl. 9, 1165–1178 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Watson, R. G. A., Baldanzi, S., Pérez-Figueroa, A., Gouws, G. & Porri, F. Morphological and epigenetic variation in mussels from contrasting environments. Mar. Biol. 165, 8 (2018).
Google Scholar 
Baldanzi, S., Watson, R., McQuaid, C. D., Gouws, G. & Porri, F. Epigenetic variation among natural populations of the South African sandhopper Talorchestia capensis. Evol. Ecol. 31, 77–91 (2017).
Google Scholar 
Ardura, A., Zaiko, A., Morán, P., Planes, S. & Garcia-Vazquez, E. Epigenetic signatures of invasive status in populations of marine invertebrates. Sci. Rep. 7, 5 (2017).
Google Scholar 
Baldanzi, S., Storch, D., Navarrete, S. A., Graeve, M. & Fernández, M. Latitudinal variation in maternal investment traits of the kelp crab Taliepus dentatus along the coast of Chile. Mar. Biol. 165, 1 (2018).
Google Scholar 
Sobarzo, M., Bravo, L., Donoso, D., Garcés-Vargas, J. & Schneider, W. Coastal upwelling and seasonal cycles that influence the water column over the continental shelf off central Chile. Prog. Oceanogr. 75, 363–382 (2007).ADS 

Google Scholar 
Letelier, J., Pizarro, O. & Nuñez, S. Seasonal variability of coastal upwelling and the upwelling front off central Chile. J. Geophys. Res. Ocean. 114, 12009 (2009).ADS 

Google Scholar 
Vargas, C. A. et al. Species-specific responses to ocean acidification should account for local adaptation and adaptive plasticity. Nat. Ecol. Evol. 1, 1–7 (2017).CAS 

Google Scholar 
Pérez, C. A. et al. Influence of climate and land use in carbon biogeochemistry in lower reaches of rivers in central southern Chile: Implications for the carbonate system in river-influenced rocky shore environments. J. Geophys. Res. Biogeosciences 120, 673–692 (2015).ADS 

Google Scholar 
Saldías, G. S. et al. Satellite-measured interannual variability of turbid river plumes off central-southern Chile: Spatial patterns and the influence of climate variability. Prog. Oceanogr. 146, 212–222 (2016).ADS 

Google Scholar 
Lara, C. et al. Coastal biophysical processes and the biogeography of rocky intertidal species along the south-eastern Pacific. J. Biogeogr. 46, 420–431 (2019).
Google Scholar 
Wieters, E. A. Upwelling control of positive interactions over mesoscales: A new link between bottom-up and top-down processes on rocky shores. Mar. Ecol. Prog. Ser. 301, 43–54 (2005).ADS 

Google Scholar 
Pérez-Matus, A., Carrasco, S. A., Gelcich, S., Fernandez, M. & Wieters, E. A. Exploring the effects of fishing pressure and upwelling intensity over subtidal kelp forest communities in Central Chile. Ecosphere 8, e01808 (2017).
Google Scholar 
Iranon, N. N. & Miller, D. L. Interactions between oxygen homeostasis, food availability, and hydrogen sulfide signaling. Front. Genet. 3, 5 (2012).
Google Scholar 
Ramajo, L., Lagos, N. A. & Duarte, C. M. Seagrass Posidonia oceanica diel pH fluctuations reduce the mortality of epiphytic forams under experimental ocean acidification. Mar. Pollut. Bull. 146, 247–254 (2019).CAS 
PubMed 

Google Scholar 
Aiken, C. & Navarrete, S. Environmental fluctuations and asymmetrical ­dispersal: Generalized stability theory for studying metapopulation persistence and marine protected areas. Mar. Ecol. Prog. Ser. 428, 77–88 (2011).ADS 

Google Scholar 
Baldanzi, S. et al. Combined effects of temperature and hypoxia shape female brooding behaviors and the early ontogeny of the Chilean kelp crab Taliepus dentatus. Mar. Ecol. Prog. Ser. 646, 93–107 (2020).ADS 
CAS 

Google Scholar 
Moran, A. L. & McAlister, J. S. Egg size as a life history character of marine invertebrates: Is it all it’s cracked up to be?. Biol. Bull. 216, 226–242 (2009).PubMed 

Google Scholar 
Doherty-Weason, D. et al. Bioenergetics of parental investment in two polychaete species with contrasting reproductive strategies: The planktotrophic Boccardia chilensis and the poecilogonic Boccardia wellingtonensis (Spionidae). Mar. Ecol. 41, 1 (2020).
Google Scholar 
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300 (1995).MathSciNet 
MATH 

Google Scholar 
Steneck, R. S. et al. Kelp forest ecosystems: Biodiversity, stability, resilience and future. Environ. Conserv. 29, 436–459 (2002).
Google Scholar 
Sayols-Baixeras, S., Irvin, M. R., Arnett, D. K., Elosua, R. & Aslibekyan, S. W. Epigenetics of lipid phenotypes. Curr. Cardiovasc. Risk Rep. 10, 1–205 (2016).
Google Scholar 
Adam, A. C. et al. Profiling DNA methylation patterns of zebrafish liver associated with parental high dietary arachidonic acid. PLoS ONE 14, e0220934 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
García-Fernández, P., García-Souto, D., Almansa, E., Morán, P. & Gestal, C. Epigenetic DNA methylation mediating Octopus vulgaris early development: Effect of essential fatty acids enriched diet. Front. Physiol. 8, 1–9 (2017).
Google Scholar 
Hearn, J., Pearson, M., Blaxter, M., Wilson, P. J. & Little, T. J. Genome-wide methylation is modified by caloric restriction in Daphnia magna. BMC Genomics 20, 1–11 (2019).
Google Scholar 
Palma, A. T., Henríquez, L. A. & Ojeda, F. P. Phytoplanktonic primary production modulated by coastal geomorphology in a highly dynamic environment of central Chile. Rev. Biol. Mar. Oceanogr. 44, 325–334 (2009).
Google Scholar 
Faúndez-Báez, P., Morales, C. E. & Arcos, D. Variabilidad espacial y temporal en la hidrografía invernal del sistema de bahías frente a la VIII región (Chile centro-sur). Rev. Chil. Hist. Nat. 74, 817–831 (2001).
Google Scholar 
Osma, N. et al. Response of phytoplankton assemblages from naturally acidic coastal ecosystems to elevated pCO2. Front. Mar. Sci. 1, 323 (2020).
Google Scholar 
Rebolledo, L. et al. Siliceous productivity changes in Gulf of Ancud sediments (42°S, 72°W), southern Chile, over the last ∼150 years. Cont. Shelf Res. 31, 356–365 (2011).ADS 

Google Scholar 
Sun, Y. et al. Genome-wide analysis of DNA methylation in five tissues of Zhikong Scallop, Chlamys farreri. PLoS ONE 9, e86232 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Bernhardt, J. R., O’Connor, M. I., Sunday, J. M. & Gonzalez, A. Life in fluctuating environments. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 375, 20190454 (2020).PubMed 
PubMed Central 

Google Scholar 
Feinberg, A. P. & Irizarry, R. A. Colloquium Paper: Stochastic epigenetic variation as a driving force of development, evolutionary adaptation, and disease. Proc. Natl. Acad. Sci. USA 107, 1757 (2010).ADS 
CAS 
PubMed 

Google Scholar 
Tapia, F. J., Largier, J. L., Castillo, M., Wieters, E. A. & Navarrete, S. A. Latitudinal discontinuity in thermal conditions along the nearshore of Central-Northern Chile. PLoS ONE 9, e110841 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Reyna-López, G. E., Simpson, J. & Ruiz-Herrera, J. Differences in DNA methylation patterns are detectable during the dimorphic transition of fungi by amplification of restriction polymorphisms. Mol. Gen. Genet. 253, 703–710 (1997).PubMed 

Google Scholar 
Pérez-Figueroa, A. msap: A tool for the statistical analysis of methylation-sensitive amplified polymorphism data. Mol. Ecol. Resour. 13, 522–527 (2013).PubMed 

Google Scholar 
Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32, 2847–2849 (2016).CAS 
PubMed 

Google Scholar 
Valladares, F., Sanches-Gomez, D. & Zavala, M. A. Quantitative estimation of phenotypic plasticity: Bridging the gap between the evolutionary concept and its ecological applications. J. Ecol. 94, 1103–1116 (2006).
Google Scholar 
Excoffier, L., Smouse, P. E. & Quattro, J. M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131, 479–491 (1992).CAS 
PubMed 
PubMed Central 

Google Scholar  More

Trisos, C. H., Merow, C. & Pigot, A. L. The projected timing of abrupt ecological disruption from climate change. Nature 580, 496–501 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Schwörer, C. et al. Holocene climate, fire and vegetation dynamics at the treeline in the Northwestern Swiss Alps. Veg. Hist. Archaeobot. 23, 479–496 (2014).Article 

Google Scholar 
Steinbauer, M. J. et al. Accelerated increase in plant species richness on mountain summits is linked to warming. Nature 556, 231–234 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Grabherr, G., Gottfried, M. & Pauli, H. Climate effects on mountain plants. Nature 369, 448 (1994).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Bennett, K. D. & Willis, K. J. Pollen. Tracking Environmental Change Using Lake Sediments (eds Smol, J. P., Birks, H. J. B., Last, W. M., Bradley, R. S. & Alverson, K.) 5–32 (Kluwer Academic Publishers, 2002).Liu, S. et al. Sedimentary ancient DNA reveals a threat of warming-induced alpine habitat loss to Tibetan Plateau plant diversity. Nat. Commun. 12, 2995 (2021).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Rijal, D. P. et al. Sedimentary ancient DNA shows terrestrial plant richness continuously increased over the Holocene in northern Fennoscandia. Sci. Adv. 7, eabf9557 (2021).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Giguet-Covex, C. et al. Long livestock farming history and human landscape shaping revealed by lake sediment DNA. Nat. Commun. 5, 3211 (2014).Article 
ADS 
PubMed 

Google Scholar 
Väre, H., Lampinen, R., Humphries, C. & Williams, P. Taxonomic diversity of vascular plants in the European alpine areas. in Alpine biodiversity in Europe (eds Nagy, L., Grabherr, G., Körner, C. & Thompson, D. B. A.) 133–148 (Springer Berlin Heidelberg, 2003).Theurillat, J.-P. & Guisan, A. Potential impact of climate change on vegetation in the European alps: A Review. Climatic Change 50, 77–109 (2001).Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 405, 907–913 (2000).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Tribsch, A. & Schönswetter, P. Patterns of endemism and comparative phylogeography confirm palaeo-environmental evidence for Pleistocene refugia in the Eastern Alps. Taxon 52, 477–497 (2003).Article 

Google Scholar 
Rudmann-Maurer, K., Weyand, A., Fischer, M. & Stöcklin, J. The role of landuse and natural determinants for grassland vegetation composition in the Swiss Alps. Basic Appl. Ecol. 9, 494–503 (2008).Article 

Google Scholar 
Walsh, K. et al. A historical ecology of the Ecrins (Southern French Alps): Archaeology and palaeoecology of the Mesolithic to the Medieval period. Quat. Int. 353, 52–73 (2014).Article 

Google Scholar 
Walsh, K. & Giguet-Covex, C. Encyclopedia of the World’s Biomes 555–573 (Elsevier, 2020).Schwörer, C., Henne, P. D. & Tinner, W. A model-data comparison of Holocene timberline changes in the Swiss Alps reveals past and future drivers of mountain forest dynamics. Glob. Chang. Biol. 20, 1512–1526 (2014).Article 
ADS 
PubMed 

Google Scholar 
Henne, P. D. et al. An empirical perspective for understanding climate change impacts in Switzerland. Reg. Environ. Change 18, 1–17 (2017).
Google Scholar 
Niedrist, G., Tasser, E., Lüth, C., Dalla Via, J. & Tappeiner, U. Plant diversity declines with recent land use changes in European Alps. Plant Ecol. 202, 195–210 (2009).Article 

Google Scholar 
Lasanta-Martínez, T., Vicente-Serrano, S. M. & Cuadrat-Prats, J. M. Mountain Mediterranean landscape evolution caused by the abandonment of traditional primary activities: A study of the Spanish Central Pyrenees. Appl. Geogr. 25, 47–65 (2005).Article 

Google Scholar 
Nautiyal, S. & Kaechele, H. Adverse impacts of pasture abandonment in Himalayan protected areas: Testing the efficiency of a Natural Resource Management Plan (NRMP). Environ. Impact Assess. Rev. 27, 109–125 (2007).Article 

Google Scholar 
Karger, D. N., Nobis, M. P. & Normand, S. CHELSA-TraCE21k v1. 0. Downscaled transient temperature and precipitation data since the last glacial maximum. Climate of the Past (2021).Landolt, E. et al. Flora indicativa: Okologische Zeigerwerte und biologische Kennzeichen zur Flora der Schweiz und der Alpen (Haupt, 2010).Heiri, O., Brooks, S. J., Birks, H. J. B. & Lotter, A. F. A 274-lake calibration data-set and inference model for chironomid-based summer air temperature reconstruction in Europe. Quat. Sci. Rev. 30, 3445–3456 (2011).Article 
ADS 

Google Scholar 
Heiri, O., Ilyashuk, B., Millet, L., Samartin, S. & Lotter, A. F. Stacking of discontinuous regional palaeoclimate records: Chironomid-based summer temperatures from the Alpine region. Holocene 25, 137–149 (2015).Article 
ADS 

Google Scholar 
Ivy-Ochs, S. et al. Latest Pleistocene and Holocene glacier variations in the European Alps. Quat. Sci. Rev. 28, 2137–2149 (2009).Article 
ADS 

Google Scholar 
Finsinger, W. & Tinner, W. Pollen and plant macrofossils at Lac de Fully (2135 m a.s.l.): Holocene forest dynamics on a highland plateau in the Valais, Switzerland. Holocene 17, 1119–1127 (2007).Article 
ADS 

Google Scholar 
Baroni, C. et al. Last Lateglacial glacier advance in the Gran Paradiso Group reveals relatively drier climatic conditions established in the Western Alps since at least the Younger Dryas. Quat. Sci. Rev. 255, 106815 (2021).Article 

Google Scholar 
Schibler, J., Elsner, J. & Schlumbaum, A. Incorporation of aurochs into a cattle herd in Neolithic Europe: Single event or breeding? Sci. Rep. 4, 5798 (2014).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Schimmelpfennig, I. et al. A chronology of Holocene and Little Ice Age glacier culminations of the Steingletscher, Central Alps, Switzerland, based on high-sensitivity beryllium-10 moraine dating. Earth Planet. Sci. Lett. 393, 220–230 (2014).Article 
ADS 
CAS 

Google Scholar 
Ilyashuk, E. A., Heiri, O., Ilyashuk, B. P., Koinig, K. A. & Psenner, R. The Little Ice Age signature in a 700-year high-resolution chironomid record of summer temperatures in the Central Eastern Alps. Clim. Dyn. 52, 1–15 (2018).
Google Scholar 
Willerslev, E. et al. Fifty thousand years of Arctic vegetation and megafaunal diet. Nature 506, 47–51 (2014).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Alsos, I. G. et al. Ancient sedimentary DNA shows rapid post-glacial colonisation of Iceland followed by relatively stable vegetation until the Norse settlement (Landnám) AD 870. Quat. Sci. Rev. 259, 106903 (2021).Article 

Google Scholar 
Pansu, J. et al. Reconstructing long-term human impacts on plant communities: an ecological approach based on lake sediment DNA. Mol. Ecol. 24, 1485–1498 (2015).Article 
PubMed 

Google Scholar 
Varotto, C. et al. A pilot study of eDNA metabarcoding to estimate plant biodiversity by an alpine glacier core (Adamello glacier, North Italy). Sci. Rep. 11, 1208 (2021).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Parducci, L. et al. Proxy comparison in ancient peat sediments: Pollen, macrofossil and plant DNA. Philos. Trans. R. Soc. Lond. B Biol. Sci. 370, 20130382 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Clarke, C. L. et al. A 24,000-year ancient DNA and pollen record from the Polar Urals reveals temporal dynamics of arctic and boreal plant communities. Quat. Sci. Rev. 247, 106564 (2020).Article 

Google Scholar 
Niemeyer, B., Epp, L. S., Stoof-Leichsenring, K. R., Pestryakova, L. A. & Herzschuh, U. A comparison of sedimentary DNA and pollen from lake sediments in recording vegetation composition at the Siberian treeline. Mol. Ecol. Resour. 17, e46–e62 (2017).Article 
CAS 
PubMed 

Google Scholar 
Wilson, J. B., Peet, R. K., Dengler, J. & Pärtel, M. Plant species richness: the world records. J. Veg. Sci. 23, 796–802 (2012).Article 

Google Scholar 
Wick, L., van Leeuwen, J. F. N., van der Knaap, W. O. & Lotter, A. F. Holocene vegetation development in the catchment of Sägistalsee (1935 m asl), a small lake in the Swiss Alps. J. Paleolimnol. 30, 261–272 (2003).Article 
ADS 

Google Scholar 
Lotter, A. F. et al. Holocene timber-line dynamics at Bachalpsee, a lake at 2265 m a.s.l. in the northern Swiss Alps. Veg. Hist. Archaeobot. 15, 295–307 (2006).Article 

Google Scholar 
Thöle, L. et al. Reconstruction of Holocene vegetation dynamics at Lac de Bretaye, a high-mountain lake in the Swiss Alps. Holocene 26, 380–396 (2016).Article 
ADS 

Google Scholar 
Heiri, O., Lotter, A. F., Hausmann, S. & Kienast, F. A chironomid-based Holocene summer air temperature reconstruction from the Swiss Alps. Holocene 13, 477–484 (2003).Article 
ADS 

Google Scholar 
Garcés-Pastor, S., Cañellas-Boltà, N., Clavaguera, A., Calero, M. A. & Vegas-Vilarrúbia, T. Vegetation shifts, human impact and peat bog development in Bassa Nera pond (Central Pyrenees) during the last millennium. Holocene 27, 553–565 (2017).Article 
ADS 

Google Scholar 
Aeschimann, D., Lauber, K., Moser, D. M. & Theurillat, J. P. Flora Alpina: Atlas des 4500 Plantes Vasculaires des Alpes (Belin, 2004).Sønstebø, J. H. et al. Using next-generation sequencing for molecular reconstruction of past Arctic vegetation and climate. Mol. Ecol. Resour. 10, 1009–1018 (2010).Article 
PubMed 

Google Scholar 
Diekmann, M. Species indicator values as an important tool in applied plant ecology—a review. Basic Appl. Ecol. 4, 493–506 (2003).Article 

Google Scholar 
Giesecke, T. et al. Postglacial change of the floristic diversity gradient in Europe. Nat. Commun. 10, 5422 (2019).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Colombaroli, D. & Tinner, W. Determining the long-term changes in biodiversity and provisioning services along a transect from Central Europe to the Mediterranean. Holocene 23, 1625–1634 (2013).Article 
ADS 

Google Scholar 
Schwörer, C., Colombaroli, D., Kaltenrieder, P., Rey, F. & Tinner, W. Early human impact (5000–3000 BC) affects mountain forest dynamics in the Alps. J. Ecol. 103, 281–295 (2015).Article 

Google Scholar 
Furtwängler, A. et al. Ancient genomes reveal social and genetic structure of Late Neolithic Switzerland. Nat. Commun. 11, 1915 (2020).Article 
ADS 
PubMed 
PubMed Central 

Google Scholar 
Gilck, F. & Poschlod, P. The origin of alpine farming: A review of archaeological, linguistic and archaeobotanical studies in the Alps. Holocene 29, 1503–1511 (2019).Article 
ADS 

Google Scholar 
Tinner, W., Nielsen, E. H. & Lotter, A. F. Mesolithic agriculture in Switzerland? A critical review of the evidence. Quat. Sci. Rev. 26, 1416–1431 (2007).Article 
ADS 

Google Scholar 
Berthel, N., Schwörer, C. & Tinner, W. Impact of Holocene climate changes on alpine and treeline vegetation at Sanetsch Pass, Bernese Alps, Switzerland. Rev. Palaeobot. Palynol. 174, 91–100 (2012).Article 

Google Scholar 
Hafner, A. & Schwörer, C. Vertical mobility around the high-alpine Schnidejoch Pass. Indications of Neolithic and Bronze Age pastoralism in the Swiss Alps from paleoecological and archaeological sources. Quat. Int. https://doi.org/10.1016/j.quaint.2016.12.049 (2017).Oveisi, M. et al. Potential for endozoochorous seed dispersal by sheep and goats: Risk of weed seed transport via animal faeces. Weed Res. 61, 1–12 (2021).Article 

Google Scholar 
Bardgett, R. D. & Wardle, D. A. Herbivore-mediated linkages between aboveground and belowground communities. Ecology 84, 2258–2268 (2003).Article 

Google Scholar 
Scherrer, D. & Körner, C. Topographically controlled thermal-habitat differentiation buffers alpine plant diversity against climate warming. J. Biogeogr. 38, 406–416 (2011).Article 

Google Scholar 
Giguet-Covex, C. et al. New insights on lake sediment DNA from the catchment: Importance of taphonomic and analytical issues on the record quality. Sci. Rep. 9, 14676 (2019).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Andres, B. Alpine settlement remains in the Bernese Alps (Switzerland) in medieval and modern times. Historical Archaeologies of Transhumance across Europe (eds Costello, E. & Svensson, E.) 155–169 (Routledge, 2018).eTopoi. Journal for Ancient Studies. 3, 279–283 (2012).Grime, J. P. Competitive exclusion in herbaceous vegetation. Nature 242, 344–347 (1973).Article 
ADS 

Google Scholar 
Yuan, Z. Y., Jiao, F., Li, Y. H. & Kallenbach, R. L. Anthropogenic disturbances are key to maintaining the biodiversity of grasslands. Sci. Rep. 6, 22132 (2016).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Spiegelberger, T., Matthies, D., Müller-Schärer, H. & Schaffner, U. Scale-dependent effects of land use on plant species richness of mountain grassland in the European Alps. Ecography 29, 541–548 (2006).Article 

Google Scholar 
Maurer, K., Weyand, A., Fischer, M. & Stöcklin, J. Old cultural traditions, in addition to land use and topography, are shaping plant diversity of grasslands in the Alps. Biol. Conserv. 130, 438–446 (2006).Article 

Google Scholar 
Kampmann, D. et al. Mountain grassland biodiversity: Impact of site conditions versus management type. J. Nat. Conserv. 16, 12–25 (2008).Article 

Google Scholar 
Pellegrini, E., Buccheri, M., Martini, F. & Boscutti, F. Agricultural land use curbs exotic invasion but sustains native plant diversity at intermediate levels. Sci. Rep. 11, 8385 (2021).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Bakker, E. S., Ritchie, M. E., Olff, H., Milchunas, D. G. & Knops, J. M. H. Herbivore impact on grassland plant diversity depends on habitat productivity and herbivore size. Ecol. Lett. 9, 780–788 (2006).Article 
PubMed 

Google Scholar 
Speed, J. D. M., Austrheim, G., Hester, A. J. & Mysterud, A. Elevational advance of alpine plant communities is buffered by herbivory. J. Veg. Sci. 23, 617–625 (2012).Article 

Google Scholar 
Filazzola, A. et al. The effects of livestock grazing on biodiversity are multi-trophic: A meta-analysis. Ecol. Lett. 23, 1298–1309 (2020).Article 
PubMed 

Google Scholar 
Evans, D. M. et al. The cascading impacts of livestock grazing in upland ecosystems: A 10-year experiment. Ecosphere 6, art42 (2015).Article 

Google Scholar 
Alexander, J. M., Diez, J. M. & Levine, J. M. Novel competitors shape species’ responses to climate change. Nature 525, 515–518 (2015).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Mathieu, J. Eine Agrargeschichte der inneren Alpen. Graubünden, Tessin, Wallis 1500–1800 (Chronos, 1992).Aerni, K, Egli, H. R & Fehn, K. Siedlungsprozesse an der Höhengrenze der Ökumene: am Beispiel der Alpen: Referate der 16 Tagung des” Arbeitskreises für genetische Siedlungsforschung in Mitteleuropa” vom 20.−23. (Siedlungsforschung: Spiez, 1991).Brugger, S. O. et al. Alpine glacier reveals ecosystem impacts of Europe’s prosperity and peril over the last millennium. Geophys. Res. Lett. 48, e2021GL095039 (2021).Merkt, J. & Streif, H. Stechrohr-Bohrgeräte für limnische und marine Lockersedimente. Geologisches Jahrbuch 88, 137–148 (1970).Lamb, A. L. Determination of organic and carbonate content in soils and sediments by loss on ignition (LOI), NERC Isotope Geosciences Laboratory Report, 197 (2004).Reimer, P. J. et al. The IntCal20 Northern Hemisphere radiocarbon age calibration curve (0–55 cal kBP). Radiocarbon https://doi.org/10.1017/RDC.2020.41 (2020).Blaauw, M. & Christen, J. A. Flexible paleoclimate age-depth models using an autoregressive gamma process. Bayesian Anal. 6, 457–474 (2011).Article 
MathSciNet 
MATH 

Google Scholar 
Brooks, S. J., Langdon, P. G. & Heiri, O. The identification and use of Palaearctic Chironomidae larvae in palaeoecology. Quat. Res. Assoc. i-vi, 1-276 (2007).Schulze, E. A Key to the Larval Chironomidae and their Instars from Austrian Danube Region Streams and Rivers with Particular Reference to a Numerical Taxonomic Approach. Part I. In: Wasser und Abwasser, Supplementband 3/93. Hrsg.: Bundesamt für Wassergüte, Wien-Kaisermühlen. Schriftenleitung: Werner Kohl. Selbstverlag, 1993, 514 S., öS 562. Acta Hydrochim. Hydrobiol. 22, 191–191 (1994).Article 

Google Scholar 
Juggins, S. C2: Software for ecological and palaeoecological data analysis and visualisation (user guide version 1.5). Newcastle upon Tyne: Newcastle University (2007). https://www.staff.ncl.ac.uk/stephen.juggins/software/code/C2.pdf.Moore, P. D., Webb, J. A. & Collison, M. E. Pollen Analysis, edn 2 (Blackwell, 1991).Stockmarr & Ja Tabletes with spores used in absolute pollen analysis. Pollen Spores 13, 615–621 (1971).
Google Scholar 
Reille, M. Pollen et spores d’Europe et d’Afrique du Nord (Laboratoire de Botanique historique et Palynologie, Marseille, 1992).van Geel, B. et al. Environmental reconstruction of a Roman Period settlement site in Uitgeest (The Netherlands), with special reference to coprophilous fungi. J. Archaeol. Sci. 30, 873–883 (2003).Article 

Google Scholar 
Bennett, K. D. Determination of the number of zones in a biostratigraphical sequence. N. Phytol. 132, 155–170 (1996).Article 
CAS 

Google Scholar 
Tinner, W. et al. Pollen and charcoal in lake sediments compared with historically documented forest fires in southern Switzerland since AD 1920. Holocene 8, 31–42 (1998).Article 
ADS 

Google Scholar 
Adolf, C. et al. The sedimentary and remote-sensing reflection of biomass burning in Europe. Glob. Ecol. Biogeogr. 27, 199–212 (2018).Article 

Google Scholar 
Tinner, W. & Hu, F. S. Size parameters, size-class distribution and area-number relationship of microscopic charcoal: Relevance for fire reconstruction. Holocene 13, 499–505 (2003).Article 
ADS 

Google Scholar 
Parducci, L. et al. Ancient plant DNA in lake sediments. N. Phytol. 214, 924–942 (2017).Article 
CAS 

Google Scholar 
Alsos, I. G. et al. The treasure vault can be opened: Large-scale genome skimming works well using herbarium and silica gel dried material. Plants 9, 432 (2020).Article 
CAS 
PubMed Central 

Google Scholar 
Taberlet, P. et al. Power and limitations of the chloroplast trnL (UAA) intron for plant DNA barcoding. Nucleic Acids Res. 35, e14 (2007).Article 
PubMed 

Google Scholar 
Voldstad, L. H. et al. A complete Holocene lake sediment ancient DNA record reveals long-standing high Arctic plant diversity hotspot in northern Svalbard. Quat. Sci. Rev. 234, 106207 (2020).Article 

Google Scholar 
Boyer, F. et al. obitools: A unix-inspired software package for DNA metabarcoding. Mol. Ecol. Resour. 16, 176–182 (2016).Article 
CAS 
PubMed 

Google Scholar 
Ficetola, G. F. et al. An in silico approach for the evaluation of DNA barcodes. BMC Genomics 11, 434 (2010).Article 
PubMed 
PubMed Central 

Google Scholar 
Soininen, E. M. et al. Highly overlapping winter diet in two sympatric lemming species revealed by DNA metabarcoding. PLoS One 10, e0115335 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
Boratyn, G. M. et al. BLAST: A more efficient report with usability improvements. Nucleic Acids Res. 41, W29–W33 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Leonard, J. A. et al. Animal DNA in PCR reagents plagues ancient DNA research. J. Archaeol. Sci. 34, 1361–1366 (2007).Article 

Google Scholar 
Deiner, K. et al. Environmental DNA metabarcoding: Transforming how we survey animal and plant communities. Mol. Ecol. 26, 5872–5895 (2017).Article 
PubMed 

Google Scholar 
Ter Braak, C. J. F. & Prentice, I. C. A theory of gradient analysis. Adv. Ecol. Res. 18, 271–317 (Elsevier, 1988).Vieira, D. C., Brustolin, M. C., Ferreira, F. C. & Fonseca, G. segRDA: Anr package for performing piecewise redundancy analysis. Methods Ecol. Evol. 10, 2189–2194 (2019).Article 

Google Scholar 
Simpson, G. L. Modelling palaeoecological time series using generalised additive models. Front. Ecol. Evol. 6, 149 (2018).Wood, S. N. Generalized Additive Models: An Introduction with R (Chapman and Hall/CRC, 2017).Lefcheck, J. S. piecewiseSEM: Piecewise structural equation modelling inr for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2016).Article 

Google Scholar 
Chen, W. & Ficetola, G. F. Numerical methods for sedimentary‐ancient‐DNA‐based study on past biodiversity and ecosystem functioning. Environ. DNA 2, 115–129 (2020).Article 

Google Scholar 
Juggins, S. Rioja: Analysis of Quaternary Science Data. R package version 0.9-26. https://cran.r-project.org/web/packages/rioja/index.html (2020).Oksanen, J. et al. vegan: Community Ecology Package. Software http://CRAN.R-project.org/package=vegan (2012).Wickham, H. ggplot2-Elegant Graphics for Data Analysis (Springer, 2016).Fox, J. & Weisberg, S. An R Companion to Applied Regression (Sage, 2019).Tinner, W. & Ammann, B. Long-term responses of mountain ecosystems to environmental changes: Resilience, adjustment, and vulnerability. In Global change and mountain regions. 133–143 (Springer, Dordrecht; 2005). More

The framework of Swin MSIWe established the first Swin Transformer-based mosquito species identification (Swin MSI) model, with the help of self-constructed image dataset and multi-adjustment-test. Gradient-weighted class activation mapping was used to visualize the identification process (Fig. 1a). The key Swin Transformer block was described on Fig. 1b. Based on practical needs, Swin MSI was additional designed to identify Culex pipiens Complex on the subspecies level (Fig. 1c) and novel mosquito (which was defined as ones beyond 17 species in our dataset) classification attribution (Fig. 1d). Detailed results are shown in the following sections.Figure 1The Framework of Swin MSI. (a)The basic architecture for mosquito features extraction and identification. Attention visualization generated by filters at each layer are shown. (b) Details for Swin Transformer block. (c) For mosquito within our dataset 17 species, output is the top 5 confidence species. (d) For mosquito beyond 17 species (defined as novel species), whether the output is a species or a genus is decided after comparing with confidence threshold.Full size imageMosquito datasetsWe established the highest-definition and most-balanced mosquito image dataset to date. The mosquito image dataset covers 7 genera and 17 species (including 3 morphologically similar subspecies in the Cx. pipiens Complex), which covers the most common and important disease-transmitting mosquitoes at the global scale, with a total of 9,900 mosquito images. The image resolution was 4464 × 2976 pixels. The specific taxonomic status and corresponding images are shown in Fig. 2. Due to the limitation of field collection, Ae. vexans, Coquillettidia ochracea, Mansonia uniformis, An. vagus and Toxorhynchites splendens only have females or only have males. In addition, each mosquito species included 300 images of both sexes, which was large enough and same number for each species, in order to balance the capacity and variety of training sets.Figure 2Taxonomic status and index of mosquito species included in this study Both male and female mosquitoes were photographed from different angles such as dorsal, left side, right side, ventral side, etc. Except for 5 species, each mosquito includes 300 images of both sexes, and the resolution of mosquito photos were 4464 × 2976. Cx. pipiens quinquefasciatus, Cx. pipiens pallens, and Cx. pipiens molestus (subspecies level, in dark gray background) were 3 subspecies in Cx. pipiens Complex (species level).Full size imageWorkflow for mosquito species identificationA three-stage flowchart of building best deep learning model for identification of mosquito species model was adopted (Fig. 3). The first learning stage was conducted by three CNNs (the Mask R-CNN, DenseNet, and YOLOv5) and three transformer models (the Detection Transformer, Vision Transformer, and Swin Transformer). Based on the performance of the first-stage model and the real mosquito labels, the second learning stage involved adjusting the model parameters of the three Swin Transformer variants (T, B, and L) to compare their performances. The third learning stage involved testing the effects of inputting differently sized images (384 × 384 and 224 × 224) to the Swin Transformer-L model; finally, we proposed a deep learning model for mosquito species identification (Swin MSI) to test the recognition effects of different mosquito species. The model was validated on different mosquito species, with a focus on the identification accuracy of three subspecies within the Cx. pipiens Complex and the detection effect of novel mosquito species.Figure 3Flowchart of testing deep learning model for intelligent identification of mosquito species.Full size imageComparison between the CNN model and Transformer model results (1st round of learning)Figure 4a shows the accuracies obtained for the six different computer vision network models tested on the mosquito picture test set. The test results show that the transformer network model had a higher mosquito species discrimination ability than the CNN.Figure 4Comparison of mosquito recognition effects of computer vision network models and variants. (a) Comparison of mosquito identification accuracy between 3 CNNs and 3 Transformer; (b) The best effect CNN (YOLOv5) training set loss curve(blue), validation set loss curve(green) and validation set accuracy curve(orange); (c) The best effect Transformer (Swin Transformer) training set loss curve, validation set loss curve and validation set accuracy curve. (d) Swin-MSI-T test result confusion matrix; (e) Swin-MSI -B test result confusion matrix; (f) Swin-MSI -L test result confusion matrix. Confusion matrix of mosquito labels in which odd numbers represent females and even numbers represent males. The small squares in the confusion matrix represent the recognition readiness rate, from red to green, the recognition readiness rate is getting higher and higher An. sinensis: 1, 2; Cx. pipiens quinquefasciatus: 3, 4; Cx. pipiens pallens: 5, 6; Cx. pipiens molestus: 7,8 Cx. modestus: 9,10; Ae. albopictus: 11, 12 Ae. aegypti: 13, 14; Cx. pallidothorax: 15, 16 Ae. galloisi: 17,18 Ae. vexans: 19, 20; Ae. koreicus: 21, 22 Armigeres subalbatus: 23, 24; Coquillettidia ochracea: 25, 26 Cx. gelidus: 27, 28 Cx. triraeniorhynchus: 29, 30 Mansonia uniformis: 31, 32 An. vagus: 33, 34 Ae. elsaie: 35,36 Toxorhynchites splendens: 37, 38.Full size imageIn the CNN training process (applied to YOLOv5), the validation accuracy requires more than 110 epochs to grow to 0.9, and the validation loss requires 110 epochs to drop to a flat interval; in contrast, during the training step, these losses represent a continuously decreasing process. These results indicate that the deep learning model derived based on the Swin Transformer algorithm was able to achieve a higher recognition accuracy in less time than the rapid convergence ability of the CNN during the iterative process (Fig. 4b).The Swin Transformer model exhibited the highest test accuracy of 96.3%. During the training process, the loss of this model could stabilize after 30 epochs, and its validation accuracy could grow to 0.9 after 20 epochs; during the validation step, the loss can drop to 0.36 after 20 epochs, after which the loss curve fluctuated but did not produce adverse effects (Fig. 4c). Based on the excellent performance of the Swin Transformer model, this model was used as the baseline to carry out the subsequent analyses.Swin Transformer model variant adjustment (2nd round of learning)Following testing performed to clarify the superior performance of the Swin Transformer algorithm, we chose different Drop_path_rate, Embed_dim and Depths parameter settings and labeled the parameter sets as the Swin Transformer-T, Swin Transformer-B, and Swin Transformer-L variants. Drop_path is an efficient regularization method, and an asymmetric Drop_path_rate is beneficial for supervised representation learning when using image classification tasks and Transformer architectures. The Embed_dim parameter represents the image dimensions obtained after the input red–green–blue (RGB) image is calculated by the Swin Transformer block in stage 1. The Depths parameter is the number of Swin Transformer blocks used in the four stages. The parameter information and test results are shown in Table 1. Due to the increase in the Swin Transformer block and Embed_dim parameters in stage 3, the recognition accuracies of the three variants were found to be 95.8%, 96.3%, and 98.2%, Correspondingly, the f1 score were 96.2%, 96.7% and 98.3%; thus, these variants could effectively improve the mosquito species identification ability in a manner similar to the CNN by increasing the number of convolutional channels to extract more features and improve the overall classification ability. In this study, the Swin Transformer-L variant, which exhibited the highest accuracy, was selected as the baseline for the next work.Table 1 Parameters and test accuracy of three variants of Swin Transformer.Full size tableBy plotting a confusion matrix of the test set results derived using the three Swin Transformer variants, we clearly obtained the proportion of correct and incorrect identifications in each category to visually reflect the mosquito species discrimination ability (Fig. 4d–f). In the matrix, the darker diagonal colors indicate higher identification rate accuracies of the corresponding mosquito categories. Among them, five mosquito species were missing because the Ae. vexans, Coquillettidia ochracea, Mansonia uniformis, An. vagus and Toxorhynchites splendens species were represented in the dataset by only females or only males. The confusion matrix shown in Panel C lists the lowest number of mosquito species identification error points and the lowest accuracy level obtained in each category, suggesting that the Swin Transformer-L model has a better classification performance than the Swin Transformer-T and Swin Transformer-B models.Effect of the input image size on the discrimination ability (3rd round of learning)To investigate the relationship between the input image size and mosquito species identification performance, in this study, we conducted a comparison test between input images with sizes of 224 × 224 and 384 × 384, based on the Swin Transformer-L model, and identified 8 categories of mosquito identification accuracy differences. These test results are shown in Table 2. When using an image size of 224 × 224 pixels, the batch_size parameter was set to 16, and when using an image size of 384 × 384 pixels, the batch_size parameter was set to 4; under these conditions, the proportion of utilized video memory accounted for 67%, as shown in Eq. 4, and this was consistent with the description of the relationship between the size of self-attentive operations during the operation of the Swin Transformer model when 384 × 384 pixels images were used. The time required for the Transformer-L model to complete all the training sessions was excessive, reaching 126 h and even exceeding the 124 h required by the YOLOv5 model, which was found to require the highest computation time during the training process in this work. Long-term training process could more fully reflect the performance differences between models. Fortunately and actually, the response speed of the model will not be affected by the training time. Compared to the accuracy of 98.2% obtained for 224 × 224 inputs, the 384 × 384 input image size derived based on the Swin Transformer-L model provided a higher mosquito species identification accuracy of 99.04%, representing an improvement of 0.84%.$$Omega ({text{W}} – {text{MSA}}) = 4{text{HWC}}^{2} + 2{text{M}}^{2} {text{HWC}}$$
(1)
Table 2 Comparison of recognition accuracy for different input image sizes.Full size tableVisualizing and understanding the Swin MSI modelsTo investigate the differences in the attentional features utilized by the Swin MSI and taxonomists for mosquito species identification, we applied the Grad-CAM method to visualize the Swin MSI attentional areas on mosquitoes at different stages. Because the Swin Transformer has different attentional ranges among its multi-head self-attention steps in different stages, different attentional weights can be found on different mosquito positions. In stage 1, the feature dimension of each patch was 4 × 4 × C, thus enabling the Swin Transformer’s multi-head self-attention mechanism to give more attention to the detailed parts of the mosquitoes, such as their legs, wings, antennae, and pronota. In stage 2, the feature dimension of each patch was 8 × 8 × 2C, enabling the Swin Transformer’s multi-head self-attention mechanism to focus on the bodies of the mosquitoes, such as their heads, thoraces, and abdomens. In stage 3, when the feature dimension of each patch was 16 × 16 × 4C, the Swin Transformer’s multi-head self-attention mechanism could focus on most regions of the mosquito, thus forming a global overall attention mechanism for each mosquito (Fig. 5). This attentional focus process is essentially the same as the process used by taxonomists when classifying mosquito morphology, changing from details to localities to the whole mosquito.Figure 5Attention visualization of representative mosquitoes of the genera Ae., Cx., An., Armigeres, Coquillettidia and Mansonia. This is a visualization for identifying the regions in the image that can explain the classification progress. Images of Toxorhynchites contain only males, with obvious differences in morphological characteristics, are not shown.Full size imageAe. aegypti is widely distributed in tropical and subtropical regions around the world and transmits Zika, dengue and yellow fever. A pair of long-stalked sickle-shaped white spots on both shoulder sides of the mesoscutum, with a pair of longitudinal stripes running through the whole mesotergum, is the most important morphological identification feature of this species. This feature was the deepest section in the attention visualization, indicating that the Swin MSI model also recognized it as the principal distinguishing feature. In addition, the abdominal tergum of A. aegypti is black and segments II-VII have lateral silvery white spots and basal white bands; the model also focused on these areas.Cx. triraeniorhynchus is the main vector of Japanese encephalitis; this mosquito has a small body size, a distinctive white ring on the proboscis (its most distinctive morphological feature), and a peppery color on its whole body. Similarly, the model constructed herein focused on both the head and abdominal regions of this species.An. sinensis is the top vector of malaria in China and has no more than three white spots on its anterior wing margin and a distinct white spot on its marginal V5.2 fringe; this feature was observed in Stage 2, at which time the modelstrongly focused on the corresponding area.The most obvious feature of Armigeres subalbatus is the lateral flattening and slightly downward curving of its proboscis; the observation of the attention visualization revealed that the constructed model focused on these regions from Stage 1 to Stage 3. The mesoscutum and abdominal tergum were not critical and were less important for identification than the proboscis, and the attention visualization results correspondingly show that the neural network focused less on these features.Coquillettidia ochracea belongs to the Coquillettidia genus and is golden yellow all over its body, with the most pronounced abdomen among the analyzed species. The model showed a consistent morphological taxonomic focus on the abdomen of this species.Mansonia uniformis is a vector of Malayan filariasis. The abdominal tergum of this species is dark brown, and its abdominal segments II-VII have yellow terminal bands and lateral white spots, which are more obvious than the dark brown feature on proboscis. Through the attention visualization, we determined that the Swin MSI model was more concerned with the abdominal region features than with the proboscis features.Subspecies-level identification tests of mosquitos in the Culex pipiens ComplexFine-grained image classification has been the focus of extensive research in the field of computer vision25,26. Based on the test set (containing 270 images) constructed herein for three subspecies of the Cx. pipiens Complex, the subspecies and sex identification accuracies were 100% when the Swin MSI model was used.The morphological characteristics of Cx. pipiens quinquefasciatus, Cx. pipiens pallens, and Cx. pipiens molestus within the Cx. pipiens Complex are almost indistinguishable, but their host preferences, self-fertility properties, breeding environments, and stagnation overwintering strategies are very different27. Among the existing features available for morphological classification, the stripes on the abdominal tergum of Cx. pipiens quinquefasciatus are usually inverted triangles and are not connected with the pleurosternums, while those of Cx. pipiens pallens are rectangular and are connected with the pleurosternums. Cx. pipiens molestus is morphologically more similar to Cx. pipiens pallens as an ecological subspecies of the Cx. pipiens Complex. However, taxonomists do not recommend using the unstable feature mentioned above as the main taxonomic feature for differentiation. By analyzing the attention visualization results of these three subspecies (last three rows on Fig. 5), we found that the neural networks of Cx. pipiens quinquefasciatus, Cx. pipiens pallens, and Cx. pipiens molestus still focused on the abdominal regions, as shown in dark red. The area of focus of these neural networks differ from that of the human eye, and the results of this study suggest that the Swin MSI model can detect finely granular features among these three mosquito subspecies that are indistinguishable to the naked human eye.Novel mosquito classification attributionAfter we performed a confidence check on the successfully identified mosquito images in the dataset, the lowest confidence value was found to be 85%. A higher confidence threshold mean stricter evaluation criteria, which can better reflect the powerful performance of the model. Therefore, 0.85 was set as the confidence threshold when judging novel mosquitoes. When identifying 10 unknown mosquito species, the highest derived species confidence level was below 85%; when the results were output to the genus level (Fig. 1d), the average probability of obtaining a correct judgment was 96.26%accuracy and 98.09% F1-score (Table 3). The images tested as novel Ae., Cx. and An. mosquito were from Minakshi and Couret et al.28,29.Table 3 Probability of correct attribution of novel species.Full size table More

Most records of N-tiger cats were from savanna environments, and it was not surprising that this vegetative formation has a key influence on the N-tiger cat range in the Amazon. The bulk of the L. t. tigrinus distribution lies in the savannas, dry forests and shrublands of the Cerrado and Caatinga biomes. These are also the areas with the vast majority of records for this lowland subspecies (Supplemental Fig. S5). Hence, L. t. tigrinus is more associated with savannas and savanna-like environments than with rainforests. In fact, more than 80% of the records in the Amazon were within 100 km of a savanna patch. Colonization of the northern savanna formations of the Amazon by the N-tiger cat likely occurred during the forest-savanna shifts of the glacial period18, and the cat currently shows a patchy distribution. Strong evidence of established biogeographic corridor connections between the savannas of the Cerrado and those of the Amazon exists, suggesting northward expansion of the former during glacial periods, perhaps predating the Last Glacial Maximum19,20,21. Further corroborating this evidence, tiger cat ‘gene flow’ niche modelling showed prior connectivity between the Guiana population and that of Central Brazil and no connectivity with the Andean population22. Additionally, Guianan tiger cat skin patterns are found in savanna and transitional savanna/Amazon areas and in the semiarid shrub-woodland of Brazil and are very distinct from the patterns of the tiger cats from the Andes of northwestern South America and Central America (Supplementary Information Fig. S6).The bioclimatic variables in the best model also supported the cat’s preference for savanna areas. The best model indicated a positive effect of precipitation in the driest month on the probability of the presence of the N-tiger cat, likely indicating the Aw/As climates of tropical savannas23. These climates are marked by seasonal variation in rainfall, with a pronounced dry season. Higher rainfall during the dry season favors the growth of vegetation, which results in some tree cover within the savannas. Thus, our results agree with previous research suggesting that tiger cats avoid open savanna formations24. Similarly, the species had a significant negative response to net primary productivity. This also supports the species’ avoidance of dense lowland rainforests, which are the most productive habitats. In the Amazon biome, the least productive areas are found in more open landscapes25.The N-tiger cat’s range considered from an ecoregion perspective12 could biogeographically explain its distribution in the Amazon. All records but 2 fell within Guiana savannas, Guiana highland forest, Guiana rainforest, part of the Uatumã-Trombetas rainforest bordering the Guianas or all of it connecting to Gurupá and Monte Alegre varzea forests, as well as Marajó varzeas, the interfluve Tocantins-Araguaia/Maranhão, and the southern block of the interfluve Xingu/Tocantins-Araguaia. There were two records from the Negro-Branco moist forest, which also includes savanna-like “campinarana” formations. The range also reaches the transitional babaçu palm forests of Maranhão and the Mato Grosso seasonal forests (Supplementary Information Fig. S7, Table S3). The N-tiger cat’s range in the Amazon was determined by combining records with species distribution modeling, also matching the ecoregion perspective.Outside the Guiana Shield and likely the savanna patches of the region of the Upper Negro River, in other parts of the Amazon, the N-tiger cat seems to be restricted to the forests of the eastern Amazon, along the arc of deforestation and to transitional areas with savanna formations. The presence and absence points at camera-trapping sites could explain the N-tiger cat’s range in the Amazon and define its distribution range in the biome. Absence points, for instance, were usually located in dense rainforest habitats throughout the Amazon biome.The species may occasionally occupy rainforests, such as those of the Guianas, where it tends to be very rare. At a site in central Suriname, after an enormous trapping effort of  > 20,000 trap days in four years by cat specialists, over an area  > 1100 km2, no records of the N-tiger cat were found (Supplementary Information Table S2), although its presence is expected in that area26. This finding attests to the inherent rarity of this felid in its limited range within the Amazon. However, could its association with the arc of deforestation be related to the replacement of forest by bushy savanna-like vegetation that succeeds abandoned pastures? The other currently recognized subspecies, L. t. pardinoides (the Andean tiger cat) and L. t. oncilla (the oncilla), and the recently split southern tiger cat L. guttulus are all associated with forested areas. Conversely, L. t. tigrinus has higher abundance and is mostly found in the nonforested habitats of the Cerrado and Caatinga domains of Brazil and only rarely in rainforests. Thus, L. t. tigrinus may be an open-habitat (sub)species. However, within savannas, N-tiger cats are restricted to denser savanna formations, with open savannas deemed unsuitable24. In the semiarid Caatinga, the N-tiger cat also prefers denser formations27,28.One of the most interesting findings was the clear relationship between the ranges of the dominant mesopredator and subordinate species. The ranges of ocelots and N-tiger cats in the Amazon were diametrically opposite (Fig. 1), a finding never recorded for felids. The reported ocelot densities and relative abundance indexes (RAIs) in the Amazon range from 0.29 to 0.95 ind/km2 and 0.07–13.2 ind/100 trap-days, respectively7,29. Thus, the expected ocelot density found using modeling that allows for N-tiger cat presence is very low (Fig. 2A). In the Rupununi, the ocelot:N-tiger cat RAI ratio was roughly 10:1, with a very low RAI and expected density for N-tiger cats (see Supplementary Material). The only other relative abundance estimate of tiger cats presented for the Amazon30 was not confirmed as an estimate of tiger cats following inspection of the original records by the authors but as an estimate of margays or ocelots. This antagonistic relationship between ocelots and all other small cat species in their area of sympatry is quite impressive. It is density-dependent, as it seems to take effect only above an ocelot density threshold of 0.12 ind./km231. The influence can range from patterns of density, distribution, and occupancy to spatial and temporal use. Conversely, such an impact was not detected when either the small cats or ocelots were compared to the larger cats31,32,33,34,35.In view of the Red List assessments and applying the limited estimates presented, the expected total population size for N-tiger cats in the Amazon would be approximately 150 and 1622 individuals, considering their AOO or EOO, respectively. Applying the IUCN’s formula for mature individuals8, these numbers would be 45 and 487 individuals for the AOO and EOO, respectively.The ocelot’s preference for very dense rainforests may explain the low probability of N-tiger cat occurrence within the Amazon biome. Notably, most tiger cat records from rainforests and all those from premontane forests came from the Guiana Shield, a region where tropical grasslands and savannas dot more forested landscapes. The Guiana Highlands and Pantepui ecoregions, which make up a considerable portion of the shield, tend to have low ocelot densities (below 0.30 ind/km2), although they do contain some rainforest. Ocelot densities reach some of their lowest values in the Guianan savanna ecoregion (mean ocelot density of 0.029 in the savanna formations), where the N-tiger cat probability of occurrence was highest. At the Karanambu site in the Rupununi, all ocelot records came from either gallery forests or forest patches embedded in the savanna. Although the data did not allow us to test further hypotheses, it is likely that spatial partitioning occurs in the Guiana Shield, with N-tiger cats favoring habitats that are more open. Conversely, areas farther west in the Amazon biome, other than the predicted area, do not have any major savanna patches and are covered mostly by lowland tropical rainforest formations, where ocelots can potentially reach densities in excess of 0.7 ind/km2. Of all Amazonian records of N-tiger cats, only one came from west of the 68th meridian: a preserved specimen from Puerto Leguizamo on the Putumayo River in Colombia. The specimen was identified as L. t. pardinoides by its collector, so it most likely represents an individual that came down from the foothills of the Andes. Alternatively, it could have been caught in the Andean foothills but labeled generally as from Puerto Leguizamo, as museum records do not always present precise locations, like most of those from our dataset; thus, they could represent a broader region, not a single collection location.The records of L. t. tigrinus in the Monte-Alegre Várzea ecoregion and Tapajós-Xingu Moist Forest ecoregion (which shares a border with the Amazon River) are actually from the small savanna patches of Terra Santa and Alter do Chão, respectively, which are imbedded within the forests of these ecoregions. Similarly, the Negro-Branco Moist Forest ecoregion includes open-canopy white sand forests with savanna-like vegetation, known as ‘campinaranas’36.Although our model predicted a high probability of N-tiger cat presence in the Marajó Várzea ecoregion, the records from the island came from savanna patches and not from flooded forests and mangroves. Hence, we did not include such large areas in the AOO for the subspecies. It is likely that the highly predicted probability of presence there is an artifact of low predicted ocelot density. Nevertheless, the environment there is not suitable for either cat. Our ocelot density model was highly significant and explained almost 50% of the variation in ocelot density. The remaining variation was related to either other variables that could not be measured via satellite imagery (such as prey availability) or the sampling design of the different studies. Nonetheless, ocelot densities predicted from our model across the Amazon were within the expected range for the species29.Why are N-tiger cats absent in camera-trapping studies in Amazonian forests throughout the biome? The most straightforward answer seems to be because they simply are not there (central and western Amazon) or, where present, their numbers are extremely low (Guianas and eastern Amazon). The lack of surveys cannot be cited as a potential reason for their apparent absence because the studies that did not detect the species were conducted throughout the Amazon biome, in all nine Amazonian countries. Some of the areas have been surveyed for several years—or decades in some cases—and have failed to record a single individual (Supplementary Information Table S2). Typically, N-tiger cats appear, even prominently, on cameras in other biomes, such as in the savannas of the Cerrado and semiarid scrub of the Caatinga domain in Brazil, including sites where ocelots are present24,27,37. Clouded tiger cats (L. t. pardinoides) have also been frequently recorded on cameras in the Andes, higher than 1500 m above sea level34,38, but not in lowland Amazonian forests. This finding indicates that the N-tiger cat is not camera-shy. In northern Brazilian savannas, its density can reach 0.25 ind/km2 24. Coincidentally, this highest density estimate of the N-tiger cat is the same as the lowest ocelot density estimate for Amazonian forests24,29.Tiger cats and margays show high similarity, making misidentifications relatively common39. However, the evaluation of  > 3000 camera trap images of small-medium felids in the Amazon revealed that only one mildly resembled a tiger cat, a finding that supports the species being absent there and does not represent a case of mistaken identity with margays or even ocelots7.The Amazonian range of L. tigrinus is very limited, and populations are expected to be very small. With the upcoming split of L. t. tigrinus and L. t. pardinoides into two different species40, this situation would have serious implications for the conservation of the former. Thus, L. t. tigrinus conservation lies outside the “Amazonian safe haven” of most other carnivore species found there7. The Brazilian drylands Cerrado and Caatinga represent such places for L. t. tigrinus populations. Unfortunately, these biomes have had  > 50% of their cover completely removed41. Very importantly, besides being extremely rare in the Amazonian savannas, this rather limited vegetative formation is also considered highly threatened and of conservation priority42. Therefore, the tiger cat could become an emblematic flagship species representing the uniqueness of this vegetative formation in dire need of protection.In short, the picture that emerges is that although the N-tiger cat uses both rainforests and deciduous forests in the Amazon, it seems to be mostly associated with savanna formations and that its distribution in the Amazon is highly influenced by the ocelot, the dominant mesopredator. The N-tiger cat’s inherent rarity, expected population size, and restricted range in the Amazon suggest that this biome does not in fact represent a safe haven for the global conservation of this small felid. In addition to shedding light on and refining the N-tiger cat distribution in the Amazon, this paper highlights the importance of including biological variables, such as the potential impacts of competitors and predators on species presence and distribution, in SDMs. More

Hirai, H. et al. Organic micropollutants in marine plastics debris from the open ocean and remote and urban beaches. Mar. Pollut. Bull. 62(8), 1683–1692. https://doi.org/10.1016/j.marpolbul.2011.06.004 (2011).Article 
CAS 
PubMed 

Google Scholar 
Masó, M., Garcés, E., Pagès, F. & Camp, J. Drifting plastic debris as a potential vector for dispersing harmful algal bloom (HAB) species. Sci. Mar. 67(1), 107–111. https://doi.org/10.3989/scimar.2003.67n1107 (2003).Article 

Google Scholar 
Pandey, D., Singh, A., Ramanathan, A. & Kumar, M. The combined exposure of microplastics and toxic contaminants in the floodplains of North India: A review. J. Environ. Manag. 279, 111557. https://doi.org/10.1016/j.jenvman.2020.111557 (2021).Article 
CAS 

Google Scholar 
Peng, L. et al. Micro- and nano-plastics in marine environment: Source, distribution and threats—A review. Sci. Total Environ. 698, 134254. https://doi.org/10.1016/j.scitotenv.2019.134254 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Rillig, M. C. & Lehmann, A. Microplastic in terrestrial ecosystems and the soil?. Environ. Sci. Technol. 46(12), 6453–6454. https://doi.org/10.1021/es302011r (2012).Article 
CAS 
PubMed 

Google Scholar 
Rochman, C. M. & Hoellein, T. The global odyssey of plastic pollution. Science 368(6496), 1184–1185. https://doi.org/10.1126/science.abc4428 (2020).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Jan Kole, P., Löhr, A. J., van Belleghem, F. G. A. J. & Ragas, A. M. J. Wear and tear of tyres: A stealthy source of microplastics in the environment. Int. J. Environ. Res. Public Health https://doi.org/10.3390/ijerph14101265 (2017).Article 

Google Scholar 
Sommer, F. et al. Tire abrasion as a major source of microplastics in the environment. Aerosol Air Qual. Res. 18(8), 2014–2028. https://doi.org/10.4209/aaqr.2018.03.0099 (2018).Article 
CAS 

Google Scholar 
Beita-Sandí, W., Selbes, M., Ersan, M. S. & Karanfil, T. Release of nitrosamines and nitrosamine precursors from scrap tires. Environ. Sci. Technol. Lett. 6(4), 251–256. https://doi.org/10.1021/acs.estlett.9b00172 (2019).Article 
CAS 

Google Scholar 
Kaminsky, W. & Mennerich, C. Pyrolysis of synthetic tire rubber in a fluidised-bed reactor to yield 1,3-butadiene, styrene and carbon black. J. Anal. Appl. Pyrolysis 58–59, 803–811. https://doi.org/10.1016/S0165-2370(00)00129-7 (2001).Article 

Google Scholar 
Sundt, P., Schulze, P. E. & Syversen, F. Sources of microplastic- pollution to the marine environment. Mepex Nor. Environ. Agency 86, 20 (2014).
Google Scholar 
White, W. C. Butadiene production process overview. Chem. Biol. Interact. 166(1–3), 10–14. https://doi.org/10.1016/j.cbi.2007.01.009 (2007).Article 
CAS 
PubMed 

Google Scholar 
Alimi, O. S., Farner Budarz, J., Hernandez, L. M. & Tufenkji, N. Microplastics and nanoplastics in aquatic environments: Aggregation, deposition, and enhanced contaminant transport. Environ. Sci. Technol. 52(4), 1704–1724. https://doi.org/10.1021/acs.est.7b05559 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Cooper, D. A. & Corcoran, P. L. Effects of mechanical and chemical processes on the degradation of plastic beach debris on the island of Kauai, Hawaii. Mar. Pollut. Bull. 60(5), 650–654. https://doi.org/10.1016/j.marpolbul.2009.12.026 (2010).Article 
CAS 
PubMed 

Google Scholar 
O’Brine, T. & Thompson, R. C. Degradation of plastic carrier bags in the marine environment. Mar. Pollut. Bull. 60(12), 2279–2283. https://doi.org/10.1016/j.marpolbul.2010.08.005 (2010).Article 
CAS 
PubMed 

Google Scholar 
Song, Y. K. et al. Combined effects of UV exposure duration and mechanical abrasion on microplastic fragmentation by polymer type. Environ. Sci Technol. 51(8), 4368–4376. https://doi.org/10.1021/acs.est.6b06155 (2017).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Adobe Inc. (2019). Adobe illustrator. Retrieved from https://www.adobe.com/Products/Illustrator.Chamas, A. et al. Degradation rates of plastics in the environment. ACS Sustain. Chem. Eng. 8(9), 3494–3511. https://doi.org/10.1021/acssuschemeng.9b06635 (2020).Article 
CAS 

Google Scholar 
Councell, T. B., Duckenfield, K. U., Landa, E. R. & Callender, E. Tire-wear particles as a source of zinc to the environment. Environ. Sci. Technol. 38(15), 4206–4214. https://doi.org/10.1021/es034631f (2004).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Awet, T. T. et al. Effects of polystyrene nanoparticles on the microbiota and functional diversity of enzymes in soil. Environ. Sci. Eur. https://doi.org/10.1186/s12302-018-0140-6 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Chung, H., Son, Y., Yoon, T. K., Kim, S. & Kim, W. The effect of multi-walled carbon nanotubes on soil microbial activity. Ecotoxicol. Environ. Saf. 74(4), 569–575. https://doi.org/10.1016/j.ecoenv.2011.01.004 (2011).Article 
CAS 
PubMed 

Google Scholar 
Huber, M., Welker, A. & Helmreich, B. Critical review of heavy metal pollution of traffic area runoff: Occurrence, influencing factors, and partitioning. Sci. Total Environ. 541, 895–919. https://doi.org/10.1016/j.scitotenv.2015.09.033 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Miandad, R., Barakat, M. A., Aburiazaiza, A. S., Rehan, M. & Nizami, A. S. Catalytic pyrolysis of plastic waste: A review. Process Saf. Environ. Prot. 102, 822–838. https://doi.org/10.1016/j.psep.2016.06.022 (2016).Article 
CAS 

Google Scholar 
Zhang, X., Li, H., Cao, Q., Jin, L. & Wang, F. Upgrading pyrolytic residue from waste tires to commercial carbon black. Waste Manag. Res. 36(5), 436–444. https://doi.org/10.1177/0734242X18764292 (2018).Article 
CAS 
PubMed 

Google Scholar 
Zhu, D., Li, G., Wang, H. T. & Duan, G. L. Effects of nano- or microplastic exposure combined with arsenic on soil bacterial, fungal, and protistan communities. Chemosphere 281, 130998. https://doi.org/10.1016/j.chemosphere.2021.130998 (2021).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Pathan, S. I. et al. Soil Pollution from micro-and nanoplastic debris: A hidden and unknown biohazard. Sustainability 12(18), 1–31. https://doi.org/10.3390/su12187255 (2020).Article 
CAS 

Google Scholar 
Rillig, M. C. & Bonkowski, M. Microplastic and soil protists: A call for research. Environ. Pollut. 241, 1128–1131. https://doi.org/10.1016/j.envpol.2018.04.147 (2018).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Zettler, E. R., Mincer, T. J. & Amaral-Zettler, L. A. Life in the “Plastisphere”: Microbial communities on plastic marine debris. Environ. Sci. Technol. 47(13), 7137–7146. https://doi.org/10.1021/es401288x (2013).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Burns, E. E. & Boxall, A. B. A. Microplastics in the aquatic environment: Evidence for or against adverse impacts and major knowledge gaps. Environ. Toxicol. Chem. 37(11), 2776–2796. https://doi.org/10.1002/etc.4268 (2018).Article 
CAS 
PubMed 

Google Scholar 
Bradney, L. et al. Particulate plastics as a vector for toxic trace-element uptake by aquatic and terrestrial organisms and human health risk. Environ. Int. 2019(131), 104937. https://doi.org/10.1016/j.envint.2019.104937 (2018).Article 
CAS 

Google Scholar 
Duis, K. & Coors, A. Microplastics in the Aquatic and Terrestrial Environment: Sources (with a Specific Focus on Personal Care Products), fate and effects. Environ. Sci. Eur. 28(1), 1–25. https://doi.org/10.1186/s12302-015-0069-y (2016).Article 
CAS 

Google Scholar 
Geyer, R., Jambeck, J. R. & Law, K. L. Production, use, and fate of all plastics ever made. Sci. Adv. 3(7), 25–29. https://doi.org/10.1126/sciadv.1700782 (2017).Article 
CAS 

Google Scholar 
Jayasiri, H. B., Purushothaman, C. S. & Vennila, A. Quantitative analysis of plastic debris on recreational beaches in Mumbai, India. Mar. Pollut. Bull. 77(1–2), 107–112. https://doi.org/10.1016/j.marpolbul.2013.10.024 (2013).Article 
CAS 
PubMed 

Google Scholar 
Lassen, C., Hansen, S. F., Magnusson, K., Hartmann, N. B., Rehne Jensen, P., Nielsen, T. G. & Brinch, A. Microplastics occurrence, effects and sources of releases (2015).Weithmann, N. et al. Organic fertilizer as a vehicle for the entry of microplastic into the environment. Sci. Adv. https://doi.org/10.1126/sciadv.aap8060 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Hidalgo-Ruz, V., Gutow, L., Thompson, R. C. & Thiel, M. Microplastics in the marine environment: A review of the methods used for identification and quantification. Environ. Sci. Technol. 46(6), 3060–3075. https://doi.org/10.1021/es2031505 (2012).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Boenigk, J., Matz, C., Jürgens, K. & Arndt, H. Confusing selective feeding with differential digestion in bacterivorous nanoflagellates. J. Eukaryot. Microbiol. 48(4), 425–432. https://doi.org/10.1111/j.1550-7408.2001.tb00175.x (2001).Article 
CAS 
PubMed 

Google Scholar 
Boenigk, J., Matz, C., Jürgens, K. & Arndt, H. Food concentration-dependent regulation of food selectivity of interception-feeding bacterivorous nanoflagellates. Aquat. Microb. Ecol. 27(2), 195–202. https://doi.org/10.3354/ame027195 (2002).Article 

Google Scholar 
Wright, S. L., Thompson, R. C. & Galloway, T. S. The physical impacts of microplastics on marine organisms: A review. Environ. Pollut. 178, 483–492. https://doi.org/10.1016/j.envpol.2013.02.031 (2013).Article 
CAS 
PubMed 

Google Scholar 
Moore, C. J. Synthetic polymers in the marine environment: A rapidly increasing, long-term threat. Environ. Res. 108(2), 131–139. https://doi.org/10.1016/j.envres.2008.07.025 (2008).Article 
CAS 
PubMed 

Google Scholar 
Fu, S. F. et al. Exposure to polystyrene nanoplastic leads to inhibition of anaerobic digestion system. Sci. Total Environ. 625, 64–70. https://doi.org/10.1016/j.scitotenv.2017.12.158 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Bock, C. et al. Factors shaping community patterns of protists and bacteria on a European scale. Environ. Microbiol. 22(6), 2243–2260. https://doi.org/10.1111/1462-2920.14992 (2020).Article 
PubMed 

Google Scholar 
Besseling, E., Wang, B., Lürling, M. & Koelmans, A. A. Nanoplastic affects growth of S. obliquus and reproduction of D. magna. Environ. Sci. Technol. 48(20), 12336–12343. https://doi.org/10.1021/es503001d (2014).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Brown, D. M., Wilson, M. R., MacNee, W., Stone, V. & Donaldson, K. Size-dependent proinflammatory effects of ultrafine polystyrene particles: A role for surface area and oxidative stress in the enhanced activity of ultrafines. Toxicol. Appl. Pharmacol. 175(3), 191–199. https://doi.org/10.1006/taap.2001.9240 (2001).Article 
CAS 
PubMed 

Google Scholar 
Jeong, C. B. et al. Microplastic size-dependent toxicity, oxidative stress induction, and p-JNK and p-P38 activation in the monogonont rotifer (Brachionus Koreanus). Environ. Sci. Technol. 50(16), 8849–8857. https://doi.org/10.1021/acs.est.6b01441 (2016).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Kang, H. C., Jeong, H. J., Jang, S. H. & Lee, K. H. Feeding by common heterotrophic protists on the phototrophic dinoflagellate Biecheleriopsis adriatica (Suessiaceae) compared to that of other suessioid dinoflagellates. Algae 34(2), 127–140. https://doi.org/10.4490/algae.2019.34.5.29 (2019).Article 
CAS 

Google Scholar 
Sjollema, S. B., Redondo-Hasselerharm, P., Leslie, H. A., Kraak, M. H. S. & Vethaak, A. D. Do plastic particles affect microalgal photosynthesis and growth?. Aquat. Toxicol. 170, 259–261. https://doi.org/10.1016/j.aquatox.2015.12.002 (2016).Article 
CAS 
PubMed 

Google Scholar 
Rossi, G., Barnoud, J. & Monticelli, L. Polystyrene nanoparticles perturb lipid membranes. J. Phys. Chem. Lett. 5(1), 241–246. https://doi.org/10.1021/jz402234c (2014).Article 
CAS 
PubMed 

Google Scholar 
Brandts, I. et al. Effects of nanoplastics on mytilus galloprovincialis after individual and combined exposure with carbamazepine. Sci. Total Environ. 643, 775–784. https://doi.org/10.1016/j.scitotenv.2018.06.257 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Ciacci, C. et al. Nanoparticle-biological interactions in a marine benthic foraminifer. Sci. Rep. https://doi.org/10.1038/s41598-019-56037-2 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Kim, J. A. et al. Low dose of amino-modified nanoparticles induces cell cycle arrest. ACS Nano 7(9), 7483–7494 (2013).Article 
CAS 
PubMed 

Google Scholar 
Mao, Y. et al. Phytoplankton response to polystyrene microplastics: Perspective from an entire growth period. Chemosphere https://doi.org/10.1016/j.chemosphere.2018.05.170 (2018).Article 
PubMed 

Google Scholar 
Wang, F. et al. Time resolved study of cell death mechanisms induced by amine-modified polystyrene nanoparticles. Nanoscale 5(22), 10868–10876. https://doi.org/10.1039/c3nr03249c (2013).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Xia, T. et al. Comparison of the abilities of ambient and manufactured nanoparticles to induce cellular toxicity according to an oxidative stress paradigm. Nano Lett. 6(8), 1794–1807. https://doi.org/10.1021/nl061025k (2006).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Lagarde, F. et al. Microplastic interactions with freshwater microalgae: Hetero-aggregation and changes in plastic density appear strongly dependent on polymer type. Environ. Pollut. 215, 331–339. https://doi.org/10.1016/j.envpol.2016.05.006 (2016).Article 
CAS 
PubMed 

Google Scholar 
Bhattacharya, P., Lin, S., Turner, J. P. & Ke, P. C. Physical adsorption of charged plastic nanoparticles affects algal photosynthesis. J. Phys. Chem. C 114(39), 16556–16561. https://doi.org/10.1021/jp1054759 (2010).Article 
CAS 

Google Scholar 
Johansen, J. L., Rønn, R. & Ekelund, F. Toxicity of cadmium and zinc to small soil protists. Environ. Pollut. 242, 1510–1517. https://doi.org/10.1016/j.envpol.2018.08.034 (2018).Article 
CAS 
PubMed 

Google Scholar 
Díaz, S., Martín-González, A. & Carlos Gutiérrez, J. Evaluation of heavy metal acute toxicity and bioaccumulation in soil ciliated protozoa. Environ. Int. 32(6), 711–717. https://doi.org/10.1016/j.envint.2006.03.004 (2006).Article 
CAS 
PubMed 

Google Scholar 
Subba, P. et al. Zinc stress induces physiological, ultra-structural and biochemical changes in mandarin orange (Citrus Reticulata Blanco) seedlings. Physiol. Mol. Biol. Plants 20(4), 461–473. https://doi.org/10.1007/s12298-014-0254-2 (2014).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Corcoll, N. et al. The effect of metals on photosynthesis processes and diatom metrics of biofilm from a metal-contaminated river: A translocation experiment. Ecol. Indic. 18, 620–631. https://doi.org/10.1016/j.ecolind.2012.01.026 (2012).Article 
CAS 

Google Scholar 
Moffett, B. F. et al. Zinc contamination decreases the bacterial diversity of agricultural soil. FEMS Microbiol. Ecol. 43(1), 13–19. https://doi.org/10.1016/S0168-6496(02)00448-8 (2003).Article 
CAS 
PubMed 

Google Scholar 
Kuperman, R. G. & Carreiro, M. M. Soil heavy metal concentrations, microbial biomass and enzyme activities in a contaminated grassland ecosystem. Soil Biol. Biochem. 29(2), 179–190. https://doi.org/10.1016/S0038-0717(96)00297-0 (1997).Article 
CAS 

Google Scholar 
Masmoudi, S. et al. Cadmium, copper, sodium and zinc effects on diatoms: From heaven to hell-a review. Cryptogam Algol 34(2), 185–225. https://doi.org/10.7872/crya.v34.iss2.2013.185 (2013).Article 

Google Scholar 
Gadd, G. M. & de Rome, L. Biosorption of copper by fungal melanin. Appl. Microbiol. Biotechnol. 29(6), 610–617. https://doi.org/10.1007/BF00260993 (1988).Article 
CAS 

Google Scholar 
Khan, M. & Scullion, J. Effects of metal (Cd, Cu, Ni, Pb or Zn) enrichment of sewage-sludge on soil micro-organisms and their activities. Appl. Soil. Ecol. 20(2), 145–155. https://doi.org/10.1016/S0929-1393(02)00018-5 (2002).Article 

Google Scholar 
Guillard, R. R. L. & Lorenzen, C. J. Yellow-green algae with chlorophyllide C12. J. Phycol. 8(1), 10–14. https://doi.org/10.1111/j.1529-8817.1972.tb03995.x (1972).Article 
CAS 

Google Scholar 
Zagata, P., Kopańska, M., Greczek-Stachura, M. & Burnecki, T. Acute toxicity of metals: Nickel and zinc to Paramecium bursaria and its endosymbionts. J. Microbiol. Biotechnol. Food Sci. 04, 128–131. https://doi.org/10.15414/jmbfs.2015.4.special2.128-131 (2015).Article 
CAS 

Google Scholar 
Lenz, R., Enders, K. & Nielsen, T. G. Microplastic exposure studies should be environmentally realistic. Proc. Natl. Acad. Sci. U. S. A. 113(29), E4121–E4122. https://doi.org/10.1073/pnas.1606615113 (2016).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Schertzinger, G., Ruchter, N. & Sures, B. Metal accumulation in sediments and amphipods downstream of combined sewer overflows. Sci. Total Environ. 616–617, 1199–1207. https://doi.org/10.1016/j.scitotenv.2017.10.199 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Erasmus, J. H. et al. Metal accumulation in riverine macroinvertebrates from a platinum mining region. Sci. Total Environ. https://doi.org/10.1016/j.scitotenv.2019.134738 (2020).Article 
PubMed 

Google Scholar 
Pradhan, S., Hedberg, J., Blomberg, E., Wold, S. & Odnevall Wallinder, I. Effect of sonication on particle dispersion, administered dose and metal release of non-functionalized, non-inert metal nanoparticles. J. Nanopart. Res. 18(9), 1–14. https://doi.org/10.1007/s11051-016-3597-5 (2016).Article 
CAS 

Google Scholar 
Taurozzi, J. S., Hackley, V. A. & Wiesner, M. R. Preparation of nanoparticle dispersions from powdered material using ultrasonic disruption. NIST Spec. Publ. 1200–2, 1–15 (2012).
Google Scholar 
Graupner, N. et al. Effects of short-term flooding on aquatic and terrestrial microeukaryotic communities: A mesocosm approach. Aquat. Microb. Ecol. 80(3), 257–272. https://doi.org/10.3354/ame01853 (2017).Article 

Google Scholar 
Strasser, R., Srivastava, A. & Tsimilli-Michael, M. The fluorescence transient as a tool to characterize and screen photosynthetic samples. In Probing Photosynthesis Mechanisms, Regulation and Adaption (eds Yanus, M. et al.) (Taylor and Francis, 2020).
Google Scholar 
Thwe, A. & Kasemsap, P. Quantification of OJIP fluorescence transient in tomato plants under acute ozone stress (2015).Amaral-Zettler, L. A., McCliment, E. A., Ducklow, H. W. & Huse, S. M. A method for studying protistan diversity using massively parallel sequencing of V9 hypervariable regions of small-subunit ribosomal RNA genes. PLoS ONE 4(7), 1–9. https://doi.org/10.1371/journal.pone.0006372 (2009).Article 
CAS 

Google Scholar 
Medlin, L., Elwood, H. J., Stickel, S. & Sogin, M. L. The characterization of enzymatically amplified eukaryotic 16S-like RRNA-coding regions. Gene 71(2), 491–499. https://doi.org/10.1016/0378-1119(88)90066-2 (1988).Article 
CAS 
PubMed 

Google Scholar 
Andrews, S. FastQC: A quality control tool for high throughput sequence data (2015).Lange, A. et al. AmpliconDuo: A split-sample filtering protocol for high-throughput amplicon sequencing of microbial communities. PLoS ONE 10(11), 1–22. https://doi.org/10.1371/journal.pone.0141590 (2015).Article 
CAS 

Google Scholar 
Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27(6), 863–864. https://doi.org/10.1093/bioinformatics/btr026 (2011).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Masella, P. A., Bartram, A. K., Truszkowski, J. M., Brow, D. G. & Neufeld, J. D. PANDAseq: Paired-end assembler for illumina sequences. BMC Bioinform. https://doi.org/10.1186/1471-2105-13-31 (2012).Article 

Google Scholar 
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27(16), 2194–2200. https://doi.org/10.1093/bioinformatics/btr381 (2011).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Mahé, F., Rognes, T., Quince, C., de Vargas, C. & Dunthorn, M. Swarm: Robust and fast clustering method for amplicon-based studies. PeerJ 2014(1), 1–13. https://doi.org/10.7717/peerj.593 (2014).Article 

Google Scholar 
Callahan, B. J. et al. DADA2: High-resolution sample inference from illumina amplicon data. Nat. Methods 13(7), 581–583. https://doi.org/10.1038/nmeth.3869 (2016).Article 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Welzel, M. et al. Natrix: A snakemake-based workflow for processing, clustering, and taxonomically assigning amplicon sequencing reads. BMC Bioinform. 21(1), 1–14. https://doi.org/10.1186/s12859-020-03852-4 (2020).Article 
CAS 

Google Scholar 
Oksanen, J. Package “vegan” Title Community Ecology Package (2022).R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.r-project.org/.Chen, W., Simpson, J. & Leveque, C. RAM: R for amplicon-sequencing-based microbial-ecology (2018).Yarza, P. et al. Uniting the classification of cultured and uncultured bacteria and archaea using 16S RRNA gene sequences. Nat. Rev. Microbiol. 12(9), 635–645. https://doi.org/10.1038/nrmicro3330 (2014).Article 
CAS 
PubMed 

Google Scholar 
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15(12), 1–21. https://doi.org/10.1186/s13059-014-0550-8 (2014).Article 
CAS 

Google Scholar 
Palarea-Albaladejo, J. & Martín-Fernández, J. A. ZCompositions—R package for multivariate imputation of left-censored data under a compositional approach. Chemom. Intell. Lab. Syst. 143, 85–96. https://doi.org/10.1016/j.chemolab.2015.02.019 (2015).Article 
CAS 

Google Scholar 
Gloor, G. B., Macklaim, J. M., Pawlowsky-Glahn, V. & Egozcue, J. J. Microbiome datasets are compositional: And this is not optional. Front. Microbiol. 8, 1–6. https://doi.org/10.3389/fmicb.2017.02224 (2017).Article 

Google Scholar 
Dusaucy, J., Gateuille, D., Perrette, Y. & Naffrechoux, E. Microplastic pollution of worldwide lakes. Environ. Pollut. 284, 117075. https://doi.org/10.1016/j.envpol.2021.117075 (2021).Article 
CAS 
PubMed 

Google Scholar 
Vardhan, K. H., Kumar, P. S. & Panda, R. C. A review on heavy metal pollution, toxicity and remedial measures: Current trends and future perspectives. J. Mol. Liq. 290, 111197. https://doi.org/10.1016/j.molliq.2019.111197 (2019).Article 
CAS 

Google Scholar 
Damare, V. S. Diversity of thraustochytrid protists isolated from brown alga, Sargassum cinereum using 18S RDNA sequencing and their morphological response to heavy metals. J. Mar. Biol. Assoc. 95(2), 265–276. https://doi.org/10.1017/S0025315414001696 (2015).Article 
CAS 

Google Scholar 
Giongo, A. et al. Adaption of Microbial communities to the hostile environment in the Doce river after the collapse of two iron ore tailing dams. Heliyon https://doi.org/10.1016/j.heliyon.2020.e04778 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Kelly, J. J., Häggblom, M. M. & Tate, R. L. Effects of heavy metal contamination and remediation on soil microbial communities in the vicinity of a zinc smelter as indicated by analysis of microbial community phospholipid fatty acid profiles. Biol. Fertil. Soils 38(2), 65–71. https://doi.org/10.1007/s00374-003-0642-1 (2003).Article 
CAS 

Google Scholar 
Baddar, Z. E., Peck, E. & Xu, X. Temporal deposition of copper and zinc in the sediments of metal removal constructed wetlands. PLoS ONE 16, 1–14. https://doi.org/10.1371/journal.pone.0255527 (2021).Article 
CAS 

Google Scholar 
Li, X., Shen, Z., Wai, O. W. H. & Li, Y. S. Chemical partitioning of heavy metal contaminants in sediments of the Pearl River Estuary. Chem. Speciat. Bioavailab. 12(1), 17–25. https://doi.org/10.3184/095422900782775607 (2000).Article 
CAS 

Google Scholar 
Müller, B. & Sigg, L. Interaction of trace metals with natural particle surfaces: Comparison between adsorption experiments and field measurements—Dedicated to Werner Stumm for his 65th birthday. Aquat. Sci. 52(1), 75–92. https://doi.org/10.1007/BF00878242 (1990).Article 

Google Scholar 
Bradl, H. B. Adsorption of heavy metal ions on soils and soils constituents. J. Colloid Interface Sci. 277(1), 1–18. https://doi.org/10.1016/j.jcis.2004.04.005 (2004).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Siegel, F. R. Environmental Geochemistry of Potentially Toxic Heavy Metals (Springer-Verlag, 2002).Book 

Google Scholar 
Vig, K., Megharaj, M., Sethunathan, N. & Naidu, R. Bioavailability and toxicity of cadmium to microorganisms and their activities in soil: A review. Adv. Environ. Res. 8(1), 121–135. https://doi.org/10.1016/S1093-0191(02)00135-1 (2003).Article 
CAS 

Google Scholar 
Nicolau, A., Mota, M. & Lima, N. Physiological responses of tetrahymena pyriformis to copper, zinc, cycloheximide and triton X-100. FEMS Microbiol. Ecol. 30(3), 209–216. https://doi.org/10.1016/S0168-6496(99)00057-4 (1999).Article 
CAS 
PubMed 

Google Scholar 
Admiraal, W. et al. Short-term toxicity of zinc to microbenthic algae and bacteria in a metal polluted stream. Water Res. 33(9), 1989–1996. https://doi.org/10.1016/S0043-1354(98)00426-6 (1999).Article 
CAS 

Google Scholar 
Bradac, P., Navarro, E., Odzak, N., Behra, R. & Sigg, L. Kinetics of cadmium accumulation in periphyton under freshwater conditions. Environ. Toxicol. Chem. 28(10), 2108–2116. https://doi.org/10.1897/08-511R1.1 (2009).Article 
CAS 
PubMed 

Google Scholar 
Collard, J. & Matagne, R. F. Cd2+ resistance in wild-type and mutant strains of Chlamydomonas reinhardtii. Environ. Exp. Bot. 34(2), 235–244 (1994).Article 
CAS 

Google Scholar 
Wright, R. J., Gibson, M. I. & Christie-Oleza, J. A. Understanding microbial community dynamics to improve optimal microbiome selection. Microbiome 7(1), 1–14. https://doi.org/10.1186/s40168-019-0702-x (2019).Article 

Google Scholar 
Buffle, J. The key role of environmental colloids/nanoparticles for the sustainability of life. Environ. Chem. 3(3), 155–158. https://doi.org/10.1071/ENv3n3_ES (2006).Article 
CAS 

Google Scholar 
Nowack, B. & Bucheli, T. D. Occurrence, behavior and effects of nanoparticles in the environment. Environ. Pollut. 150(1), 5–22. https://doi.org/10.1016/j.envpol.2007.06.006 (2007).Article 
CAS 
PubMed 

Google Scholar 
Fetzer, I. et al. The extent of functional redundancy changes as species’ roles shift in different environments. Proc. Natl. Acad. Sci. U. S. A. 112(48), 14888–14893. https://doi.org/10.1073/pnas.1505587112 (2015).Article 
ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Biggs, C. R. et al. Does functional redundancy affect ecological stability and resilience? A review and meta-analysis. Ecosphere https://doi.org/10.1002/ecs2.3184 (2020).Article 

Google Scholar 
Fleeger, J. W. How do indirect effects of contaminants inform ecotoxicology? A review. Processes https://doi.org/10.3390/pr8121659 (2020).Article 

Google Scholar 
Oriekhova, O. & Stoll, S. Heteroaggregation of nanoplastic particles in the presence of inorganic colloids and natural organic matter. Environ. Sci. Nano. 5(3), 792–799. https://doi.org/10.1039/c7en01119a (2018).Article 
CAS 

Google Scholar 
Rowenczyk, L. et al. Heteroaggregates of polystyrene nanospheres and organic matter: Preparation, characterization and evaluation of their toxicity to algae in environmentally relevant conditions. Nanomaterials 11(2), 1–15. https://doi.org/10.3390/nano11020482 (2021).Article 
CAS 

Google Scholar 
Saavedra, J., Stoll, S. & Slaveykova, V. I. Influence of nanoplastic surface charge on eco-corona formation, aggregation and toxicity to freshwater zooplankton. Environ. Pollut. 252, 715–722. https://doi.org/10.1016/j.envpol.2019.05.135 (2019).Article 
CAS 
PubMed 

Google Scholar 
Bižic-Ionescu, M., Ionescu, D. & Grossart, H. P. Organic particles: Heterogeneous hubs for microbial interactions in aquatic ecosystems. Front. Microbiol. 9, 1–15. https://doi.org/10.3389/fmicb.2018.02569 (2018).Article 

Google Scholar 
Lespes, G., Faucher, S. & Slaveykova, V. I. natural nanoparticles, anthropogenic nanoparticles, where is the Frontier?. Front. Environ. Sci. 8, 1–5. https://doi.org/10.3389/fenvs.2020.00071 (2020).Article 

Google Scholar 
Stabnikova, O. et al. Microbial life on the surface of microplastics in natural waters. Appl. Sci. 11(24), 1–19. https://doi.org/10.3390/app112411692 (2021).Article 
CAS 

Google Scholar 
Suominen, S., Doorenspleet, K., Sinninghe Damsté, J. S. & Villanueva, L. Microbial community development on model particles in the deep sulfidic waters of the Black Sea. Environ. Microbiol. 23(6), 2729–2746. https://doi.org/10.1111/1462-2920.15024 (2021).Article 
CAS 
PubMed 

Google Scholar 
Wagner, S., Gondikas, A., Neubauer, E., Hofmann, T. & von der Kammer, F. Spot the difference: Engineered and natural nanoparticles in the environment-release, behavior, and fate. Angew. Chem. Int. Ed. 53(46), 12398–12419. https://doi.org/10.1002/anie.201405050 (2014).Article 
CAS 

Google Scholar 
Amelia, T. S. et al. Marine microplastics as vectors of major ocean pollutants and its hazards to the marine ecosystem and humans. Prog. Earth Planet. Sci. https://doi.org/10.1186/s40645-020-00405-4 (2021).Article 

Google Scholar 
Liu, J., Huang, J. & Che, F. Microalgae as feedstocks for biodiesel production. In Biodiesel—Feedstocks and Processing Technologies (ed. Stoytcheva, M.) (InTech, 2011). https://doi.org/10.5772/25600.Chapter 

Google Scholar 
Takamura, N., Kasai, F. & Watanabe, M. M. Effects of Cu, Cd and Zn on photosynthesis of freshwater benthic algae. J. Appl. Phycol. 1(1), 39–52. https://doi.org/10.1007/BF00003534 (1989).Article 
CAS 

Google Scholar 
Brembu, T., Jørstad, M., Winge, P., Valle, K. C. & Bones, A. M. Genome-wide profiling of responses to cadmium in the diatom Phaeodactylum tricornutum. Environ. Sci. Technol. 45(18), 7640–7647. https://doi.org/10.1021/es2002259 (2011).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Fernandez, J. C. & Henriques, F. S. Biochemical, physiological and structural effects of excess copper in plants. Bot. Rev. 57(3), 246–273 (1991).Article 

Google Scholar 
Haq, R. U., Rehman, A. & Shakoori, A. R. Effect of dichromate on population and growth of various protozoa isolated from industrial effluents. Folia Microbiol. 45(3), 275–278. https://doi.org/10.1007/bf02908959 (2000).Article 
CAS 

Google Scholar 
Rehman, A., Shakoori, F. R. & Shakoori, A. R. Heavy metal resistant freshwater ciliate, Euplotes mutabilis, isolated from industrial effluents has potential to decontaminate wastewater of toxic metals. Bioresour. Technol. 99(9), 3890–3895. https://doi.org/10.1016/j.biortech.2007.08.007 (2008).Article 
CAS 
PubMed 

Google Scholar 
Rehman, A., Ashraf, S., Qazi, J. I. & Shakoori, A. R. Uptake of lead by a ciliate, stylonychia mytilus, isolated from industrial effluents: Potential use in bioremediation of wastewater. Bull. Environ. Contam. Toxicol. 75(2), 290–296. https://doi.org/10.1007/s00128-005-0751-7 (2005).Article 
CAS 
PubMed 

Google Scholar 
Shakoori, A. R., Rehman, A. & ul-Haq, R. Multiple metal resistance in the ciliate protozoan, vorticella microstoma, isolated from industrial effluents and its potential in bioremediation of toxic wastes. Bull. Environ. Contam. Toxicol. 72(5), 1046–1051. https://doi.org/10.1007/s00128-004-0349-5 (2004).Article 
CAS 
PubMed 

Google Scholar 
Falasco, E. et al. Morphological abnormalities of diatom silica walls in relation to heavy metal contamination and artificial growth conditions. Water SA 35(5), 595–606. https://doi.org/10.4314/wsa.v35i5.49185 (2009).Article 
CAS 

Google Scholar 
Tadros, M. G., Mbuthia, P. & Smith, W. Differential response of marine diatoms to trace metals. Bull. Environ. Contam. Toxicol. 44(6), 826–831. https://doi.org/10.1007/BF01702170 (1990).Article 
CAS 
PubMed 

Google Scholar 
Wanner, M. et al. Soil testate amoebae and diatoms as bioindicators of an old heavy metal contaminated floodplain in Japan. Microb. Ecol. 79(1), 123–133. https://doi.org/10.1007/s00248-019-01383-x (2020).Article 
CAS 
PubMed 

Google Scholar 
Shi, J., Podola, B. & Melkonian, M. Application of a prototype-scale twin-layer photobioreactor for effective N and P removal from different process stages of municipal wastewater by immobilized microalgae. Bioresour. Technol. 154, 260–266. https://doi.org/10.1016/j.biortech.2013.11.100 (2014).Article 
CAS 
PubMed 

Google Scholar 
Li, T., Lin, G., Podola, B. & Melkonian, M. Continuous removal of zinc from wastewater and mine dump leachate by a microalgal biofilm PSBR. J. Hazard. Mater. 297, 112–118. https://doi.org/10.1016/j.jhazmat.2015.04.080 (2015).Article 
CAS 
PubMed 

Google Scholar 
Bruning, K. Infection of the diatom Asterionella by a chytrid. I. Effects of light on reproduction and infectivity of the parasite. J. Plankton Res. 13(1), 103–117. https://doi.org/10.1093/plankt/13.1.103 (1991).Article 

Google Scholar 
Carney, L. T. & Lane, T. W. Parasites in algae mass culture. Front. Microbiol. 5, 1–8. https://doi.org/10.3389/fmicb.2014.00278 (2014).Article 

Google Scholar 
Hanic, L. A., Sekimoto, S. & Bates, S. S. Oomycete and chytrid infections of the marine diatom Pseudo-nitzschia pungens (Bacillariophyceae) from Prince Edward Island. Botany 87(11), 1096–1105. https://doi.org/10.1139/B09-070 (2009).Article 
CAS 

Google Scholar 
Sun, A. et al. Fertilization alters protistan consumers and parasites in crop-associated microbiomes. Environ. Microbiol. 23(4), 2169–2183. https://doi.org/10.1111/1462-2920.15385 (2021).Article 
CAS 
PubMed 

Google Scholar 
Scholz, B., Guillou, L., Marano, A. V., Neuhauser, S. & Brooke, K. Europe PMC funders group zoosporic parasites infecting marine diatoms—A black box that needs to be opened. Fungal Ecol. https://doi.org/10.1016/j.funeco.2015.09.002.Zoosporic (2017).Article 

Google Scholar 
Peacock, E. E., Olson, R. J. & Sosik, H. M. Parasitic infection of the diatom Guinardia delicatula, a recurrent and ecologically important phenomenon on the New England Shelf. Mar. Ecol. Prog. Ser. 503, 1–10. https://doi.org/10.3354/meps10784 (2014).Article 
ADS 

Google Scholar 
Duarte, S., Pascoal, C. & Cássio, F. Effects of zinc on leaf decomposition by fungi in streams: Studies in microcosms. Microb. Ecol. 48(3), 366–374. https://doi.org/10.1007/s00248-003-2032-5 (2004).Article 
CAS 
PubMed 

Google Scholar 
Kammerlander, B. et al. High diversity of protistan plankton communities in remote high mountain lakes in the European Alps and the Himalayan Mountains. FEMS Microbiol. Ecol. 91(4), 1–10. https://doi.org/10.1093/femsec/fiv010 (2015).Article 
CAS 

Google Scholar 
Sieber, G., Beisser, D., Bock, C. & Boenigk, J. Protistan and fungal diversity in soils and freshwater lakes are substantially different. Sci. Rep. 10(1), 1–11. https://doi.org/10.1038/s41598-020-77045-7 (2020).Article 
CAS 

Google Scholar 
Gunaalan, K., Fabbri, E. & Capolupo, M. The hidden threat of plastic leachates: A critical review on their impacts on aquatic organisms. Water Res. https://doi.org/10.1016/j.watres.2020.116170 (2020).Article 
PubMed 

Google Scholar 
Tetu, S. G., Sarker, I., Schrameyer, V., Pickford, R., Elbourne, L. D., Moore, L.R. & Paulsen, I.T. Plastic leachates impair growth and oxygen production in Prochlorococcus, the ocean’s most abundant photosynthetic bacteria. Commun. Biol. 2(1), 1–9. https://doi.org/10.1038/s42003-019-0410-x (2019).Gouin, T., Roche, N., Lohmann, R. & Hodges, G. A Thermodynamic approach for assessing the environmental exposure of chemicals absorbed to microplastic. Environ. Sci. Technol. 45(4), 1466–1472. https://doi.org/10.1021/es1032025 (2011).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Lohmann, R. Microplastics are not important for the cycling and bioaccumulation of organic pollutants in the oceans—But should microplastics be considered POPs themselves?. Integr. Environ. Assess. Manag. 13(3), 460–465. https://doi.org/10.1002/ieam.1914 (2017).Article 
CAS 
PubMed 

Google Scholar 
Sukkasem, C. & Laehlah, S. An economical upflow bio-filter circuit (UBFC): A biocatalyst microbial fuel cell for sulfate-sulfide rich wastewater treatment. Environ. Sci. 1(2), 161–168. https://doi.org/10.1039/c4ew00028e (2015).Article 
CAS 

Google Scholar 
Abatenh, E., Gizaw, B., Tsegaye, Z. & Wassie, M. The role of microorganisms in bioremediation-A review. Open J. Environ. Biol. 2(1), 38–46. https://doi.org/10.17352/ojeb (2017).Article 

Google Scholar 
Zrimec, J., Kokina, M., Jonasson, S., Zorrilla, F. & Zelezniak, A. Plastic-degrading potential across the global microbiome correlates with recent pollution trends. MBio https://doi.org/10.1128/mBio (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Siver, P. A. Synurophyte algae. In Freshwater Algae of North America. Ecology and classification (eds Wehr, J. D. & Sheath, R. G.) 523–558 (Elsevier, 2003).Chapter 

Google Scholar 
Andersen, R. A. Molecular systematics of the chrysophyceae and synurophyceae. In Unravelling the Algae: The Past, Present, and Future of Algal Systematics (eds Brodie, J. & Lewis, J.) 285–314 (CRC Press, Boca Raton, 2007).Chapter 

Google Scholar 
Engin, I. K., Cekmecelioglu, D., Yücel, A. M. & Oktem, H. A. Evaluation of heterotrophic and mixotrophic cultivation of novel Micractinium Sp. ME05 on vinasse and its scale up for biodiesel production. Bioresour. Technol. 251, 128–134. https://doi.org/10.1016/j.biortech.2017.12.023 (2018).Article 
CAS 
PubMed 

Google Scholar 
Patrick, R. Ecology of freshwater diatoms and diatom communities. In The Biology of Diatoms (ed. Werner, D.) 284–332 (University of California Press, 1977).
Google Scholar 
Findenig, B. M., Chatzinotas, A. & Boenigk, J. Taxonomic and ecological characterization of stomatocysts of spumella-like flagellates (Chrysophyceae). J. Phycol. 46(5), 868–881. https://doi.org/10.1111/j.1529-8817.2010.00892.x (2010).Article 

Google Scholar 
Perez-Garcia, O., Escalante, F. M. E., de-Bashan, L. E. & Bashan, Y. Heterotrophic cultures of microalgae: Metabolism and potential products. Water Res. 45(1), 11–36. https://doi.org/10.1016/j.watres.2010.08.037 (2011).Article 
CAS 
PubMed 

Google Scholar 
Preisig, H. R. & Hibberd, D. J. Ultrastructure and taxonomy of Paraphysomonas (Chrysophyceae) and related genera 3. Nord. J. Bot. 3(6), 695–723. https://doi.org/10.1111/j.1756-1051.1983.tb01481.x (1983).Article 

Google Scholar 
Atkins, M. S. et al. Tolerance of flagellated protists to high sulfide and metal concentrations potentially encountered at deep-sea hydrothermal vents. Mar. Ecol. Prog. Ser. 226, 63–75. https://doi.org/10.3354/meps226063 (2002).Article 
ADS 
CAS 

Google Scholar 
Manru, G., Weisong, F. & Yunfen, S. Ecological study on protozoa in the sediment of the three-gorges area of the Changjiang River. Chin. J. Oceanol. Limnol. 6(3), 272–280. https://doi.org/10.1007/BF02846505 (1988).Article 

Google Scholar 
Tomilina, I. I., Gremyachikh, V. A., Myl’Nikov, A. P. & Komov, V. T. The effect of metal oxide nanoparticles (CeO2, TiO2, and ZnO) on biological parameters of freshwater nanoflagellates and crustaceans. Dokl. Biol. Sci. 436(1), 53–55. https://doi.org/10.1134/S0012496611010169 (2011).Article 
CAS 
PubMed 

Google Scholar 
Schampera, C. et al. Exposure to nanoplastics affects the outcome of infectious disease in phytoplankton. Environ. Pollut. https://doi.org/10.1016/j.envpol.2021.116781 (2021).Article 
PubMed 

Google Scholar 
Gonçalves, J. M., Sousa, V. S., Teixeira, M. R. & Bebianno, M. J. Chronic toxicity of polystyrene nanoparticles in the marine mussel Mytilus galloprovincialis. Chemosphere https://doi.org/10.1016/j.chemosphere.2021.132356 (2021).Article 
PubMed 

Google Scholar 
Kelpsiene, E., Torstensson, O., Ekvall, M. T., Hansson, L. A. & Cedervall, T. Long-term exposure to nanoplastics reduces life-time in Daphnia magna. Sci. Rep. 10(1), 1–7. https://doi.org/10.1038/s41598-020-63028-1 (2020).Article 
CAS 

Google Scholar 
Amin, N. M. Techniques for assessment of heavy metal toxicity using Acanthamoeba Sp, a small, naked and free-living amoeba. Funct. Ecosyst. https://doi.org/10.5772/36008 (2012).Article 

Google Scholar 
Amin, N. M., Azhar, N. & Shazili, M. Cytotoxic effects of mercury, cadmium, lead and zinc on Acanthamoeba Castellanii (2006).Gnecco, I., Berretta, C., Lanza, L. G. & la Barbera, P. Storm water pollution in the urban environment of Genoa, Italy. Atmos. Res. 77, 60–73. https://doi.org/10.1016/j.atmosres.2004.10.017 (2005).Article 
CAS 

Google Scholar 
Heim, R. R. An overview of weather and climate extremes—Products and trends. Weather Clim. Extrem. 10, 1–9. https://doi.org/10.1016/j.wace.2015.11.001 (2015).Article 

Google Scholar 
Saiki, M. K., Castleberry, D. T., May, T. W., Martin, B. A. & Bullard, F. N. Copper, cadmium, and zinc concentrations in aquatic food chains from the upper Sacramento River (California) and selected tributaries. Arch. Environ. Contam. Toxicol. 29(4), 484–491. https://doi.org/10.1007/BF00208378 (1995).Article 
CAS 

Google Scholar 
Wagner, S. et al. Tire wear particles in the aquatic environment—A review on generation, analysis, occurrence, fate and effects. Water Res. 139, 83–100. https://doi.org/10.1016/j.watres.2018.03.051 (2018).Article 
CAS 
PubMed 

Google Scholar 
Zhang, L., Zhao, B., Xu, G. & Guan, Y. Characterizing fluvial heavy metal pollutions under different rainfall conditions: Implication for aquatic environment protection. Sci. Total Environ. 635, 1495–1506. https://doi.org/10.1016/j.scitotenv.2018.04.211 (2018).Article 
ADS 
CAS 
PubMed 

Google Scholar 
Zhao, B. et al. Characterization of nitrosamines and nitrosamine precursors as non-point source pollutants during heavy rainfall events in an urban water environment. J. Hazard. Mater. 424, 127552. https://doi.org/10.1016/j.jhazmat.2021.127552 (2022).Article 
CAS 
PubMed 

Google Scholar 
Hüffer, T., Wagner, S., Reemtsma, T. & Hofmann, T. Sorption of organic substances to tire wear materials: Similarities and differences with other types of microplastic. TrAC Trends Anal. Chem. 113, 392–401. https://doi.org/10.1016/j.trac.2018.11.029 (2019).Article 
CAS 

Google Scholar 
Tamis, J. E. et al. Environmental risks of car tire microplastic particles and other road runoff pollutants. Microplastics Nanoplastics 1(1), 1–17. https://doi.org/10.1186/s43591-021-00008-w (2021).Article 

Google Scholar 
Chèvre, N. et al. Substance flow analysis as a tool for urban water management. Water Sci. Technol. 63(7), 1341–1348. https://doi.org/10.2166/wst.2011.132 (2011).Article 
PubMed 

Google Scholar 
Šourková, M., Adamcová, D. & Vaverková, M. D. The influence of microplastics from ground tyres on the acute, subchronical toxicity and microbial respiration of soil. Environ. MDPI 8(11), 1–14. https://doi.org/10.3390/environments8110128 (2021).Article 

Google Scholar 
Ye, G., Zhang, X., Yan, C., Lin, Y. & Huang, Q. Polystyrene microplastics induce microbial dysbiosis and dysfunction in surrounding seawater. Environ. Int. 156, 106724. https://doi.org/10.1016/j.envint.2021.106724 (2021).Article 
CAS 
PubMed 

Google Scholar  More

Sites overview and characteristicsThis study focused on three regions located in subantarctic, arctic, and subarctic latitudes. The respective latitudinal and longitudinal ranges covered in this study were: 54.95 to 52.08 °S, and 72.03 to 67.34 °W in Patagonia; 67.44 to 67.54 °N, and 86.59 to 86.71 °E in Siberia; 63.21 to 68.63 °N, and −150.79 to −145.98 °W in Alaska (Figs. 1 and 2). The exact coordinates for each sample were included in the submitted dataset. The field campaigns were conducted in 2016, during the summer for each respective region: January-February in Chilean Patagonia, June-July in Alaska and July-August in Siberia.Fig. 1Location of the three areas included in this study (panel a). The permafrost state and the number of sites and samples per region is indicated for each area. General views of 5 sites are provided as examples (b–f). Panel B provides a large view of the ecosystem surrounding the wetland ALP2 (Alaska, exact location indicated by the white circle). Lake PCL1 (panel c) is representative of the lakes on Navarino island (Chilean Patagonia). The glacial lake SIL2 is shown in panel d. At site SIP5, the hollow at first plan is surrounded by palsa (hummock, second plan), characterized by dark organic matter and lichen vegetation (panel e). The PPP3 peatland shown in panel f is dominated by Sphagnum magellanicum, like most peatlands in the area.Full size imageFig. 2Maps of sampling sites in Patagonia, Alaska and Siberia, indicating the ecosystem type (lake, wetland, soil). The tables show the complete- (in white) and the partial- (in grey) characterization sites. The exact coordinates of each sample are provided in the data record (See data records section).Full size imageFor every site included in the present study, a set of nine qualitative environmental and/or ecological site-scale descriptors was selected and adapted from ENVO Environment Ontology40, which included for example permafrost state, biome, environmental feature and vegetation type (Table 1, Fig. 3). Permafrost state was obtained from the NSIDC permafrost map41. The biome, large-scale descriptor based on climate and vegetation criteria, was derived from Olson et al.42. Temperate forest, boreal forest, and tundra biomes were included. The environmental features that were representative for the three regions were considered: lakes, wetlands, broadleaf/coniferous/mixed forest soils, grassland, tundra, and palsa. All the metadata was included in the submitted dataset. Table 2 summarizes the main types of sampled ecosystems and their main characteristics in the three regions, while Supplementary Table S1 provides the details of each sampling site.Table 1 Overview of the dataset contained in Mimarks sheet.Full size tableFig. 3Description of the qualitative environmental/ecological descriptors used to describe every sample, derived from ENVO Environment Ontology40.Full size imageTable 2 Main types of sampled ecosystems in the three studied regions.Full size tableIn Alaska, the studied area ranged from the Alaska Range and Fairbanks area (interior, continental climate, 63–65°N, discontinuous permafrost) up to Toolik Field Station (North Slope, arctic climate, 66–69°N, continuous permafrost; Fig. 2). The physiochemistry and CH4 emissions of lakes ALL1 (Killarney lake), ALL2 (Otto lake), ALL3 (Nutella lake), and ALL4 (Goldstream lake) were previously characterized35. A number of heterogeneous soil and wetland samples were collected around the studied Alaskan lakes and/or from monitored sites, as detailed in Supplementary Table S1. In the Alaska Range and Fairbanks area, soils were mostly covered by mixed or taiga forests, alpine tundra, and bogs or fens wetlands. In the norther Brooks Ranges mountain system, the landscape was piedmont hills with a predominant soil of porous organic peat underlain by silt and glacial till, all in a permafrost state, characterized mainly by Sphagnum and Eriophorum vegetation, as well as dwarf shrubs.In Siberia, the studied area was located in the discontinuous permafrost region surrounding Igarka, on the eastern bank of the Yenisei River (Fig. 2). This region was mainly covered by forest, dominated by larch (Larix Siberica), birch (Betula Pendula), and Siberian pine (Pinus Siberica), and palsa landscapes (frozen peat mounts), the latter being dominated by moss, lichens, Labrador tea and dwarf birch. In degraded areas, thermokarst bogs were dominated by Sphagnum spp. and Eriophorum spp. Land cover was an indicator of permafrost status, since forested areas reflected a deep permafrost table ( >2 m) associated with Pleistocene permafrost, while palsa-dominated landscapes were indicative of the presence of near-surface ( More