Ecology

A single long day induces the photoperiodic molecular responseFigure 1c shows results from the experiment 1, as evidenced from the qPCR measurement of mRNA expression of genes of known biological functions in the blood and hypothalamus. Clearly, the exposure to extended light period induced a molecular response by hour 18 of the first long day, as shown by change in mRNA levels of candidate genes in both central (hypothalamus) peripheral (blood) tissues of photosensitive buntings. Blood mRNA levels of peroxiredoxin 4 (prdx4) were significantly lower at hour 18 mimicking a long 18 h photoperiod than those at hour 10 mimicking a short 10 h photoperiod (p = 0.002, t = 5.18, n = 4/time point). Paradoxically, this indicated a reduced cellular response against oxidative stress in the otherwise photo stimulated birds on the first long day. We speculate that prdx4 expression pattern would be inversed (i.e. increased prdx4 mRNA levels) after several long days when birds show photoperiodically stimulated hyperphagia (increased food intake) and lipogenesis (fat accumulation). Intriguingly, however, blood mRNA levels of gpx1 (p = 0.399, t = 0.91, n = 4/time point) and sod1 (p = 0.845, t = 0.20, n = 4/time point) genes were not different between hours 10 and 18 (Student’s t-test, Fig. 1c(a–c)). Taken together differences in the expression pattern of these enzymes, we speculate differential activation of the enzymatic pathways that are probably involved in the oxidative cellular response when migratory birds are exposed to an acute change in their photoperiodic environment.On the other hand, blood il1β mRNA levels were significantly higher at hour 18 than the hour 10 (p = 0.041, t = 2.58, n = 4/time point; Student’s t-test, Fig. 1c(d)). It is consistent with the known role of il1β-encoded interleukin 1β, as a crucial mediator of the inflammation and a marker of the innate immune system22,23. Increased il1β mRNA expression on the first long day is consistent with the idea of parallel photoperiodic induction of multiple biological processes, including those associated with the innate immune response, body fattening and gonadal maturation in migratory songbirds28; however, the possibility that an upregulated interleukin was an indicative a stress response cannot be excluded at this time.Changes in hypothalamic gene expressions further confirm a rapid molecular response to the extended light period when it surpasses the threshold photoperiod, i.e. acts as the stimulatory long day. Reciprocal switching of genes involved in the thyroid hormone responsive pathway at hour 18 particularly evidences this. Hypothalamic mRNA levels of tshβ (p = 0.033, t = 2.75, n = 4/time point) and dio2 (p = 0.0004, t = 7.14, n = 4/time point) genes were higher, and that of dio3 gene expression was lower at hour 18 than the hour 10 (p = 0.036, t = 2.68, n = 4/time point). This is also in agreement with the rapid photoperiodic response found on the first long day in plasma LH secretion, and in hypothalamic expressions of Fos-immunoreactivity and thyroid hormone responsive genes in blackheaded buntings14,33 and other photoperiodic birds15,17,19,32,34,35,36,37,38. However, gnrh mRNA levels were not found significantly different between hours 10 and 18 of the first long day (p = 0.324, t = 1.07, n = 4/time point; Student’s t-test, Fig. 1c(e–h) indicating that hour 18 was probably too early a time for an upregulated gnrh expression on the first long day37,38,39.RNA-Seq reveals differences in time course of the photoperiodic responseTable S2 summarizes the primary statistics used for RNA-Seq results. Using only transcripts with non-zero abundance, we compared the time course of transcriptome-wide response in the hypothalamus both as the function of time (within photosensitive or photorefractory state) and LHS (photosensitive vs. photorefractory state; n = 2/time point/state except at hour 22 in photorefractory state which had n = 1 sample size). Further, to show a functional linkage of differentially expressed genes (DEGs), we performed STRING analysis that predicts the protein–protein interaction (see methods for details).Results on hypothalamic gene expressions suggest that buntings react to the acute photoperiodic change in photorefractory state almost as they do in the photosensitive state. However, the comparison of the overall RNASeq data from both states revealed LHS-dependent pattern in the time course of transcriptional response, with differences in the number and functions of DEGs and associated physiological pathways.Within state differences in time course of transcriptional responseWe examined the time course of response on the first long day, by comparing gene expressions at the hours 14, 18 and 22 of the extended light period that mimicked 14 h, 18 h and 22 h long photoperiods, respectively, with those at hour 10 that mimicked a 10 h short photoperiod.Photosensitive stateAt hour 14, we found 10 differentially expressed genes (DEGs) with 4 upregulated and 6 downregulated genes (Figs. 2a, 3a, Table S3). Of the 10 DEGs, atp6v1e1, atp6v1b2, uqcrc1 and pgam1 genes enriched the oxidative phosphorylation, metabolic pathways, phagosome and mTOR signalling pathways (Table 1). The oxidative phosphorylation and metabolic pathways were upregulated at hour 10, while the phagosome and mTOR signalling pathways were enriched by two genes that were opposite in the expression trend: atp6v1e1 was upregulated while atp6v1b2 was downregulated at hour 14. The STRING analysis showed a significant interaction of atp6v1e1 and atp6v1b2 encoded proteins (ATP6V1E1 and ATP6V1B2). These proteins are the components of vacuolar ATPase enzyme that mediates the acidification of eukaryotic intracellular organelles necessary for protein sorting and zymogen activation. Further, at hour 14, ttr gene that codes for transthyretin (a preferential T3 binder) and pomc gene that codes for the proopiomelanocortin receptor had significantly lower expressions. Whereas, low ttr gene expression, as in photostimulated redheaded buntings40, might indicate a reduced trafficking of thyroid hormones via ttr-encoded transthyretins in the photosensitive state, the low pomc gene expression might suggest the removal of inhibitory effects of the opioids (e.g. β-endorphin, a pomc-encoded proopiomelanocortin product) on hypothalamic GnRH and, in turn, pituitary LH secretion41,42.Figure 2Top panel: Volcano plots showing results of differential gene expression analysis (− log10 padj. vs. log2 fold change values) in the hypothalamus within the photosensitive (a–c) and photorefractory states (e–g). The comparison protocol is shown on the left. In each state, the comparisons were done with respect to the hour 10 value (akin to short day control). Venn diagram shows common and unique DEGs in photosensitive (d) and photorefractory states (h). Bottom panel: Volcano plots showing results of differential gene expression analysis (− log10 padj. vs. log2 fold change values) between the photosensitive and photorefractory states. The pairwise comparisons were made at all the four time points (hours 10 (i), 14 (j), 18 (k) and 22 (l)). Venn diagram shows common and unique DEGs between states at hours 10, 14, 18 and 22 (m). Genes in a volcano plot with log2 fold change  > 2 are marked by green colour, and those with log2 fold change  > 2 and p value (padj.)  More

1.Rothschild, L. & Mancinelli, R. Life in extreme environments. Nature 409, 1092–1101 (2001).ADS 
CAS 
PubMed 
Article 

Google Scholar 
2.Cavicchioli, R. Extremophiles and the search for extraterrestrial life. Astrobiology 2, 281–292 (2002).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
3.Parro, V. et al. A microbial oasis in the hypersaline atacama subsurface discovered by a life detector chip: Implications for the search for life on mars. Astrobiology 11, 969–996 (2011).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
4.Lee, C. J. D. et al. NaCl-saturated brines are thermodynamically moderate, rather than extreme, microbial habitats. FEMS Microbiol. Rev. 42, 672–693 (2018).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
5.Coleine, C. et al. Specific adaptations are selected in opposite sun exposed Antarctic cryptoendolithic communities as revealed by untargeted metabolomics. PLoS ONE 15, 1–17 (2020).Article 
CAS 

Google Scholar 
6.Rathour, R. et al. A comparative metagenomic study reveals microbial diversity and their role in the biogeochemical cycling of Pangong lake. Sci. Total Environ. 731, 139074 (2020).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
7.Suosaari, E. P. et al. New multi-scale perspectives on the stromatolites of Shark Bay, Western Australia. Sci. Rep. 6, 1–13 (2016).Article 
CAS 

Google Scholar 
8.Suosaari, E. P. et al. Stromatolite provinces of Hamelin pool: Physiographic controls on stromatolites and associated lithofacies. J. Sediment. Res. 89, 207–226 (2019).ADS 
Article 

Google Scholar 
9.Wong, H., Ahmed-Cox, A. & Burns, B. Molecular ecology of hypersaline microbial mats: Current insights and new directions. Microorganisms 4, 6 (2016).PubMed Central 
Article 
CAS 

Google Scholar 
10.Grotzinger, J. R. & Knoll, A. H. Stromatolites in Precambrian carbonates: Evolutionary mileposts or environmental dipsticks?. Annu. Rev. Earth Planet. Sci. 27, 313–358 (1999).ADS 
CAS 
PubMed 
Article 

Google Scholar 
11.Grotzinger, J. P. & James, N. P. Precambrian carbonates: Evolution of understanding. In Carbonate Sedimentation and Diagenesis in the Evolving Precambrian World (eds Grotzinger, J. P. & James, N. P.) 3–20 (Society for Sedimentary Geology, 2000).Chapter 

Google Scholar 
12.Demergasso, C. et al. Microbial mats from the Llamara salt flat, northern Chile. Rev. Chil. Hist. Nat. 76, 485–499 (2003).Article 

Google Scholar 
13.Demergasso, C. et al. Distribution of prokaryotic genetic diversity in athalassohaline lakes of the Atacama Desert, Northern Chile. FEMS Mirobiol. Ecol. 48, 57–69 (2004).CAS 
Article 

Google Scholar 
14.Saghaï, A. et al. Unveiling microbial interactions in stratified mat communities from a warm saline shallow pond. Environ. Microbiol. 19, 2405–2421 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
15.Shen, J., Zerkle, A. L., Stueeken, E. & Claire, M. W. Nitrates as a potential N supply for microbial ecosystems in a hyperarid mars analog system. Life 9, 79 (2019).PubMed Central 
Article 
CAS 
PubMed 

Google Scholar 
16.Albarracín, V. H., Galván, F. S. & Farías, M. E. Extreme microbiology at Laguna Socompa: A high-altitude Andean lake (3570 m a.s.l.) in Salta, Argentina. In Microbial Ecosystems in Central Andes Extreme Environments: Biofilms, Microbial Mats, Microbialites and Endoevaporites (ed. Farías, M. E.) 205–220 (Springer, 2020).Chapter 

Google Scholar 
17.Aszalós, J. M. et al. Bacterial diversity of a high-altitude permafrost thaw pond located on Ojos del Salado (Dry Andes, Altiplano-Atacama Region). Astrobiology 20, 754–765 (2020).ADS 
PubMed 
Article 
CAS 

Google Scholar 
18.Boidi, F. J., Mlewski, E. C., Gomez, F. J. & Gérard, E. Characterization of microbialites and microbial mats of the Laguna Negra hypersaline lake (Puna of Catamarca, Argentina). In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 183–203 (Springer, 2020).Chapter 

Google Scholar 
19.Farías, M. E. & Saona Acuña, L. A. Modern microbial mats and endoevaporite systems in Andean lakes: A general approach. In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 21–33 (Springer, 2020).Chapter 

Google Scholar 
20.Farías, M. E., Villafañe, P. G. & Lencina, A. I. Integral propsection of andean microbial ecosystem project. In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 245–260 (Springer, 2020).Chapter 

Google Scholar 
21.Gomez, F. J., Boidi, F. J., Mlewski, E. C. & Gérard, E. The carbonate system in Hypersaline Lakes: The case of Laguna Negra (in the Puna Region of Catamarca, Argentina). In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 231–242 (Springer, 2020).Chapter 

Google Scholar 
22.Otálora, F. et al. Hydrochemical and mineralogical evolution through evaporitic processes in Salar de Llamara Brines (Atacama, Chile). ACS Earth Sp. Chem. 4, 882–896 (2020).Article 
CAS 

Google Scholar 
23.Rasuk, M. C., Visscher, P. T., Leiva, M. C. & Farías, M. E. Mats and microbialites from Laguna La Brava. In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 221–230 (Springer, 2020).Chapter 

Google Scholar 
24.Demergasso, C. et al. Novelty and spatio-temporal heterogeneity in the bacterial diversity of hypersaline Lake Tebenquiche (Salar de Atacama). Extremophiles 12, 491–504 (2008).CAS 
PubMed 
Article 

Google Scholar 
25.del Rocío Mora-Ruiz, M. & Díaz-Gil, C. Microbial diversity in athalassohaline Argentinean Salterns. In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 165–179 (Springer, 2020).Chapter 

Google Scholar 
26.Vignale, F. A. et al. Lithifying and non-lithifying microbial ecosystems in the wetlands and salt flats of the central Andes. Microb. Ecol. https://doi.org/10.1007/s00248-021-01725-8 (2021).Article 
PubMed 

Google Scholar 
27.Stivaletta, N., Barbieri, R., Cevenini, F. & López-García, P. Physicochemical conditions and microbial diversity associated with the evaporite deposits in the Laguna de la Piedra (Salar de Atacama, Chile). Geomicrobiol. J. 28, 83–95 (2011).CAS 
Article 

Google Scholar 
28.Farías, M. E. et al. Characterization of bacterial diversity associated with microbial mats, gypsum evaporites and carbonate microbialites in thalassic wetlands: Tebenquiche and La Brava, Salar de Atacama, Chile. Extremophiles 18, 311–329 (2014).PubMed 
Article 
CAS 

Google Scholar 
29.Fernandez, A. B. et al. Microbial diversity in sediment ecosystems (evaporites domes, microbial mats, and crusts) of Hypersaline Laguna Tebenquiche, Salar de Atacama, Chile. Front. Microbiol. 7, 1–18 (2016).Article 

Google Scholar 
30.Rasuk, M. C. et al. Bacterial diversity in microbial mats and sediments from the Atacama Desert. Microb. Ecol. 71, 44–56 (2016).CAS 
PubMed 
Article 

Google Scholar 
31.Farias, M. E. et al. Prokaryotic diversity and biogeochemical characteristics of benthic microbial ecosystems at La Brava, a hypersaline lake at Salar de Atacama, Chile. PLoS ONE 12, 1–25 (2017).Article 
CAS 

Google Scholar 
32.Gutiérrez-Preciado, A. et al. Functional shifts in microbial mats recapitulate early Earth metabolic transitions. Nat. Ecol. Evol. 2, 1700–1708 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
33.Escudero, L. et al. A thiotrophic microbial community in an acidic brine lake in Northern Chile. Antonie Van Leeuwenhoek 111, 1403–1419 (2018).PubMed 
Article 
PubMed Central 

Google Scholar 
34.Chong-Díaz, G. Die Salare in Nordchile—Geologie, Struktur und Geochemie. Geotekton. Forsch. 67, 1–146 (1984).
Google Scholar 
35.Risacher, F. & Fritz, B. Geochemistry of Bolivian salars, Lipez, southern Altiplano: Origin of solutes and brine evolution. Geochim. Cosmochim. Acta 55, 687–705 (1991).ADS 
CAS 
Article 

Google Scholar 
36.Pueyo, J. J., Chong, G. & Jensen, A. Neogene evaporites in Desert volcanic environments: Atacama Desert, northern Chile. Sedimentology 48, 1411–1431 (2001).ADS 
CAS 
Article 

Google Scholar 
37.Simicic Hernández, Y. P. Thickness Distribution of the Oligo-Neogenous Sedimentary Cover of the Tamarugal Pampas, Northern Chile (20 ° 45 ’to 21 ° 30’S) (Universidad de Chile, 2015).
Google Scholar 
38.Cabrera, S., Bozzo, S. & Fuenzalida, H. Variations in UV radiation in Chile. J. Photochem. Photobiol. B Biol. 28, 137–142 (1995).CAS 
Article 

Google Scholar 
39.Cabrol, N. A. et al. Life in the Atacama: Searching for life with rovers (science overview). J. Geophys. Res. Biogeosci. https://doi.org/10.1029/2006JG000298 (2007).Article 

Google Scholar 
40.Solari, M. The unexplored geobiological heritage of Chile: Key to understand the past and future. In XIV Congr. Geológico Chil. 1–5 (2015).41.Rasuk, M. C. et al. Microbial characterization of microbial ecosystems associated to evaporites domes of gypsum in Salar de Llamara in Atacama Desert. Microb. Ecol. 68, 483–494 (2014).CAS 
PubMed 
Article 

Google Scholar 
42.Surma, J., Assonov, S., Herwartz, D., Voigt, C. & Staubwasser, M. The evolution of 17O-excess in surface water of the arid environment during recharge and evaporation. Sci. Rep. 8, 4972 (2018).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
43.Rasuk, M. C., Leiva, M. C., Kurth, D. & Farías, M. E. Complete characterization of stratified ecosystems of the Salar de Llamara (Atacama Desert). In Microbial Ecosystems in Central Andes Extreme Environments (ed. Farías, M. E.) 153–164 (Springer, 2020).Chapter 

Google Scholar 
44.Kiefer, E., Dorr, M., Ibbeken, H. & Gotze, H. Gravity-based mass balance of an alluvial fan giant: The Arcas Fan, Pampa del Tamarugal, Northern Chile. Rev. Geol. Chile 24, 165–185 (1997).
Google Scholar 
45.Houston, J. & Hartley, A. J. The central Andean west-slope rainshadow and its potential contribution to the origin of hyper-aridity in the Atacama Desert. Int. J. Climatol. 23, 1453–1464 (2003).Article 

Google Scholar 
46.Dunai, T. J., López, G. A. G. & Juez-Larré, J. Oligocene-Miocene age of aridity in the Atacama Desert revealed by exposure dating of erosion-sensitive landforms. Geology 33, 321–324 (2005).ADS 
Article 

Google Scholar 
47.Hartley, A. J., Chong, G., Houston, J. & Mather, A. 150 million years of climatic stability: Evidence from the Atacama Desert, northern Chile. J. Geol. Soc. Lond. 162, 421–424 (2005).Article 

Google Scholar 
48.Clarke, J. D. A. Antiquity of aridity in the Chilean Atacama Desert. Geomorphology 73, 101–114 (2006).ADS 
Article 

Google Scholar 
49.Houston, J. Evaporation in the Atacama Desert: An empirical study of spatio-temporal variations and their causes. J. Hydrol. 330, 402–412 (2006).ADS 
Article 

Google Scholar 
50.Fuenzalida, H. & Rutllant, J. Estudio Sobre el Origen del Vapor de agua que Precipita en el Invierno Altiplánico (1986).51.Grosjean, M., Geyh, M. A., Messerli, B. & Schotterer, U. Late-glacial and early Holocene lake sediments, ground-water formation and climate in the Atacama Altiplano 22–24°S. J. Paleolimnol. 14, 241–252 (1995).ADS 
Article 

Google Scholar 
52.Garreaud, R. Multiscale analysis of the summertime precipitation over the central Andes. Mon. Weather Rev. 127, 901–921 (1999).ADS 
Article 

Google Scholar 
53.Houston, J. Groundwater recharge through an alluvial fan in the Atacama Desert, northern Chile: Mechanisms, magnitudes and causes. Hydrol. Process. 16, 3019–3035 (2002).ADS 
Article 

Google Scholar 
54.Marazuela, M. A., Vázquez-Suñé, E., Ayora, C., García-Gil, A. & Palma, T. Hydrodynamics of salt flat basins: The Salar de Atacama example. Sci. Total Environ. 651, 668–683 (2019).ADS 
CAS 
PubMed 
Article 

Google Scholar 
55.Cereceda, P., Larrain, H., Osses, P., Farías, M. & Egaña, I. The spatial and temporal variability of fog and its relation to fog oases in the Atacama Desert, Chile. Atmos. Res. 87, 312–323 (2008).Article 

Google Scholar 
56.del Río, C. et al. The role of topography in the spatial distribution of the low stratocumulus cloud and fog in the Peruvian coastal Desert. In AGU Fall Meeting Abstracts 2018, A31J-2979 (2018).57.Hasler, K., Jaque, I., Pucheu, A. & Ortiz, C. Análisis de la Información Histórica de la Operación de la Medida de Mitigación. Estudio de Impacto Ambiental: Modification parcial del Sistema del Sistema de Reinyección en los puquios de Llamara, Elaborado por Geobiota (2020).58.Ordoñez, R., Hasler, K., Pucheu, A. & Ortiz, C. Modelo Numérico Hidrogeológico Acuífero Salar de Llamara. Estudio de Impacto Ambiental, Modificación Parcial del Sistema de reinyección en los Puquios de Llamara, elaborado por Geobiota (2020).59.Babel, M. Models for evaporite, selenite and gypsum microbialite deposition in ancient saline basins. Acta Geol. Pol. 54, 219-U6 (2004).
Google Scholar 
60.Rumrich, U., Lange-Bertalot, H. & Rumrich, M. Diatoms of the Andes. Annotated diatom monographs. Iconogr. Diatomol. 9, 671 (2000).
Google Scholar 
61.Lowe, R. L. Keeled and canalled raphid diatoms. In Freshwater Algae of North America (ed. Lowe, R. L.) 669–684 (Elsevier, 2003).Chapter 

Google Scholar 
62.Whitton, B. A. & Kelly, M. G. Use of algae and other plants for monitoring rivers. Aust. J. Ecol. 20, 45–56 (1995).Article 

Google Scholar 
63.Burow, L. C. et al. Identification of Desulfobacterales as primary hydrogenotrophs in a complex microbial mat community. Geobiology 12, 221–230 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
64.Garcés, I. et al. Características geoquímicas generales del sistema salino del Salar de Llamara (Chile). Estud. Geol. 52, 23–35 (1996).Article 

Google Scholar 
65.López, P. L., Auqué, L. F., Garcés, I. & Chong, G. Características geoquímicas y pautas de evolución de las salmueras superficiales del Salar de Llamara, Chile. Rev. Geol. de Chile 26, 89–108 (1999).Article 

Google Scholar 
66.Kampf, S. K. & Tyler, S. W. Spatial characterization of land surface energy fluxes and uncertainty estimation at the Salar de Atacama, Northern Chile. Adv. Water Resour. 29, 336–354 (2006).ADS 
Article 

Google Scholar 
67.Des-Marais, D. J. The Biogeochemistry of Hypersaline Microbial Mats. In Advances in Microbial Ecology (ed. Jones, J. G.) (Springer, 1995).
Google Scholar 
68.Vogel, M. B. et al. The role of biofilms in the sedimentology of actively forming gypsum deposits at Guerrero Negro, Mexico. Astrobiology 9, 875–893 (2009).ADS 
CAS 
PubMed 
Article 

Google Scholar 
69.Vogel, M. B. et al. Biological influences on modern sulfates: Textures and composition of gypsum deposits from Guerrero Negro, Baja California Sur, Mexico. Sediment. Geol. 223, 265–280 (2010).ADS 
CAS 
Article 

Google Scholar 
70.Ali-Bik, M. W., Metwally, H. I. M., Wali, A. M. A. & Kamel, M. G. Facies and geochemistry of non-marine gypsum, EMISAL, Egypt. Geol. Acta 11, 409–420 (2013).CAS 

Google Scholar 
71.Taher, A. G. Formation and calcification of modern gypsum-dominated stromatolites, EMISAL, Fayium, Egypt. Facies 60, 721–735 (2014).Article 

Google Scholar 
72.Handford, C. Sedimentology and evaporite genesis in a Holocene continental-sabkha playa basin—Bristol Dry Lake, California. Sedimentology 29, 239–253 (1982).ADS 
Article 

Google Scholar 
73.Gerdes, G., Krumbein, W. E. & Holtkamp, E. Salinity and water activity related zonation of microbial communities and potential stromatolites of the Gavish Sabkha. In Hypersaline Ecosystems. Ecological Studies (Analysis and Synthesis) (eds Friedman, G. M. & Krumbein, W. E.) 238–236 (Springer, 1985).Chapter 

Google Scholar 
74.Davie, A. W., Mitrovic, S. M. & Lim, R. P. Succession and accrual of benthic algae on cobbles of an upland river following scouring. Inl. Waters 2, 89–100 (2012).Article 

Google Scholar 
75.Cohen, Y., Jørgensen, B. B., Padan, E. T. & Shilo, M. Sulphide-dependent anoxygenic photosynthesis in the cyanobacterium Oscillatoria limnetica. Nature 257, 489–492 (1975).ADS 
CAS 
Article 

Google Scholar 
76.Oren, A. & Shilo, M. Anaerobic heterotrophic dark metabolism in the cyanobacterium Oscillatoria limnetica: Sulfur respiration and lactate fermentation. Arch. Microbiol. 122, 77–84 (1979).CAS 
Article 

Google Scholar 
77.Muñoz, J., Amat, F., Green, A. J., Figuerola, J. & Gómez, A. Bird migratory flyways influence the phylogeography of the invasive brine shrimp Artemia franciscana in its native American range. PeerJ 2013, 1–28 (2013).
Google Scholar 
78.Clegg, J. S. & Trotman, C. N. A. Physiological and biochemical aspects of Artemia ecology. In Artemia: Basic and Applied Biology (eds Abatzopoulos, T. J. et al.) 129–170 (Springer, 2002).Chapter 

Google Scholar 
79.Collado, G. A., Valladares, M. A. & Méndez, M. A. Hidden diversity in spring snails from the andean altiplano, the second highest plateau on earth, and the Atacama Desert, the driest place in the world. Zool. Stud. 52, 1–13 (2013).Article 

Google Scholar 
80.Herbst, D. B., Conte, F. P. & Brookes, V. J. Osmoregulation in an alkaline salt lake insect, Ephydra (Hydropyrus) hians Say (Diptera: Ephydridae) in relation to water chemistry. J. Insect Physiol. 34, 903–909 (1988).CAS 
Article 

Google Scholar 
81.Cycil, L. M. et al. Metagenomic insights into the diversity of halophilic microorganisms indigenous to the Karak Salt Mine, Pakistan. Front. Microbiol. 11, 1567 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
82.Dillon, J. G., Carlin, M., Gutierrez, A., Nguyen, V. & McLain, N. Patterns of microbial diversity along a salinity gradient in the Guerrero Negro solar saltern, Baja CA Sur, Mexico. Front. Microbiol. https://doi.org/10.3389/fmicb.2013.00399 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
83.Benlloch, S. et al. Prokaryotic genetic diversity throughout the salinity gradient of a coastal solar saltern. Environ. Microbiol. 4, 349–360 (2002).PubMed 
Article 

Google Scholar 
84.Casamayor, E. O. et al. Changes in Archaeal, bacterial and Eukaryal assemblages along a salinity gradient by comparison of genetic fingerprinting methods in a multipond solar saltern. Environ. Microbiol. 4, 338–348 (2002).PubMed 
Article 

Google Scholar 
85.Gorrasi, S. et al. Spatio-temporal variation of the bacterial communities along a salinity gradient within a thalassohaline environment (Saline di Tarquinia Salterns, Italy). Molecules 26, 1338 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
86.Campbell, B. J. & Kirchman, D. L. Bacterial diversity, community structure and potential growth rates along an estuarine salinity gradient. ISME J. 7, 210–220 (2013).CAS 
PubMed 
Article 

Google Scholar 
87.Gorbushina, A. A. Life on the rocks. Environ. Microbiol. 9, 1613–1631 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
88.Wierzchos, J. et al. Adaptation strategies of endolithic chlorophototrophs to survive the hyperarid and extreme solar radiation environment of the Atacama Desert. Front. Microbiol. 6, 1–17 (2015).Article 

Google Scholar 
89.Cody, R. D. & Cody, A. M. Gypsum nucleation and crystal morphology in analog saline terrestrial environments. J. Sediment. Res. 58, 247–255 (1988).CAS 

Google Scholar 
90.Arp, G., Thiel, V., Reimer, A., Michaelis, W. & Reitner, J. Biofilm exopolymers control microbialite formation at thermal springs discharging into the alkaline Pyramid Lake, Nevada, USA. Sediment. Geol. 126, 159–176 (1999).ADS 
CAS 
Article 

Google Scholar 
91.Dupraz, C. et al. Processes of carbonate precipitation in modern microbial mats. Earth-Sci. Rev. 96, 141–162 (2009).ADS 
CAS 
Article 

Google Scholar 
92.Cabestrero, Ó. & Sanz-Montero, M. E. Brine evolution in two inland evaporative environments: Influence of microbial mats in mineral precipitation. J. Paleolimnol. 59, 139–157 (2016).Article 

Google Scholar 
93.Farías, M. E. Microbial Ecosystems in Central Andes Extreme Environments (Springer, 2020).Book 

Google Scholar 
94.Oren, A. Thermodynamic limits to microbial life at high salt concentrations. Environ. Microbiol. 13, 1908–1923 (2011).CAS 
PubMed 
Article 

Google Scholar 
95.Rice, E. W., Baird, R. B., Eaton, A. D. & Clesceri, L. S. Standard Methods for the Examination of Water and Wastewater, 22nd Edition. (APHA American Public Health Association, 2012).96.Díaz, C. & Maidana, N. I. Diatomeas de los Salares Atacama y Punta Negra II Región-Chile (Centro de Ecología Aplicada Ltda. & Minera Escondida Ltda, 2005).
Google Scholar 
97.Patrick, R. Results of Research in the Antofagasta Ranges of Chile and Bolivia. II. Diatoms (Bacillariophyceae) from the Alimentary tract of Phoenicoparrus jamesi (1961).98.Frenguelli, J. Diatomeas del Río de la Plata. Rev. del Mus. la Plata Sección Bot. 3, 213–334 (1941).
Google Scholar 
99.Parra, O., González, M., Dellarossa, V., Rivera, P. & Orellana, M. Taxonomic Manual of Phytoplankton of Continental Waters with Special Reference to the Phytoplankton of Chile: Chlorophyceae. Part III: Cryptophyceae, Dinophyceae, Euglenophyceae (1982).100.Parra, O. & González, M. Taxonomic Manual of Phytoplankton of Continental Waters with Special Reference to the Phytoplankton of Chile: Chlorophyceae. Part I: Volvocales, Tetrasporales, Chlorococcales and Ulothricales (1983).101.Seeligmann, C. & Maidana, N. I. Diatomeas (Bacillariophyceae) en ambientes de altura de la provincia de Catamarca (Argentina). Boletín Soc. Argentina Bot. 38, 39–50 (2003).
Google Scholar 
102.Seeligmann, C., Maidana, N. I. & Morales, M. Diatoms (Bacillariophyceae) of high altitude wetlands in the Province of Jujuy-Argentina. Boletín Soc. Argentina Bot. 43, 1–17 (2008).
Google Scholar 
103.Maidana, N. I. & Seeligmann, C. Diatomeas (Bacillariophyceae) de ambientes acuáticos de altura de la Provincia de Catamarca, Argentina II. Boletín la Soc. Argentina Bot. 41, 1–13 (2006).
Google Scholar 
104.Álvarez Blanco, I., Cejudo Figueiras, C., Godos, I. F., Múñoz Torre, R. & White Lance, S. The diatoms of the salt flats of the Bolivian Altiplano: Floristic singularities. Bull. R. Span. Soc. Nat. Hist. 105, 67–82 (2011).
Google Scholar 
105.Maidana, N. I. & Seeligmann, C. T. Diatoms (Bacillariophyceae) in high-altitude wetlands of Catamarca Province (Argentina). III. Bol. LA Soc. Argentina Bot. 50, 447–466 (2015).Article 

Google Scholar 
106.Woelfl, S., Caputo, L., García-Chicote, J. & de Los Ríos, P. Manuales Para la Bioindicación: Zooplancton Vol. 1 (Manuales Sociedad Chilena de Limnología, 2008).
Google Scholar 
107.De los Rios-Escalante, P. & Salgado, I. Artemia (Crustacea, Anostraca) in Chile: A review of basic and applied biology. Lat. Am. J. Aquat. Res. 40, 487–496 (2017).Article 

Google Scholar 
108.Araya, J. M. & Zúñiga, L. R. Taxonomic manual of the lacustrine zooplankton of Chile. Limnol. Bull. Univ. Austral Chile 8, 1–69 (1985).
Google Scholar 
109.Fernández, H. R. & Domínguez, E. Guide for the Determination of South American Benthic Arthropods, Entomotropica16(3), 219 (2001).110.Crespo, J. E. & Baessolo, L. A. Biogeografia y taxonomia del género Artemis (Crustacea, Anostraca) en Chile: una revisión. Hist. Nat. I(1), 17–21 (2002).
Google Scholar 
111.Villalobos, A., Ness, J. E., Gustafsson, C., Minshull, J. & Govindarajan, S. Gene Designer: A synthetic biology tool for constructing artificial DNA segments. BMC Bioinform. 7, 285 (2006).Article 
CAS 

Google Scholar 
112.Palma, A. T., Schwarz, A. O. & Fariña, J. M. Experimental evidence of the tolerance to chlorate of the aquatic macrophyte Egeria densa in a Ramsar wetland in southern Chile. Wetlands 33, 129–140 (2013).Article 

Google Scholar 
113.Echeverría-Vega, A. et al. Watershed-induced limnological and microbial status in two oligotrophic Andean Lakes exposed to the same climatic scenario. Front. Microbiol. 9, 357 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
114.Lane, D. J. 16S/23S rRNA sequencing. In Nucleic Acid Techniques in Bacterial Systematic (eds Stackebrandt, E. & Goodfellow, M.) 115–175 (Wiley, 1991).
Google Scholar 
115.Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
116.Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C. & Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
117.Clarke, K. R. & Gorley, R. N. Getting Started with PRIMER v7 20 (Plymouth Marine Laboratory, 2015).
Google Scholar 
118.Langille, M. G. I. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821 (2014).Article 
CAS 

Google Scholar 
119.Clesceri, L. S., Greenberg, A. E. & Eaton, A. D. Standard Methods for the Examination of Water and Wastewater, 20th Edition. (APHA American Public Health Association, 1998). More

1.Alfaro, M. E. et al. Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. PNAS 106, 134–14 (2009).Article 

Google Scholar 
2.Seehausen, O. Process and pattern in cichlid radiations – inferences for understanding unusually high rates of evolutionary diversification. N. Phytol. 207, 304–312 (2015).Article 

Google Scholar 
3.Tank, D. C. et al. Nested radiations and the pulse of angiosperm diversification: Increased diversification rates often follow whole genome duplications. N. Phytol. 207, 454–467 (2015).Article 

Google Scholar 
4.Neubauer, T. A., Harzhauser, M., Georgopoulou, E., Kroh, A. & Mandic, O. Tectonics, climate, and the rise and demise of continental aquatic species richness hotspots. PNAS 112, 11478–11483 (2015).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
5.Schluter, D. The Ecology of Adaptive Radiation (Oxford University Press, 2000).6.Schenk, J. J., Rowe, K. C. & Steppan, S. J. Ecological opportunity and incumbency in the diversification of repeated continental colonizations by muroid rodents. Syst. Biol. 62, 837–864 (2013).PubMed 
Article 
PubMed Central 

Google Scholar 
7.Jenkins, N. C., Pimm, S. L. & Joppa, L. N. Global vertebrate diversity and conservation. PNAS 110, E2602–E2610 (2013).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
8.Kier, G. et al. A global assessment of endemism and species richness across island and mainland regions. PNAS 106, 9322–9327 (2009).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
9.Reyjol, Y. et al. Patterns in species richness and endemism of European freshwater fish. Glob. Ecol. Biogeogr. 16, 65–75 (2007).Article 

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

Google Scholar 
11.Albrecht, C., Trajanovski, S., Kuhn, K., Streit, B. & Wilke, T. Rapid evolution of an ancient lake species flock: freshwater limpets (Gastropoda: Ancylidae) in the Balkan lake Ohrid. Org. Divers. Evol. 6, 294–307 (2006).Article 

Google Scholar 
12.Vonlanthen, P. et al. Eutrophication causes speciation reversal in whitefish adaptive radiations. Nature 482, 357–362 (2012).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
13.Renema, W. et al. Hopping hotspots: Global shifts in marine biodiversity. Science 321, 654–657 (2008).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
14.Gargani, J. & Rigollet, C. Mediterranean Sea level variations during the Messinian salinity crisis. Geophys. Res. Lett. 34, 1–5 (2007).Article 

Google Scholar 
15.Culver, D. C. et al. The mid-latitude biodiversity ridge in terrestrial cave fauna. Ecography 29, 120–128 (2006).Article 

Google Scholar 
16.Zagmajster, M. et al. Geographic variation in range size and beta diversity of groundwater crustaceans: Insights from habitats with low thermal seasonality. Glob. Ecol. Biogeogr. 23, 1135–1145 (2014).Article 

Google Scholar 
17.Trontelj, P., Blejec, A. & Fišer, C. Ecomorphological convergence of cave communities. Evolution 66, 3852–3865 (2012).PubMed 
Article 
PubMed Central 

Google Scholar 
18.Trontelj, P., Borko, Š. & Delić, T. Testing the uniqueness of deep terrestrial life. Sci. Rep. 9, 15188 (2019).ADS 
PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
19.Morvan, C. et al. Timetree of Aselloidea reveals species diversification dynamics in groundwater. Syst. Biol. 62, 512–522 (2013).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
20.Eme, D. et al. Do cryptic species matter in macroecology? Sequencing European groundwater crustaceans yields smaller ranges but does not challenge biodiversity determinants. Ecography 41, 424–436 (2017).Article 

Google Scholar 
21.Lukić, M., Delić, T., Pavlek, M., Deharveng, L. & Zagmajster, M. Distribution pattern and radiation of the European subterranean genus Verhoeffiella (Collembola, Entomobryidae). Zool. Scr. 49, 86–100 (2019).Article 

Google Scholar 
22.Väinölä, R. et al. Global diversity of amphipods (Amphipoda; Crustacea) in freshwater. Hydrobiologia 595, 241–255 (2008).Article 

Google Scholar 
23.Fišer, C., Delić, T., Luštik, R., Zagmajster, M. & Altermatt, F. Niches within a niche: ecological differentiation of subterranean amphipods across Europe’s interstitial waters. Ecography 42, 1212–1223 (2019).Article 

Google Scholar 
24.Stroud, J. T. & Losos, J. B. Ecological opportunity and adaptive radiation. Annu. Rev. Ecol. Evol. Syst. 47, 507–532 (2016).Article 

Google Scholar 
25.McInerney, C. E. et al. The ancient Britons: groundwater fauna survived extreme climate change over tens of millions of years across NW Europe. Mol. Ecol. 23, 1153–1166 (2014).PubMed 
Article 
PubMed Central 

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

Google Scholar 
27.Culver, D. C. & Pipan, T. The Biology of Caves and Other Subterranean Habitats 2nd edn (OUP, 2019).28.Kralj-Fišer, S. et al. The interplay between habitat use, morphology and locomotion in subterranean crustaceans of the genus. Niphargus. Zool. 139, 125742 (2020).Article 

Google Scholar 
29.Delić, T., Trontelj, P., Rendoš, M. & Fišer, C. The importance of naming cryptic species and the conservation of endemic subterranean amphipods. Sci. Rep. 7, 1–12 (2017).Article 
CAS 

Google Scholar 
30.Harmon, L. J., Schulte, J. A., Losos, J. B. & Larson, A. Tempo and mode of evolutionary radiation in iguanian lizards. Science 301, 961–964 (2003).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
31.Murrell, D. J. A global envelope test to detect non‐random bursts of trait evolution. Methods Ecol. Evol. 9, 1739–1748 (2018).Article 

Google Scholar 
32.Freckleton, R. P. & Harvey, P. H. Detecting non-Brownian trait evolution in adaptive radiations. PLoS Biol. 4, e373 (2006).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
33.Clavel, J., Escarguel, G. & Merceron, G. mvMORPH: an R package for fitting multivariate evolutionary models to morphometric data. Methods Ecol. Evol. 6, 1311–1319 (2015).Article 

Google Scholar 
34.Kováč, M. et al. The central paratethyspalaeoceanography: a water circulation model based on microfossilproxies, climate, and changes of depositional environment. Acta Geol. Slov. 9, 75–114 (2017).
Google Scholar 
35.Kováč, M. et al. Towards better correlation of the Central Paratethys regional time scale with the standard geological time scale of the Miocene Epoch. Geol. Carpath. 69, 283–300 (2018).ADS 
Article 

Google Scholar 
36.Mahler, D. L., Ingram, T., Revell, L. J. & Losos, J. B. Exceptional convergence on the macroevolutionary landscape in island lizard radiations. Science 341, 6143 (2013).Article 
CAS 

Google Scholar 
37.Ingram, T. & Mahler, D. SURFACE: detecting convergent evolution from comparative data by fitting Ornstein‐Uhlenbeck models with stepwise Akaike Information Criterion. Methods Ecol. Evol. 4, 416–425 (2013).Article 

Google Scholar 
38.Hansen, T. F. Stabilizing selection and the comparative analysis of adaptation. Evolution 51, 1341 (1997).Article 

Google Scholar 
39.Popov, S. V., Rögl, F. & Rozanov, A. Y. Lithological-Paleogeographic Maps of Paratethys: 10 Maps Late Eocene to Pliocene (Schweizerbart’sche Verlagsbuchhandlung, 2004).40.Barrier, E., Vrielynck, B., Brouillet, J. F. & Brunet, M. F. Paleotectonic Reconstruction of the Central Tethyan Realm. Tectonono-Sedimentary-Palinspastic Maps from Late Permian to Pliocene (CCGM/CGMW, 2018).41.Handy, M. R., Ustaszewski, K. & Kissling, E. Reconstructing the Alps–Carpathians–Dinarides as a key to understanding switches in subduction polarity, slab gaps and surface motion. Int. J. Earth. Sci. 104, 1–26 (2015).CAS 
Article 

Google Scholar 
42.Esmaeili-Rineh, S., Sari, A., Delić, T., Moškrič, A. & Fišer, C. Molecular phylogeny of the subterranean genus Niphargus (Crustacea: Amphipoda) in the Middle East: a comparison with European Niphargids. Zool. J. Linn. Soc. 175, 812–826 (2015).Article 

Google Scholar 
43.Pimm, S. L. et al. The biodiversity of species and their rates of extinction, distribution, and protection. Science 344, 6187 (2014).Article 
CAS 

Google Scholar 
44.Hou, Z. & Sket, B. A review of Gammaridae (Crustacea: Amphipoda): the family extent, its evolutionary history, and taxonomic redefinition of genera. Zool. J. Linn. Soc. 176, 323–348 (2016).Article 

Google Scholar 
45.Corrigan, L. J., Horton, T., Fotherby, H., White, T. A. & Hoelzel, A. R. Adaptive evolution of deep‐sea amphipods from the superfamily lysiassanoidea in the North Atlantic. Evol. Biol. 41, 154–165 (2014).Article 

Google Scholar 
46.Clarke, A. & Johnston, I. A. Evolution and adaptive radiation of Antarctic fishes. Trends Ecol. Evol. 11, 212–218 (1996).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
47.Macdonald, K. S. 3rd, Yampolsky, L. & Duffy, J. E. Molecular and morphological evolution of the amphipod radiation of Lake Baikal. Mol. Phylogenet. Evol. 35, 323–343 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
48.Elmer, K. R. et al. Parallel evolution of Nicaraguan crater lake cichlid fishes via non-parallel routes. Nat. Commun. 5, 5168 (2014).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
49.Esquerré, D. & Keogh, J. S. Parallel selective pressures drive convergent diversification of phenotypes in pythons and boas. Ecol. Lett. 19, 800–809 (2016).PubMed 
Article 
PubMed Central 

Google Scholar 
50.von Saltzwedel, H., Scheu, S. & Schaefer, I. Founder events and pre-glacial divergences shape the genetic structure of European Collembola species. BMC Evol. Biol. 16, 148 (2016).Article 
CAS 

Google Scholar 
51.Mammola, S. et al. Scientists’ warning on the conservation of subterranean ecosystems. BioScience 69, 641–650 (2019).Article 

Google Scholar 
52.Lefébure, T., Douady, C. J., Malard, F. & Gibert, J. Testing dispersal and cryptic diversity in a widely distributed groundwater amphipod (Niphargus rhenorhodanensis). Mol. Phylogenet. Evol. 42, 676–686 (2007).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
53.Zagmajster, M., Turjak, M., & Sket, B. Database on subterranean biodiversity of the Dinarides and neighboring regions – SubBioDatabase. In 21st International Conference on Subterranean Biology, 2–7 September, 2012, Košice, Slovakia, Abstract book (ed. Kováč, Ĺ., et al.) 116–117 https://doi.org/10.13140/2.1.4518.0487 (Pavol Jozef Šafárik University, Košice, 2012).54.Clark, K., Karsch-Mizrachi, I., Lipman, D. J., Ostell, J. & Sayers, E. W. GenBank. Nucleic Acids Res. 44, D67–D72 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
55.Fišer, C. et al. Translating Niphargus barcodes from Switzerland into taxonomy with a description of two new species (Amphipoda, Niphargidae). ZooKeys 760, 113–141 (2018).Article 

Google Scholar 
56.Jurado-Rivera, J. A. et al. Molecular systematics of Haploginglymus, a genus of subterranean amphipods endemic to the Iberian Peninsula (Amphipoda: Niphargidae). Contrib. Zool. 86, 239–260 (2017).Article 

Google Scholar 
57.Copilaş-Ciocianu, D., Borko, Š. & Fišer, C. The late blooming amphipods: global change promoted post-Jurassic ecological radiation despite Palaeozoic origin. Mol. Phylogenet. Evol. 143, 106664 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
58.Horton, T. et al. World Register of Marine Species. https://www.marinespecies.org. Accessed 6 Mar 2020 (2020).59.Fišer, C., Trontelj, P., Luštrik, R. & Sket, B. Toward a unified taxonomy of Niphargus (Crustacea: Amphipoda): a review of morphological variability. Zootaxa 2061, 1–22 (2009).Article 

Google Scholar 
60.Katoh, K. & Standley, D. M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
61.Talavera, G. & Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 56, 564–577 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
62.Guindon, S. et al. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst. Biol. 59, 307–321 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
63.Lanfear, R., Calcott, B., Ho, S. Y. & Guindon, S. PartitionFinder: combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol. Biol. Evol. 29, 1695–1701 (2012).64.Ronquist, F. et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 61, 539–542 (2012).PubMed 
PubMed Central 
Article 

Google Scholar 
65.Nguyen, L. T., Schmidt, H. A., von Haeseler, A. & Minh, B. Q. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum likelihood phylogenies. Mol. Biol. Evol. 32, 268–274 (2015).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
66.Rambaut, A., Drummond, A. J., Xie, D., Baele, G. & Suchard, M. A. Tracer v1.7. http://tree.bio.ed.ac.uk/software/tracer/ (2018).67.Minh, B. Q., Nguyen, M. A. T. & von Haeseler, A. Ultrafast approximation for phylogenetic bootstrap. Mol. Biol. Evol. 30, 1188–1195 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
68.Bouckaert, R. et al. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput. Biol. 10, e1003537 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
69.Jażdżewskii, K. & Kupryjanowicz, J. One more fossil Niphargid (Malacostraca: Amphipoda) from Baltic Amber. J. Crustac. Biol. 30, 413–416 (2010).Article 

Google Scholar 
70.Brikiatis, L. Late Mesozic North Atlantic land bridges. Earth-Sci. Rev. 159, 47–57 (2016).ADS 
Article 

Google Scholar 
71.Allegrucci, G., Trucchi, E. & Sbordoni, V. Tempo and mode of species diversification in Dolichopoda cave crickets (Orthoptera, Rhaphidophoridae). Mol. Phylogenet. Evol. 60, 108–121 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
72.R Core Team. R: A Language and Environment for Statistical Computing. https://www.R-project.org/ (R Foundation for Statistical Computing, 2018).73.Pybus, O. G. & Harvey, P. H. Testing macro-evolutionary models using incomplete molecular phylogenies. Proc. R. Soc. B. 267, 2267–2272 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
74.Etienne, R. S. & Haegeman, B. A conceptual and statistical framework for adaptive radiations with a key role for diversity dependence. Am. Nat. 180, E75–E89 (2012).PubMed 
Article 
PubMed Central 

Google Scholar 
75.Bollback, J. P. SIMMAP: Stochastic character mapping of discrete traits on phylogenies. BMC Bioinformatics 7, 88 (2006).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
76.Pagel, M., Meade, A. & Barker, D. Bayesian estimation of ancestral character states on phylogenies. Syst. Biol. 53, 673–684 (2004).PubMed 
Article 
PubMed Central 

Google Scholar 
77.QGIS.org. QGIS Geographic Information System. http://www.qgis.org (QGIS Association, 2021).78.Esri & U.S. National Park Service. “Physical” [basemap]. Scale Not Given. “World Physical Map”. https://www.arcgis.com/home/item.html?id=c4ec722a1cd34cf0a23904aadf8923a0. Accessed 12 Dec 2019 (2019).79.Borko, Š., Trontelj, P., Seehausen, O., Moškrič, A., Fišer, C. Supplementary data and code: a subterranean adaptive radiation of amphipods in Europe [Data set]. Zenodo https://doi.org/10.5281/zenodo.4779097 (2021). More

Ethics and vertebrate animalsThe field surveys and collections were conducted on accessible public areas or private residential areas with property owners’ permission. The study did not involve human participants, or endangered or protected species. Laboratory mice were used as a blood source for mosquitoes. All experimental protocols were approved by the Institutional Animal Care and Use Committee (IACUC) of the University of California, Irvine (UCI) (IACUC protocol number: AUP-19-165). All methods were carried out in accordance with relevant IACUC guidelines and regulations.Study sites and mosquito larval habitat surveillanceThe study was carried out in Orange County, California, USA. Orange County is a highly urbanized county with an estimated population density of approximately 1470 people/km2 according to U.S. Census Bureau, an average annual low/high temperature range of 13–25 °C, 65% relative humidity, and annual precipitation of about 350 mm according to U.S. Climate Data. Annual rainfall was 261 mm, 311 mm, 198 mm and 475 mm for 2016, 2017, 2018 and 2019, respectively. A major drought event occurred in December 2017 and February 2018 when the total rainfall in the 3-month period was 20.6% of the 30-year average. Both Ae. aegypti and Ae. albopictus were discovered in the county in 20158. Culex quinquefasciatus is the most abundant mosquito in the county and breeds readily in a variety of residential, commercial and USDS water sources, and is the primary vector of West Nile virus in southern California18.Larval mosquito surveillance in Orange County was conducted from 2016 to 2019 by the Orange County Mosquito and Vector Control District (OCMVCD) through its routine mosquito surveillance and treatment program, following the recommendations of the California Department of Public Health and the Mosquito and Vector Control Association of California19. Briefly, OCMVCD staff conducted routine inspection for aquatic habitats in randomly selected public areas, and performed door-to-door mosquito larval and adult sampling on residential or commercial premises upon the request of the residents or business owners while distributing public education materials for vector control and personal protection. Arial photography was used to examine the presence of abandoned swimming pools in residential areas. In addition to surface aquatic habitats, subsurface habitats (e.g., catch basins, underground drains, manhole chambers, and public utility vaults) were examined for larval abundance of all mosquito species. In 2019, OCMVCD completed 5,622 mosquito service requests, and conducted 11,813 inspection and treatments on routine sites using a variety of public health-approved adulticides and larvicides. A total of 38,099 underground drains and catch basins and 6925 km of flood channels were treated. In addition, a total of 17,783 km of gutters and 3562 neglected swimming pools were inspected and treated. The larval distribution data reported here were based on this extensive field sampling effort20.Larval sampling used standard mosquito dippers or pipettes, and specialized modifications of these to sample hard to reach areas. Mosquito larvae from each source were collected, transferred into a uniquely-numbered vial with isopropyl alcohol (70%), and submitted to the laboratory for identification; if present, live pupae were collected and held in site-specific labelled rearing chambers (BioQuip Products, Inc., Rancho Dominguez, CA) until emergence. Third and fourth instar mosquito larvae (1–100, depending on sample size) and emerged adults were identified to species using a stereo microscope (40–50x) and morphological features described in taxonomic keys21,22. Results were uploaded to OCMVCD’s data management system, along with collection date, GPS location, and habitat type for each sample site. For this study, larval habitats were classified into six types: small container, underground system, ornamental water features, marsh, pools/spas, and creek (Table S1). The container classification included flowerpots/vases, saucers, tires, bowls, boxes, buckets, dishes, tree holes, etc. Underground storm drain system referred to larval habitats such as catch basins, manhole chambers, underground drains, and public utility vaults that were below the ground. Water feature included flood control channels, ponds, fountains, birdbaths, street gutters and small reservoirs, etc. Marsh included both fresh and salt water marshes.Mosquito strains and water source for laboratory studiesWe examined the effect of USDS water on oviposition substrate preference and larval development in microcosms in an insectary with climate control (27 ± 1 °C, 70 ± 10% relative humidity, and 12 h light/12 h dark photoperiod) at UCI. To minimize potential bias on behavior and ecology from mosquito colonization, this study did not use previously established laboratory mosquito colonies. Instead, we used Ae. aegypti and Ae. albopictus adults reared from field-collected eggs using ovicups in residential areas of Orange and Los Angeles Counties, California, respectively. Culex quinquefasciatus were also reared from eggs of field-collected, blood-engorged adult mosquitoes using gravid traps in Orange County23.All experiments reported here used two types of habitat water: (1) USDS water collected from seven manhole chambers or catch basins (33°47′01.9″N, 117°53′19.0″W, Orange City, manhole; 33°52′25.0″N, 117°57′02.6″W, Fullerton City, manhole; 33°44′44.4″N, 118°06′24.2″W, Seal Beach City, manhole; 33°55′38.9″N, 117°56′51.4″W, La Habra City, manhole; 33°52′48.9″N, 117°55′21.4″W, Fullerton City, catch basin; 33°54′35.2″N, 117°56′02.5″W, Fullerton City, catch basin; 33°52′25.0″N, 117°57′02.6″W, Fullerton City, catch basin); and 2) flowerpot water from vases of three cemeteries in Orange County (33°50′29.0″N, 117°53′57.9″W; 33°46′21.5″N, 117°50′35.8″W; 33°46′12.3″N, 117°50′21.4″W). Water (including sediments) from each breeding source was collected with mosquito dippers and mixed together by habitat type into 18.9 L (five-gallon) Nalgene™ containers. The containers were transported to the laboratory in shaded ice containers, and stored overnight in a refrigerator at 4 °C. The experiments described below were conducted on the field-collected water for the two habitat types. We selected flowerpot water as the comparison substrate because flowerpot containers showed the highest larval positivity rate in the study area.Oviposition preference testTo examine whether USDS water attracts or repels egg laying by Ae. aegypti and Ae. albopictus mosquitoes, a two-choice oviposition preference test was conducted. Briefly, this experiment used two ovicups placed within a mosquito cage (1 × 0.5 × 0.5 m3), one ovicup with 200 ml USDS water and another with 200 ml flowerpot water. Adult mosquitoes were bloodfed on mice; fully engorged females 3-days post-bloodfeeding were used for oviposition preference tests. Ten gravid Ae. aegypti females were released into a cage and allowed to lay eggs for three days, and the number of eggs in each ovicup were counted. Five replicates were used. The same experiment was conducted for Ae. albopictus.To evaluate whether the presence of Cx. quinquefasciatus larvae has any impact on the egg laying behavior of invasive Aedes mosquitoes, the two-choice oviposition preference test described above was used. One ovicup contained 200 ml USDS water and ten first-instar Cx. quinquefasciatus larvae, while the second ovicup contained 200 ml USDS water only. Ten gravid Ae. aegypti or Ae. albopictus females were released into a cage and allowed to lay eggs for three days. Five replicates were used. We also conducted this experiment using flowerpot water with the same design and same number of replicates to determine whether the impact of Cx. quinquefasciatus larvae on Aedes mosquito egg laying behavior was similar across different water substrate types.Egg hatchingTo investigate the effects of different habitat water sources on egg hatching, 50 Ae. aegypti or Ae. albopictus eggs on separate filter papers were introduced into ovicups with 200 ml USDS water or flowerpot water. Deoxygenized distilled water that we routinely use in laboratory mosquito colony maintenance was used as a positive control. The experiment was conducted in an insectary with climate control (27 ± 1 °C). The number of larvae hatched were counted daily for six days continuously. Five replicates were used.Larval survivorshipA life table study was conducted on Ae. aegypti and Ae. albopictus larvae to determine the effect of USDS water and flowerpot water on larval development and survivorship. Twenty-five newly hatched Ae. aegypti or Ae. albopictus larvae were introduced into a microcosm that contained 200 ml USDS or field flowerpot water. The number of dead and surviving larvae was recorded daily until they pupated. Pupae were counted, and removed to different paper cups for emergence to adults. Four replicates were used for each type of habitat water per species. We included Cx. quinquefasciatus in the larval life table study for method validation purposes because the larvae of this species were known to successfully develop into pupae and adults in USDS water in southern California10.Larval survivorship experiments were conducted in two different seasons. The first was in the summer (August–September) 2019 when the density of invasive Aedes species peaked19, and also insecticide runoff from mosquito and residential/agricultural pest control applications were at the highest levels in southern California24. The second was in the winter (December) 2019 when there was little insecticide treatment for mosquito and pest control. This design enabled us to examine seasonality in larval survivorship and the impact of environmental insecticide runoff in USDS water. To determine whether USDS water’s nutritional deficiency plays a major role in limiting Aedes larval development, we repeated the larval survival experiment by adding 0.1 g Tetramin Tropical Flakes, the standard larval mosquito diet in insectaries, to the microcosms every 2 days. The number of dead and surviving larvae, pupae, and emergent adults was recorded daily.Data analysisAll aquatic habitats that were positive or negative for the larvae of Ae. aegypti, Ae. albopictus and Cx. quinquefasciatus (the predominant species), were mapped using ArcGIS 10.7.1. The proportion of aquatic habitats positive for Ae. aegypti and Cx. quinquefasciatus was calculated for each habitat type from 2016 to 2019. To examine variation in Aedes and Culex larval positivity rate among different groups of larval habitats within the USDS, larval positivity rates for Ae. aegypti and Cx. quinquefasciatus were calculated for underground water retention vaults, underground catch basins/manholes, and underground pipelines/tunnels. The Chi-square test was used to examine the statistical significance. Culex quinquefasciatus was analyzed because it was the most common species, whereas Ae. albopictus was not included in the analysis due to insufficient number of Ae. albopictus positive habitats. To determine whether USDS water attracted or repelled oviposition of invasive Aedes mosquitoes, a pairwise t test was used to compare egg number in USDS water ovicups to flowerpot water ovicups for each Aedes species. Similarly, a pairwise t-test was used to test the effect of Cx. quinquefasciatus larvae on Aedes mosquito oviposition choice.To examine the effect of water sources on egg hatching, the t-test was used to analyze the egg hatching rate. The analysis of larval life table study data focused on pupation rates and larval-to-pupal development times. The pupation rate was calculated as the proportion of first-instar larvae that molted into pupae. The effect of water sources and larval food supplementation on pupation rate was analyzed using non-parametric Wilcoxon test. The t-test was used to analyze the duration of larval-to-pupal development. Kaplan–Meier survival analysis was used to determine the effects of food supplementation and water source on larval development for each species, and the log-rank test was conducted to determine their statistical significance. All statistical analyses were performed using JMP software (JMP 14.2, SAS Institute Inc.). More

