Philippa Kaur

Schneider-Crease, I. et al. Identifying wildlife reservoirs of neglected taeniid tapeworms: Non-invasive diagnosis of endemic Taenia serialis infection in a wild primate population. PLoS Negl Trop Dis 11, e0005709 (2017).Article 

Google Scholar 
Schneider-Crease, I. et al. Ecology eclipses phylogeny as a major driver of nematode parasite community structure in a graminivorous primate. Funct. Ecol. 34, 1898–1906 (2020).Article 

Google Scholar 
Gillespie, T. R. Noninvasive assessment of gastrointestinal parasite infections in free-ranging primates. Int. J. Primatol. 27, 1129 (2006).Article 

Google Scholar 
Budischak, S. A., Hoberg, E. P., Abrams, A., Jolles, A. E. & Ezenwa, V. O. A combined parasitological molecular approach for noninvasive characterization of parasitic nematode communities in wild hosts. Mol. Ecol. Resour. 15, 1112–1119 (2015).Article 

Google Scholar 
Ghalehnoei, H., Bagheri, A., Fakhar, M. & Mishan, M. A. Circulatory microRNAs: promising non-invasive prognostic and diagnostic biomarkers for parasitic infections. Eur. J. Clin. Microbiol. Infect. Dis. 39, 395–402 (2020).Article 
CAS 

Google Scholar 
Hing, S., Narayan, E. J., Andrew Thompson, R. C. & Godfrey, S. S. The relationship between physiological stress and wildlife disease: consequences for health and conservation. Wildl Res. 43, 51–60 (2016).Article 

Google Scholar 
Kersey, D. C. & Dehnhard, M. The use of noninvasive and minimally invasive methods in endocrinology for threatened mammalian species conservation. Gen. Comp. Endocrinol. 203, 296–306 (2014).Article 
CAS 

Google Scholar 
Behringer, V. & Deschner, T. Non-invasive monitoring of physiological markers in primates. Horm. Behav. 91, 3–18 (2017).Article 
CAS 

Google Scholar 
Higham, J. P., Stahl-Hennig, C. & Heistermann, M. Urinary suPAR: A non-invasive biomarker of infection and tissue inflammation for use in studies of large free-ranging mammals. R Soc Open Sci. 7, 191825 (2020).Article 
ADS 
CAS 

Google Scholar 
Heistermann, M. & Higham, J. P. Urinary neopterin, a non-invasive marker of mammalian cellular immune activation, is highly stable under field conditions. Sci Rep. 5, 16308 (2015).Article 
ADS 
CAS 

Google Scholar 
Higham, J. P. et al. Evaluating noninvasive markers of nonhuman primate immune activation and inflammation. Am. J. Phys. Anthropol. 158, 673–684 (2015).Article 

Google Scholar 
Behringer, V. et al. Elevated neopterin levels in wild, healthy chimpanzees indicate constant investment in unspecific immune system. BMC Zool. 4, 1–7 (2019).Article 

Google Scholar 
Dibakou, S. E., Basset, D., Souza, A., Charpentier, M. & Huchard, E. Determinants of variations in fecal neopterin in free-ranging mandrills. Front. Ecol. Evol. https://doi.org/10.3389/fevo.2019.00368 (2019).Article 

Google Scholar 
Löhrich, T., Behringer, V., Wittig, R. M., Deschner, T. & Leendertz, F. H. The use of neopterin as a noninvasive marker in monitoring diseases in wild chimpanzees. EcoHealth 15, 792–803 (2018).Article 

Google Scholar 
Behringer, V., Stevens, J. M. G., Leendertz, F. H., Hohmann, G. & Deschner, T. Validation of a method for the assessment of urinary neopterin levels to monitor health status in nonhuman primate species. Front Physiol. 8, 51 (2017).Article 

Google Scholar 
Negrey, J. D., Behringer, V., Langergraber, K. E. & Deschner, T. Urinary neopterin of wild chimpanzees indicates that cell-mediated immune activity varies by age, sex, and female reproductive status. Sci. Rep. 11, 9298 (2021).Article 
ADS 
CAS 

Google Scholar 
Behringer, V. et al. Cell-mediated immune ontogeny is affected by sex but not environmental context in a long-lived primate species. Front. Ecol. Evol. 9, 272 (2021).Article 

Google Scholar 
Müller, N., Heistermann, M., Strube, C., Schülke, O. & Ostner, J. Age, but not anthelmintic treatment, is associated with urinary neopterin levels in semi-free ranging Barbary macaques. Sci. Rep. 7, 41973 (2017).Article 
ADS 

Google Scholar 
Dibakou, S. E. et al. Ecological, parasitological and individual determinants of plasma neopterin levels in a natural mandrill population. Int. J. Parasitol. Parasites Wildl. 11, 198–206 (2020).Article 

Google Scholar 
Eisenhut, M. Neopterin in diagnosis and monitoring of infectious diseases. J. Biomark. 2013, 196432 (2013).Article 

Google Scholar 
Franceschi, C. & Campisi, J. Chronic inflammation (inflammaging) and its potential contribution to age-associated diseases. J Gerontol. A Biol. Sci. Med Sci. 69(Suppl 1), S4–9 (2014).Article 

Google Scholar 
Basha, S., Surendran, N. & Pichichero, M. Immune responses in neonates. Expert. Rev. Clin. Immunol. 10, 1171–1184 (2014).Article 
CAS 

Google Scholar 
Werner, E. R. et al. Determination of neopterin in serum and urine. Clin. Chem. 33, 62–66 (1987).Article 
CAS 

Google Scholar 
Lucore, J. M., Marshall, A. J., Brosnan, S. F. & Benítez, M. E. Validating urinary neopterin as a biomarker of immune response in captive and wild capuchin monkeys. Front. Vet. Sci. 9, 918036. https://doi.org/10.3389/fvets.2022.918036 (2022).Article 

Google Scholar 
Berdowska, A. & Zwirska-Korczala, K. Neopterin measurement in clinical diagnosis. J. Clin. Pharm. Ther. 26, 319–329 (2001).Article 
CAS 

Google Scholar 
Denz, H. et al. Value of urinary neopterin in the differential diagnosis of bacterial and viral infections. Klin. Wochenschr. 68, 218–222 (1990).Article 
CAS 

Google Scholar 
Reibnegger, G. et al. Urinary neopterin reflects clinical activity in patients with rheumatoid arthritis. Arthritis Rheum. 29, 1063–1070 (1986).Article 
CAS 

Google Scholar 
Huber, C. et al. Immune response-associated production of neopterin. Release from macrophages primarily under control of interferon-gamma. J. Exp. Med. 160, 310–316 (1984).Article 
CAS 

Google Scholar 
Horak, E., Gassner, I., Sölder, B., Wachter, H. & Fuchs, D. Neopterin levels and pulmonary tuberculosis in infants. Lung 176, 337–344 (1998).Article 
CAS 

Google Scholar 
Fendrich, C. et al. Urinary neopterin concentrations in rhesus monkeys after infection with simian immunodeficiency virus (SIVmac 251). AIDS 3, 305–307 (1989).Article 
CAS 

Google Scholar 
Chan, C. P. Y. et al. Detection of serum neopterin for early assessment of dengue virus infection. J. Infect. 53, 152–158 (2006).Article 

Google Scholar 
Wu, D. F., Behringer, V., Wittig, R. M., Leendertz, F. H. & Deschner, T. Urinary neopterin levels increase and predict survival during a respiratory outbreak in wild chimpanzees (Taï National Park, Côte d’Ivoire). Sci. Rep. 8, 13346 (2018).Article 
ADS 

Google Scholar 
Maizels, R. M. & McSorley, H. J. Regulation of the host immune system by helminth parasites. J. Allergy Clin. Immunol. 138, 666–675 (2016).Article 
CAS 

Google Scholar 
Maizels, R. M. & Yazdanbakhsh, M. Immune regulation by helminth parasites: Cellular and molecular mechanisms. Nat. Rev. Immunol. 3, 733–744 (2003).Article 
CAS 

Google Scholar 
Yazdanbakhsh, M., Kremsner, P. G. & van Ree, R. Allergy, parasites, and the hygiene hypothesis. Science 296, 490–494 (2002).Article 
ADS 
CAS 

Google Scholar 
Faz-López, B., Morales-Montor, J. & Terrazas, L. I. Role of macrophages in the repair process during the tissue migrating and resident helminth infections. Biomed. Res. Int. 2016, 8634603 (2016).Article 

Google Scholar 
Garcia, H. H., Rodriguez, S., Friedland, J. S. Cysticercosis Working Group in Peru. Immunology of Taenia solium taeniasis and human cysticercosis. Parasite Immunol. 36, 388–396. https://doi.org/10.1111/pim.12126 (2014)Article 
CAS 

Google Scholar 
Thaiss, C. A., Zmora, N., Levy, M. & Elinav, E. The microbiome and innate immunity. Nature 535, 65–74 (2016).Article 
ADS 
CAS 

Google Scholar 
Schneider-Crease, I. A., Griffin, R. H., Gomery, M. A., Bergman, T. J. & Beehner, J. C. High mortality associated with tapeworm parasitism in geladas (Theropithecus gelada) in the Simien Mountains National Park, Ethiopia. Am. J. Primatol. https://doi.org/10.1002/ajp.22684 (2017).Article 

Google Scholar 
Nguyen, N. et al. Fitness impacts of tapeworm parasitism on wild gelada monkeys at Guassa, Ethiopia. Am. J. Primatol. 77, 579–594 (2015).Article 

Google Scholar 
Schneider-Crease, I. A. et al. Helminth infection is associated with dampened cytokine responses to viral and bacterial stimulations in Tsimane forager-horticulturalists. Evol. Med. Public Health 9, 349–359 (2021).Article 

Google Scholar 
Roberts, E. K., Lu, A., Bergman, T. J. & Beehner, J. C. Female reproductive parameters in wild geladas (Theropithecus gelada). Int. J. Primatol. 38, 1–20 (2017).Article 

Google Scholar 
Beehner, J. C. et al. Corrigendum to “Testosterone related to age and life-history stages in male baboons and geladas” [Horm. Behav. 56/4 (2009) 472-480]. Horm Behav. 80, 149 (2016).Article 

Google Scholar 
Erb, R. E., Tillson, S. A., Hodgen, G. D. & Plotka, E. D. Urinary creatinine as an index compound for estimating rate of excretion of steroids in the domestic sow. J. Anim. Sci. 30, 79–85 (1970).Article 
CAS 

Google Scholar 
Tinsley Johnson, E., Snyder-Mackler, N., Lu, A., Bergman, T. J. & Beehner, J. C. Social and ecological drivers of reproductive seasonality in geladas. Behav. Ecol. 29, 574–588 (2018).Article 

Google Scholar 
Kaushik, S. & Kaur, J. Effect of chronic cold stress on intestinal epithelial cell proliferation and inflammation in rats. Stress 8, 191–197 (2005).Article 
CAS 

Google Scholar 
Jarvey, J. C., Low, B. S., Pappano, D. J. & Bergman, T. J. Graminivory and fallback foods: annual diet profile of geladas (Theropithecus gelada) living in the Simien Mountains National Park, Ethiopia. Int. J. Primatol. https://doi.org/10.1007/s10764-018-0018-x (2018).Article 

Google Scholar 
Gowda, C., Hadley, C. & Aiello, A. E. The association between food insecurity and inflammation in the US adult population. Am. J. Public Health. 102, 1579–1586 (2012).Article 

Google Scholar 
Becker, D. J. et al. Macroimmunology: The drivers and consequences of spatial patterns in wildlife immune defence. J. Anim. Ecol. 89, 972–995 (2020).Article 

Google Scholar 
Bates, D., Maechler, M., Bolker, B. & Walker, S. lme4: Linear mixed-effects models using Eigen and S4. R package version 1, 1–7 (2014).R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2021. https://www.R-project.org/.Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest Package: Tests in linear mixed effects models. J. Stat. Softw. 82, 1–26 (2017).Article 

Google Scholar 
Simon, A. K., Hollander, G. A. & McMichael, A. Evolution of the immune system in humans from infancy to old age. Proc. Biol. Sci. 282, 20143085 (2015).
Google Scholar 
Heinonen, S. et al. Infant immune response to respiratory viral infections. Immunol. Allergy Clin. North Am. 39, 361–376 (2019).Article 

Google Scholar 
Teran, R. et al. Immune system development during early childhood in tropical Latin America: Evidence for the age-dependent down regulation of the innate immune response. Clin. Immunol. 138, 299–310 (2011).Article 
CAS 

Google Scholar 
van de Pol, M. & Verhulst, S. Age-dependent traits: a new statistical model to separate within- and between-individual effects. Am. Nat. 167, 766–773 (2006).Article 

Google Scholar 
Furman, D. et al. Chronic inflammation in the etiology of disease across the life span. Nat. Med. 25, 1822–1832 (2019).Article 
CAS 

Google Scholar 
Petrovsky, N., McNair, P. & Harrison, L. C. Diurnal rhythms of pro-inflammatory cytokines: regulation by plasma cortisol and therapeutic implications. Cytokine 10, 307–312 (1998).Article 
CAS 

Google Scholar 
Lasselin, J., Rehman, J.-U., Åkerstedt, T., Lekander, M. & Axelsson, J. Effect of long-term sleep restriction and subsequent recovery sleep on the diurnal rhythms of white blood cell subpopulations. Brain Behav. Immun. 47, 93–99 (2015).Article 

Google Scholar 
Auzéby, A., Bogdan, A., Krosi, Z. & Touitou, Y. Time-dependence of urinary neopterin, a marker of cellular immune activity. Clin Chem. 34, 1866–1867 (1988).Article 

Google Scholar 
Ansari, A. & Williams, J. F. The eosinophilic response of the rat to infection with Taenia taeniaeformis. J. Parasitol. 62, 728–736 (1976).Article 
CAS 

Google Scholar 
Schneider-Crease, I. A., Snyder-Mackler, N., Jarvey, J. C. & Bergman, T. J. Molecular identification of Taenia serialis coenurosis in a wild Ethiopian gelada (Theropithecus gelada). Vet. Parasitol. 198, 240–243 (2013).Article 
CAS 

Google Scholar 
Terrazas, L. I., Bojalil, R., Govezensky, T. & Larralde, C. Shift from an early protective Th1-type immune response to a late permissive Th2-type response in murine cysticercosis (Taenia crassiceps). J. Parasitol. 84, 74–81 (1998).Article 
CAS 

Google Scholar 
Toenjes, S. A., Spolski, R. J., Mooney, K. A. & Kuhn, R. E. The systemic immune response of BALB/c mice infected with larval Taenia crassiceps is a mixed Th1/Th2-type response. Parasitology 118(Pt 6), 623–633 (1999).Article 
CAS 

