A Multidisciplinary Evaluation of Spatial Heterogeneity in the Puquios of the Salar de Llamara, Atacama Desert, Northern Chile

The Atacama Desert, Central Andes of Northern Chile, is an extreme environment characterized by high UV radiation, wide temperature variation, minimum precipitation and is reputed as the driest desert in the world. Scarce lagoons associated with salt ats (salars) in this desert are the surface expression of shallow groundwater which serve as refugia for life, and often host microbial communities associated with evaporitic mineral deposition. Recent investigations of the Puquios of the Salar de Llamara in the Atacama Desert based on multidisciplinary eld studies provide unprecedented detail regarding the spatial heterogeneity of physical, chemical, and biological characteristics of such saline lake environments. Four main lagoons (‘Puquios’) and more than 400 smaller ponds, occur in an area less than 5 km 2 , are characterized by high variability in electrical conductivity, benthic and planktonic biota and microbiota, lagoon bottom type, and style of mineral deposition. The heterogeneity of system parameters as observed spatially in the Puquios is likely to be expanded with temporal observations incorporating seasonality. Results provide new insight into the complexity of these Andean ecosystems, which may be key to resilience in extreme environments at the edge of habitability. seconds and 72° C for 1 minute, after which nal elongation step at 72°C for 5 minutes. PCR products were visualized in 2% agarose gel to determine the success of amplication and the relative intensity of the bands. Multiple equivalent samples of PCR products were pooled together in equal proportions based on their molecular weight and DNA concentrations. Pooled samples were puried using calibrated Ampure XP beadsto prepare a DNA library by following the Illumina TruSeq DNA library preparation protocol. Sequencing was performed at MR DNA (www.mrdnalab.com, Shallowater, TX, USA) on a MiSeq, following the manufacturer’s guidelines. The sequence data were processed using a proprietary analysis pipeline (MR DNA) with the following steps: the barcodes were eliminated from sequences, then sequences < 150 bp and those with ambiguous base calls were removed. Then, sequences were denoised, OTUs generated and chimeras removed. The reads were added to the Sequence Read Archive (SRA) National Center for Biotechnology Information (NCBI), (http://trace.ncbi.nlm.nih.gov/Traces/sra/). (Puquios of the Salar de Llamara Raw sequence read); BioSample: (A multidisciplinary evaluation of spatial heterogeneity in the Puquios of the Salar de Llamara, Atacama Desert, Northern Chile); SRA: SRR12536357 Desert). The sequences were processed further by means of the QIIME (Quantitative Insights Into Microbial Ecology) pipeline 106 : Raw sequences were ltered by base quality score, average base content read distribution in the reads.


Introduction
Life on Earth can thrive in almost every ecological niche, and many extremophiles prosper in the harshest of environments [1][2][3][4][5][6] . Examples of extreme environments, where biologic communities persist at the edge of habitability, are common throughout some Andean ecosystems of South America 3,7,8 . Studies of more than 40 microbial ecosystems in Chile, Argentina, and Bolivia document a wide range of geological, chemical, biological and mineralogical processes occurring in these extreme environments [8][9][10][11][12][13][14][15][16][17][18][19][20][21][22] . Previous studies of these microbial ecosystems are predominantly focused on the deposition of carbonate minerals 12-14, 19,23 , while the other studies have focused on the deposition of evaporitic minerals including gypsum, halite, and anhydrite, which dominate many of these systems 15,18,19,24−30 . One such example of a gypsum-depositing system is the series of lagoons locally known as Puquios, in the Salar de Llamara of the Atacama Desert. Recent studies of the Puquios suggest that multidisciplinary studies of the strong environmental gradients and a variety of microbial communities 17,24,31,32 , observed in this polyextreme environment will provide new insight into conditions that promote habitability at the dry limits of life.
The Salar de Llamara is located in one of the main endorheic basins in the Pampa del Tamarugal, Central Depression, of the Norte Grande de Chile [33][34][35] . The Puquios are variably sized depressions lled with ground brines. The Pampa del Tamarugal is the largest and last base level before the Paci c Ocean for the sedimentary ll and water input provided by the High and Preandean Cordilleras. The sedimentary ll corresponds to detrital sediments and volcanic rocks that can reach up to 1300 m thickness 36 . The Puquios in this arid to hyperarid environment of the Atacama Desert, are in uenced by a complex interplay of physical, chemical, geological, and biological processes. Known environmental pressures include high uxes of UV irradiation, extreme aridity, and signi cant uctuations of temperature and salinity 27,37,38 that contribute to the characterization as the dry limit of life 39 , of which the Puquios are a modern example. Despite these harsh environmental conditions, the Puquios support the development of a variety of pelagic and benthic communities, including microbial mats 17,19,24,29−31 , and exhibit diverse styles of gypsum deposition.
Previous studies in the Puquios have focused primarily on microscopy and molecular analyses of microbial communities in subaerial crusts and domes 15,18,20,24,27,29,30 but detailed investigations of the relationship between physical, chemical, geological, ecological, and microbiological processes that may in uence mineral deposition in these lagoons are limited. Initial studies in the Puquios of brines 32,40 have begun to quantify strati cation and lateral gradients of brine chemistry within a portion of the Puquios system in the Salar de Llamara. In this paper, we expand on this foundation by incorporating two additional areas within the Puquios system that have yet to be described. We incorporate a multidisciplinary approach to document the spatial heterogeneity of water chemistry, planktonic and benthic biota, community structure of the microbial ecosystems, lagoon bottom types, and mineral phases and morphologies in the Puquios. We base our observations on more than 1,285 in situ measurements of the brines, detailed ecological and chemical descriptions of each lagoon, as well as microbiological and morphological analysis of representative bottom types. The variability in spatial distributions observed throughout the system provide unprecedented insight into the complexity of this intriguing gypsum depositional ecosystem and highlight the need for spatial and temporal multi-disciplinary studies to understand ecosystem function and mineral deposition in these extreme environments. Results also lead to speculation that heterogeneity may be integral to the resilience of ecosystems that thrive at the edge of habitability.

Environmental Setting
The Puquios (Fig. 1) are located in the eastern part of the Salar de Llamara (21°23´S -69°37´W), an actively depositing salar (salt at) in the Pampa del Tamarugal, Central Depression 24 . This salar is located in the distal part of the Arcas Alluvial Fan that provides most of the basin input 41 . The hyperaridity of the Central Depression has fostered the deposition of evaporitic sediments since at least the Middle Miocene 35 , and is controlled by the geologic and oceanographic setting 35,42−44 . Rainfall is extremely low, estimated to be less than 1 mm in the Atacama Desert core 45 as a result of the rainshadows cast by both the Coastal Cordillera and the High Cordillera. Furthermore, evaporation rates signi cantly exceed rainfall, with up to 90% of rainfall in the Western Cordillera lost via evaporation and runoff 46 . The permanent surfacing of sub-terrain waters is recharged by a combination of summer monsoons (Invierno Altiplánico or Invierno Boliviano) in the High Cordillera 47-49 that continuously drain through the alluvial fan systems punctuating the eastern margin of the basin 35,50,51 with occurrences of coastal dripping fog 52,53 . The combination of these hydrological processes associated with the alluvial fan systems and the geomorphological characteristics at the boundary between the Coastal Cordillera and the Pampa del Tamarugal allows the shallow ground water levels to create the small lagoons called puquios with their biological communities. Brines in the Puquios are managed by SQM Ltd.; environmental protection measures were enacted in 2012 to preserve brine chemistry and water levels in the salar environment 32,40,54,55 .
The Puquios system is characterized by four main lagoons (Figs. 1b, 2), plus a series of smaller peripheral ponds (ranging from sub-meter to tens of meters in diameter). Puquio 1 consists of a main shallow lagoon with low-turbidity waters, an area of 1760 m 2 and a maximum depth of 50 cm. A multitude of small ponds exist between the main lagoons of Puquio 1 and 2 that are collectively referred to as "the Transition Zone". The Transition Zone contains nearly 400 small lagoons of various water turbidities, sizes, and depths (maximum depth recorded 80 cm). Puquio 2 is largest lagoon in the Salar de Llamara with very low-turbidity brines, an area of 4650 m 2 and a depth less than 1 meter. The main lagoon of the Puquio 3 system has an area of 1660 m 2 , elongated north -south, surrounded by a number of smaller, peripheral ponds. Water in Puquio 3 is clear, whereas the shallow peripheral ponds often have milky-white to peach-colored water. Puquio 4 consists of one main lagoon (1485 m 2 ) elongated northsouth, surrounded by numerous peripheral ponds, mostly to the east. The main lagoon of Puquio 4 appears turquoise in color, and the water is clear, whereas the surrounding ponds also often have milkywhite to peach-colored water similar to those adjacent to Puquio 3.

Electrical Conductivity
In the austral summer 2017-18, more than 1285 in situ measurements of electrical conductivity (EC, Fig. 2, Supplemental Table S1) were collected from both the surface brines and depths of the lagoons in the Puquios to characterize environmental gradients, and is interpreted to be a proxy for the salinity of the lagoon water. Chemical gradients on multiple scales were revealed; in situ measurements document spatial heterogeneity not only within the Puquios system as a whole, but also de ne a high degree of variability within a single lagoon.
Within the system as a whole, EC exhibited a high degree of spatial heterogeneity in the Puquios (Table 1, Table S1). The lowest EC value was observed in the surface brine of a peripheral lagoon of Puquio 3 (6.2 mS/cm, Figs. 2c, 2d), and the highest EC values were observed in the surface brine of the main lagoon of Puquio 4 (164.9 mS/cm, Figs. 2e, 2f). Spatial trends in EC were also observed within the Puquios 1 and 2 system, where EC values grade from lower values in the main lagoon of Puquio 1 to signi cantly higher EC values in the main lagoon of Puquio 2 in the surface brine (Table 1, Figs. 2a, 2b). Observation of such variability in the Puquios system suggests that hydrological connections between each of the four Puquios may not be direct, and that higher EC values re ect higher degrees of evaporation and/or lower rates of groundwater ux into the brine pools.
Spatial heterogeneity was also found to be a signi cant characteristic of surface brines in the Puquios, even within a single lagoon. EC values within a single lagoon were found to vary through space, with ranges of more than 89.5, 89.6, and 138.0 mS/cm in the surface brines of the main lagoons in Puquios 2, 3 and 4 respectively (Figs. 2a, 2e). Comparison of surface brine EC measurements between the peripheral lagoons also demonstrated a large range in measured values across short distances, suggesting that the peripheral ponds can be isolated from the main lagoon water and may develop distinct EC values and degrees of spatial heterogeneity as a result. EC values of surface brines in the Transition Zone between Puquios 1 and 2 exhibit a range of more than 130.9 mS/cm (Fig. 2a).
Gradients in EC measurements were also a common characteristic of bottom brines in the Puquios. Like the surface brine measurements, the highest bottom brine EC value (161.4 mS/cm) was observed in Puquio 2, while the lowest EC measurement (21.5 mS/cm) was recorded in the bottom brines of Puquio 1 (Fig. 2b). Spatial heterogeneity of EC values in the bottom brines was reduced relative to the range observed in surface brines (Figs. 2a, 2b).  Table 1. Notably, strati cation also exhibited spatial heterogeneity between lagoons: Puquios 1, 2, and 4 exhibited normal strati cation, whereas measurements in Puquio 3 were nearly uniform, suggesting a well-mixed water column (Table 1) at the time of measurement.

Free-Living Biota
The free-living biota that thrives in the aquatic component of the Puquios system is comprised of four major groups, namely: phytobenthos (microalgae that exist on the bottom), phytoplankton (microalgae in the water column), zoobenthos (invertebrates found on the bottom or buried within the rst few centimeters of substrate) and zooplankton (invertebrates in the water column). The rst two groups, which are photosynthetic primary producers, are represented mainly by Diatoms (51 species) and by the less diverse Cyanobacteria (7 species). The primary and secondary consumers present on the bottom and in the water column are represented by 8 taxa: arachnids (spiders), annelids (worms), branchiopods (brine shrimps), coleopterans (beetles), copepods, dipterans ( ies), hexapods (springtails), gastropods (snails) and odonates (dragon ies). Several of these taxa have complex life cycles and thus occur in different ontogenetic stages (i.e. from larval stages to adults). Some, like dipterans and odonates, only spend their larval phase in the water while some, like the brine shrimp Artemia franciscana, are holoplanktonic.
In the austral summer 2017-18, all four assemblages of biota exhibited ample variability in several ecological descriptors (e.g. number of taxa and Shannon diversity index), both within and between Puquios. The abundance of many species within each assemblage exhibited notable differences between Puquios (Fig. 3). The abundance of phytobenthos ranged from a few cells/mm 2 to ca. 700 cells/mm 2 and was highest in Puquio 1, which also exhibited important differences among sampling points, followed by Puquio 2. Phytobenthos was least abundant in Puquios 3 and 4 (Fig. 3a). Phytoplankton exhibited an even more pronounced pattern of abundance variation between and within Puquios. The highest abundance of phytoplankton was recorded in the southern-most sampling point of Puquio 1, more than twice of that found on its northern-most point. Puquios 2, 3 and 4 shared a uniform lower abundance in the range of 1,000 cells/L (Fig. 3b). The abundance of zoobenthos was also greater, and highly variable in Puquio 1, with up to 90,000 individuals/m 2 in its northern section and around 20,000 individuals/m 2 along its southern edge. Puquio 2, followed by 3, displayed the lowest abundance, and Puquio 4 exhibited a varied and intermediate level of abundance (Fig. 3c). Zooplankton displayed a consistently low abundance in Puquios 1, 2 and 3, in contrast to Puquio 4, which was characterized by a relatively higher abundance of zooplankton, particularly in its southern-most section (Fig. 3d).
Besides abundance, other community-level descriptors, such as number of taxa (S) and Shannon diversity index (H'), quanti ed for each of these four assemblages, varied amply ( Table 2)

Bacterial Community Analysis
The microbial communities were assessed in March 2017 from subaerial gypsum structures (Fig. 4a) and brines. Bacterial community composition and predicted metabolisms reveal a heterogenous distribution of assemblages that is largely driven by local environmental factors. The microbial mat communities inhabiting the subaerial gypsum structures around the margins of the Puquios exhibit a stacking pattern at each sampling location, with a gradient of colors changing from green (E1) to brown (E2, E3) to colorless (E4) with increasing depth (Fig. 4b).
Samples of the E1 bacterial communities from twelve gypsum structures and four brines were analyzed by the 16S rRNA gene bacterial amplicon sequencing. The bacterial diversity measured by H' ranged from 2.5 to 5.8 in the E1 layer and 2.6-3.9 in the brines. The evenness index was less than 0.82 for the E1 layer and less than 0.61 for brine samples. The similarity between the E1 samples was > 40% at the phylum level and > 25% at the genus level. Bacteria in the E1 layer and brines a liate mostly to high-rank prokaryotic taxa, Proteobacteria, Planctomycetes, Verrucomicrobia, Cyanobacteria, and Actinobacteria in important abundances. Within the Proteobacteria phylum, the alpha (orders Rhodobacterales, Rhizobiales) and gamma (orders Alteromonadales, Chromatiales) classes were predominant. The similarity among the brine samples was > 55% at the phylum and > 30% and the genus level. Bacteria in these groups a liate mostly to the high-rank prokaryotic taxa Proteobacteria and Bacteroidetes as previously reported (Demergasso et al., 2007). Classes alpha (order Sphingomonadales) and gamma (orders Oceanospirillales and Vibrionales) were predominant within the Proteobacteria phylum.
The differences between the bacterial community pro les across the Puquios is driven largely by brine chemistry. The non-Metric Multi-dimensional scaling (nMDS) analysis grouped the Puquio 1 E1 samples separately from the other three Puquios, displaying a strong correlation to EC, pH, and Mg concentrations in the E1 layer (Pearson correlation > 0.8) (Fig. 4c). The grouping of the brine bacterial communities was also related to brine chemistry (EC, SO 4 2− and Mg 2+ ) (Pearson correlation > 0.9), grouping Puquio 1 away from the other higher salinity locations (Fig. 4d). E1 samples from Puquio 1 contained higher abundances (Pearson correlation > 0.85) of sulfate reducing bacteria (Desulfobacterales, Thermoanerobacterales) and some anoxygenic phototrophic and sulfur oxidizing Alphaproteobacteria (Rodobiaceae). Puquios 2, 3, and 4 were characterized by higher EC brine values and higher abundances of the halophilic and halotolerant genera of Cyanobacteria, Synechocystis, Halomicronema, Cyanotece and Euhalothece. The variability in the microbial communities of the brine samples was related to several genera of anoxygenic phototrophic bacteria from Alpha and Gammaproteobacteria and Chloro exi groups.
The metabolisms from the predicted metagenomes by PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) suggests that there are not only key metabolic differences between E1 and brines, but also between E1 layers of structures related to the salinity characteristics at each sample location. The occurrence of the sulfur oxidation pathway appears to be dominant in both brines and subaerial structures (Fig. 4e). Phototrophy is more abundant in E1 samples, whereas heterotrophic and aerobic metabolisms are more important in the brine samples (Fig. 4e). There is also a clear differentiation between the subaerial structures. Anoxygenic and oxygenic photosynthesis are more abundant in E1 from the most saline ponds. Meanwhile, nitrate and sulfate reduction and the use of CO and H 2 as electron donors are more represented in E1 from the lowest saline pond.
Desulfobacterales was observed among the sulfate reducers which has been reported previously in the subaqueous structures 29  suggesting microbial in uence on mineral precipitation in this lagoon (Table 3).
Bottom types in the Transition Zone were highly variable, grading from occulent sediment to a variety of orange, brown, and black microbial mats, with various morphologies. SEM analyses of critically dried samples of mat types throughout the Transition Zone (Figs. 6a-6h) showed intimate mixtures of EPS, microbes and diatoms with ne-grained minerals (primarily gypsum), often with distinct lamination of the mineral substrate. Gypsum crystals within the microbial mats exhibited various habits including bipyramidal, acicular, platy, lenticular, etc., and small crystals often aggregated into larger crystals within organic matrices (Figs. 6c, 6d, 6g). Accumulations of minerals and microbes were present in most samples examined. Another striking characteristic of samples collected from ponds in the transitional zone was the common appearance of small acicular crystals accruing along lamentous cyanobacterial sheaths (Fig. 6h). Additional mineral accumulations found in abundance throughout the Transition Zone included various carbonate species, sodium sulfates, magnesium clays, and manganese oxides. As in Puquio 1, mineral deposition in the Transition Zone was largely microbially in uenced (Table 3).
In the Transition Zone to the southwest of Puquio 2 and within Puquio 2 the bottom type was characterized by euhedral gypsum spar (Supplementary Figs. S2aa, S2ab, Fig. 6a). Samples displayed cm-scale selenite crystals and lacked organic matrices. These large crystals had well-developed crystal faces and were sometimes coated with manganese oxides (Fig. 6m). The bottom of Puquio 2 was blanketed by gypsum spar crystals and lacked dense bacterial communities. The side walls of Puquio 2 were comprised mainly of gypsum botryoids, which often merged into plates. Diatoms were prevalent, often in great abundance. The large crystals, lack of agglutinated grains, and paucity of observed microbes and EPS suggests that microbial in uence was limited, and mineral precipitation was largely physiochemical ( Table 3).
The main lagoon of Puquio 3 was dominated by pinnacles (Fig. 5c, Supplemental Figs. S2v, S2w, and S3c) along the northwestern margin of Puquio 3, and orange/brown microbial mats in the southeast (Fig. 5c, Supplemental Figs. S2t, S3c). SEM/EDS analyses revealed granular and laminated gypsum accumulations closely associated with microbes (Supplemental Table S2, Figs. 6q-6t). A sample collected in the south exhibited larger crystals and aggregates, whereas a sample collected in the north of Puquio 3 (Fig. 6r) was largely microcrystalline. Mineral precipitation appears both microbially and physicochemically mediated ( Table 3).

Discussion
Results above allow us to characterize the Puquios in terms of heterogeneity. Heterogeneity in turn may foster ecosystem resilience, enabling habitability within the polyextreme conditions that characterize the hyper arid environment of the Atacama Desert (sensu 57 ).
Heterogeneity of biota, microbial communities, lagoon bottom types and type of mineral precipitation within the Puquios as documented above, appears to be driven by spatial variability in brine chemistry. Variations in water level and brine chemistry may drive changes in microbial communities as their abundances and functions respond to changing environmental conditions through time. Functional redundancy in the system is hypothesized to contribute to a high degree of resilience in response to changes in brine chemistry and lifestyle allowing these microbial communities to persist in harsh polyextreme conditions. Where EC is relatively low, as in Puquio 1 and the Transition Zone, microbial mats are common (Supplemental Fig. S4a -S4c) and the microscopic analysis revealed a heterogeneous distribution of lamentous and unicellular cyanobacteria, cell morphologies compatible with chromatium and diatoms (Supplemental Fig. S5). Abundances of phytobenthos and zoobenthos are relatively high, and mineral deposition, primarily gypsum, is largely microbially driven. In contrast, when brines have high EC values, such as those observed in Puquio 2 (Fig. 2c, F, Supplemental Fig. S1) the occurrence of microbial mats is rare and the bottom is covered by gypsum crystals, water column biota are present in lower proportions, and gypsum precipitation appears to be largely physicochemical (Supplemental Fig.   S4d). Although more work is needed to delineate physicochemical versus microbial controls on gypsum precipitation, mineral precipitation in Puquio 3 appears in uenced by both microbial and physicochemical controls, whereas Puquio 4 appears dominated by physicochemical controls. In both sites microbial mats are covered by different gypsum morphologies (Supplemental Fig. S4e, S4f, respectively).
Heterogeneity observed in all studied aspects of the Puquios is hypothesized to generate a system that contains diverse ecological niches within the brines, porewaters, subaerial structures, and microbial mats that characterize this system. Furthermore, seasonal changes in brine chemistry and water level in the system, are likely to impact all of the ecological niches to varying degrees. Dynamic interplay between physical, chemical, geological and biological processes is thought to generate unique communities that foster different styles of life, including both the community structure and function. For example, the Invierno Altiplanico that occurred in January -February 2017 is hypothesized to have impacted both water levels and microbial communities in Puquio 1. Higher water levels after the Invierno Altiplanico likely lled the pore spaces of the subaerially exposed structures, reduced available oxygen in the environment that fosters endoevaporitic microbial communities, and produced the observed enrichment in the anaerobic microbial population in Puquio 1 samples during the March 2017 eld campaign. The surface water level therefore is an important determinant of the shift from a subaerial to a subaqueous lifestyle as was previously observed 58 . Considering that EC has an indirect correlation with surface water level in the Puquios (R 2 = 0.50 to 0.83), seasonally variable water levels in Puquio 1 may drive the relationship between EC and microbial community structure in Puquio 1, distinguishing it from similar endoevaporitic communities near Puquios 2, 3 and 4. Similar seasonal differences in Anaerobic Deltaproteobacteria abundances were observed between summer and winter samples in a previous study 17 and between the depths analyzed in subaqueous niches in March 2012 29 .
Gradients in environmental parameters have previously been described in the Puquios 32,59,60 , as well as in a variety of Andean microbial ecosystems, such as Lake Tebenquiche in Chile 25  Beyond Andean microbial ecosystems, extreme environments have been observed in other modern and ancient hypersaline lagoons, lakes in arid settings, sabkha environments, and anthropogenically-driven evaporite depositional settings, such as the salt works at Guerrero Negro 69-71 and the EMISAL salt works in Egypt 72,73 . Such extreme conditions typically subject resident biota to wide uctuations in temperature, salinity, pH, dissolved oxygen, total dissolved solids and redox potential 74,75 . In these environments, high light intensity is often coupled with low-redox potential and low oxygen concentration; heavy metals and nutrients may also uctuate (sensu 76 ). As such, observations of horizontal and vertical heterogeneity in brine chemistry, in particular electrical conductivity, and the corresponding diversity of the biological and mineralogical components of the Puquios appear to be characteristic of extreme environments. Similarly, shallow lagoons within these settings may be subject to long periods of desiccation during seasonal changes. Moreover, seasonal and interannual variability in climate, aridity, water activity, UV, and temperature will likely produce even more dramatic environmental gradients than those observed during the November 2017 eld campaign.
Heterogeneity fosters resilience to environmental change as a result of multiple factors at various levels of biological organization 79 . At the landscape level, spatial heterogeneity affects localized responses to perturbations 82 , providing a greater range of resources and microenvironments that can act as buffers to inhabitants 78,83−85 . In the Puquios, spatial heterogeneity of brine chemistry is a dominant driver of microbial diversity and the dominant metabolic pathways between brine communities and E1 layer of the subaerial structures. Most notable was the diversity within the E1 layer communities given they reside above the brines and are largely phototrophic, yet they are impacted by the surrounding brines. The strong brine chemical gradient across the Puquio system provides juxtaposed habitats that can provide nutrients and resources to support heterotrophically dominated aqueous communities as well as halotolerant communities on the extreme end, and autotrophic communities in the subaerial structures.
The overall diversity appears to be created by a geochemical architecture where mineral precipitation creates multiple niches where metabolically exible communities can persist in a polyextreme environment. Indeed, studies have shown that community diversity increases the capacity of biota to survive and/or recover from perturbations (e.g. 86-88 ).
Spatial heterogeneity coupled with complexity promotes species coexistence by providing a wide range of niche environments, which can enhance species diversity 89 . Spatially variable brine chemistry in the Puquios was shown to generate a diversity of biological communities in both planktic and benthic communities of eukaryotes and prokaryotes, a gradient in the role of biological in uence on mineral deposition, and a variety of mineral assemblages. These intra-and interspeci c variations allow redundancy, combined functioning and interaction of species, as well as differing responses, tolerances, and adaptability (sensu 79 ), which is critical for the ecosystem to thrive in response to environmental perturbations, leading to resilience at multiple scales.
On the scale of microns to hundreds of meters, these multi-disciplinary observations of the Puquios provide a case study of the inherent heterogeneity unique to these extreme environments. We argue that the importance of variability observed in the Puquios at multiple scales can likely be extended to encompass the wider complex of Andean microbial ecosystems in all of northern Chile, Argentina, and Bolivia 8 . Spanning several orders of magnitude, this spatial heterogeneity of physical, chemical, biological and geological processes observed is hypothesized to be the cornerstone for the resilience of these ecosystems that persist in some of the harshest conditions on Earth.

Conclusions
The Puquios are a highly diverse ecosystem, and just one example of the high degree of heterogeneity across multiple scales of observation within Andean polyextreme environments that host microbial ecosystems. Each lagoon has distinct spatial trends in electrical conductivity, biological diversity of water column and benthic primary producers and consumers, heterogeneous microbial communities, and a variety of sedimentary structures and textures and mineral phases. While primary producers like phytoplankton and phytobenthos thrive at comparable high diversities throughout the Puquios, their total abundance is greater in the Puquio with the lowest EC. On the other hand, consumers (zooplankton and zoobenthos) exhibit lower diversity throughout, mainly represented by few very specialized taxa (e.g. Artemia franciscana), with varied abundance patterns that do not necessarily relate to EC ranges. Microbial communities in the Puquios also exhibit spatially complex trends that appear to be related to a variety of ecological niches (brines, subaerial structures, submerged microbial mats). Differences in lifestyle are likely produced by the niche's inherent biogeochemical response to environmental change, leading to the observed diversity, function, and structure of the microbial communities that inhabit them. The ability of a community to shift between lifestyles, as observed in the endoevaporitic communities near the margins of the Puquios, should be considered a source of diversity and resilience. However, for a better understanding of the communities' response to typical environmental changes experienced on a seasonal time scale (sensu 32,54,55 ), a multi-disciplinary surveying effort across broader spatio-temporal scales is needed. Bottom types and coupled mineral precipitation revealed a wide depositional spectrum ranging from ne-grained granular gypsum precipitation with microbial in uences to large sparry accumulations of selenite with little or no microbial in uences. Moreover, it is likely that these biological, physical, chemical, and geological observations from a single point in time do not capture the full range of heterogeneity through time and at multiple scales. Seasonal changes in precipitation, evaporation, temperature, and climatic conditions drive temporal variability at least as extreme as the spatial heterogeneity documented in this study 32,54,55 . Initial characterization of heterogeneity across multiple scales and disciplines at a single point in time provides important baseline data to allow comparisons with other active systems and to better understand relationships between heterogeneity and resilience, including the capacity of an ecological community to reorganize after disturbances. Documenting the heterogeneity of the Puquio ecosystem across multiple disciplines provides unprecedented insight into the complexity of an ecosystem at the edge of habitability.

Quanti cation of benthic and pelagic biota
Biotic survey of water column and benthic organisms was performed at 5 sites within each of the 4 Puquios. At each site 3 replicate samples were obtained. For zooplankton and phytoplankton 36 L of water was passed through a 35 µm and 60 µm sieve respectively by means of an electrical immersion pump powered by a portable 12 V battery. Each sample was summarized and introduced into a 50 mL Falcon tube and xed using 2% Lugol for phytoplankton and 5% diluted formalin for zooplankton. For

Characterization of Microbial Communities
Bacterial communities were collected from subaerial gypsum structures at 3 sites from each of the four Puquios (n = 12) using a 20 cm sediment core device. A sample of the microbial mat inhabiting the gypsum structure was collected from the E1 'green layer'. A 1 L brine sample was collected from the water surrounding each subaerial gypsum structure for comparison between the planktonic microbial community and the microbial mat community. The E1 from each core sample of the subaerial gypsum structures was isolated in the lab for DNA extraction and the brines were ltered by a nitrocellulose membrane with a 0. The V1-V3 variable region of the 16S rRNA gene were ampli ed in a PCR thermocycler using the primers 27f (5'-AGAGTTTGATCCTGGCTCAG-3') 105 and 534r (5'-ATTACCGCGGCTGCTGG-3') 22 for bacteria with a barcode on the forward primer. PCR was performed using the HotStarTaq Plus Master Mix Kit (Qiagen, USA) under the following conditions: 94°C for 3 minutes, followed by 28 cycles of 94°C for 30 seconds, 53°C for 40 seconds and 72° C for 1 minute, after which we performed a nal elongation step at 72°C for 5 minutes. PCR products were visualized in 2% agarose gel to determine the success of ampli The sequences were processed further by means of the QIIME (Quantitative Insights Into Microbial Ecology) pipeline 106 : Raw sequences were ltered by base quality score, average base content per read and GC distribution in the reads. Reads that did not cluster with other sequences, i.e. singletons (abundance < 2) were removed. Chimeras were also removed using the UCHIME program 107 . The preprocessed sequences were nally grouped into operational taxonomic units (OTUs) using the clustering program UCLUST at a similarity threshold of 0.97 107 . The pre-processed reads were used to identify the OTUs using QIIME and aligned the representative sequences against the Greengenes core set reference database using PyNAST 106 . A representative sequence for each OTU was classi ed using RDP classi er and the Greengenes OTU database. The alpha-rarefaction was then calculated by means of the "core_diversity_analysis" command and standardized the number of sequences to the smaller sample size by means of Chao 1 (2599 sequences). The rarefacted data using the Primer-7 (Primer-E) software 108 was used to determine the beta diversity and plot the main coordinates graphs.

In Situ Electrical Conductivity Measurements
In

Bottom type mapping and microbe-mineral interactions
To distinguish physio-chemical substrates in the Puquios from substrates that are products of benthic organisms, we adopted a multiscale approach. At the macroscale, broad facies patterns were mapped using aerial imagery; these bottom type maps were re ned using comprehensive ground truthing. This mapping approach identi ed principal bottom types based on upper few centimeters of the substrate surface across the Puquios system (Supplementary Fig. S2 and S3). At the mesoscale, hand samples of the different bottom types were described (Supplemental Table S2). At the microscale, crystal morphologies in context with surrounding matrices were analyzed in the lab using scanning electron microscopy (SEM) with electron dispersive spectroscopy (EDS). This comprehensive approach de ned mapping units used to document spatial distribution of bottom types and distinguished between abiotic mineral structures and those mineral deposits in uenced and/or induced by microbial mats.
High resolution drone images of the Puquios, collected in November 2017, were loaded onto an iPad mini into Global Mapper Mapper Mobile, a powerful GIS data viewing and eld data collection application for iOS. The software application uses the iPad's GPS capability to provide location information for remote mapping projects. The iPad mini was used in the eld during the January 2018 eld campaign where drone imagery and shape les were loaded into the application and subsequently, multiple locations throughout and around each Puquio, including all accessory ponds were investigated. At each location, a photograph was taken, and the bottom type was identi ed and described.