1.Bezemer, T. M. & van Dam, N. M. Linking aboveground and belowground interactions via induced plant defenses. Trends Ecol. Evol. 20, 617–624. https://doi.org/10.1016/j.tree.2005.08.006 (2005).Article 
PubMed 

Google Scholar 
2.Kaplan, I. et al. Physiological integration of roots and shoots in plant defense strategies links above-and belowground herbivory. Ecol. Lett. 11, 841–851. https://doi.org/10.1111/j.1461-0248.2008.01200.x (2008).Article 
PubMed 

Google Scholar 
3.Huang, W., Siemann, E., Carrillo, J. & Ding, J. Below-ground herbivory limits induction of extrafloral nectar by above-ground herbivores. Ann. Bot. 115, 841–846. https://doi.org/10.1093/aob/mcv011 (2015).Article 
PubMed 
PubMed Central 

Google Scholar 
4.Omer, A. D., Thaler, J. S., Granett, J. & Karban, R. Jasmonic acid induced resistance in grapevines to a root and leaf feeder. J. Econ. Entomol. 93, 840–845. https://doi.org/10.1603/0022-0493-93.3.840 (2000).CAS 
Article 
PubMed 

Google Scholar 
5.Yamawo, A., Ohsaki, H. & Cahill, J. F. Jr. Damage to leaf veins suppresses root foraging precision. Am. J. Bot. 106, 1126–1130. https://doi.org/10.1002/ajb2.1338 (2019).Article 
PubMed 