Google Scholar 
Gaze, S. et al. Characterising the mucosal and systemic immune responses to experimental human hookworm infection. PLoS Pathog. 8, e1002520 (2012).Article 
CAS 

Google Scholar 
Johnston, M. J. G., MacDonald, J. A. & McKay, D. M. Parasitic helminths: a pharmacopeia of anti-inflammatory molecules. Parasitology 136, 125–147 (2009).Article 
CAS 

Google Scholar 
Cortés, A., Muñoz-Antoli, C., Esteban, J. G. & Toledo, R. Th2 and Th1 responses: Clear and hidden sides of immunity against intestinal helminths. Trends Parasitol. 33, 678–693 (2017).Article 

Google Scholar 
White, M. P. J., McManus, C. M. & Maizels, R. M. Regulatory T-cells in helminth infection: induction, function and therapeutic potential. Immunology 160, 248–260 (2020).Article 
CAS 

Google Scholar 
Maizels, R. M. & Holland, M. J. Parasite immunity: Pathways for expelling intestinal helminths. Curr Biol. 8, R711–R714 (1998).Article 
CAS 

Google Scholar 
Zhang, D. & Frenette, P. S. Cross talk between neutrophils and the microbiota. Blood 133, 2168–2177 (2019).Article 
CAS 

Google Scholar 
Wang, J., Chen, W.-D. & Wang, Y.-D. The relationship between gut microbiota and inflammatory diseases: The role of macrophages. Front. Microbiol. 11, 1065 (2020).Article 

Google Scholar 
Pallikkuth, S. et al. Age associated microbiome and microbial metabolites modulation and its association with systemic inflammation in a rhesus macaque model. Front. Immunol. 12, 748397 (2021).Article 
CAS 

Google Scholar 
Pierce, Z. et al. The infant gut microbiome is associated with stool markers of macrophage and neutrophil activity. FASEB J. 30, 668–9 (2016).
Google Scholar 
Levast, B. et al. Impact on the gut microbiota of intensive and prolonged antimicrobial therapy in patients with bone and joint infection. Front. Med. https://doi.org/10.3389/fmed.2021.586875 (2021).Article 

Google Scholar 
Hooper, L. V., Littman, D. R. & Macpherson, A. J. Interactions between the microbiota and the immune system. Science 336, 1268–1273 (2012).Article 
ADS 
CAS 

Google Scholar 
Round, J. L. & Mazmanian, S. K. The gut microbiota shapes intestinal immune responses during health and disease. Nat. Rev. Immunol. 9, 313–323 (2009).Article 
CAS 

Google Scholar 
Libertucci, J. & Young, V. B. The role of the microbiota in infectious diseases. Nat. Microbiol. 4, 35–45 (2019).Article 
CAS 

Google Scholar 
Huttenhower, C. et al. Structure, function and diversity of the healthy human microbiome. Nature 486, 207–214 (2012).Article 
ADS 
CAS 

Google Scholar 
Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K. & Knight, R. Diversity, stability and resilience of the human gut microbiota. Nature 489, 220–230 (2012).Article 
ADS 
CAS 

Google Scholar 
Blaser, M. J. & Falkow, S. What are the consequences of the disappearing human microbiota?. Nat. Rev. Microbiol. 7, 887–894 (2009).Article 
CAS 

Google Scholar 
Brown, E. M., Kenny, D. J. & Xavier, R. J. Gut microbiota regulation of T cells during inflammation and autoimmunity. Annu. Rev. Immunol. 37, 599–624 (2019).Article 
CAS 

Google Scholar 
Gollwitzer, E. S. & Marsland, B. J. Impact of early-life exposures on immune maturation and susceptibility to disease. Trends Immunol. 36, 684–696 (2015).Article 
CAS 

Google Scholar 
Ravi, A. et al. Loss of microbial diversity and pathogen domination of the gut microbiota in critically ill patients. Microbial. Genomics https://doi.org/10.1099/mgen.0.000293 (2019).McLaren, M. R. & Callahan, B. J. Pathogen resistance may be the principal evolutionary advantage provided by the microbiome. Philos. Trans. R. Soc. Lond. B Biol. Sci. 375, 20190592 (2020).Article 

Google Scholar 
Ragonnaud, E. & Biragyn, A. Gut microbiota as the key controllers of “healthy” aging of elderly people. Immun Ageing 18, 2 (2021).Article 

Google Scholar 
Wilmanski, T. et al. Gut microbiome pattern reflects healthy ageing and predicts survival in humans. Nat. Metab. 3, 274–286 (2021).Article 

Google Scholar 
Cattadori, I. M. et al. Impact of helminth infections and nutritional constraints on the small intestine microbiota. PLoS ONE 11, e0159770 (2016).Article 

Google Scholar 
Houlden, A. et al. Chronic Trichuris muris infection in C57BL/6 mice causes significant changes in host microbiota and metabolome: Effects reversed by pathogen clearance. PLoS ONE 10, e0125945 (2015).Article 

Google Scholar 
Holm, J. B. et al. Chronic Trichuris muris infection decreases diversity of the intestinal microbiota and concomitantly increases the abundance of Lactobacilli. PLoS ONE 10, e0125495 (2015).Article 

Google Scholar 
Peachey, L. E., Jenkins, T. P. & Cantacessi, C. This gut ain’t big enough for both of us. Or is it? Helminth–microbiota interactions in veterinary species. Trends Parasitol. 33, 619–632 (2017).Article 

Google Scholar  More

Franchini, M. & Mannucci, P. M. Air pollution and cardiovascular disease. Thromb. Res. 129, 230–234 (2012).Article 
CAS 

Google Scholar 
Langrish, J. P. et al. Reducing personal exposure to particulate air pollution improves cardiovascular health in patients with coronary heart disease. Environ. Health Perspect. 120, 367–372 (2012).Article 
CAS 

Google Scholar 
Tanwar, V., Katapadi, A., Adelstein, J. M., Grimmer, J. A. & Wold, L. E. Cardiac pathophysiology in response to environmental stress: A current review. Curr. Opin. Physiol. 1, 198–205 (2018).Article 

Google Scholar 
Franchini, M. & Mannucci, P. M. Particulate air pollution and cardiovascular risk: short-term and long-term effects. in Seminars in Thrombosis and Hemostasis. Vol. 35. 665–670 (© Thieme Medical Publishers, 2009).Shah, A. S. V. et al. Global association of air pollution and heart failure: A systematic review and meta-analysis. Lancet 382, 1039–1048 (2013).Article 
CAS 

Google Scholar 
Cesaroni, G. et al. Long term exposure to ambient air pollution and incidence of acute coronary events: Prospective cohort study and meta-analysis in 11 European cohorts from the ESCAPE Project. BMJ 348, f7412 (2014).Article 

Google Scholar 
Brook, R. D. et al. Particulate matter air pollution and cardiovascular disease: An update to the scientific statement from the American Heart Association. Circulation 121, 2331–2378 (2010).Article 
CAS 

Google Scholar 
Héroux, M.-E. et al. Quantifying the health impacts of ambient air pollutants: Recommendations of a WHO/Europe project. Int. J. Public Health 60, 619–627 (2015).Article 

Google Scholar 
An, Z., Jin, Y., Li, J., Li, W. & Wu, W. Impact of particulate air pollution on cardiovascular health. Curr. Allergy Asthma Rep. 18, 1–7 (2018).Article 
CAS 

Google Scholar 
Zanobetti, A., Baccarelli, A. & Schwartz, J. Gene-air pollution interaction and cardiovascular disease: A review. Prog. Cardiovasc. Dis. 53, 344–352 (2011).Article 
CAS 

Google Scholar 
Liang, R. et al. Effect of exposure to PM2.5 on blood pressure: A systematic review and meta-analysis. J. Hypertens. 32, 2130–2141 (2014).Article 
CAS 

Google Scholar 
Jerrett, M. et al. Traffic-related air pollution and obesity formation in children: A longitudinal, multilevel analysis. Environ. Heal. 13, 49 (2014).Article 

Google Scholar 
McConnell, R. et al. A longitudinal cohort study of body mass index and childhood exposure to secondhand tobacco smoke and air pollution: The Southern California Children’s Health Study. Environ. Health Perspect. 123, 360–366 (2015).Article 

Google Scholar 
Renzi, M. et al. Air pollution and occurrence of type 2 diabetes in a large cohort study. Environ. Int. 112, 68–76 (2018).Article 
CAS 

Google Scholar 
Jomova, K. et al. Arsenic: Toxicity, oxidative stress and human disease. J. Appl. Toxicol. 31, 95–107 (2011).CAS 

Google Scholar 
Al-Kindi, S. G., Brook, R. D., Biswal, S. & Rajagopalan, S. Environmental determinants of cardiovascular disease: Lessons learned from air pollution. Nat. Rev. Cardiol. 17, 656–672 (2020).Article 

Google Scholar 
Noroozian, M. The elderly population in iran: An ever growing concern in the health system. Iran. J. Psychiatry Behav. Sci. 6, 1 (2012).
Google Scholar 
Chokshi, D. A. & Farley, T. A. The cost-effectiveness of environmental approaches to disease prevention. N. Engl. J. Med. 367, 295–297 (2012).Article 
CAS 

Google Scholar 
Nieuwenhuijsen, M. J. Influence of urban and transport planning and the city environment on cardiovascular disease. Nat. Rev. Cardiol. 15, 432–438 (2018).Article 

Google Scholar 
Barnett, A. G. et al. The effects of air pollution on hospitalizations for cardiovascular disease in elderly people in Australian and New Zealand cities. Environ. Health Perspect. 114, 1018–1023 (2006).Article 
CAS 

Google Scholar 
Koken, P. J. M. et al. Temperature, air pollution, and hospitalization for cardiovascular diseases among elderly people in Denver. Environ. Health Perspect. 111, 1312–1317 (2003).Article 

Google Scholar 
Institute for Health Metrics and Evaluation. GBD 2019. (University of Washington, 2022).IHME. GBD 2019 Data and Tools Overview. (University of Washington, 2020).Dicker, D. et al. Global, regional, and national age-sex-specific mortality and life expectancy, 1950–2017: A systematic analysis for the Global Burden of Disease Study 2017. Lancet 392, 1684–1735 (2018).Article 

Google Scholar 
Miller, D. C. & Salkind, N. J. Handbook of Research Design and Social Measurement (Sage, 2002).Book 

Google Scholar 
Rosenthal, F. S., Carney, J. P. & Olinger, M. L. Out-of-hospital cardiac arrest and airborne fine particulate matter: A case–crossover analysis of emergency medical services data in Indianapolis, Indiana. Environ. Health Perspect. 116, 631–636 (2008).Article 

Google Scholar 
Ensor, K. B., Raun, L. H. & Persse, D. A case-crossover analysis of out-of-hospital cardiac arrest and air pollution. Circulation 127, 1192–1199 (2013).Article 
CAS 

Google Scholar 
Forastiere, F. et al. A case-crossover analysis of out-of-hospital coronary deaths and air pollution in Rome, Italy. Am. J. Respir. Crit. Care Med. 172, 1549–1555 (2005).Article 

Google Scholar 
Levy, D. et al. A case-crossover analysis of particulate matter air pollution and out-of-hospital primary cardiac arrest. Epidemiology 12, 193–199 (2001).Article 
CAS 

Google Scholar 
Silverman, R. A. et al. Association of ambient fine particles with out-of-hospital cardiac arrests in New York City. Am. J. Epidemiol. 172, 917–923 (2010).Article 

Google Scholar 
Naghavi, M. et al. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980–2016: A systematic analysis for the Global Burden of Disease Study 2016. Lancet 390, 1151–1210 (2017).Article 

Google Scholar 
Berend, N. Contribution of air pollution to COPD and small airway dysfunction. Respirology 21, 237–244 (2016).Article 

Google Scholar 
Vahedian, M., Khanjani, N., Mirzaee, M. & Koolivand, A. Associations of short-term exposure to air pollution with respiratory hospital admissions in Arak, Iran. J. Environ. Health Sci. Eng. 15, 17 (2017).Yaser, H. S., Narges, K., Yaser, S. & Rasoul, M. Air pollution and cardiovascular mortality in Kerman from 2006 to 2011. Am. J. Cardiovasc. Dis. Res. 2, 27–30 (2014).
Google Scholar 
Khaefi, M. et al. Association of particulate matter impact on prevalence of chronic obstructive pulmonary disease in Ahvaz, southwest Iran during 2009–2013. Aerosol Air Qual. Res. 17, 230–237 (2017).Article 
CAS 

Google Scholar 
Khaniabadi, Y. O. et al. Exposure to PM10, NO2, and O3 and impacts on human health. Environ. Sci. Pollut. Res. Int. 24, 2781–2789 (2017).Article 
CAS 

Google Scholar 
Momtazan, M. et al. An investigation of particulate matter and relevant cardiovascular risks in Abadan and Khorramshahr in 2014–2016. Toxin Rev. 38, 1–8 (2018).Martinelli, N., Olivieri, O. & Girelli, D. Air particulate matter and cardiovascular disease: A narrative review. Eur. J. Intern. Med. 24, 295–302 (2013).Article 
CAS 

Google Scholar 
Khaniabadi, Y. O. et al. Mortality and morbidity due to ambient air pollution in Iran. Clin. Epidemiol. Glob. Health 7, 222–227 (2019).Article 

Google Scholar 
Almeida-Silva, M. et al. Exposure and dose assessment to particle components among an elderly population. Atmos. Environ. 102, 156–166 (2015).Article 
ADS 
CAS 

Google Scholar 
Suh, H. H., Zanobetti, A., Schwartz, J. & Coull, B. A. Chemical properties of air pollutants and cause-specific hospital admissions among the elderly in Atlanta, Georgia. Environ. Health Perspect. 119, 1421–1428 (2011).Article 

Google Scholar 
Chien, L.-C., Yang, C.-H. & Yu, H.-L. Estimated effects of Asian dust storms on spatiotemporal distributions of clinic visits for respiratory diseases in Taipei children (Taiwan). Environ. Health Perspect. 120, 1215–1220 (2012).Article 

Google Scholar 
Khaniabadi, Y. O. et al. Chronic obstructive pulmonary diseases related to outdoor PM10, O3, SO2, and NO2 in a heavily polluted megacity of Iran. Environ. Sci. Pollut. Res. 25, 17726–17734 (2018).Article 
CAS 

Google Scholar 
Omidi Khaniabadi, Y. et al. Air quality modeling for health risk assessment of ambient PM10, PM2.5 and SO2 in Iran. Hum. Ecol. Risk Assess. Int. J. 25, 1298–1310 (2019).Newell, K., Kartsonaki, C., Lam, K. B. H. & Kurmi, O. P. Cardiorespiratory health effects of particulate ambient air pollution exposure in low-income and middle-income countries: A systematic review and meta-analysis. Lancet Planet. Health 1, e368–e380 (2017).Article 