Google Scholar 
6.Heil, M. & Karban, R. Explaining evolution of plant communication by airborne signals. Trends Ecol. Evol. 25, 137–144. https://doi.org/10.1016/j.tree.2009.09.010 (2010).Article 
PubMed 

Google Scholar 
7.Karban, R., Yang, L. J. & Edwards, K. F. Volatile communication between plants that affects herbivory: A meta-analysis. Ecol. Lett. 17, 44–52. https://doi.org/10.1111/ele.12205 (2014).Article 
PubMed 

Google Scholar 
8.Yoneya, K. & Takabayashi, J. Plant-plant communication mediated by airborne signals: Ecological and plant physiological perspectives. Plant Biotechnol. 31, 409–416. https://doi.org/10.5511/plantbiotechnology.14.0827a (2014).CAS 
Article 

Google Scholar 
9.Morrell, K. & Kessler, A. Plant communication in a widespread goldenrod: Keeping herbivores on the move. Funct. Ecol. 31, 1049–1061. https://doi.org/10.1111/1365-2435.12793 (2017).Article 

Google Scholar 
10.Arimura, G. et al. Herbivory-induced volatiles elicit defence genes in lima bean leaves. Nature 406, 512–515 (2000).ADS 
CAS 
Article 

Google Scholar 
11.Shiojiri, K. et al. Weeding volatiles reduce leaf and seed damage to field-grown soybeans and increase seed isoflavones. Sci. Rep. 7, 1–8. https://doi.org/10.1038/srep41508 (2017).CAS 
Article 