Google Scholar 
Dominici, F. et al. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA 295, 1127–1134 (2006).Article 
CAS 

Google Scholar 
Qiu, H. et al. Inflammatory and oxidative stress responses of healthy elders to solar-assisted large-scale cleaning system (SALSCS) and changes in ambient air pollution: A quasi-interventional study in Xi’an, China. Sci. Total Environ. 806, 151217 (2022).Article 
ADS 
CAS 

Google Scholar 
Fiordelisi, A. et al. The mechanisms of air pollution and particulate matter in cardiovascular diseases. Heart Fail. Rev. 22, 337–347 (2017).Article 
CAS 

Google Scholar 
Yang, D., Yang, X., Deng, F. & Guo, X. Ambient air pollution and biomarkers of health effect. Ambient Air Pollut. Health Impact China 1017, 59–102 (2017).Lim, S. S. et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet 380, 2224–2260 (2012).Article 

Google Scholar 
Newby, D. E. et al. Expert position paper on air pollution and cardiovascular disease. Eur. Heart J. 36, 83–93 (2015).Article 
CAS 

Google Scholar 
Brook, R. D., Newby, D. E. & Rajagopalan, S. Air pollution and cardiometabolic disease: An update and call for clinical trials. Am. J. Hypertens. 31, 1–10 (2018).Article 
CAS 

Google Scholar 
Kang, S.-H. et al. Ambient air pollution and out-of-hospital cardiac arrest. Int. J. Cardiol. 203, 1086–1092 (2016).Article 

Google Scholar 
Thurston, G. D. et al. Ambient particulate matter air pollution exposure and mortality in the NIH-AARP diet and health cohort. Environ. Health Perspect. 124, 484–490 (2016).Article 

Google Scholar 
Gallagher, L. G. et al. Applying a moving total mortality count to the cities in the NMMAPS database to estimate the mortality effects of particulate matter air pollution. Circulation 172, 872–879 (2010).
Google Scholar 
Rodopoulou, S., Samoli, E., Chalbot, M.-C.G. & Kavouras, I. G. Air pollution and cardiovascular and respiratory emergency visits in Central Arkansas: A time-series analysis. Sci. Total Environ. 536, 872–879 (2015).Article 
ADS 
CAS 

Google Scholar 
Teng, T.-H.K. et al. A systematic review of air pollution and incidence of out-of-hospital cardiac arrest. J. Epidemiol. Commun. Health 68, 37–43 (2014).Article 

Google Scholar 
Brook, R. D. et al. Air pollution and cardiovascular disease: A statement for healthcare professionals from the expert panel on population and prevention science of the American Heart Association. Circulation 109, 2655–2671 (2004).Article 

Google Scholar 
Raza, A. et al. Short-term effects of air pollution on out-of-hospital cardiac arrest in Stockholm. Eur. Heart J. 35, 861–868 (2014).Article 
CAS 

Google Scholar 
Baccarelli, A. et al. Effects of exposure to air pollution on blood coagulation. J. Thromb. Haemost. 5, 252–260 (2007).Article 
CAS 

Google Scholar 
Franchini, M. & Mannucci, P. M. Thrombogenicity and cardiovascular effects of ambient air pollution. Blood 118, 2405–2412 (2011).Article 
CAS 

Google Scholar 
Yin, F. et al. Diesel exhaust induces systemic lipid peroxidation and development of dysfunctional pro-oxidant and pro-inflammatory high-density lipoprotein. Arterioscler. Thromb. Vasc. Biol. 33, 1153–1161 (2013).Article 
CAS 

Google Scholar 
Chirinos, J. A. et al. Elevation of endothelial microparticles, platelets, and leukocyte activation in patients with venous thromboembolism. J. Am. Coll. Cardiol. 45, 1467–1471 (2005).Article 
CAS 

Google Scholar 
Adar, S. D. et al. Fine particulate air pollution and the progression of carotid intima-medial thickness: A prospective cohort study from the multi-ethnic study of atherosclerosis and air pollution. PLoS Med. 10, e1001430 (2013).Article 

Google Scholar 
Kampfrath, T. et al. Chronic fine particulate matter exposure induces systemic vascular dysfunction via NADPH oxidase and TLR4 pathways. Circ. Res. 108, 716–726 (2011).Article 
CAS 

Google Scholar 
Sun, Q. et al. Long-term air pollution exposure and acceleration of atherosclerosis and vascular inflammation in an animal model. JAMA 294, 3003–3010 (2005).Article 
CAS 

Google Scholar 
Dennekamp, M. et al. Outdoor air pollution as a trigger for out-of-hospital cardiac arrests. Epidemiology 21, 494–500 (2010).Straney, L. et al. Evaluating the impact of air pollution on the incidence of out-of-hospital cardiac arrest in the Perth Metropolitan Region: 2000–2010. J. Epidemiol. Commun. Health 68, 6–12 (2014).Article 

Google Scholar 
Sullivan, J. et al. Exposure to ambient fine particulate matter and primary cardiac arrest among persons with and without clinically recognized heart disease. Am. J. Epidemiol. 157, 501–509 (2003).Article 
CAS 

Google Scholar 
Barton, T. J. et al. Traditional cardiovascular risk factors strongly underestimate the 5-year occurrence of cardiovascular morbidity and mortality in spinal cord injured individuals. Arch. Phys. Med. Rehabil. 102, 27–34 (2021).Article 

Google Scholar 
Burg, M. M. et al. Risk for incident hypertension associated with PTSD in military veterans, and the effect of PTSD treatment. Psychosom. Med. 79, 181 (2017).Article 

Google Scholar 
Hinojosa, R. Veterans’ likelihood of reporting cardiovascular disease. J. Am. Board Fam. Med. 32, 50–57 (2019).Article 

Google Scholar 
Rush, T., LeardMann, C. A. & Crum-Cianflone, N. F. Obesity and associated adverse health outcomes among US military members and veterans: Findings from the millennium cohort study. Obesity 24, 1582–1589 (2016).Article 

Google Scholar 
Stefanovics, E. A., Potenza, M. N. & Pietrzak, R. H. Smoking, obesity, and their co-occurrence in the US military veterans: Results from the national health and resilience in veterans study. J. Affect. Disord. 274, 354–362 (2020).Article 

Google Scholar 
Brook, R. D. et al. Insights into the mechanisms and mediators of the effects of air pollution exposure on blood pressure and vascular function in healthy humans. Hypertension 54, 659–667 (2009).Article 
CAS 

Google Scholar 
Rajagopalan, S. & Brook, R. D. Air pollution and type 2 diabetes: Mechanistic insights. Diabetes 61, 3037–3045 (2012).Article 
CAS 

Google Scholar 
Franklin, S. S. & Wong, N. D. Hypertension and cardiovascular disease: Contributions of the Framingham Heart Study. Glob. Heart 8, 49–57 (2013).Article 

Google Scholar 
Gu, D. et al. Blood pressure and risk of cardiovascular disease in Chinese men and women. Am. J. Hypertens. 21, 265–272 (2008).Article 

Google Scholar 
Wang, H. et al. Blood pressure, body mass index and risk of cardiovascular disease in Chinese men and women. BMC Public Health 10, 189 (2010).Article 

Google Scholar 
O’Brien, E. The Lancet Commission on hypertension: Addressing the global burden of raised blood pressure on current and future generations. J. Clin. Hypertens. 19, 564–568 (2017).Article 

Google Scholar 
Cai, Y. et al. Associations of short-term and long-term exposure to ambient air pollutants with hypertension: A systematic review and meta-analysis. Hypertension 68, 62–70 (2016).Article 
CAS 

Google Scholar 
Zhang, Z., Laden, F., Forman, J. P. & Hart, J. E. Long-term exposure to particulate matter and self-reported hypertension: A prospective analysis in the Nurses’ Health Study. Environ. Health Perspect. 124, 1414–1420 (2016).Article 

Google Scholar 
Cosselman, K. E., Navas-Acien, A. & Kaufman, J. D. Environmental factors in cardiovascular disease. Nat. Rev. Cardiol. 12, 627 (2015).Article 
CAS 

Google Scholar 
Baccarelli, A. et al. Effects of particulate air pollution on blood pressure in a highly exposed population in Beijing, China: A repeated-measure study. Environ. Heal. 10, 108 (2011).Article 

Google Scholar 
Mordukhovich, I. et al. Black carbon exposure, oxidative stress genes, and blood pressure in a repeated-measures study. Environ. Health Perspect. 117, 1767–1772 (2009).Article 
CAS 

Google Scholar 
Chen, H. et al. Spatial association between ambient fine particulate matter and incident hypertension. Circulation 129, 562–569 (2014).Article 
CAS 

Google Scholar 
Honjo, K. et al. The effects of smoking and smoking cessation on mortality from cardiovascular disease among Japanese: Pooled analysis of three large-scale cohort studies in Japan. Tob. Control 19, 50–57 (2010).Article 

Google Scholar 
Lawlor, D. A., Song, Y.-M., Sung, J., Ebrahim, S. & Smith, G. D. The association of smoking and cardiovascular disease in a population with low cholesterol levels: A study of 648 346 men from the Korean national health system prospective cohort study. Stroke 39, 760–767 (2008).Article 
CAS 

Google Scholar 
Wold, L. E. et al. Cardiovascular remodeling in response to long-term exposure to fine particulate matter air pollution. Circ. Hear. Fail. 5, 452–461 (2012).Article 
CAS 

Google Scholar 
Zoeller, R. T. et al. Endocrine-disrupting chemicals and public health protection: A statement of principles from The Endocrine Society. Endocrinology 153, 4097–4110 (2012).Article 
CAS 

Google Scholar 
Ruiz, D., Becerra, M., Jagai, J. S., Ard, K. & Sargis, R. M. Disparities in environmental exposures to endocrine-disrupting chemicals and diabetes risk in vulnerable populations. Diabetes Care 41, 193–205 (2018).Article 
CAS 

Google Scholar 
Taylor, D. Toxic Communities: Environmental Racism, Industrial Pollution, and Residential Mobility (NYU Press, 2014).
Google Scholar 
Newbold, R. R., Padilla-Banks, E. & Jefferson, W. N. Environmental estrogens and obesity. Mol. Cell. Endocrinol. 304, 84–89 (2009).Article 
CAS 

Google Scholar 
Szyszkowicz, M., Rowe, B. H. & Brook, R. D. Even low levels of ambient air pollutants are associated with increased emergency department visits for hypertension. Can. J. Cardiol. 28, 360–366 (2012).Article 
CAS 

Google Scholar 
van den Hooven, E. H. et al. Air pollution, blood pressure, and the risk of hypertensive complications during pregnancy: The generation R study. Hypertension 57, 406–412 (2011).Article 

Google Scholar 
Vali, M. et al. Effect of meteorological factors and Air Quality Index on the COVID-19 epidemiological characteristics: An ecological study among 210 countries. Environ. Sci. Pollut. Res. 38, 1–11 (2021).Kiani, B. et al. Association between heavy metals and colon cancer: An ecological study based on geographical information systems in North-Eastern Iran. BMC Cancer 21, 1–12 (2021).Article 

Google Scholar 
Cyranoski, D. China tests giant air cleaner to combat smog. Nature 555, 152–154 (2018).Article 
ADS 
CAS 

Google Scholar  More