Google Scholar 
12.Karban, R., Shiojiri, K., Huntzinger, M. & McCall, A. C. Damage-induced resistance in sagebrush: Volatiles are key to intra- and interplant communication. Ecology 87, 922–930. https://doi.org/10.1890/0012-9658(2006)87[922:drisva]2.0.co;2 (2006).Article 
PubMed 

Google Scholar 
13.Kikuta, Y. et al. Specific regulation of pyrethrin biosynthesis in Chrysanthemum cinerariaefolium by a blend of volatiles emitted from artificially damaged conspecific plants. Plant Cell Physiol. 52, 588–596. https://doi.org/10.1093/pcp/pcr017 (2011).CAS 
Article 
PubMed 

Google Scholar 
14.Karban, R., Baldwin, I. T., Baxter, K. J., Laue, G. & Felton, G. W. Communication between plants: Induced resistance in wild tobacco plants following clipping of neighboring sagebrush. Oecologia 125, 66–71. https://doi.org/10.1007/PL00008892 (2000).ADS 
CAS 
Article 
PubMed 

Google Scholar 
15.Karban, R., Huntzinger, M. & McCall, A. C. The specificity of eavesdropping on sagebrush by other plants. Ecology 85, 1845–1852. https://doi.org/10.1890/03-0593 (2004).Article 