Study areaPenang is situated in the northern part of Peninsular Malaysia and lies within the latitudes 5°12’N to 5°30′ N and longitudes 100°09’E to 100°26’E (Fig. 7). Penang with a land area of 295 Km2, has an estimated population of 720,000 and is regarded as the most populated island in Malaysia. Penang shares the same border on the north and east with Kedah State and the south with Perak State. There are two main parts of Penang State: Penang Island and the mainland which is also regarded as Seberang Perai. These two parts of the State are connected by the two bridges. The eastern part of Penang Island is the most urbanized area comprising industries, commercial centres and residential buildings. However, the western part is less developed comprising mainly hilly terrain and forests22. This study is focused on the Island part of Penang. This island is endowed with a yearly equatorial climate (hot and humid). It has a mean annual temperature ranging between 27 and 30 °C while the mean annual relative humidity ranges between 70 and 90%. Also, the mean annual rainfall is about 267–624 cm.Figure 7The map of Penang State showing the Penang Island (created by the authors using ArcMap 10.8 software).Full size imageData acquisitionThe flow chart of the methodology is presented in Fig. 8. Landsat satellite images were used for the assessment of changes in land use covering a period of 2010–2021 (11 years).Figure 8The flow chart of the methodology.Full size imageThese images were gotten from the website of the United State Geological Survey (https://earthexplorer.usgs.gov). The Landsat images include the Landsat 5 TM (thematic mapper) and Landsat 8 OLI / TIRS (operational land imager / thermal infrared sensor). These were downloaded from the Landsat level 1 dataset (Table 6) with additional criteria which reduced the.Table 6 The characteristics of the satellite data used.Full size tableDetermination of LST and NDVI for Landsat 5 and 8Band 6 of Landsat 5 and band 10 of Landsat 8 were used for the determination of the land surface temperature (LST). The LST and normalized difference vegetation index were determined using the following steps:Conversion of top of atmosphere (TOA) radianceUsing the radiance rescaling factor, thermal infra-red digital numbers were converted to TOA spectral radiance using the equation below29: (frac{Red – NIR}{{Red + NIR}}) (frac{Red – NIR}{{Red + NIR}}) For Landsat 8,$$ {text{L}}lambda = left( {{text{ML}} times {text{ Qcal}}} right) + left( {{text{AL}} – {text{Oi}}} right) $$
(1)
For Landsat 5,$$ {text{L}}lambda = left( {{ }frac{{{text{LMax}}lambda – {text{LMin}}lambda }}{{{text{QcalMax}} – {text{QcalMin}}}}} right) times left( {left( {{text{Qcal }} – {text{QcalMin}}} right) + {text{LMin}}lambda } right) $$
(2)
where Lλ is TOA spectral radiance, ML is radiance multiplicative band Number, AL is radiance add band number, Qcal is quantized and calibrated standard product pixel values (DN for band 6 or band 10), Oi is the correction value for the respective bands, LMaxλ is spectral radiance scaled to QcalMax, LMinλ is spectral radiance scaled to QcalMin, QcalMax is maximum quantized calibrated pixel value, and QcalMin is minimum quantized calibrated pixel value.Conversion to TOA brightness temperature (BT)Spectral radiance data were converted to TOA brightness temperature using the thermal constant values in the Metadata file29.Kelvin (K) to Celcius (°C) degrees$$ BT = {raise0.5exhbox{$scriptstyle {K2}$} kern-0.1em/kern-0.15em lower0.25exhbox{$scriptstyle {ln left( {frac{K1}{{{text{L}}lambda { } + { }1}}} right)}$}} – 273.15 $$
(3)
where BT is the Top of atmosphere brightness temperature (°C), Lλ is TOA spectral radiance (W.m−2 .sr−1 .µm−1)), K1 is the K1 constant band number, and K2 is the K2 constant band number. For Landsat 5, K1 is 607.76, and K2 is 1260.56.Normalized difference vegetation index (NDVI)The Normalized Difference Vegetation Index (NDVI) is a standardized vegetation index which reveals the intensity of greenness and surface radiant temperature of the area30,31. The index value of NDVI usually ranges from − 1 to 1. The higher NDVI value indicates that the vegetation of the area is denser and healthier. This shows that the NDVI values of normal healthy vegetation range from 0.1– 0.75, while it is almost zero for rock and soil, and negative value for water bodies24. The NDVI is calculated using the followings:$$ {text{NDVI }} = frac{{left( {{text{NIR }}{-}{text{ RED}}} right){ }}}{{left( {{text{NIR }} + {text{ RED}}} right)}} $$
(4)
In Landsat 4–7$$ {text{NDVI }} = , left( {{text{Band 4 }}{-}{text{ Band 3}}} right) , / , left( {{text{Band 4 }} + {text{ Band 3}}} right) $$In Landsat 8$$ {text{NDVI }} = , left( {{text{Band 5 }}{-}{text{ Band 4}}} right) , / , left( {{text{Band 5 }} + {text{ Band 4}}} right) $$where: RED = DN values from the RED band, and NIR = DN values from the Near Infra-red band.Land Surface Emissivity (LSE)Land Surface Emissivity is the average emissivity of an element on the surface of the earth calculated from NDVI values.$$ {text{PV }} = left{ {frac{{left( {{text{NDVI }} – {text{ NDVImin}}} right)}}{{left( {{text{NDVImax }} – {text{ NDVImin}}} right)}}} right}^{2} $$
(5)
where PV is the Proportion of vegetation, NDVI is the DN value from the NDVI image, NDVImin is the minimum DN value from the NDVI image, and NDVImax is the maximum DN value from the NDVI image.$$ {text{E }} = left( {0.004{ } times {text{PV}}} right) + 0.986 $$
(6)
where E is land surface emissivity, PV is the Proportion of vegetation, 0.986 corresponds to a correction value of the equation.Land Surface Temperature (LST)Land Surface Temperature (LST) is the radiative temperature which is calculated using top of atmosphere brightness temperature, the wavelength of emitted radiance and land surface emissivity.$$ {text{LST}} = {raise0.5exhbox{$scriptstyle {BT}$} kern-0.1em/kern-0.15em lower0.25exhbox{$scriptstyle {left( {1 + left( {{lambda } times { }frac{{{text{BT}}}}{{{text{c}}2}}} right) times {text{ln}}left( {text{E}} right)} right)}$}} $$
(7)
Here c2 is 14388. The value of λ for Landsat 5 (Band 6) is 11.5 µm and Landsat 8 (Band 10) is 10.8 µm.Where BT is the top of atmosphere brightness temperature, λ is the wavelength of emitted radiance, and E is land surface emissivity.c2 = h*c/s (1.4388*10–2 mK = 14388 mK), h is Planck’s constant (6.626*1034 Js), s is Boltzmann constant (1.38*1023 JK), c is velocity of light (2.998*108 m/s).Determination of land use and land cover (LULC) of the study areaThe Landsat images were pre-screened and subjected to clipping and classification32. The boundary shape file of Penang was used to clip out the area of study.Image classificationThe unsupervised method involving a random assignment of sample training points and supervised methods of satellite image classification was employed in this study for determining the LULC types. This mixture of image classification methods has been reported as vital in achieving a high accuracy level33. Bands 5, 4 and 3 were used to classify Landsat 8 while bands 4, 3 and 2 were used for classifying Landsat 5. We used the extraction by mask in the spatial analyst tool of ArcMap 10.2.1 software to extract the study area from the selected bands of the Landsat satellite images. A widely used supervised image classification method was adopted for classifying the Landsat bands in this study32,34. The principle of operation of this method involves the identification of known sample training points which are then used to classify other unknown points with related spectral signatures35. The three monochromatic satellite bands were combined to produce the false colour composite (FCC) using the data management tool36. This involves drawing polygons on the LULC type to select the training points. The LULC types adopted for this study include urbanized areas, agricultural land, rocks, forests, bare surfaces, and water bodies. These were modified LULC types from IPOC Good Practice Guidance37. To achieve this, a minimum of 40 sample points were selected randomly for each category of LULC type36. Having prior knowledge of the study area assisted in the selection of the training points38.The multivariate maximum likelihood classification (MLC) technique was used for transforming the images. Other image transformation techniques have been used by researchers. These include the fuzzy set classifier, neural networks (NN) classifier, extraction and classification of homogenous objects (ECHO) classifier, per-field classifier, sub-pixel classifier, decision trees (DTs), support vector machines (SVMs), minimum distance classifier (MDC) and so-on39. The adoption of any of these techniques is dependent on the knowledge of the area of study, band selection, accessibility of data, the complexity of the landscape, the classification algorithm, and the proficiency of the analyst39. We preferred MLC to other techniques in this study due to its reported high level of accuracy in tropical regions32,34. Another reason for choosing MLC is that it is readily incorporated in many widely used GIS software packages. This MLC algorithm operates based on assigning pixels to the highest probability class and establishing the class ownership of such pixels. It is also regarded as a parametric classifier whose data follows almost a normal distribution39. We ensured the accuracy of this classifier by assigning a large number of training sample points using our prior knowledge of the study area.Description of the LULC categoriesThe urbanized area is the developed part of the study area. This includes houses, roads, railways, and industries. This is also known to be a settlement in other literature40. Agricultural land is the part of the study area dominated by agricultural activities and herbaceous plants and grasses. Agricultural land is generally a product of deforestation36. Rocks are part of the study area comprising solid mineral materials (rocks). Bare land is the bare soil which is either made open by natural or human activities.Forests are parts of the study area dominated by trees. They can be primary or secondary forests depending on the rate of disturbances. According to41, forest land is an area having more than 0.5 ha of flora comprising trees (height is above 5 m) with a canopy greater than 10%. The forests in Penang are generally both primary and secondary42. Water bodies are parts of the study area covered by water seasonally or permanently. These include seas, rivers, lakes, ponds, streams, or reservoirs40.Determination of change in the LULCThe rate and extent of change in the LULC of Penang within the periods under consideration were determined following the formula below43:$$ {text{Changed area }}left( {{text{C}}_{{text{a}}} } right) , = {text{ T}}_{{text{a}}} left( {text{year 2}} right) , {-}{text{ T}}_{{text{a}}} left( {text{year 1}} right) $$
(8)
$$ {text{Changed extent }}left( {{text{C}}_{{text{e}}} } right) , = {text{ C}}_{{text{a}}} /{text{ T}}_{{text{a}}} left( {text{year 1}} right) $$
(9)
$$ {text{Percentage of change }} = {text{ C}}_{{text{e}}} {text{x 1}}00 $$
(10)
where Ta means the total area.Determination of relationship between LST and NDVIThe values of LST and NDVI at 20 random points of each LULC class were used. The relationship between the LST and NDVI across all the LULC classes in each year was determined using the bivariate linear regression analysis. This was done in Paleontological Statistical (PAST) package 3.0.Classification accuracy assessmentThe classification accuracy was assessed by taking ground truth coordinate data of the LULC of the study area using a geographical positioning system (GPS) device (Garmin Etrex 10). These data were compared with the LULC classified in this study32. Consequently, an error matrix was generated. This normally uses ground truth data to explain the accuracy of the classified LULC. The error matrix comprises the user’s accuracy, the producer’s accuracy, overall accuracy and the Kappa index32.The producer’s accuracy (omission error) represents the probability of the correctly classified reference pixel and it is determined using this formula below:$${text{Producer’s accuracy }}left( % right) , = { 1}00% , – {text{ error of omission}} $$
(11)
Also, the user’s accuracy (commission error) represents the probability that the classified pixel matches the one on the ground36 and it is determined using the formula below:$$ {text{User’s accuracy }}left( % right) , = { 1}00% , – {text{ error of commission}} $$
(12)
The statistical accuracy of the matrix was determined using the Kappa coefficient44. This Kappa coefficient ranges from − 1 to + 145. Therefore, the overall accuracy of the classification was determined by dividing the total number of correctly classified pixels by the total number of sampled ground truth data40. More

Field parameters and major ionsThe results of the hydrochemistry and environmental isotopes for the 40 lakes are presented spatially in Figs. S1–S11 and are located in Tables S1 and S2.The lake waters are oxic (8.6–12.6 mg l−1) and range from slightly acidic (pH 6.0) to slightly alkaline (pH 9.2). Lake water temperatures are generally highest for lakes along the west coast (greater than 10 °C, Table S2). Phosphate concentrations are below detection level (0.1 mg l−1) for all lakes and nitrate was low ranging from below detection limit ( More

Chemotaxonomic characteristicsThe predominant respiratory quinone for all novel strains was ubiquinone 8 (Q-8), consistent with other Ramlibacter species. C16:0 and summed feature 3 (consisting of C16:1 ω7c and/or C16:1 ω6c) were identified as the common major fatty acids ( > 10%) of the novel strains USB13T, AW1T, GTP1T, and HM2T. Other than the aforementioned fatty acids, strain USB13T had C10:0 3-OH additionally as its major fatty acid, whereas strains AW1T and HM2T shared C17:0 cyclo and summed feature 8 (consisting of C18:1 ω7c and/or C18: 1 ω6c) as its additional fatty acids. Detailed comparisons of the fatty acid profiles of the novel strains and their reference strains are summarized in Table S1.Strains USB13T, AW1T, GTP1T, and HM2T shared major polar lipids diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE), which was consistent with the major polar lipids of the reference strains. Additionally, the polar lipid profile of USB13T consisted of one unidentified phosphoaminolipid, two unidentified phosphoglycoaminolipids, and six unidentified polar lipids while the polar lipid profile of AW1T had one unidentified lipid, one unidentified phosphoglycolipid, and three unidentified glycolipids in addition. The polar lipid profile of strain GTP1T additionally consisted of two unidentified phosphoaminolipids, and that of strain HM2T additionally had one unidentified phosphoaminolipid, one unidentified phosphoglycolipid, one unidentified phosphoglycoaminolipid, and two unidentified phospholipids. Polar lipid profiles of the novel strains USB13T, AW1T, GTP1T, and HM2T are shown in Figure S1.Physiological, morphological characteristics, and screening of cellulose-degrading strainsWhen grown on R2A agar, strain USB13T produced reddish white and flat colonies while strain AW1T produced orange, convex colonies, strain GTP1T produced white, convex colonies, and strain HM2T produced cream-colored, flat, transparent colonies. Under TEM, monotrichous flagella were observed only in strain HM2T, and when tested for motility, strain USB13T and AW1T showed gliding motility, whereas strain GTP1T was non-motile. Strains USB13T and HM2T showed positive results for both catalase and oxidase activities; strain AW1T showed positive results for catalase and negative results for oxidase activity, and strain GTP1T showed negative results for catalase and positive results for oxidase activity. All strains were identified to be strictly aerobic, while showing negative results for urea, gelatin, starch, chitin, and DNA hydrolysis and positive results for hydrolysis of Tween 80. In addition, strain USB13T was the only strain to produce iron-chelating siderophores. When tested for NaCl tolerance, growth of strain USB13T was observed in NaCl concentrations of 0–7% (w/v), possibly due to the fact the strain was isolated from a marine environment. A detailed comparison of physiological and morphological characteristics between the novel species and its closely related Ramlibacter strains is presented in Table 1, while TEM images of the novel strains are shown in Figure S2. Results of the reference strains in Table 1 coincided with the data from the original literature1,3,4,5,7,8.Table 1 Characteristics differentiating strains USB13T, AW1T, GTP1T, and HM2T from closely related strains of the genus Ramlibacter.Full size tableStrains: 1, USB13T; 2, AW1T; 3, GTP1T; 4, HM2T; R. monticola KACC 19175T; 6, R. alkalitolerans KACC 19305T; 7, R. ginsenosidimutans KACC 17527T; 8, R. humi KCTC 52922T; 9, R. henchirensis KACC 11925T; 10, R. tataouinensis KACC 11924T; 11, R. rhizophilus KCTC 52083T. All strains are positive for esterase lipase (C8), while all strains are negative for chitin hydrolysis. All data were obtained from this study unless indicated otherwise. + , Positive; w + , weakly positive; -, negative.R2A agar plates supplemented with 1% (w/v) CMC were stained with Congo red dye after 7 days of incubation. Clear zones only formed around colonies of strain USB13T, indicating that strain USB13T solely possessed CMC-hydrolyzing activity among the four novel strains. When inoculated in basal salt medium, filter paper from the USB13T sample underwent degradation, whereas samples containing strains AW1T, GTP1T, and HM2T did not show any signs of degradation.Phylogenetic and genomic analysesEzBioCloud search results and BLASTn searches revealed that the novel strains belonged to the family Comamonadaceae and genus Ramlibacter. Using BLASTn, 16S rRNA gene sequence similarities were determined where strain USB13T was closest to strain GTP1T (98.5%), followed by strain HM2T (98.1%) and strain AW1T (97.1%). Strain AW1T shared the highest similarity with strain GTP1T (97.3%), followed by strain HM2T (97.1%), while strain GTP1T shared a similarity of 98.2% with strain HM2T. Phylogenetic analysis based on the MP method (Fig. 1) showed the clustering of the novel strains USB13T, AW1T, GTP1T, and HM2T with strains such as R. monticola G-3-2T, R. ginsenosidimutans BXN5-27T, R. alkalitolerans CJ661T, and R. rhizophilus YS3.2.7T. Similar topologies were observed in trees reconstructed by ML (Figure S3) and MP methods. The UBCG phylogenomic tree (Fig. 2), which was reconstructed using whole genome sequences, also showed close clustering of the selected reference strains and novel strains.Figure 1Maximum-parsimony (MP) tree reconstructed based on 16S rRNA gene sequences, showing the relationship between strains USB13T, AW1T, GTP1T, and HM2T and other closely related type strains. Bootstrap values based on 1000 replications are listed as percentages at branching points. Only bootstrap values exceeding 50% are shown. Bar, 50 substitutions per nucleotide position.Full size imageFigure 2Phylogenomic tree of strains USB13T, AW1T, GTP1T, and HM2T and their closely related taxa was reconstructed based on core genomes using UBCG version 3.0 pipeline42. NCBI GenBank accession numbers are shown in parentheses. Bootstrap analysis was carried out using 1000 replications. Percentage bootstrap values ( > 50%) are given at branching points. Bar, 0.050 substitution per position.Full size imageDraft genome sequences of the novel strains USB13T, AW1T, GTP1T, and HM2T were deposited in the GenBank database under the accession numbers JACORT000000000, JAEQNA000000000, JACORU000000000, and JADDIV000000000, respectively. In addition, the draft genome sequences of R. monticola KACC 19175T, R. alkalitolerans KACC 19305T, and R. ginsenosidimutans KACC 17527T were also deposited in GenBenk under the accession numbers JAEQNE000000000, JAEQND000000000, and JAEPWM000000000, respectively. The assembled genome size of the novel strains USB13T, AW1T, GTP1T, and HM2T was 5.53 Mbp, 5.11 Mbp, 6.15 Mbp, 4.31 Mbp, respectively. G + C content ranged from 67.9% to 69.9%, which was similar to those of the reference strains. The genomic features of the novel strains and their closely related Ramlibacter strains are presented in Table S2. CheckM analysis showed the following estimations for each strain: USB13T, had a 99.84% completeness and 0.68% contamination; AWIT, had a 99.84% completeness and 0.86% contamination; GTP1T, had a 99.38% completeness and 1.32% contamination; HM2T, had a 97.51% completeness and 0.16% contamination. These results indicated that the draft genome results for all strains were reliable. ANI values between the novel strains and reference strains ranged from 76.5–83.4% while dDDH values ranged from 20.7–26.7%, and AAI values ranged from 65.7–80.4%. All values were below the threshold for delineation of a new species54. ANI values between the novel strains and their reference strains are presented in Fig. 3, while a detailed comparison of GGDC and AAI values are shown in Table 2.Figure 3Heatmap of strains USB13T, AW1T, GTP1T, and HM2T and other closely related strains within the genus Ramlibacter, generated with OrthoANI values calculated using OAT software45. Bacterial strains and accession numbers are indentical to those of Fig. 2.Full size imageTable 2 Average amino acid identity (AAI) and digital DNA-DNA hybridization (dDDH) value comparisons between the closely related Ramlibacter type species and the novel strains, USB13T, AW1T, GTP1T, and HM2T. AAI values were calculated by two-way AAI, while dDDH values were calculated based on formula 246.Full size tableBased on NCBI PGAP annotation and CAZyme prediction results, strain USB13T, which was the only strain to show cellulolytic activity, possessed a total of four protein CDs encoding CAZymes, namely, two GH15 proteins, one glycosyl hydrolase protein, and one GH99-like domain-containing protein. Despite not showing any cellulolytic activity, strain AW1T possessed eight CAZyme CDs; the most amount among the novel strains. The enzymes include, two GH2 proteins, one GH5 protein, three GH15 proteins, one glycoside hydrolase protein, and one cellulase family glycosyl hydrolase. Strain GTP1T possessed two CDs encoding one GH15 protein and one GH16 protein; strain HM2T possessed three CDs encoding one GH2, one GH15, and one GH18 protein. All strains possessed GH15, which is known for its glucoamylase activity in fungi55. A detailed summary of the novel strains CAZymes are presented in Table S3 and a comparison of CAZyme numbers between strains USB13T, AW1T, GTP1T, and HM2T is summarized in Table S4. The presence of these genes may suggest the cellulolytic activity of strain USB13T, while it is uncertain why GH families responsible for endoglucanase (GH 5–8, 12, 16, 44, 45, 48, 51, 64, 71, 74, 81, 87, 124, and 128), exoglucanase (GH 5–7, and 48), and β-glucosidase (GH 1, 3, 4, 17, 30, and 116) were not present in the genome11.COG predictions (Fig. 4) revealed that the majority of the core genes of the four novel strains accounted for genes belonging to the functional categories C (energy production and conversion), E (amino acid transport and metabolism), I (lipid transport and metabolism), T (signal transduction mechanisms), and K (transcription). Meanwhile, the number of core genes belonging in category G, carbohydrate transport and metabolism, was the highest for strain USB13T (258), followed by GTP1T (230), HM2T (212), and AW1T (181). The high number of genes in strain USB13T may be a contributing factor in the strain’s cellulolytic activity. A comparison of COG gene count distribution of the novel strains is presented in Table S5.Figure 4Comparison of total number of matched genes of strains USB13T, AW1T, GTP1T, and HM2T according to functional classes based on Cluster of Orthologous Groups of proteins (COG) predictions48.Full size imageAntiSMASH analysis results showed four gene clusters within the genome of strain USB13T: ribosomally synthesized and post-translationally modified peptides (RIPP)-like cluster (989,516–1,000,916 nt; JACORT010000001), terpene synthesis (8,622–30,347 nt; JACORT010000003), RIPP precursor peptide recognition element (RRE)-containing cluster (311,469–333,619 nt; JACORT010000004), and redox-cofactor (281,860–303,948 nt; JACORT010000007). Among the clusters, the RRE-containing cluster showed 11% similarity to streptobactin, a tricatechol-type siderophore isolated from Streptomyces sp. YM5-79956. Strain AW1T had a total of eight gene clusters which encoded for: arylpolyene (165,946–207,130 nt), terpene (618,322–640,854 nt), RIPP-like proteins (804,411–819,137 nt), non-ribosomal peptide synthetase cluster (NRPS)-like (61,798–104,764 nt), betalactone (323,399–348,739 nt), N-acetylglutaminylglutamine amide (NAGGN; 106,834–121,648 nt), type I polyketide synthase (T1PKS; 56,584–107,578 nt), and heterocyst glycolipid synthase-like polyketide synthase (hglE-KS; 75,419–113,566 nt). Strain GTP1T possessed four gene clusters that encoded for RRE-containing cluster (175,155–199,102 nt), homoserine lactone (110,293–130,892 nt), a signaling molecule known for its involvement in bacterial quorum sensing, the RIPP-like cluster (38,002–48,856 nt), and terpene synthesis (47,942–69,701 nt). Strain HM2T had two gene clusters that encoded for resorcinol (403,967–445,901 nt), an organic compound known for its antiseptic properties, and terpene (697,660–721,242 nt), which showed 100% similarity for carotenoid synthesis. BRIG analysis results showed that a majority of the regions within the four analyzed genomes were conserved with at least 70% similarity (Figure S4).Cellulolytic potential and FE-SEM analysis of strain USB13T
A USB13T-inoculated basal salt medium sample containing degraded filter paper was examined under FE-SEM to observe the morphological interactions between cellulose fibers and USB13T cells. Images in Fig. 5 show individual rod cells of strain USB13T surrounding filter paper fibers, indicating bacterial adherence.Figure 5Field emission-scanning electron microscopy (FE-SEM) images of adhesion of strain USB13T to degraded filter paper fibers. Arrows indicate filter paper fibers. (A) low magnification (5000(times)) and (B), high magnification (20,000(times)) images of strain USB13T surrounding filter paper fibers.Full size imageThe enzymatic assay results showed endoglucanase, exoglucanase, β-glucosidase, and filter paper cellulase (FPCase) activities of strain USB13T, wherein activities for endoglucanase was the highest and β-glucosidase was the lowest in all experiments. As seen in Fig. 6A, enzyme activity for all cellulolytic enzymes increased along with its cultivation time. In addition, enzyme activities showed the highest results when tested on buffer solutions of pH 6.0 (Fig. 6B), indicating the enzymes’ resistance to moderately acidic conditions. The pH of the buffer solution seemed to be an important factor in enzyme activity, as activity of endoglucanase, exoglucanase, and FPCase drastically decreased when the pH was altered from pH 6.0 to pH 7.0. Meanwhile, β-glucosidase activity was relatively resistant to pH change as its activity decreased less than 50%. On day 7, enzyme activities were measured as 1.91 IU/mL for endoglucanase, 1.77 IU/mL for exoglucanase, 0.76 IU/mL for β-glucosidase, and 1.12 IU/mL for FPCase at pH 6.0. When measured at pH 8.0, where enzyme activity was the lowest, enzyme activities were measured as 0.51 IU/mL for endoglucanase, 0.25 IU/mL for exoglucanase, 0.45 IU/mL for β-glucosidase, and 0.23 IU/mL for FPCase; all values were less than half of the measured activity at pH 6.0. The results of strain USB13T are comparable to FPCase results of other species such as Mucilaginibacter polytrichastri RG4-7T (0.98 U/mL) isolated from the moss Polytrichastrum formosum14, Paenibacillus lautus BHU3 (2.9 U/mL) isolated from a landfill site57, and Serratia rubidaea DBT4 (0.5 U/mL) isolated from the gastrointestinal tract of a black Bengal goat58.Figure 6Cellulolytic enzyme activity of strain USB13T. Enzyme activity was defined in international units (IU); one unit of enzymatic activity was defined as the amount of enzyme that releases 1 μmol of glucose per mL per 1 min of reaction. (A) cellulase activity results under different cultivation time; (B) cellulase activity under different buffer solution pH. Values in the figure are mean values of triplicate data with standard deviation.Full size imageDespite the absence of the main three cellulolytic enzymes, endoglucanase, exoglucanase, and β-glucosidase, the cellulolytic activity of strain USB13T was confirmed through SEM images, CMC agar screening, and enzymatic assay results. However, because PGAP annotation results showed that other non-cellulolytic strains also possessed CAZymes, in some cases more than strain USB13T, further research is necessary to understand the mechanics of how CAZymes and other cellulases interact to degrade cellulose, and how these genes are expressed under certain conditions. Furthermore, the cellulolytic activity of strain USB13T can be further optimized for commercial use by adjusting growth conditions such as pH, temperature, and growth media.While cellulolytic bacteria are known to inhabit animal intestinal tracts, the rumen, and soil, they can be found almost everywhere, such as ocean floors, municipal landfills, and even extreme environments such as hot springs59. In these habitats, cellulolytic bacteria utilize cellulose while cohabiting with non-cellulolytic bacteria. There have been many studies suggesting the synergistic role non-cellulolytic bacteria play in cellulose degradation, where non-cellulolytic bacteria aid cellulose degradation by neutralizing pH or removing harmful metabolites60,61,62.Bacterial cellulases have shown immense value in various industries such as animal feed processing, food and brewery production, and agriculture, not to mention biofuel synthesis through biomass utilization11. Due to the versatile uses of bacterial cellulases, the cellulolytic strain USB13T has the potential to become an invaluable resource. However, further research of the novel strain’s cellulose-degradation mechanisms is necessary to develop and commercially make use of its bacterial cellulases in the future. In addition, research regarding co-culturing non-cellulolytic bacteria and strain USB13T may also help in developing effective methods to use an otherwise underutilized bioresource.Taxonomy of novel Ramlibacter speciesWhile phylogenetic analyses indicated that the novel strains USB13T, AW1T, GTP1T, and HM2T should be assigned to the genus Ramlibacter, differences in fatty acid compositions, polar lipid profiles, and physiological characteristics suggested that the four novel strains are noticeably distinct from other validly published species of the genus. Additionally, genomic characteristics such as ANI, dDDH, and AAI values further supported the novel strains’ position as a distinct species within the genus Ramlibacter. Therefore, we propose that the strains USB13T, AW1T, GTP1T, and HM2T represent novel species within the genus Ramlibacter.Description of the novel Ramlibacter speciesThe descriptions of the novel species are given according to the standards of the Judicial Commission of the International Committee on Systematic Bacteriology63.Description of Ramlibacter cellulosilyticus sp. nov
Ramlibacter cellulosilyticus (cel.lu.lo.si.ly’ti.cus. N.L. n. cellulosum, cellulose; N.L. adj. lyticus from Gr. lytikos, dissolving; N.L. masc. adj. cellulosilyticus, cellulose-dissolving).Cells of strain USB13T are Gram-negative, rod-shaped, non-flagellated and motile by gliding. The strain is positive for both oxidase and catalase activity, while cells have a width of 0.3–0.5 μm and length of 2.0–2.4 μm. When observed on R2A agar, colonies are reddish white, flat with entire margins, and have a diameter of 1–2 mm. Growth of strain USB13T is observed at 7–50 °C (optimum, 28–30 °C), at pH 5.0–10.0 (optimum, pH 6.0), and at NaCl concentrations of 0–7% (optimum, 0–3%). The strain is unable to grow in anaerobic conditions. Produces siderophores and hydrolyzes Tween 20, Tween 80, CMC, and esculin. According to the API ZYM results, the strain showed positive results for alkaline phosphatase, esterase lipase (C8), leucine arylamidase, acid phosphatase, β-galactosidase, α-glucosidase, and β-glucosidase. In the API 20NE assay, strain USB13T showed positive results only for β-galactosidase. The predominant respiratory quinone is ubiquinone 8 (Q-8). The major fatty acids are C16:0, C10:0 3-OH, and summed feature 3 (consisting of C16:1 ω7c and/or C16:1 ω6c). The polar lipid profile consists of diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), phosphatidylethanolamine (PE), one unidentified phosphoaminolipid, two unidentified phosphoglycoaminolipids, and six unidentified polar lipids. The G + C content is 69.7%. The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA gene sequence and the assembled genome sequence of strain USB13T are MN603953 and JACORT000000000, respectively.The type strain USB13T (= KACC 21656T = NBRC 114839T) was isolated from shallow coastal water at Haeundae Beach, Busan, Republic of Korea.Description of Ramlibacter aurantiacus sp. nov
Ramlibacter aurantiacus (au.ran.ti’a.cus. L. masc. adj. aurantiacus, orange-colored, referring to the orange colonies of the strain).Cells of strain AW1T are Gram-negative, coccoid to short rod-shaped, non-flagellated, and motile by gliding. The strain is negative for oxidase activity, and positive for catalase activity. When observed on R2A agar, colonies are orange, convex, with entire margins, and 0.5–1.0 mm in diameter. Under TEM cells have and approximate width of 0.3–0.5 μm and length of 0.6–0.8 μm. Growth of strain AW1T can be observed at 7–45 °C (optimum, 30 °C), at pH 7.0–10.0 (optimum, 7.0–8.0), and at NaCl concentrations of 0–3% (optimum, 0–1%). The strain does not grow under anaerobic conditions but is able to hydrolyze Tween 80. In addition, AW1T is not able to produce siderophores. In the API ZYM assay, positive for alkaline phosphatase, esterase (C4), esterase lipase (C8), leucine arylamidase, and β-glucosidase. In the API 20NE assay, positive for esculin hydrolysis. The predominant respiratory quinone is ubiquinone 8 (Q-8). The major fatty acids are C16:0, C17:0 cyclo, summed feature 3 (consisting of C16:1 ω7c and/or C16:1 ω6c), and summed feature 8 (consisting of C18:1 ω7c and/or C18:1 ω6c). The polar lipid profile consists of diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), phosphatidylethanolamine (PE), one unidentified phosphoglycolipid, one unidentified lipid, and three unidentified glycolipids. The G + C content is 68.6%. The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA gene sequence and the assembled genome sequence of strain AW1T are MN498045 and JAEQNA000000000, respectively.The type strain AW1T (= KACC 21544T = NBRC 114862T) was isolated from soil at Aewol, Jeju Island, Republic of Korea.Description of Ramlibacter albus sp. nov
Ramlibacter albus (al’bus. L. masc. adj. albus, white, referring to the white colonies of the strain).Strain GTP1T is non-motile, Gram-negative, strictly aerobic, positive for oxidase activity, and negative for catalase activity. When observed on R2A, colonies are white, convex, with entire margins, and 1–2 mm in diameter. Under TEM, cells lack flagella, are rod-shaped, and have a width of 0.7–0.8 μm and length of 1.6–1.9 μm. Growth of strain GTP1T can be observed at 10–45 °C (optimum, 30 °C), at pH 5.0–8.0 (optimum, pH 7.0), and at NaCl concentrations of 0–2% (optimum, 0%). The strain shows positive results for Tween 20 and Tween 80 hydrolysis. GTP1T does not produce siderophores when tested on CAS-blue agar. According to API ZYM results, strain GTP1T is positive for alkaline phosphatase, esterase (C4), esterase lipase (C8), and leucine arylamidase, while the API 20NE assay results show negative results for all substrates. The predominant respiratory quinone is ubiquinone 8 (Q-8). The major fatty acids are C16:0 and summed feature 3 (consisting of C16:1 ω7c and/or C16:1 ω6c). The polar lipid profile consists of diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), phosphatidylethanolamine (PE), and two unidentified phosphoaminolipids. The predominant respiratory quinone is ubiquinone 8 (Q-8). The major fatty acids are C16:0, C17:0 cyclo, summed feature 3 (consisting of C16:1 ω7c and/or C16:1 ω6c), and summed feature 8 (consisting of C18:1 ω7c and/or C18:1 ω6c). The polar lipid profile consists of diphosphatidylglycerol (DPG), phosphatidylglycerol (PG), phosphatidylethanolamine (PE), one unidentified phosphoaminolipid, one unidentified phosphoglycolipid, one unidentified phosphoglycoaminolipid, and two unidentified polar lipids. The G + C content is 67.9%. The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA gene sequence and the assembled genome sequence of strain GTP1T are MN498046 and JACORU000000000, respectively.The type strain GTP1T (= KACC 21702T = NBRC 114488T) was isolated from soil at Seogwipo, Jeju Island, Republic of Korea.Description of Ramlibacter pallidus sp. nov
Ramlibacter pallidus (pal’li.dus. L. masc. adj. pallidus, pale, referring to the color of the colonies).Cells of strain HM2T are Gram-negative, and positive for both oxidase and catalase activities. When observed on R2A agar, colonies are cream-colored, transparent, 1.0–2.5 mm in diameter, and flat with entire margins. Under TEM, monotrichous flagella are observed, and cells are rod-shaped with a width of 0.4–0.78 μm and length of 1.7–1.8 μm. The strain shows the fastest growth at a temperature range of 25–35 °C and at pH values between 8.0 and 9.0. When NaCl is present, growth is observed at concentrations of 0–3% (w/v), with optimal growth was observed at concentrations of 0–1% (w/v). The strain is not able to tolerate anaerobic conditions. Strain HM2T hydrolyzes Tween 80 and weakly hydrolyzes casein. However, siderophore production cannot be observed when tested on CAS-blue agar. According to API ZYM tests, strain HM2T shows positive results for alkaline phosphatase, esterase (C4), esterase lipase (C8), leucine arylamidase, valine arylamidase, acid phosphatase, and naphthol-AS-BI-phosphohydrolase. In addition, API 20NE tests show positive results for nitrate (NO3) to nitrite (NO2-) reduction and esculin hydrolysis. The G + C content is 69.9%. The GenBank/EMBL/DDBJ accession numbers for the 16S rRNA gene sequence and the assembled genome sequence of strain HM2T are MN498047 and JADDIV000000000, respectively.The type strain HM2T (= KCTC 82557T = NBRC 114489T) was isolated from soil at Seopjikoji, Jeju Island, Republic of Korea. More