Google Scholar 
16.Shiojiri, K. et al. Exposure to artificially damaged goldenrod volatiles increases saponins in seeds of field-grown soybean plants. Phytochem. Lett. 36, 7–10. https://doi.org/10.1016/j.phytol.2020.01.014 (2020).CAS 
Article 

Google Scholar 
17.Glinwood, R., Ninkovic, V. & Pettersson, J. Chemical interaction between undamaged plants—Effects on herbivores and natural enemies. Phytochemistry 72, 1683–1689. https://doi.org/10.1016/j.phytochem.2011.02.010 (2011).CAS 
Article 
PubMed 

Google Scholar 
18.Lokesh, R. A. V. I., Manasvi, V. & Lakshmi, B. P. Antibacterial and antioxidant activity of saponin from Abutilon indicum leaves. Asian J. Pharm. Clin. Res. 9, 344–347. https://doi.org/10.22159/ajpcr.2016.v9s3.15064 (2016).CAS 
Article 

Google Scholar 
19.Raji, P., Samrot, A. V., Keerthana, D. & Karishma, S. Antibacterial activity of alkaloids, flavonoids, saponins and tannins mediated green synthesised silver nanoparticles against Pseudomonas aeruginosa and Bacillus subtilis. J. Cluster Sci. 30, 881–895. https://doi.org/10.1007/s10876-019-01547-2 (2019).CAS 
Article 

Google Scholar 
20.Toth, R., Toth, D. & Starke, D. Vesicular-arbuscular mycorrhizal colonization in Zea mays affected by breeding for resistance to fungal pathogens. Can. J. Bot. 68, 1039–1044. https://doi.org/10.1139/b90-131 (1990).Article 

Google Scholar 
21.Matyssek, R. et al. The plant’s capacity in regulating resource demand. Plant Biol. 7, 560–580. https://doi.org/10.1055/s-2005-872981 (2005).CAS 
Article 
PubMed 

Google Scholar 
22.Brundrett, M. C. Coevolution of roots and mycorrhizas of land plants. New Phytol. 154, 275–304. https://doi.org/10.1046/j.1469-8137.2002.00397.x (2002).Article 
PubMed 

Google Scholar 
23.Kiers, E. T., Rousseau, R. A., West, S. A. & Denison, R. F. Host sanctions and the legume-rhizobium mutualism. Nature 425, 78–81. https://doi.org/10.1038/nature01931 (2003).ADS 
CAS 
Article 

Google Scholar 
24.Walters, D. R. Plant Defense: Warding Off Attack by Pathogens, Herbivores, and Parasitic Plants (Blackwell Publishing, 2011).
Google Scholar 
25.Szakiel, A., Pączkowski, C. & Henry, M. Influence of environmental biotic factors on the content of saponins in plants. Phytochem. Rev. 10, 493–502. https://doi.org/10.1007/s11101-010-9164-2 (2011).CAS 
Article 

Google Scholar 
26.Singh, B. & Kaur, A. Control of insect pests in crop plants and stored food grains using plant saponins: A review. LWT 87, 93–101. https://doi.org/10.1016/j.lwt.2017.08.077 (2018).CAS 
Article 

Google Scholar 
27.Barton, K. E. & Koricheva, J. The ontogeny of plant defense and herbivory: Characterizing general patterns using meta-analysis. Am. Nat. 175, 481–493. https://doi.org/10.1086/650722 (2010).Article 
PubMed 

Google Scholar 
28.Feeny PP. Seasonal changes in oak leaf tannins and nutrients as a cause of spring feeding by winter moth caterpillars. Ecology 51, 565–581 https://doi.org/10.2307/1934037 (1970).Article 

Google Scholar 
29.Dudt, J. F. & Shure, D. J. The influence of light and nutrients on foliar phenolics and insect herbivory. Ecology 75, 86–98. https://doi.org/10.2307/1939385 (1994).Article 

Google Scholar 
30.Julkunen-Tiitto, R. Phenolic constituents in the leaves of northern willows: Methods for the analysis of certain phenolics. Asian J. Pharm. Clin. Res. 33, 213–217. https://doi.org/10.1021/jf00062a013 (1985).CAS 
Article 

Google Scholar 
31.Folin, O. & Denis, W. A colorimetric method for the determination of phenols (and phenol derivatives) in urine. J. Biol. Chem. 22, 305–308 (1915).CAS 
Article 

Google Scholar 
32.Huang, D., Ou, B. & Prior, R. L. The chemistry behind antioxidant capacity assays. J. Agric. Food Chem. 53, 1841–1856. https://doi.org/10.1021/jf030723c (2005).CAS 
Article 
PubMed 

Google Scholar 
33.Mukai, T., Horie, H. & Goto, T. A simple method for determining saponin in tea seed. Chagyo Kenkyu Hokoku 75, 29–31 (1992) (Japanese with English abstract).CAS 
Article 

Google Scholar 
34.Dubois, M., Gilles, K. A., Hamilton, J. K., Rebers, P. A. & Smith, F. A colorimetric method for the determination of sugars. Nature 168, 167–167 (1951).ADS 
CAS 
Article 

Google Scholar 
35.Dubois, M., Gilles, K. A., Hamilton, J. K., Rebers, P. T. & Smith, F. Colorimetric method for determination of sugars and related substances. Anal. Chem. 28, 350–356. https://doi.org/10.1021/ac60111a017 (1956).CAS 
Article 

Google Scholar 
36.Harad, Y. Cation and anion exchange capacity of soil background and methods. 302 Jpn. J. Soil Scie. Plant Nutr. 55, 273–283. 303 (1984) (in Japanese).
37.R Development Core Team. R A language and environment for statistical computing (R Foundation for Statistical Computing, 2020).
Google Scholar 
38.De Geyter, E., Lambert, E., Geelen, D. & Smagghe, G. Novel advances with plant saponins as natural insecticides to control pest insects. Pest Technol. 1, 96–105 (2007).
Google Scholar 
39.Pankhurst, C. E. & Sprent, J. I. Effects of water stress on the respiratory and nitrogen-fixing activity of soybean root nodules. J. Exp. Bot. 26, 287–304. https://doi.org/10.1093/jxb/26.2.287 (1975).CAS 
Article 

Google Scholar 
40.Sparg, S., Light, M. E. & Van Staden, J. Biological activities and distribution of plant saponins. J. Ethnopharmacol. 94, 219–243. https://doi.org/10.1016/j.jep.2004.05.016 (2004).CAS 
Article 
PubMed 

Google Scholar 
41.Saha, S., Walia, S., Kumar, J., Dhingra, S. & Parmar, B. S. Screening for feeding deterrent and insect growth regulatory activity of triterpenic saponins from Diploknema butyracea and Sapindus mukorossi. J. Agric. Food Chem. 58, 434–440. https://doi.org/10.1021/jf902439m (2010).CAS 
Article 
PubMed 

Google Scholar 
42.Hoagland, R. E., Zablotowicz, R. M. & Oleszek, W. A. Effects of alfalfa saponins on in vitro physiological activity of soil and rhizosphere bacteria. J. Crop. Prod. 4, 349–361. https://doi.org/10.1300/J144v04n02_16 (2001).CAS 
Article 

Google Scholar 
43.Killeen, G. F. et al. Antimicrobial saponins of Yucca schidigera and the implications of their in vitro properties for their in vivo impact. J. Agric. Food Chem. 46, 3178–3186. https://doi.org/10.1021/jf970928j (1998).CAS 
Article 

Google Scholar 
44.Sugiyama, A. The soybean rhizosphere: Metabolites, microbes, and beyond—A review. J. Adv. Res. 19, 67–73. https://doi.org/10.1016/j.jare.2019.03.005 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
45.Lucas-Barbosa, D. Integration studies on plant-pollinator and plant–herbivore interactions. Trends Plant Sci. 21, 125–133. https://doi.org/10.1016/j.tplants.2015.10.013 (2016).CAS 
Article 
PubMed 

Google Scholar 
46.De Deyn, G. B., Raaijmakers, C. E. & Van der Putten, W. H. Plant community development is affected by nutrients and soil biota. J. Ecol. 92(5), 824–834. https://doi.org/10.1111/j.0022-0477.2004.00924.x (2004).Article 

Google Scholar 
47.Knelman, J. E. et al. Nutrient addition dramatically accelerates microbial community succession. PLoS ONE 9(7), e102609. https://doi.org/10.1371/journal.pone.0102609 (2014).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
48.Yoneya, K. & Takeshi, M. Co-evolution of foraging behaviour in herbivores and their natural enemies predicts multifunctionality of herbivore-induced plant volatiles. Funct. Ecol. 29, 451–461. https://doi.org/10.1111/1365-2435.12398 (2015).Article 

Google Scholar 
49.Karban, R. Plant communication increases heterogeneity in plant phenotypes and herbivore movement. Funct. Ecol. 31, 990–991. https://doi.org/10.1111/1365-2435.12806 (2017).Article 

Google Scholar  More

1.May, R. M. Thresholds and breakpoints in ecosystems with a multiplicity of stable states. Nature 269, 471–477 (1977).Article 

Google Scholar 
2.Scheffer, M., Carpenter, S., Foley, J., Folke, C. & Walker, B. Catastrophic shifts in ecosystems. Nature 413, 591–596 (2001).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
3.Archer, S., Schimel, D. S. & Holland, E. A. Mechanisms of shrubland expansion: land use, climate, or carbon dioxide. Clim. Change 29, 91–99 (1995).Article 

Google Scholar 
4.Maher, E. L. & Germino, M. J. Microsite variation among conifer species during seedling establishment at alpine treeline. Ecoscience 13, 334–341 (2006).Article 

Google Scholar 
5.Knapp, A. K. et al. Shrub encroachment in North American grasslands: shifts in growth form dominance alters control of ecosystem carbon inputs. Glob. Change Biol. 14, 615–623 (2008).Article 

Google Scholar 
6.McKee, K. L. & Rooth, J. E. Where temperate meets tropical: multi-factorial effects of elevated CO2, nitrogen enrichment, and competition on a mangrove-salt marsh community. Glob. Change Biol. 14, 971–984 (2008).Article 

Google Scholar 
7.Huang, H., Zinnert, J. C., Wood, L. K., Young, D. R. & D’Odorico, P. Non-linear shift from grassland to shrubland in temperate barrier islands. Ecology 99, 1671–1681 (2018).PubMed 
Article 
PubMed Central 

Google Scholar 
8.Huang, H., Anderegg, L. D. L., Dawson, T. E., Mote, S. & D’Odorico, P. Critical transition to woody plant dominance through microclimate feedbacks in North American coastal ecosystems. Ecology 101, e03107 (2020).PubMed 
Article 
PubMed Central 

Google Scholar 
9.Huenneke, L. F., Anderson, J. P., Remmenga, M. & Schlesinger, W. H. Desertification alters patterns of aboveground net primary production in Chihuahuan ecosystems. Glob. Change Biol. 8, 247–264 (2002).Article 

Google Scholar 
10.Li, J., Okin, G. S., Hartman, L. J. & Epstein, H. E. Quantitative assessment of wind erosion and soil nutrient loss in desert grasslands of southern New Mexico, USA. Biogeochemistry 85, 317–332 (2007).Article 

Google Scholar 
11.D’Odorico, P., Okin, G. S. & Bestelmeyer, B. T. A synthetic review of feedbacks and drivers of shrub encroachment in arid grasslands. Ecohydrology 5, 520–530 (2012).Article 

Google Scholar 
12.Van Auken, O. Shrub invasions of North American semiarid grasslands. Annu. Rev. Ecol. Evol. Syst. 31, 197–215 (2000).Article 

Google Scholar 
13.Sankaran, M. et al. Determinants of woody cover in African savannas. Nature 438, 846–849 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
14.Gehrig-Fasel, J., Guisan, A. & Zimmermann, N. E. Tree line shifts in the Swiss Alps: climate change or land abandonment? J. Veg. Sci. 18, 571–582 (2007).Article 

Google Scholar 
15.Cavanaugh, K. C. et al. Poleward expansion of mangroves is a threshold response to decreased frequency of extreme cold events. Proc. Natl Acad. Sci. USA 111, 723–727 (2014).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
16.D’Odorico, P. et al. Vegetation–microclimate feedbacks in woodland–grassland ecotones. Glob. Ecol. Biogeogr. 22, 364–379 (2013).Article 

Google Scholar 
17.He, Y., D’Odorico, P. & De Wekker, S. F. The relative importance of climate change and shrub encroachment on nocturnal warming in the Southwestern United States. Int. J. Climatol. 35, 475–480 (2014).Article 

Google Scholar 
18.He, Y., D’Odorico, P., De Wekker, S. F., Fuentes, J. D. & Litvak, M. On the impact of shrub encroachment on microclimate conditions in the northern Chihuahuan desert. J. Geophys. Res. Atmos. 115, D21120 (2010).Article 

Google Scholar 
19.Ives, A. R. & Carpenter, S. R. Stability and diversity of ecosystems. Science 317, 58–62 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
20.Scheffer, M. et al. Early-warning signals for critical transitions. Nature 461, 53–59 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
21.Ratajczak, Z. et al. Abrupt change in ecological systems: inference and diagnosis. Trends Ecol. Evol. 33, 513–526 (2018).PubMed 
Article 
PubMed Central 

Google Scholar 
22.Sugihara, G. & May, R. M. Applications of fractals in ecology. Trends Ecol. Evol. 5, 79–86 (1990).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
23.Rodriguez-Iturbe, I. & Rinaldo, A. Fractal River Basins: Chance and Self- Organization. (Cambridge University Press, New York, 1997).24.Majumder, S., Tamma, K., Ramaswamy, S. & Guttal, V. Inferring critical thresholds of ecosystem transitions from spatial data. Ecology 100, e02722 (2019).PubMed 
Article 