Lenoir, J. et al. Species better track climate warming in the oceans than on land. Nat. Ecol. Evol. 4(8), 1044–1059 (2020).Article 

Google Scholar 
Hobbs, R. J., Valentine, L. E., Standish, R. J. & Jackson, S. T. Movers and stayers: Novel assemblages in changing environments. Trends Ecol. Evol. 33, 116–128 (2017).Article 

Google Scholar 
Gilman, S. E., Urban, M. C., Tewksbury, J., Gilchrist, G. W. & Holt, R. D. A framework for community interactions under climate change. Trends Ecol. Evol. 25, 325–331 (2010).Article 

Google Scholar 
Ockendon, N. et al. Mechanisms underpinning climatic impacts on natural populations: Altered species interactions are more important than direct effects. Glob. Change Biol. 20, 2221–2229 (2014).Article 
ADS 

Google Scholar 
Gómez-Aparicio, L., García-Valdés, R., Ruíz-Benito, P. & Zavala, M. A. Disentangling the relative importance of climate, size and competition on tree growth in Iberian forests: Implications for forest management under global change. Glob. Change Biol. 17, 2400–2414 (2011).Article 
ADS 

Google Scholar 
Pecl, G. T. et al. Biodiversity redistribution under climate change: Impacts on ecosystems and human well-being. Science https://doi.org/10.1126/science.aai9214 (2017).Article 

Google Scholar 
Scheffers, B. R. et al. The broad footprint of climate change from genes to biomes to people. Science 354, aaf7671. https://doi.org/10.1126/science.aaf7671 (2016).Article 
CAS 

Google Scholar 
Vergés, A. et al. The tropicalization of temperate marine ecosystems: Climate-mediated changes in herbivory and community phase shifts. Proc. R. Soc. B-Biol. Sci. 281, 20140846. https://doi.org/10.1098/rspb.2014.0846 (2014).Article 

Google Scholar 
Poore, A. G. B. et al. Global patterns in the impact of marine herbivores on benthic primary producers. Ecol. Lett. 15, 912–922. https://doi.org/10.1111/j.1461-0248.2012.01804.x (2012).Article 

Google Scholar 
Bennett, S., Wernberg, T., Harvey, E. S., Santana-Garcon, J. & Saunders, B. J. Tropical herbivores provide resilience to a climate-mediated phase shift on temperate reefs. Ecol. Lett. 18, 714–723 (2015).Article 

Google Scholar 
Vergés, A. et al. Long-term empirical evidence of ocean warming leading to tropicalization of fish communities, increased herbivory and loss of kelp. Proc. Natl. Acad. Sci. 113(48), 13791–13796 (2016).Article 
ADS 

Google Scholar 
Vergés, A. et al. Tropical rabbitfish and the deforestation of a warming temperate sea. J. Ecol. 102, 1518–1527. https://doi.org/10.1111/1365-2745.12324 (2014).Article 

Google Scholar 
Kumagai, N. H. et al. Ocean currents and herbivory drive macroalgae-to-coral community shift under climate warming. Proc. Natl. Acad. Sci. 115, 8990–8995 (2018).Article 
ADS 
CAS 

Google Scholar 
Demko, A. M. et al. Declines in plant palatability from polar to tropical latitudes depend on herbivore and plant identity. Ecology 98, 2312–2321. https://doi.org/10.1002/ecy.1918 (2017).Article 

Google Scholar 
Floeter, S. R., Behrens, M. D., Ferreira, C. E. L., Paddack, M. J. & Horn, M. H. Geographical gradients of marine herbivorous fishes: Patterns and processes. Mar Biol 147, 1435–1447 (2005).Article 

Google Scholar 
Longo, G. O., Hay, M. E., Ferreira, C. E. L. & Floeter, S. R. Trophic interactions across 61 degrees of latitude in the Western Atlantic. Glob. Ecol. Biogeogr. 28, 107–117. https://doi.org/10.1111/geb.12806 (2019).Article 

Google Scholar 
Bolser, R. & Hay, M. Are tropical plants better defended? Palatability and defenses of temperate versus tropical seaweeds. Ecology 77, 2269–2286 (1996).Article 

Google Scholar 
Borer, E. T. et al. Global biogeography of autotroph chemistry: is insolation a driving force?. Oikos 122, 1121–1130. https://doi.org/10.1111/j.1600-0706.2013.00465.x (2013).Article 
CAS 

Google Scholar 
Miranda, T. et al. Convictfish on the move: Variation in growth and trophic niche space along a latitudinal gradient. ICES J. Mar. Sci. https://doi.org/10.1093/icesjms/fsz098%JICESJournalofMarineScience (2019).Article 

Google Scholar 
Linton, S. M. The structure and function of cellulase (endo-β-1, 4-glucanase) and hemicellulase (β-1, 3-glucanase and endo-β-1, 4-mannase) enzymes in invertebrates that consume materials ranging from microbes, algae to leaf litter. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 240, 110354 (2020).Article 
CAS 

Google Scholar 
Poloczanska, E. S. et al. Global imprint of climate change on marine life. Nat. Clim. Change 3, 919–925. https://doi.org/10.1038/nclimate1958 (2013).Article 
ADS 

Google Scholar 
Nakamura, Y., Feary, D. A., Kanda, M. & Yamaoka, K. Tropical fishes dominate temperate reef fish communities within western Japan. PLoS ONE 8, e81107 (2013).Article 
ADS 

Google Scholar 
Tanaka, K., Taino, S., Haraguchi, H., Prendergast, G. & Hiraoka, M. Warming off southwestern Japan linked to distributional shifts of subtidal canopy-forming seaweeds. Ecol. Evol. 2, 2854–2865. https://doi.org/10.1002/ece3.391 (2012).Article 

Google Scholar 
Pessarrodona, A. et al. Homogenization and miniaturization of habitat structure in temperate marine forests. Glob. Change Biol. 27, 5262–5275 (2021).Article 
CAS 

Google Scholar 
Yamano, H., Sugihara, K. & Nomura, K. Rapid poleward range expansion of tropical reef corals in response to rising sea surface temperatures. Geophys. Res. Lett. 38, L04601. https://doi.org/10.1029/2010gl046474 (2011).Article 
ADS 

Google Scholar 
Mezaki, T. & Kubota, S. Changes of hermatypic coral community in coastal sea area of Kochi, high-latitude Japan. Aquabiology 201, 332–337 (2012).
Google Scholar 
Serisawa, Y., Imoto, Z., Ishikawa, T. & Ohno, M. Decline of the Ecklonia cava population associated with increased seawater temperatures in Tosa Bay, southern Japan. Fish Sci 70, 189–191. https://doi.org/10.1111/j.0919-9268.2004.00788.x (2004).Article 
CAS 

Google Scholar 
Kiriyama, T., Mitsunaga, N., Yasumoto, S., Fujii, A. & Yotsui, T. Undergrown phenomenon of brown alga, Hizikia fusiformis, thought to be caused by grazing of herbivores at Tsutsuura, Tsushima Islands [Japan]. Bulletin of Nagasaki Prefectural Institute of Fisheries (Japan) (1999).Kiriyama, T., Fujii, A. & Fujita, Y. Feeding and characteristic bite marks on Sargassum fusiforme by several herbivorous fishes. Aquac. Sci. 53, 355–365 (2005).
Google Scholar 
Yatsuya, K., Kiriyama, T., Kiyomoto, S., Taneda, T. & Yoshimura, T. On the deterioration process of Ecklonia and Eisenia beds observed in 2013 at Gounoura, Iki Island, Nagasaki Prefecture, Japan.-Initiation of the bed degradation due to high water temperature in summer and subsequent cascading effect by the grazing of herbivorous fish in autumn. Algal Resour. 7, 79–94 (2014).
Google Scholar 
Noda, M., Ohara, H., Murase, N., Ikeda, I. & Yamamoto, K. The grazing of Eisenia bicyclis and several species of Sargassaceous and Cystoseiraceous seaweeds by Siganus fuscescens in relation to the differences of species composition of their seaweed beds. Nippon Suisan Gakkaishi 80, 201–213 (2014).Article 

Google Scholar 
Noda, M., Kinoshita, J., Tanada, N. & Murase, N. Characteristics of bite scars observed in kelp forests of Lessoniaceae denuded by short-term foraging damages of the herbivorous fish Siganus fuscecens. J. Natl. Fish. Univ. 66, 111–122 (2018).
Google Scholar 
Wernberg, T. et al. Seaweed communities in retreat from ocean warming. Curr. Biol. 21, 1828–1832. https://doi.org/10.1016/j.cub.2011.09.028 (2011).Article 
CAS 

Google Scholar 
Terazono, Y., Nakamura, Y., Imoto, Z. & Hiraoka, M. Fish response to expanding tropical Sargassum beds on the temperate coasts of Japan. Mar. Ecol. Prog. Ser. 464, 209–220. https://doi.org/10.3354/meps09873 (2012).Article 
ADS 

Google Scholar 
Duffy, J. E. & Hay, M. E. Seaweed adaptations to herbivory – chemical, structural, and morphological defenses are often adjusted to spatial or temporal patterns of attack. Bioscience 40, 368–375 (1990).Article 

Google Scholar 
Endo, H., Suehiro, K., Kinoshita, J. & Agatsuma, Y. Combined effects of temperature and nutrient enrichment on palatability of the brown alga Sargassum yezoense (Yamada) Yoshida & T. Konno. Am. J. Plant Sci. 6, 275 (2015).Article 
CAS 

Google Scholar 
Clements, K. D., German, D. P., Piché, J., Tribollet, A. & Choat, J. H. Integrating ecological roles and trophic diversification on coral reefs: Multiple lines of evidence identify parrotfishes as microphages. Biol. J. Linn. Soc. 120, 729–751. https://doi.org/10.1111/bij.12914 (2017).Article 

Google Scholar 
Wang, Y., Naumann, U., Wright, S. T. & Warton, D. I. mvabund–an R package for model-based analysis of multivariate abundance data. Methods Ecol. Evol. 3, 471–474 (2012).Article 

Google Scholar 
Wilson, S. K., Bellwood, D. R., Choat, J. H. & Furnas, M. J. Detritus in the epilithic algal matrix and its use by coral reef fishes. Oceanogr. Mar. Biol. Annu. Rev. 41, 279–309 (2003).
Google Scholar 
Helfman, G. S. in The Behaviour of Teleost Fishes 366–387 (Springer, 1986).Prince, J., LeBlanc, W. & Maciá, S. Design and analysis of multiple choice feeding preference data. Oecologia 138, 1–4 (2004).Article 
ADS 