Google Scholar 
25.Staver, A. C., Asner, G. P., Rodriguez-Iturbe, I., Levin, S. A. & Smit, I. Spatial patterning among savanna trees in high resolution, spatially extensive data. Proc. Natl Acad. Sci. USA 116, 10685 (2019).Article 
CAS 

Google Scholar 
26.Kéfi, S. et al. Spatial vegetation patterns and imminent desertification in Mediterranean arid ecosystems. Nature 449, 213–217 (2007).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
27.Kéfi, S. et al. Robust scaling in ecosystems and the meltdown of patch size distributions before extinction. Ecol. Lett. 14, 29–35 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
28.Weissmann, H., Kent, R., Michael, Y. & Shnerb, N. M. Empirical analysis of vegetation dynamics and the possibility of a catastrophic desertification transition. PLoS ONE 12, e0189058 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
29.Sankaran, S., Majumder, S., Viswanathan, A. & Guttal, V. Clustering and correlations: inferring resilience from spatial patterns in ecosystems. Methods Ecol. Evol. 10, 2079–2089 (2019).Article 

Google Scholar 
30.Zinnert, J. C. et al. Spatial–temporal dynamics in barrier island upland vegetation: the overlooked coastal landscape. Ecosystems 19, 685–697 (2016).CAS 
Article 

Google Scholar 
31.Thompson, J. A., Zinnert, J. C. & Young, D. R. Immediate effects of microclimate modification enhance native shrub encroachment. Ecosphere 8, e01687 (2017).Article 

Google Scholar 
32.Wood, L. K., Hays, S. & Zinnert, J. C. Decreased temperature variance associated with biotic composition enhances coastal shrub encroachment. Sci. Rep. 10, 8210 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
33.Stauffer, D. & Aharony, A. Introduction to Percolation Theory (Taylor and Francis, London, 1985).34.Taubert, F. et al. Global patterns of tropical forest fragmentation. Nature 554, 519–522 (2018).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
35.Scanlon, T. M., Caylor, K. K., Levin, S. A. & Rodriguez-Iturbe, I. Positive feedbacks promote power-law clustering of Kalahari vegetation. Nature 449, 209–212 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
36.Kéfi, S., Rietkerk, M., van Baalen, M. & Loreau, M. Local facilitation, bistability and transitions in arid ecosystems. Theor. Popul. Biol. 71, 367–379 (2007).PubMed 
Article 
PubMed Central 

Google Scholar 
37.Berdugo, M., Kéfi, S., Soliveres, S. & Maestre, F. T. Plant spatial patterns identify alternative ecosystem multifunctionality states in global drylands. Nat. Ecol. Evol. 1, 3 (2017).PubMed 
Article 
PubMed Central 

Google Scholar 
38.Dakos, V., Kéfi, S., Rietkerk, M., van Nes, E. H. & Scheffer, M. Slowing down in spatially patterned ecosystems at the brink of collapse. Am. Nat. 177, E153–E166 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
39.Weerman, E. J. et al. Changes in diatom patch-size distribution and degradation in a spatially self-organized intertidal mudflat ecosystem. Ecology 93, 608–618 (2012).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
40.Kéfi, S., Dakos, V., Scheffer, M., Van Nes, E. H. & Rietkerk, M. Early warning signals also precede non-catastrophic transitions. Oikos 122, 641–648 (2012).Article 

Google Scholar 
41.van Belzen, J. et al. Vegetation recovery in tidal marshes reveals critical slowing down under increased inundation. Nat. Commun. 8, 15811 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
42.USACE-TEC & JALBTCX. Hyperspectral imagery for Hog Island, VA, 2013 ver 7. Environmental Data Initiative. https://doi.org/10.6073/pasta/6a5cc305e93c2baf9283facee688c504 (2018).43.Young, D. R. et al. Cross-scale patterns in shrub thicket dynamics in the Virginia barrier complex. Ecosystems 10, 854–863 (2007).Article 

Google Scholar 
44.Young, D. R., Shao, G. & Porter, J. H. Spatial and temporal growth dynamics of barrier island shrub thickets. Am. J. Bot. 82, 638–645 (1995).Article 

Google Scholar 
45.McGarigal, K., Cushman, S. A. & Ene, E. FRAGSTATS v4: Spatial Pattern Analysis Program for Categorical and Continuous Maps (University of Massachusetts at Amherst, MA, 2012).46.Clauset, A., Shalizi, C. R. & Newman, M. E. J. Power-law distributions in empirical data. SIAM Rev. 51, 661–703 (2009). 2009.Article 

Google Scholar 
47.Hayden, B. P. Ecosystem feedbacks on climate at the landscape scale. Philos. Trans. R. Soc. B, Biol. Sci. 353, 5–18 (1998).Article 

Google Scholar  More

We untangled how the direct and indirect paths of forest cover and water quality variables interact and shape anopheline assemblages in two seasons. Although previous studies determined how environmental variables at different spatial extents affected anopheline distributions in Amazônia, most studies focused on a single effect of an environmental variable or focused on single habitat types (terrestrial or aquatic)22,23,41,42. Our most important finding is that seasonality modulates the direct and indirect effects of forest cover on Amazônian anopheline larval distributions. In particular, we found that forest cover had stronger direct and indirect influence on larval anopheline assemblage composition in the rainy season than the dry season.The different paths and strengths of forest cover influences on anopheline assemblages during the rainy and dry seasons can be associated with the responses of adults and larvae to forest characteristics. Forest cover influences water quality variables of ponds by shading, organic matter inputs and erosion processes43. These effects have consequences for pond water quality44 and favor the establishment of different culicid species45. We showed that during the rainy season, forest cover directly and indirectly influenced site water quality. Greater forest cover in the rainy season directly and indirectly affected A. nimbus and the secondary malaria vectors A. triannulatus and A. braziliensis positively. In the dry season, greater forest cover positively but marginally affected A. peryassui, A. nuneztovari and A. albitarsis, but only indirectly through water quality. Some species like A. triannulatus, A. nuneztovari and A. braziliensis coexist with the malaria vector, A. darlingi, in breeding sites46, and these species have been positively associated with pH, dissolved oxygen and total suspended solids in natural and artificial habitats20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47, which are environmental conditions favored by greater forest cover. The marginal indirect effect of forest cover on anopheline assemblage in the dry season suggests that we need caution in the interpretation of this result and long-term temporal data is required to confirm if this effect is corroborated.Forest conditions influence mosquito vectors and their hosts. For example, some mosquitoes are zoophiles that feed on the blood of birds, reptiles, and mammals48, which are often more abundant in conserved areas. Other species are anthropophilic and prefer to feed on human blood49 and altered environments can force these species to migrate and, consequently, to change hosts48. In our study, A. triannulatus and A. minbus were more abundant in sites with more natural characteristics, whereas A. darlingi and A. nuneztovari were more abundant in altered landscapes. In addition, urbanization and deforestation increase the proximities of humans and domestic animals to mosquito vectors and their hosts, thereby maintaining and increasing transmission cycles50.Forest conditions influence anopheline diversity by different paths, which may alter the strength of their seasonal effects. During the dry season, mosquito survival is also affected by altered microclimate (e.g., lower humidity)51 and lentic habitats contain less water, increased nutrient concentrations and decreased abundance and richness of mosquitoes52,53. We observed that rainfall plays an important role in the larval abundance of Anopheles in artificial larval habitats in Manaus. In addition, climatic factors such as rainfall and river levels are strongly associated with vector abundance and malaria cases in the region54,55. During the rainy season, increased water volume in artificial habitats provides more areas for distribution and development of mosquito species56 and we detected a significant increase in abundance of A. triannulatus, A. darlingi and A. nuneztovari. These observations may partially explain why we found a direct effect of forest cover on mosquitoes only during the rainy season.Our results add more evidence that managing and conserving forest cover is important to control anophelines, thereby decreasing the contact of potential vectors (e.g., A. darlingi) with humans. In general, our results support the idea that mosquitoes are directly affected by the loss of native forest cover57 in the rainy season. Mosquitoes associated with serious human diseases (e.g., malaria, yellow fever, dengue, leishmaniasis) are more abundant in areas with low levels of native forest cover14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55,56,57,58. This is a critically important finding because recent studies have shown that forest cover plays an important role in the vector dynamics of mosquitoes and forest conservation keeps pathogens within the forest, avoiding spillover to human settlements59. On the other hand, deforestation provides favorable conditions for these vectors, thereby increasing malária cases and decreasing scores of the Human Development Index60. In addition, there is a positive correlation between mosquito abundance in fragmented forests and the prevalence of Plasmodium, the protozoan that causes malaria61.Artificial larval habitats promote conditions for malaria vectors in Amazônia62,63. Therefore, the best way to develop control techniques would be to understand larval ecology in these habitats, where they are more sensitive to infections by pathogens, parasites, predation, larvicides and growth regulators64. This information is necessary to minimize failures in programs to control or eradicate the vector and the disease. Under this perspective, our study adds a new piece in the puzzle of mosquito control in Amazônia. For example, during the rainy season when forest cover directly and indirectly influences larval habitats, control programs can strengthen the control of key limnological variables, habitat structure, and entomological aspects, intensifying the environmental filter, particularly in areas with little forest cover and greater human concentrations near those habitats. The limnological study of Anopheles larval habitats is still far from complete, as each case has peculiarities inherent to them. Despite attempts, anophelines demonstrate versatility in relation to abiotic parameters20,21,22,23,24,25,26,65,66. However, we can use approaches that modify the larval environments. For example, more efficient management of water levels in fish farming ponds could decrease larval numbers and anopheline reproduction, Similarly, greater rationing of fish feed would decrease the supply of food resources for mosquito larvae. It is also worth mentioning that some variables are related to the efficiency of others. Regarding biological control via entomopathogenic bacteria, environmental factors (solar radiation) and water quality (amounts of total suspended solids and organic matter), can interfere with the effectiveness of the formulated Bacillus sphaericus applied in habitats for vector control62,63,64,65,66,67. Furthermore, eutrophication decreased the assemblages of aquatic invertebrates predating mosquito larvae.Another alternative is the use of physical control (removal of grasses and macrophytes from the edge of habitats), helping to reduce microhabitats that provide larval refuges. Also, increased light and water temperature at the edges favor natural predation and biological control processes from potential fish and macroinvertebrates. The conservation of natural enemies and the use of biotic agents in the population control of vector mosquitoes have been recommended in small and medium-sized natural and artificial breeding sites19,20,21,22,23,24,25,26,27,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53. A combination of techniques that shape the important environmental variables for the establishment of these species are essential for vector control.The analytical approach used here opens some windows of opportunity for improvements that are important to be recognized. First, our model did not incorporate important complexity of natural systems, such as ecological interactions among vectors and hosts, including human behavior. Agent-based models, including different host behavior, could provide important insights in this way. Second, our study is very limited in terms of temporal climatic variability. Additional information is needed to better understand the effects of long-term changes in land-use, water quality and climate and their interactions with mosquito assemblages in the region, particularly considering an ecological-evolutionary perspective. Third, it is important to highlight that the magnitude of effects of the estimated drivers were not the same in the rainy and dry seasons. Also, they may not remain constant in coming decades, especially considering potential regional process on mosquito assemblages, such as spillover effects, mass effects and host changes. Fourth, our study was carried out in an area of Amazonia that has experienced, a relatively old land use conversion from forest to urban areas (urban expansion rate of around 12% per year for the past 34 years)68. Beginning in the 1970s, human population increased at a rate of around 23% per decade and 25% in Manaus11. Therefore, the region we studied is very relevant in terms of historical interactions among human populations, mosquitoes and land use changes. However, understanding the effect of these changes on mosquito assemblages in areas with different land-use change dynamics, provides us with important information69, particularly those with very rapid urbanization processes, such as in the Arch of Deforestation70. Lastly, we need studies that consider the nexus among climate and land use changes, human and animal population health, economic conditions, and ecosystem services provided by these forest-urban transitional regions. Such information would facilitate including mosquito information in land use planning and climate mitigation programs based on forest management in and around cities.Therefore, identifying ecological factors and paths that affect the composition of species of epidemiological importance are essential because they inform vector integrated management strategies. We emphasize that larval control in lentic habitats requires knowledge about larval ecology and the effects of biotic and abiotic variables on larvae, especially when it comes to biological controls. The application of integrated pest management can be conducted in both dry and rainy seasons. However, we recommend focusing on the dry season when larval habitats are more limited, in smaller volumes and more accessible for entry and application of vector control techniques. These are critically important considerations because over 2 million people live in Amazonas state11 and anophelines transmitted over 59,637 malaria cases in the Amazon region in the first half of 2020, and about 44.4% came from the state of Amazonas71. More

1.Aitken, S. N., Yeaman, S., Holliday, J. A., Wang, T. & Curtis-McLane, S. Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol. Appl. 1, 95–111 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
2.Renner, S. S. & Zohner, C. M. Climate change and phenological mismatch in trophic interactions among plants, insects, and vertebrates. Annu. Rev. Ecol. Evol. Syst. 49, 165–182 (2018).Article 

Google Scholar 
3.Piao, S. et al. Plant phenology and global climate change: current progresses and challenges. Glob. Change Biol. 25, 1922–1940 (2019).Article 

Google Scholar 
4.Kooyers, N. J., Greenlee, A. B., Coloicchio, J. M., Oh, M. & Blackman, B. K. Replicate altitudinal clines reveal that evolutionary flexibility underlies adaptation to drought stress in annual Mimulus guttatus. New Phytol. 206, 152–165 (2015).5.Evans, L. M. et al. Population genomics of Populus trichocarpa identifies signature of selection and adaptive trait associations. Nat. Genet. 46, 1089–1096 (2016).Article 
CAS 

Google Scholar 
6.Wadgymar, S. M., Daws, S. C. & Anderson, J. T. Integrating viability and fecundity selection to illuminate the adaptive nature of genetic clines. Evol. Lett. 1, 26–39 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
7.Nord, E. A. & Lynch, J. P. Plant phenology: a critical controller of soil resource acquisition. J. Exp. Bot. 60, 1927–1937 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
8.Polgar, C. A. & Primack, R. B. Tansley review: leaf-out phenology of temperate woody plants: from trees to ecosystems. New Phytol. 191, 926–941 (2011).PubMed 
Article 
PubMed Central 

Google Scholar 
9.Parmesan, C. & Yohe, G. A globally coherent fingerprint of climate change impacts across natural systems. Nature 421, 37–42 (2003).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
10.Root, T. L. et al. Fingerprints of global warming on wild animals and plants. Nature 421, 57–60 (2003).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
11.Stephens, P. A. et al. Consistent response of bird population to climate change on two continents. Science 352, 84–87 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
12.Wagner, M. R. et al. Natural soil microbes alter flowering phenology and the intensity of selection on flowering time in a wild Arabidopsis relative. Ecol. Lett. 17, 717–726 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
13.Lu, T. et al. Rhizosphere microorganisms can influence the timing of plant flowering. Microbiome 6, 231 (2018).14.Friesen, M. et al. Microbially mediated plant functional traits. Ann. Rev. Ecol. Evol. Syst. 42, 23–46 (2011).Article 

Google Scholar 
15.Mendes, R., Garbeva, P. & Raaijmakers, J. M. The rhizosphere microbiome: significance of plant beneficial, plant pathogenic, and human pathogenic microorganisms. FEMS Microbiol. Rev. 37, 634–663 (2013).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
16.Compant, S., Samad, A., Faist, H. & Sessitsch, A. A review on the plant microbiome: Ecology, functions, and emerging trends in microbial application. J. Adv. Res. 19, 29–37 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
17.Panke-Buisse, K., Poole, A., Goodrich, J., Ley, R. & Kao-Kniffin, J. Selection on soil microbiomes reveals reproducible impacts on plant function. ISME J. 9, 980–989 (2015).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
18.Fitzpatrick, C. R., Mustafa, Z. & Viliunas, J. Soil microbes alter plant fitness under competition and drought. J. Evol. Biol. https://doi.org/10.1111/jeb.13426 (2019).19.Lau, J. A. & Lennon, J. T. Rapid responses of soil microorganisms improve plant fitness in novel environments. Proc. Natl Acad. Sci. USA 109, 14058–14062 (2012).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
20.Van Nuland, M. E., Ware, I. M., Bailey, J. K. & Schweitzer, J. A. Ecosystem feedbacks contribute to geographic variation in the plant-soil evolutionary dynamics across fertility gradient. Funct. Ecol. 33, 95–106 (2019).21.Zolla, G., Badri, D. V., Bakker, M. G., Manter, D. K. & Vivanco, J. M. Soil microbiome vary in their ability to confer drought tolerance to Arabidopsis. Appl. Soil Ecol. 68, 1–9 (2013).Article 

Google Scholar 
22.Gehring, C. A., Sthultz, C. M., Flores-Renteria, L., Whipple, A. V. & Whitham, T. G. Tree genetics defines fungal partner communities that may confer drought tolerance. Proc. Natl Acad. Sci. USA 114, 11169–11174 (2017).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
23.Woolbright, S. A., Whitham, T. G., Gehring, C. A., Allan, G. J. & Bailey, J. K. Climate relicts and their associated communities as natural ecology and evolution laboratories. Trends Ecol. Evol. 29, 406–416 (2014).PubMed 
Article 
PubMed Central 

Google Scholar 
24.Lankau, R. A., Zhu, K. & Ordonez, A. Mycorrhizal strategies of tree species correlate with trailing range edge responses to current and past climate change. Ecology 96, 1451–1458 (2015).Article 

Google Scholar 
25.Ware, I. M. et al. Climate-driven reduction of genetic variation in plant phenology alters soil communities and nutrient pools. Glob. Change Biol. 25, 1514–1528 (2019).Article 

Google Scholar 
26.Nguyen, N. H. et al. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fung. Ecol. 20, 241–248 (2016).Article 

Google Scholar 
27.Nottingham, A. T. et al. Microbes follow Humboldt: temperature drives plant and soil Microbial diversity patterns from the Amazon to the Andes. Ecology 99, 2455–2466 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
28.Martiny, J. B. et al. Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–111 (2006).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
29.Fierer, N., Strickland, M. S., Liptzin, D., Bradford, M. A. & Cleveland, C. C. Global patterns of belowground communities. Ecol. Lett. 12, 1238–1249 (2009).PubMed 
Article 
PubMed Central 

Google Scholar 
30.Fierer, N. et al. Comparative metagenomic, phylogenetic, and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 6, 1007–1017 (2012).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
31.Waldrop, M. P. et al. The interacting roles of climate, soils, and plant production, on soil microbial communities at a continental scale. Ecology 98, 1957–1967 (2017).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
32.Nelson, M. B., Martiny, A. C. & Martiny, J. B. Global biogeography of microbial nitrogen-cycling traits in soil. Proc. Natl Acad. Sci. USA 113, 8033–8040 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
33.Schweitzer, J. A. et al. Plant-soil-microorganism interactions: heritable relationship between plant genotype and associated soil microorganisms. Ecology 89, 773–781 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
34.de Vries, F. T. et al. Abiotic drivers and plant traits explain landscape-scale patterns in soil microbial communities. Ecol. Lett. 15, 1230–1239 (2012).PubMed 
Article 
PubMed Central 

Google Scholar 
35.Prober, S. M. et al. Plant diversity predicts beta but not alpha diversity of soil microbes across grasslands worldwide. Ecol. Lett. 18, 85–95 (2015).PubMed 
Article 
PubMed Central 

Google Scholar 
36.Delgado-Baquerizo, M. et al. Plant attributes explain the distribution of soil microbial communities in two contrasting regions of the globe. New Phytol. 219, 574–587 (2018).PubMed 
Article 
PubMed Central 

Google Scholar 
37.Menzel, A. Trends in phenological phases in Europe between 1951 and 1996. Int. J. Biometeorol. 44, 76–81 (2000).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
38.Morin, X. et al. Leaf phenology in 22 North American tree species during the 21st century. Glob. Change Biol. 15, 961–975 (2010).Article 

Google Scholar 
39.Loarie, S. R. et al. The velocity of climate change. Nature 462, 1052–1055 (2009).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
40.Renwick, K. M. & Rocca, M. E. Temporal context affects the observed rate of climate-driven range shifts in tree species. Glob. Ecol. Biogeog. 24, 44–51 (2015).Article 

Google Scholar 
41.Elena, S. F. & Lenski, R. E. Evolution experiments with microorganisms: the dynamics and genetic bases of adaptation. Nat. Rev. Genet. 4, 457–469 (2003).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
42.Finlay, B. J. Global dispersal of free-living microbial eukaryote species. Science 296, 1061–1063 (2002).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
43.Kivlin, S. N., Emery, S. M. & Rudgers, J. A. Fungal symbionts alter plant responses to global change. Am. J. Bot. 100, 1445–1457 (2013).PubMed 
Article 
PubMed Central 

Google Scholar 
44.Fernandez, C. W. & Kennedy, P. G. Revisiting the ‘Gadgil effect’: do interguild fungal interactions control carbon cycling in forest soils. New Phytol. 209, 1382–1394 (2016).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
45.Fisher, D. G. et al. Plant genetic effects on soils under climate change. Plant Soil 379, 1–19 (2014).46.van der Wal, A., Geyden, T. D., Kuyper, T. W. & de Boer, W. A thready affair: linking fungal diversity and community dynamics to terrestrial decomposition processes. FEMS Microbiol. Rev. 37, 477–494 (2013).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
47.Perez-Izquierdo, L. et al. Plant intraspecific variation modulates nutrient cycling through its belowground rhizospheric microbiome. J. Ecol. 107, 1594–1605 (2019).Article 

Google Scholar 
48.Crowther, T. W. et al. The global soil community and its influence on biogeochemistry. Science 365, https://doi.org/10.1126/science.aav0550 (2019).49.Steidinger, B. S. et al. Climatic controls of decomposition drive the global biogeography of forest-tree symbioses. Nature 569, 404–408 (2019).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
50.van der Putten, W. H., Bradford, M. A., Brinkman, E. P., van de Voorde, T. F. J. & Veen, G. F. Where, when and how plant-soil feedback matters in a changing world. Funct. Ecol. 30, 1109–1121 (2016).Article 

Google Scholar 
51.Van Nuland, M. E. et al. Plant-soil feedbacks: connecting ecosystem ecology and evolution. Funct. Ecol. 30, 1032–1042 (2016).Article 

Google Scholar 
52.Ware, I. M. et al. Feedbacks link ecosystem ecology and evolution across spatial and temporal scales: empirical evidence and future directions. Funct. Ecol. 33, 31–42 (2019).Article 

Google Scholar 
53.Cooke, J. E. K. & Rood, S. B. Trees of the people: the growing science of poplars in Canada and worldwide. Can. J. Bot. 85, 1103–1110 (2007).Article 

Google Scholar 
54.Evans, L. M., Allan, G. J., Meneses, N., Max, T. L. & Whitham, T. G. Herbivore host-associated genetic differentiation depends on the scale of plant genetic variation examined. Evol. Ecol. 27, 65–81 (2013).Article 

Google Scholar 
55.Hijmans, R. J., Cameron, S. E., Para, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978 (2005).Article 

Google Scholar 
56.Castro, H. F., Classen, A. T., Austin, E. E., Norby, R. J. & Schadt, C. W. Soil microbial community responses to multiple experimental climate change drivers. Appl. Environ. Microbiol. 76, 999–1007 (2010).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
57.Wilson, R. M. et al. Stability of peatland carbon to rising temperatures. Nat. Commun. 7, 13723 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
58.Van Nuland, M. E., Bailey, J. B. & Schweitzer, J. A. Divergent plant-soil feedbacks could alter future elevation ranges and ecosystem dynamics. Nature Ecol. Evol. https://doi.org/10.1038/s41559-017-0150 (2017).59.Richardson, A. D. et al. Influence of spring phenology on seasonal and annual carbon balance in two contrasting New England forests. Tree Phys. 29, 321–331 (2009).CAS 
Article 

Google Scholar 
60.Richardson, A. D. et al. Influence of spring and autumn phenological transitions on forest ecosystem productivity. Philos. Trans. Roy. Soc. B Biol. Sci. 365, 3227–3246 (2010).Article 

Google Scholar 
61.Callahan, B. J. et al. DADA2: High resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
62.Tremblay, J. et al. Primer and platform effects on 16S rRNA tag sequencing. Front. Microbiol. 6, 771 (2015).PubMed 
PubMed Central 

Google Scholar 
63.Chao, A., Chiu, C. H. & Jost, L. Unifying species diversity, phylogenetic diversity, functional diversity, and related similarity and differentiation measures through hill numbers. Ann. Rev. Ecol. Evol. Syst. 45, 297–324 (2014).Article 

Google Scholar 
64.Chao, A. et al. Rarefaction and extrapolation with Hill numbers: a framework for sampling and estimation in species diversity studies. Ecol. Monogr. 84, 45–67 (2014).Article 

Google Scholar 
65.Ma, Z. Measuring microbiome diversity and similarity with Hill numbers. Metagenomics https://doi.org/10.1016/B978-0-08-102268-9.00008-2 (2018).66.Ferrier, S., Manion, G., Elith, J. & Richardson, K. Using generalized dissimilarity modelling to analyse and predict patterns of beta diversity in regional biodiversity assessment. Diversity Distrib. 13, 252–264 (2007).Article 

Google Scholar 
67.Fitzpatrick, M. C. et al. Environmental and historical imprints on beta diversity: insights from variation in rates of species turnover along gradients. Proc. Biol. Sci. 280, 20131201 (2013).PubMed 
PubMed Central 

Google Scholar 
68.R Core Team. R: A Language and Environment for Statistical Computing http://www.R-project.org (R Foundation for Statistical Computing, 2016). More