Google Scholar 
Hartig, F. DHARMa: residual diagnostics for hierarchical (multi-level/mixed) regression models. R package version 0.3 3 (2020).Ohno, M. & Ishikawa, M. Physiological ecology of brown alga, Ecklonia on coast of Tosa Bay, southern Japan. I. Seasonal variation of Ecklonia bed. Rep. USA Marine Biol. Inst. Kochi Univ. 4, 59–73 (1982).
Google Scholar 
Agostini, S. et al. Simplification, not “tropicalization”, of temperate marine ecosystems under ocean warming and acidification. Glob. Change Biol. 27, 4771–4784 (2021).Article 
CAS 

Google Scholar 
Clements, K. & Choat, J. Influence of season, ontogeny and tide on the diet of the temperate marine herbivorous fish Odax pullus (Odacidae). Mar. Biol. 117, 213–220 (1993).Article 

Google Scholar 
Mizuta, H., Hayasaki, J. & Yamamoto, H. Relationship between nitrogen content and sorus formation in the brown alga Laminaria japonica cultivated in southern Hokkaido, Japan. Fish. Sci. 64, 909–913 (1998).Article 
CAS 

Google Scholar 
Kumura, T., Yasui, H. & Mizuta, H. Nutrient requirement for zoospore formation in two alariaceous plants Undaria pinnatifida (Harvey) Suringar and Alaria crassifolia Kjellman (Phaeophyceae: Laminariales). Fish. Sci. 72, 860–869 (2006).Article 
CAS 

Google Scholar 
Qiu, Z. et al. Future climate change is predicted to affect the microbiome and condition of habitat-forming kelp. Proc. R. Soc. B 286, 20181887 (2019).Article 

Google Scholar 
Hoey, A. S. & Bellwood, D. R. Limited functional redundancy in a high diversity system: Single species dominates key ecological process on coral reefs. Ecosystems 12, 1316–1328. https://doi.org/10.1007/s10021-009-9291-z (2009).Article 

Google Scholar 
Streit, R. P., Hoey, A. S. & Bellwood, D. R. Feeding characteristics reveal functional distinctions among browsing herbivorous fishes on coral reefs. Coral Reefs 34, 1037–1047 (2015).Article 
ADS 

Google Scholar 
Van Alstyne, K. L. & Paul, V. J. The biogeography of polyphenolic compounds in marine macroalgae – Temperate brown algal defenses deter feeding by tropical herbivorous fishes. Oecologia 84, 158–163 (1990).Article 
ADS 

Google Scholar 
Targett, N. M., Boettcher, A. A., Targett, T. E. & Vrolijk, N. H. Tropical marine herbivore assimilation of phenolic-rich plants. Oecologia 103, 170–179 (1995).Article 
ADS 

Google Scholar 
Prado, P. & Heck, K. L. Seagrass selection by omnivorous and herbivorous consumers: Determining factors. Mar. Ecol. Prog. Ser. 429, 45–55. https://doi.org/10.3354/meps09076 (2011).Article 
ADS 

Google Scholar 
Montgomery, W. L. & Gerking, S. D. Marine macroalgae as foods for fishes: an evaluation of potential food quality. Environ. Biol. Fish. 5, 143–153 (1980).Article 

Google Scholar 
Duffy, J. & Paul & V.J.,. Prey nutritional quality and the effectiveness of chemical defenses against tropical reef fishes. Oecologia 90, 333–339 (1992).Article 
ADS 
CAS 

Google Scholar 
Michael, P. J., Hyndes, G. A., Vanderklift, M. A. & Vergés, A. Identity and behaviour of herbivorous fish influence large-scale spatial patterns of macroalgal herbivory in a coral reef. Mar. Ecol. Prog. Ser. 482, 227–240 (2013).Article 
ADS 

Google Scholar 
Bennett, S. & Bellwood, D. R. Latitudinal variation in macroalgal consumption by fishes on the Great Barrier Reef. Mar. Ecol. Prog. Ser. 426, 241–252 (2011).Article 
ADS 

Google Scholar 
Zarco-Perello, S., Wernberg, T., Langlois, T. J. & Vanderklift, M. A. Tropicalization strengthens consumer pressure on habitat-forming seaweeds. Sci. Rep. 7, 820. https://doi.org/10.1038/s41598-017-00991-2 (2017).Article 
ADS 
CAS 

Google Scholar 
Smith, S. M. et al. Tropicalisation and kelp loss shift trophic composition and lead to more winners than losers in fish communities. Glob. Change Biol. 27(11), 2537–2548 (2021).Article 
ADS 
CAS 

Google Scholar 
Zarco-Perello, S. et al. Range-extending tropical herbivores increase diversity, intensity and extent of herbivory functions in temperate marine ecosystems. Funct. Ecol. 34, 2411–2421. https://doi.org/10.1111/1365-2435.13662 (2020).Article 

Google Scholar  More

Coyne, J. A. & Orr, H. A. Speciation 83–178 (Sinauer, 2004).Dieckmann, U., Doebeli, M., Metz, J. A. & Tautz, D. Adaptive Speciation (Cambridge University Press, 2004).Butlin, R. K., Galindo, J. & Grahame, J. W. Sympatric, parapatric or allopatric: the most important way to classify speciation? Philos. T. Roy. Soc. B 363, 2997–3007 (2008).Article 

Google Scholar 
Smadja, C. M. & Butlin, R. K. A framework for comparing processes of speciation in the presence of gene flow. Mol. Ecol. 20, 5123–5140 (2011).Article 

Google Scholar 
Seehausen, O. et al. Genomics and the origin of species. Nat. Rev. Genet. 15, 176–192 (2014).Article 
CAS 

Google Scholar 
Kulmuni, J., Butlin, R. K., Lucek, K., Savolainen, V. & Westram, A. M. Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers. Philos. T. Roy. Soc. B 375, 20190528 (2020).Article 

Google Scholar 
Hofreiter, M. & Stewart, J. Ecological change, range fluctuations and population dynamics during the Pleistocene. Curr. Biol. 19, R584–R594 (2009).Article 
CAS 

Google Scholar 
Longman, J., Mills, B. J. W., Manners, H. R., Gernon, T. M. & Palmer, M. R. Late Ordovician climate change and extinctions driven by elevated volcanic nutrient supply. Nat. Geosci. 14, 924–929 (2021).Article 
ADS 
CAS 

Google Scholar 
Thomson, R. C., Spink, P. Q. & Shaffer, H. B. A global phylogeny of turtles reveals a burst of climate-associated diversification on continental margins. Proc. Natl Acad. Sci. USA 118, e2012215118 (2021).Article 
CAS 

Google Scholar 
Chaboureau, A. C., Sepulchre, P., Donnadieu, Y. & Franc, A. Tectonic-driven climate change and the diversification of angiosperms. Proc. Natl Acad. Sci. USA 111, 14066–14070 (2014).Article 
ADS 
CAS 

Google Scholar 
Hewitt, G. The genetic legacy of the Quaternary ice ages. Nature 405, 907–913 (2000).Article 
ADS 
CAS 

Google Scholar 
Schmitt, T. Molecular biogeography of Europe: Pleistocene cycles and postglacial trends. Front. Zool. 4, 11 (2007).Article 

Google Scholar 
Haffer, J. Speciation in Amazonian forest birds. Science 165, 131–137 (1969).Article 
ADS 
CAS 

Google Scholar 
Ebdon, S. et al. The Pleistocene species pump past its prime: evidence from European butterfly sister species. Mol. Ecol. 30, 3575–3589 (2021).Article 

Google Scholar 
Excoffier, L., Foll, M. & Petit, R. J. Genetic consequences of range expansions. Annu. Rev. Ecol. Evol. Syst. 40, 481–501 (2009).Article 

Google Scholar 
Baker, H. G. Self-compatibility and establishment after ‘long-distance’ dispersal. Evolution 9, 347–349 (1955).
Google Scholar 
Fisher, R. The Genetical Theory of Natural Selection 125–129 (Oxford University Press, 1930).Endler, J. A. Geographic Variation, Speciation, and Clines. Monographs in Population Biology Vol. 10, 53–65, 142–150 (Princeton University Press, 1977).Doebeli, M. & Dieckmann, U. Speciation along environmental gradients. Nature 421, 259–264 (2003).Article 
ADS 
CAS 

Google Scholar 
Ispolatov, J. & Doebeli, M. Diversification along environmental gradients in spatially structured populations. Evol. Ecol. Res. 11, 295–304 (2009).
Google Scholar 
Rettelbach, A., Servedio, M. R. & Hermisson, J. Speciation in peripheral populations: effects of drift load and mating systems. J. Evol. Biol. 29, 1073–1090 (2016).Article 
CAS 

Google Scholar 
Wright, S. I., Kalisz, S. & Slotte, T. Evolutionary consequences of self-fertilization in plants. Proc. R. Soc. Lond. Ser. B 280, 20130133 (2013).
Google Scholar 
Hu, X.-S. Mating system as a barrier to gene flow. Evolution 69, 1158–1177 (2015).Article 
CAS 

Google Scholar 
Glémin, S. How are deleterious mutations purged? Drift versus nonrandom mating. Evolution 57, 2678–2687 (2003).
Google Scholar 
Warwick, S. I., Francis, A. & Al-Shehbaz, I. A. Brassicaceae: species checklist and database on CD-Rom. Plant Syst. Evol. 259, 249–258 (2006).Article 

Google Scholar 
Warwick, S. I., Al-Shehbaz, I. A. & Sauder, C. A. Phylogenetic position of Arabis arenicola and generic limits of Aphragmus and Eutrema (Brassicaceae) based on sequences of nuclear ribosomal DNA. Can. J. Bot. 84, 269–281 (2006).Article 
CAS 

Google Scholar 
Hohmann, N. et al. Taming the wild: resolving the gene pools of non-model Arabidopsis lineages. BMC Evol. Biol. 14, e224 (2014).Article 

Google Scholar 
Novikova, P. Y. et al. Sequencing of the genus Arabidopsis identifies a complex history of nonbifurcating speciation and abundant trans-specific polymorphism. Nat. Genet. 48, 1077–1082 (2016).Article 
CAS 

Google Scholar 
Perrier, A. & Willi, Y. Intraspecific variation in reproductive barriers between two recently-diverged, allopatric Arabidopsis species. J. Evol. Biol. https://doi.org/10.1111/jeb.14122 (2022). (in press).Griffin, P. C. & Willi, Y. Evolutionary shifts to self-fertilisation restricted to geographic range margins in North American Arabidopsis lyrata. Ecol. Lett. 17, 484–490 (2014).Article 
CAS 

Google Scholar 
Willi, Y., Fracassetti, M., Zoller, S. & Van Buskirk, J. Accumulation of mutational load at the edges of a species range. Mol. Biol. Evol. 35, 781–791 (2018).Article 
CAS 

Google Scholar 
Schmickl, R., Jørgensen, M. H., Brysting, A. K. & Koch, M. A. The evolutionary history of the Arabidopsis lyrata complex: a hybrid in the Amphi-Beringian area closes a large distribution gap and builds up a genetic barrier. BMC Evol. Biol. 10, e98 (2010).Article 

Google Scholar 
Pyhäjärvi, T., Aalto, E. & Savolainen, O. Time scales of divergence and speciation among natural populations and subspecies of Arabidopsis lyrata (Brassicaceae). Am. J. Bot. 99, 1314–1322 (2012).Article 

Google Scholar 
Dyke, A. S. in Quaternary Glaciations – Extent and Chronology, Part II: North America (Elsevier, Amsterdam, 2004).Kirkpatrick, M. & Ravigné, V. Speciation by natural and sexual selection: models and experiments. Am. Nat. 159, S22–S35 (2002).Article 

Google Scholar 
Igic, B., Lande, R. & Kohn, J. R. Loss of self‐incompatibility and its evolutionary consequences. Int. J. Plant Sci. 169, 93–104 (2008).Article 

Google Scholar 
Willi, Y. & Määttänen, K. Evolutionary dynamics of mating system shifts in Arabidopsis lyrata. J. Evol. Biol. 23, 2123–2131 (2010).Article 
CAS 

Google Scholar 
Lucek, K. & Willi, Y. Drivers of linkage disequilibrium across a species’ geographic range. PLoS Genet. 17, e1009477 (2021).Article 
CAS 

Google Scholar 
Pironon, S. et al. Geographic variation in genetic and demographic performance: new insights from an old biogeographical paradigm: the centre-periphery hypothesis. Biol. Rev. 92, 1877–1909 (2017).Article 

Google Scholar 
Encinas-Viso, F., Young, A. G. & Pannell, J. R. The loss of self-incompatibility in a range expansion. J. Evol. Biol. 33, 1235–1244 (2020).Article 

Google Scholar 
Jarne, P. & Auld, J. R. Animals mix it up too: the distribution of self-fertilization among hermaphroditic animals. Evolution 60, 1816–1824 (2006).
Google Scholar 
Foxe, J. P. et al. Reconstructing origins of loss of self-incompatibility and selfing in North American Arabidopsis lyrata: a population genetic context. Evolution 64, 3495–3510 (2010).Article 

Google Scholar 
Koski, M. H., Layman, N. C., Prior, C. J., Busch, J. W. & Galloway, L. F. Selfing ability and drift load evolve with range expansion. Evol. Lett. 3, 500–512 (2019).Article 

Google Scholar 
Prior, C. J. & Busch, J. W. Selfing rate variation within species is unrelated to life‐history traits or geographic range position. Am. J. Bot. 108, 2294–2308 (2021).Article 

Google Scholar 
Skeels, A. & Cardillo, M. Reconstructing the geography of speciation from contemporary biodiversity data. Am. Nat. 193, 240–254 (2019).Article 

Google Scholar 
Sánchez-Castro, D., Perrier, A. & Willi, Y. Reduced climate adaptation at range edges in North American Arabidopsis lyrata. Glob. Ecol. Biogeogr. 31, 1066–1077 (2022).Article 

Google Scholar 
Roessler, K. et al. The genome-wide dynamics of purging during selfing in maize. Nat. Plants 5, 980–990 (2019).Article 
CAS 

Google Scholar 
Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. Preprint at https://arxiv.org/abs/1303.3997 (2013).Hu, T. T. et al. The Arabidopsis lyrata genome sequence and the basis of rapid genome size change. Nat. Genet. 43, 476–481 (2011).Article 

Google Scholar 
Li, H. et al. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).Article 

Google Scholar 
McKenna, A. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).Article 
CAS 

Google Scholar 
Cingolani, P. et al. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 6, 80–92 (2012).Article 
CAS 

Google Scholar 
Alexander, D. H., Novembre, J. & Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664 (2009).Article 
CAS 

Google Scholar 
Danecek, P. et al. The variant call format and VCFtools. Bioinformatics 27, 2156–2158 (2011).Article 
CAS 

Google Scholar 
Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313 (2014).Article 
CAS 

Google Scholar 
Pickrell, J. K. & Pritchard, J. K. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 8, e1002967 (2012).Article 
CAS 

Google Scholar 
Excoffier, L., Dupanloup, I., Huerta-Sánchez, E., Sousa, V. C. & Foll, M. Robust demographic inference from genomic and SNP data. PLoS Genet. 9, e1003905 (2013).Article 

Google Scholar 
Marchi, N. et al. The genomic origins of the world’s first farmers. Cell 185, 1842–1859 (2022).Article 
CAS 

Google Scholar 
Li, H. & Durbin, R. Inference of human population history from individual whole-genome sequences. Nature 475, 493–496 (2011).Article 
CAS 

Google Scholar 
Genete, M., Castric, V. & Vekemans, X. Genotyping and de novo discovery of allelic variants at the Brassicaceae self-incompatibility locus from short-read sequencing data. Mol. Biol. Evol. 7, 1193–1201 (2020).Article 

Google Scholar 
Lynch, M. et al. Genome-wide linkage-disequilibrium profiles from single individuals. Genetics 198, 269–281 (2014).Article 

Google Scholar 
R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria (2021).Paradis, E. & Schliep, K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 35, 526–528 (2019).Article 
CAS 

Google Scholar 
Nychka, D., Furrer, R., Paige, J. & Sain, S. fields: tools for spatial data. R package version 14.1 https://github.com/dnychka/fieldsRPackage (2021).Asquith, W. lmomco—L-moments, censored L-moments, trimmed L-moments, L-comoments, and many distributions. R package version 2.4.7 (2022).Revell, L. J. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol. Evol. 3, 217–223 (2012).Article 

Google Scholar 
Lemon, J. Plotrix: a package in the red light district of R. R. N. 6, 8–12 (2006).
Google Scholar 
Pebesma, E. J. & Bivand, R. S. Classes and methods for spatial data in R. R. N. 5, 9–13 (2005).
Google Scholar 
Bivand, R. S., Pebesma, E. & Gomez-Rubio, V. Applied Spatial Data Analysis with R Second edition (Springer, 2013). More

Daszak, P. Emerging infectious diseases of wildlife-threats to biodiversity and human health. Science 287, 443–449 (2000).Article 
ADS 
CAS 

Google Scholar 
Jones, K. E. et al. Global trends in emerging infectious diseases. Nature 451, 990–993 (2008).Article 
ADS 
CAS 

Google Scholar 
Altizer, S., Ostfeld, R. S., Johnson, P. T. J., Kutz, S. & Harvell, C. D. Climate change and infectious diseases: From evidence to a predictive framework. Science 1979(341), 514–519 (2013).Article 
ADS 

Google Scholar 
Kilpatrick, A. M., Briggs, C. J. & Daszak, P. The ecology and impact of chytridiomycosis: An emerging disease of amphibians. Trends Ecol. Evol. 25, 109–118 (2010).Article 

Google Scholar 
Blehert, D. S. et al. Bat white-nose syndrome: An emerging fungal pathogen?. Science 1979(323), 227–227 (2009).Article 

Google Scholar 
Wilfert, L. et al. Deformed wing virus is a recent global epidemic in honeybees driven by Varroa mites. Science 1979(351), 594–597 (2016).Article 
ADS 

Google Scholar 
Garamszegi, L. Z. Climate change increases the risk of malaria in birds. Glob. Change Biol. 17, 1751–1759 (2011).Article 
ADS 

Google Scholar 
Zamora-Vilchis, I., Williams, S. E. & Johnson, C. N. Environmental temperature affects prevalence of blood parasites of birds on an elevation gradient: Implications for disease in a warming climate. PLoS ONE 7, e39208 (2012).Article 
ADS 
CAS 

Google Scholar 
Harvell, D., Altizer, S., Cattadori, I. M., Harrington, L. & Weil, E. Climate change and wildlife diseases: When does the host matter the most?. Ecology 90, 912–920 (2009).Article 

Google Scholar 
Burge, C. A. et al. Climate change influences on marine infectious diseases: Implications for management and society. Ann. Rev. Mar. Sci. 6, 249–277 (2014).Article 

Google Scholar 
Tracy, A. M., Pielmeier, M. L., Yoshioka, R. M., Heron, S. F. & Harvell, C. D. Increases and decreases in marine disease reports in an era of global change. Proc. R. Soc. B Biol. Sci. 286, 20191718 (2019).Article 

Google Scholar 
Lejeusne, C., Chevaldonné, P., Pergent-Martini, C., Boudouresque, C. F. & Pérez, T. Climate change effects on a miniature ocean: The highly diverse, highly impacted Mediterranean Sea. Trends Ecol. Evol. 25, 250–260 (2010).Article 

Google Scholar 
Basso, L. et al. The Pen Shell, Pinna nobilis: A review of population status and recommended research priorities in the Mediterranean Sea. Adv. Mar. Biol. 71, 109–160 (2015).Article 

Google Scholar 
Catanese, G. et al. Haplosporidium pinnae sp. nov., a haplosporidan parasite associated with mass mortalities of the fan mussel, Pinna nobilis, in the Western Mediterranean Sea. J. Invertebr. Pathol. 157, 9–24 (2018).Article 
CAS 

Google Scholar 
Vázquez-Luis, M. et al. S.O.S. Pinna nobilis: A mass mortality event in western Mediterranean sea. Front. Mar. Sci. 4, 220 (2017).Article 

Google Scholar 
García-March, J. R. et al. Can we save a marine species affected by a highly infective, highly lethal, waterborne disease from extinction?. Biol. Conserv. 243, 108498 (2020).Article 

Google Scholar 
Prado, P. et al. Pinna nobilis in suboptimal environments are more tolerant to disease but more vulnerable to severe weather phenomena. Mar. Environ. Res. 163, 105220 (2021).Article 
CAS 

Google Scholar 
Cabanellas-Reboredo, M. et al. Tracking a mass mortality outbreak of pen shell Pinna nobilis populations: A collaborative effort of scientists and citizens. Sci. Rep. 9, 13355 (2019).Article 
ADS 

Google Scholar 
Kersting, D. K. et al. Recruitment disruption and the role of unaffected populations for potential recovery after the Pinna nobilis mass mortality event. Front. Mar. Sci. 7, 1–11 (2020).Article 
ADS 

Google Scholar 
Box, A. et al. Reduced antioxidant response of the fan mussel Pinna nobilis related to the presence of haplosporidium pinnae. Pathogens 9, 1–14 (2020).Article 

Google Scholar 
Peyran, C., Morage, T., Nebot-Colomer, E., Iwankow, G. & Planes, S. Unexpected residual habitats raise hope for the survival of the fan mussel Pinna nobilis along the Occitan coast (Northwest Mediterranean Sea). Endanger Species Res. 48, 123–137 (2022).Article 

Google Scholar 
Rosa, R. D. et al. A hemocyte gene expression signature correlated with predictive capacity of oysters to survive Vibrio infections. BMC Genomics 13, 1–12 (2012).Article 

Google Scholar 
van de Vijver, M. J. et al. A gene-expression signature as a predictor of survival in breast cancer. N. Engl. J. Med. 347, 1999–2009 (2002).Article 

Google Scholar 
Seppey, M., Manni, M. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness. In Methods in Molecular Biology 227–245 https://doi.org/10.1007/978-1-4939-9173-0_14 (2019).Smith-Unna, R., Boursnell, C., Patro, R., Hibberd, J. M. & Kelly, S. TransRate: Reference-free quality assessment of de novo transcriptome assemblies. Genome Res. 26, 1134–1144 (2016).Article 
CAS 

Google Scholar 
Guo, X. & Ford, S. E. Infectious diseases of marine mollusks and host responses as revealed by genomic tools. Philos. Trans. R. Soc. B Biol. Sci. https://doi.org/10.1098/rstb.2015.0206 (2016).Article 

Google Scholar 
Pauletto, M. et al. Deep transcriptome sequencing of Pecten maximus hemocytes: A genomic resource for bivalve immunology. Fish Shellfish Immunol. 37, 154–165 (2014).Article 
CAS 

Google Scholar 
Caurcel, C. et al. MolluscDB: A genome and transcriptome database for molluscs. Philos. Trans. R. Soc. Lond. B Biol. Sci. 376, 20200157 (2021).Article 
CAS 

Google Scholar 
de Oliveira, A. L. et al. Comparative transcriptomics enlarges the toolkit of known developmental genes in mollusks. BMC Genomics 17, 1–23 (2016).Article 

Google Scholar 
Richardson, M. F. & De Sherman, C. D. H. De novo assembly and characterization of the invasive Northern Pacific Seastar transcriptome. PLoS ONE 10, e0142003 (2015).Article 

Google Scholar 
Zhang, D., Wang, F., Dong, S. & Lu, Y. D. De novo assembly and transcriptome analysis of osmoregulation in Litopenaeus vannamei under three cultivated conditions with different salinities. Gene 578, 185–193 (2016).Article 
CAS 

Google Scholar 
Werner, G. D. A., Gemmell, P., Grosser, S., Hamer, R. & Shimeld, S. M. Analysis of a deep transcriptome from the mantle tissue of Patella vulgata Linnaeus (Mollusca: Gastropoda: Patellidae) reveals candidate biomineralising genes. Mar. Biotechnol. 15, 230–243 (2013).Article 
CAS 

Google Scholar 
Ding, J. et al. Transcriptome sequencing and characterization of Japanese scallop Patinopecten yessoensis from different shell color lines. PLoS ONE 10, e0116406 (2015).Article 

Google Scholar 
Harney, E. et al. De novo assembly and annotation of the European abalone Haliotis tuberculata transcriptome. Mar Genomics 28, 11–16 (2016).Article 

Google Scholar 
Khalturin, K., Hemmrich, G., Fraune, S., Augustin, R. & Bosch, T. C. G. More than just orphans: Are taxonomically-restricted genes important in evolution?. Trends Genet. 25, 404–413. https://doi.org/10.1016/j.tig.2009.07.006 (2009).Article 
CAS 

Google Scholar 
Gibson, A. K., Smith, Z., Fuqua, C., Clay, K. & Colbourne, J. K. Why so many unknown genes? Partitioning orphans from a representative transcriptome of the lone star tick Amblyomma americanum. BMC Genomics 14, 135 (2013).Article 
CAS 

Google Scholar 
Albertin, C. B. et al. The octopus genome and the evolution of cephalopod neural and morphological novelties. Nature 524, 220–224 (2015).Article 
ADS 
CAS 

Google Scholar 
Vogeler, S., Galloway, T. S., Lyons, B. P. & Bean, T. P. The nuclear receptor gene family in the Pacific oyster, Crassostrea gigas, contains a novel subfamily group. BMC Genomics 15, 369 (2014).Article 

Google Scholar 
Allam, B. & Raftos, D. Immune responses to infectious diseases in bivalves. J. Invertebr. Pathol. 131, 121–136. https://doi.org/10.1016/j.jip.2015.05.005 (2015).Article 
CAS 

Google Scholar 
Allam, B. & Pales Espinosa, E. Bivalve immunity and response to infections: Are we looking at the right place?. Fish Shellfish Immunol. 53, 4–12. https://doi.org/10.1016/j.fsi.2016.03.037 (2016).Article 
CAS 

Google Scholar 
Qiu, L., Song, L., Xu, W., Ni, D. & Yu, Y. Molecular cloning and expression of a Toll receptor gene homologue from Zhikong Scallop, Chlamys farreri. Fish Shellfish Immunol. 22, 451–466 (2007).Article 
CAS 

Google Scholar 
Zhang, L., Li, L., Zhu, Y., Zhang, G. & Guo, X. Transcriptome analysis reveals a rich gene set related to innate immunity in the eastern oyster (Crassostrea virginica). Mar. Biotechnol. 16, 17–33 (2014).Article 

Google Scholar 
Moreira, R. et al. Transcriptomics of in vitro immune-stimulated hemocytes from the Manila clam Ruditapes philippinarum using high-throughput sequencing. PLoS ONE 7, e35009 (2012).Article 
ADS 
CAS 

Google Scholar 
Toubiana, M. et al. Toll-like receptors and MyD88 adaptors in Mytilus: Complete cds and gene expression levels. Dev. Comp. Immunol. 40, 158–166 (2013).Article 
CAS 

Google Scholar 
He, Y. et al. Transcriptome analysis reveals strong and complex antiviral response in a mollusc. Fish Shellfish Immunol. 46, 131–144 (2015).Article 
CAS 

Google Scholar 
Zhang, L. et al. Massive expansion and functional divergence of innate immune genes in a protostome. Sci. Rep. 5, 8693 (2015).Article 
CAS 

Google Scholar 
Casadevall, A. & Pirofski, L. A. Host–pathogen interactions: The attributes of virulence. J. Infect. Dis. 184, 337–344. https://doi.org/10.1086/322044 (2001).Article 
CAS 

Google Scholar 
Jones, B., Pascopella, L. & Falkow, S. Entry of microbes into the host: Using M cells to break the mucosal barrier. Curr. Opin. Immunol. 7, 474–478 (1995).Article 
CAS 

Google Scholar 
Liévin-Le Moal, V. & Servin, A. L. The front line of enteric host defense against unwelcome intrusion of harmful microorganisms: Mucins, antimicrobial peptides, and microbiota. Clin. Microbiol. Rev. 19, 315–337. https://doi.org/10.1128/CMR.19.2.315-337.2006 (2006).Article 
CAS 

Google Scholar 
Trigos, S., Vicente, N., Prado, P. & Espinós, F. J. Adult spawning and early larval development of the endangered bivalve Pinna nobilis. Aquaculture 483, 102–110 (2018).Article 

Google Scholar 
Vázquez-Luis, M., Nebot-Colomer, E., Deudero, S., Planes, S. & Boissin, E. Natural hybridization between pen shell species: Pinna rudis and the critically endangered Pinna nobilis may explain parasite resistance in P. nobilis. Mol. Biol. Rep. https://doi.org/10.1007/s11033-020-06063-5 (2021).Article 

Google Scholar 
Katsares, V., Tsiora, A., Galinou-Mitsoudi, S. & Imsiridou, A. Genetic structure of the endangered species Pinna nobilis (Mollusca: Bivalvia) inferred from mtDNA sequences. Biologia 63, 412–417 (2008).Article 
CAS 

Google Scholar 
Gonzalez-Wanguemert, M. et al. Highly polymorphic microsatellite markers for the Mediterranean endemic fan mussel Pinna nobilis. Mediterr. Mar. Sci. 16, 31 (2014).Article 

Google Scholar 
Peyran, C., Planes, S., Tolou, N., Iwankow, G. & Boissin, E. Development of 26 highly polymorphic microsatellite markers for the highly endangered fan mussel Pinna nobilis and cross-species amplification. Mol. Biol. Rep. 47, 2551–2559 (2020).Article 
CAS 

Google Scholar 
Peakall, R. & Smouse, P. E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—an update. Bioinformatics 28, 2537–2539 (2012).Article 
CAS 

Google Scholar  More