Functional Green-Tuned Proteorhodopsin from Modern Stromatolites

RESEARCH ARTICLE Functional Green-Tuned Proteorhodopsin from Modern Stromatolites

Virginia Helena Albarracín1,2*, Ivana Kraiselburd3, Christian Bamann4, Phillip G. Wood4, Ernst Bamberg4, María Eugenia Farias4, Wolfgang Gärtner5*

1 Planta Piloto de Procesos Industriales y Microbiológicos (PROIMI), CCT, CONICET. Av. Belgrano y Pasaje Caseros. 4000- S. M. de Tucumán, Argentina, 2 Facultad de Ciencias Naturales e Instituto Miguel Lillo, Universidad Nacional de Tucumán, 4000, S. M. de Tucumán, Argentina, 3 Instituto de Biología Molecular y Celular de Rosario (IBR - CONICET), Facultad de Ciencias Bioquímicas y Farmacéuticas (FBIOYF - UNR), Suipacha 590, 2000, Rosario, Santa Fe, Argentina, 4 Max-Planck-Institute for Biophysics, Max-von-Laue-Straße 3, D-60438 Frankfurt am Main, Germany, 5 Max-Planck-Institute for Chemical Energy a11111 Conversion, Stiftstrasse 34–36, D-45470 Mülheim, Germany

* [email protected] (VHA); [email protected] (WG)

Abstract

OPEN ACCESS The sequenced genome of the poly-extremophile Exiguobacterium sp. S17, isolated from modern stromatolites at Laguna Socompa (3,570 m), a High-Altitude Andean Lake (HAAL) Citation: Albarracín VH, Kraiselburd I, Bamann C, Wood PG, Bamberg E, Farias ME, et al. (2016) in Argentinean Puna revealed a putative proteorhodopsin-encoding gene. The HAAL area Functional Green-Tuned Proteorhodopsin from is exposed to the highest UV irradiation on Earth, making the microbial community living in Modern Stromatolites. PLoS ONE 11(5): e0154962. the stromatolites test cases for survival strategies under extreme conditions. The heterolo- doi:10.1371/journal.pone.0154962 gous expressed protein E17R from Exiguobacterium (248 amino acids, 85% sequence Editor: Erika Kothe, Friedrich Schiller University, identity to its ortholog ESR from E. sibiricum) was assembled with retinal displaying an GERMANY absorbance maximum at 524 nm, which makes it a member of the green-absorbing PR-sub- Received: February 9, 2016 family. Titration down to low pH values (eventually causing partial protein denaturation) indi- Accepted: April 21, 2016 cated a pK value between two and three. Global fitting of data from laser flash-induced

Published: May 17, 2016 absorption changes gave evidence for an early red-shifted intermediate (its formation being below the experimental resolution) that decayed (τ1 =3.5μs) into another red-shifted inter- Copyright: © 2016 Albarracín et al. This is an open τ μ τ access article distributed under the terms of the mediate. This species decayed in a two-step process ( 2 =84 s, 3 = 11 ms), to which the Creative Commons Attribution License, which permits initial state of E17-PR was reformed with a kinetics of 2 ms. Proton transport capability of unrestricted use, distribution, and reproduction in any the HAAL protein was determined by BLM measurements. Additional blue light irradiation medium, provided the original author and source are reduced the proton current, clearly identifying a blue light absorbing, M-like intermediate. credited. The apparent absence of this intermediate is explained by closely matching formation and Data Availability Statement: All relevant data are decay kinetics. within the paper and its Supporting Information files.

Funding: This work was supported by the PICT 2010 2991 Project (FONCyT, Argentina), PIP-CONICET 2013-0519 Project, the Proalar Agreement (MINCYT- DAAD) DA/13/01, The Max-Planck-Society (MPI for Biophysics and MPI for Chemical Energy Conversion, Introduction Germany). A Marie Curie FP7-People-2010-IIF Microbial rhodopsins are considered to be the simplest energy-harvesting photoreceptors con- EXTREMOPHIL (273831) in Germany and its return phase in Argentina (PIIFR-GA-2010-910831- sisting of a single, membrane-embedded protein covalently binding a retinal (Vitamin-A alde- EXTREMOPHIL) supported research work of VHA. hyde) chromophore [1]. The study of the microbial archetype, bacteriorhodopsin (BR) from The funders had no role in study design, data the archaeon Halobacterium salinarum [2], led to the suggestion that BR evolved as a

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 1/18 proteorhodopsin from Modern Stromatolites collection and analysis, decision to publish, or consequence of exposure to the extreme halophilic, oxygen-deprived environment of H. sali- preparation of the manuscript. narum. However, further surveys demonstrated the worldwide existence of many microbial Competing Interests: The authors have declared rhodopsins with a spectrum of functions in ecological niches, thus being probably the most that no competing interests exist. abundant phototrophic system on Earth [3]. Out of this versatile family of proteins, the proteorhodopsins (PRs), bacterial ion pumps, represent the largest group [3]. They were mostly found in microbial communities living in environments where solar exposure is high i.e. the upper layer of oceans, freshwater, brackish and hypersaline aquatic ecosystems [4–18], sea-ice matrix [19], non-marine permafrost [20, 21] and phyllosphere [22]. Moreover, proteorhodopsin genes are widely distributed among divergent bacterial taxa included alpha, beta, delta and gamma proteobacteria, Bacteroidetes, Flavobacteria, Firmicutes, Deinococcus-Thermus, Chloroflexi, Planctomycetes and Actinobac- teria [23–28]. Detailed molecular studies on various PRs and their three-dimensional structure [21, 29] (in addition to the enormous information on structure and function of BR) have generated a precise picture of their light-induced reactions [3]. Key amino acids, already identified as instrumental for the function of BR have been allocated also to the PRs [30–32]. A further evolutionary aspect of PRs was detected from sequence alignments, together with spectroscopic measurements. Two subfamilies of PRs were identified, distinguished by their absorbance maxima [33–35]. The group coined GPRs (green-absorbing) showed absorbance maxima around 525 nm, whereas the other group near 490 nm (BPRs, blue-absorbing). This switch in absorbance maxima is accomplished by a single amino acid exchange and is understood as an adaptation to the prevailing spectral conditions at different depths in a water column with green light dominating close to the surface and solely blue light remaining at greater water depths [33–35]. Microbes might profit from PRs in a variety of ways; planktonic ones may gain a competi- tive advantage by using these photoproteins to generate a light-driven proton gradient for the synthesis of ATP [36, 37]. In extreme environments though, PRs may even be paramount to counteract starvation and other hostile conditions thanks to the potential for harvesting extra energy by means of phototrophy [16, 19]. A characteristic example of a highly irradiated, extreme environment is the High-Altitude Andean Lakes (HAAL) area in the Central Andes region in South America [38, 39]. These shallow lakes and salterns dispersed at altitudes above 3,000 m are exposed to a unique combination of severe conditions including high solar global and UV irradiation, hypersalinity, wide fluctuations in daily temperatures, desiccation, high pH, and high concentrations of toxic ele- ments including arsenic, which have brought forth a distinguished microbiodiversity of extremophiles [38–42,62,63]. Among microbialites of the HAAL [41, 43–47], modern stromatolites were first described as thriving in the shore of Laguna Socompa, a hypersaline and moderate alkaline lake located at 3,570 m at the base of the active Socompa volcano. These stromatolites (41) exhibit an atypi- cal microbial community with abundant representatives of Deinococcus-Thermus, Rhodobac- teraceae, Desulfobacterales and Spirochaetes. Their metagenome revealed a high proportion of sequences that could not be classified at phylum level, depicting less than 80% identity to the best hit in the NCBI database, suggesting the presence of novel distant lineages [41]. Thus, the study of microbes associated with Socompa stromatolites and their environment could provide further understanding of their adaptation and the function of proteins working under extremely harsh conditions. Several bacterial strains from Socompa stromatolites were isolated, physiologically characterized and subjected to genome pyrosequencing [18, 48, 49]. Screening of the genomes of these strains i.e. Exiguobacterium sp. S17, Sphingomonas sp. S17 and Salinivibrio spp. S10B and

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 2/18 proteorhodopsin from Modern Stromatolites

S34, yielded genes putatively encoding essential traits for survival under multiple environmental extreme conditions, e.g., high levels of UV radiation, elevated salinity, and the presence of poisoning arsenic concentrations [18, 45, 48–50]. Likewise, genes putatively coding for the bacteriorhodopsin protein family were found showing sequential features indicating the presence of proteorhodopsins and xantorhodopsins [18, 45–47]. The former findings called for a more detailed, molecular investigation of the mechanisms involved in the resistance of these strains to extreme but common impairing factors in its original environment. We herein present the first functional characterization of E17R, a PR from Exiguobacterium sp. S17, a halotolerant, highly arsenic resistant extremophile isolated from a modern stromatolite, located at an altitude of more than 3,500 m above sea level.

Materials and Methods Strain, media and culture conditions The extremophile Exiguobacterium sp. S17 was previously isolated from modern stromatolites located at the shore of Lake Socompa (3,570 m; Fig 1)[50], and is currently maintained in the culture collection of Laboratory of Microbial Research on Andean Lakes (National System of Biological Data). For preculture and DNA extraction, S17 was grown in LB media or on LB- agar plates (15.0 g/L).

DNA isolation Genomic DNA from Exiguobacterium sp. S17 was purified from cells grown on LB broth for 24 h at 30°C and harvested by centrifugation (3,000 g for 10 min at 4°C). The pellets were washed twice with distilled water. Total genomic DNA was extracted with the DNeasy Blood and Tissue Kit (Qiagen) using the specific protocol for Gram-positive bacteria following the manufacture’s recommendations.

Gene cloning E. coli strain XL-1 Blue (Stratagene) was used for all DNA manipulations. To amplify the putative PR gene from Exiguobacterium sp. S17, two gene-specific primers PR-S17F 50 GATCGAAT TCGATGGACGAAGAAGTCAATTTAC and PR-S17R 50 GATTCTCGAGCGTTTTAATTTGTTTG AGCGTCGCGACGCTCAAACAAATTAAAACG were designed according to the published genomic DNA sequence of S17 [48]. The gene start codon (ATG) is shown in bold-type in the forward primer. Restriction sites for EcoRI and XhoI were introduced in the forward and reverse primers (underlined). E17R -gene was PCR amplified and cloned into the pET-26b(+) vector (Novagene). The vector was transformed into Escherichia coli XL1 blue cells for further analysis.

Phylogenetic analysis of protein sequences The sequence of the E17R protein from Exiguobacterium sp. S17 is available in GenBank (WP_016509804) together with the draft genome obtained by pyrosequencing (ASXD01). All sequences of PR-related proteins with significant homology (S1 File) were retrieved from pub- lic databases via the National Center for Biotechnology Information web site (www.ncbi.nlm. nih.gov). Multiple sequence alignments and phylogenetic analyses were carried out using MUSCLE [51] and Phylogeny.fr [52, 53]. The evolutionary history was inferred using the Max- imum likelihood method [54]. The tree was drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The percentages of replicate trees in which the associated taxa cluster together in the bootstrap test (1,000

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 3/18 proteorhodopsin from Modern Stromatolites

Fig 1. A. Socompa Stromatolites at Laguna Socompa, the original habitat of Exiguobacterium sp. S17. B. Macroscopic appearance of S17 colonies. C. Scanning electron micrograph of S17 cells. Bar: 1 μm. doi:10.1371/journal.pone.0154962.g001

replicates) are shown next to the branches. The analysis involved 22 amino acid sequences of PR-genes of Exiguobacterium strains. All positions containing gaps and missing data were eliminated.

Genomic environment analysis of proteorhodopsin Genes for microbial rhodopsin and its genomic environment were identified using the RAST annotation platform [55] and the IMG analysis system [56]. The search and comparison was performed with all genomes available for the Exiguobacterium genus.

Protein expression and purification Heterologous expression was performed with E. coli BL21 (DE3) cells transformed with the recombinant vector pET26-E17R. Cells from an overnight culture were used to inoculate a 200 mL culture, which then was used to inoculate 6 L of LB medium. Cells were grown with 40 μg/

mL kanamycin, until an OD600 of 0.6 was reached. At that point, expression was induced with

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 4/18 proteorhodopsin from Modern Stromatolites

0.5 mM IPTG. Cells were harvested by centrifugation after 12 h of incubation at 30°C. Cells μ were resuspended in 50 mM Tris, 5 mM MgCl2 pH 8.5 added with 1 mM PMSF, 5 M retinal (final concentrations) and DNase. Cells were then stirred for 1 h at 4°C and then lysed using a microfluidizer (Microfluidics Corporation, M110LA). The cell lysate was centrifuged for 1 h at 186,000 g at 4°C. The pellet was resuspended in 20 mM HEPES, 100 mM NaCl, pH 7.4, added with 5 μM retinal and 1% (w/v) n-dodecyl-β-maltoside (DDM) (final concentrations). Follow- ing a second centrifugation (1 h at 186,000 g, 4°C), the supernatant was filtered and proteorhodopsin was obtained by IMAC purification using a His-tag (ÄKTA, GE Healthcare).

Steady state absorbance and fluorescence spectroscopy Steady state spectra were recorded at room temperature either on a Hitachi U3000 spectropho- tometer or on a Hitachi F-4500 spectrofluorometer. The extinction coefficient of E17R was determined from the difference spectrum of a sample before and after treatment with 10 mM hydroxylamine (final concentration) and illumination (>480 nm) at room temperature. The value for E17R can be quantified from the known extinction coefficient of the retinal oxime, ε ’ -1 -1 max =33600 M cm [57].

pH-titration – For pH-tritation experiments, proteorhodopsin samples were adjusted to A540 = 0.1 0.5 in a buffer containing 20 mM HEPES, 100 mM NaCl, and 0.1% (w/v) DDM at pH 7.4. The titration was carried out in a 1 × 1 cm cuvette with continuous stirring. The pH values were adjusted with 1 M NaOH or 1 M HCl, before the spectra were recorded (Cary50, Varian, USA). Revers- ibility was checked at the end of the titration. Data were fitted using the Henderson-Hassel- balch equation.

Flash photolysis Samples were prepared with identical buffer conditions as for the pH titration with pH = 7.4. Transient absorbance changes were recorded on a home-built flash-photolysis setup. The reac- tion was started by a 10 ns laser pulse from an excimer laser pumped dye laser (Coumarin 307, λ = 503 nm, 5–10 mJ cm-2). Light from 75 W XBO lamp was filtered by narrow-band interfer- ence filters and passed through the sample and a monochromator, before it was detected by a photodiode. Absorbance changes were followed over time on two oscilloscopes for the following wavelengths: 380, 400, 420, 440, 460, 480, 500, 517, 540, 562, 580, 600, 620, 645 nm. The number of data points was reduced by a logarithmic interpolation procedure before the data was analyzed by a global fitting routine. To this end, data blocks of centuries are averaged over 2(n-1) points per nth century yielding 880 data points for the initial 45k points of the raw data for each of the oscilloscope with sampling rates of 200 ns per point and 20 μs per point, respectively. The two datasets were combined yielding 1520 point covering a range up to 1 second.

Electrophysiology of reconstituted E17R in liposomes Proteorhodopsin was reconstituted into liposomes. Liposomes were formed from E. coli polar lipids, resuspended in 20 mM HEPES, 100 mM NaCl, and 2% (w/v) n-octyl-β-glucopyranoside (OG) at pH 7.4, to which 7.5 μM protein was added. The resulting mixture was incubated on ice for 15 min. Detergent-adsorbing beads (BioBeads, BioRad, Germany) were added to the solubilized protein and the solution was stirred overnight at 4°C. − Black Lipid Membranes (BMLs) were formed across an orifice (area ~0.5×10 3 cm2) between two compartments of a Teflon cuvette filled with a buffer containing 50 mM Tris,

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 5/18 proteorhodopsin from Modern Stromatolites

5 mM MgCl2, pH 8.5. The membrane-forming solutions contained 1.5% (w/v) diphytanoyl- phosphatidylcholine (Avanti Biochemicals, Birmingham, AL) and 0.025% (w/v) octadecyla- mine (Riedel-de-Haën, Hannover, Germany) in n-decane to obtain a positively charged membrane. Proteorhodopsin-containing proteoliposomes were added to the compartment under gentle stirring. A combination of ionophores (2 μM monensin and FCCP) was added, which effectively permeabilizes the compound membrane system (final conductance 50–100 nS). A 100 W Osram HBO mercury arc lamp with a combination of bandpass filters removing UV and IR irradiation (<380 nm and >750 nm) and a cut-off filter removing irradiation <495 nm was used for sample illumination. Blue light effect was generated by additional irradiation with light longer than 360 nm. The measurements were carried out at room temperature. The photocurrents were recorded with pClamp9.0 software via a MiniDigi 1A digitizer (Axon Instruments).

Results A green proteorhodopsin from Exiguobacterium sp. S17 E17R, a gene predicted as belonging to the bacteriorhodopsin family (L479_RS10205/ WP_016509804), was found in the genome of Exiguobacterium sp. S17. It encodes a protein of 248 amino acid residues (744 bp), E17R, with features typical of a microbial rhodopsin. All essential amino acid residues of the energy transducing rhodopsins are conserved (Figs 2 and 3). E17R has 85% sequence identity with the green-PR of the psycrophilic Exiguobacterium sibiricum (ESR) for which the crystal structure was recently solved (PDB: 4HYJ) [21]. Align- ment of these two proteins sequences with that of the blue-light absorbing proteorhodopsin (BPR) of uncultured Gamma-proteobacterium Hot 75m4 (http://www.uniprot.org/taxonomy/ 245185) clearly showed the instrumental difference in an amino acid position associated with spectral tuning. Both Exiguobacterium proteins contain a leucine residue at position 94 (E17R, position 93 at ESR), which fine-tunes the PR absorption peak towards green light (absorption maximum 524 nm, Fig 4). This is in accordance with the shallow nature of the stromatolite formation in the shore of L. Socompa, the original environment of S17 (Fig 1). In turn, the corresponding position in BPR (106) is occupied by a glutamine, which is known to yield an absorption maximum near 490 nm [5, 32], a feature for proteorhodopsin-bearing microbes from deep-sea in coincidence with the light quality available at depths greater than 50 m. A BLASTP search in all publicly available genomes of the Exiguobacterium genus was performed. Twenty-one genomes displayed a sequence with significant similarity (higher than 80%) to that of S17 proteorhodopsin (S1 File). E17R shares the highest sequence identity with Exiguobacterium sp. 8-11-1 (98%), Exiguobacterium sp. AB2 (95%) and Exiguobacterium sp. AT1b (94%). All PR sequences display Leu in position 106 (following the numbering of the BPR Gamma-proteobacterium Hot 75m4), assigning them to the “green-tuned” PR subfamily. A phylogenetic tree (Fig 3) constructed with the 22 related rhodopsins (S1 File) revealed two divergent clades, one including E17R (mostly form by alkaliphiles/thermophiles) and the other one mostly formed by psycrophiles/psycrotolerants including two branches, the one of Exiguo- bacterium acetylicum and the one of E. sibiricum ESR. Thus, in view of the phylogenetic dis- tance from ESR, further biochemical and spectroscopic studies of the novel proteorhodopsin E17R were performed.

Biochemical and steady-state spectral properties of E17R Gene E17R (L479_RS10205/WP_016509804) encodes for a protein with 248 amino acids (MW: 27.4 kDa) that showed typical features of a microbial rhodopsin. Secondary structure prediction identified seven transmembrane α-helices, a lysine residue serving as putative

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 6/18 proteorhodopsin from Modern Stromatolites

Fig 2. Multiple protein alignment of PR from Exiguobacterium sp. S17 (E17R), green-light absorbing proteorhodopsin from Exiguobacterium sibiricum (ESR) and blue-light absorbing proteorhodopsin from the uncultured gamma-proteobacterium “Hot 75m4” (BPR). Residues shared between all the protein variants are marked with asterisks. Single amino acid residue at position 106 (BPR numbering) that functions as a spectral tuning switch and accounts for most of the spectral difference between the two pigment families is highlighted in light grey. Primary proton acceptor and donor are highlighted in dark grey (D86) and with a frame (K97), respectively. The Schiff base (K232 for ESR, K226 for E17R) is indicated by a diamond. Residues differing between both green-PRs are indicated with arrows. The seven transmembrane α- helices are indicated with blue bubbles. doi:10.1371/journal.pone.0154962.g002

retinal binding site at the seventh helix and proton donor and acceptor residues that identify this protein as a light-driven proton pump (Fig 2). Expression and purification of E17R yielded λ after purification and retinal addition a protein with max = 524 nm (solubilized in 0.03% [w/v] DDM, Fig 4). The chromophore extinction coefficient was determined as 50,000 M-1 cm-1 from the difference spectrum after bleaching in the presence of hydroxylamine, which formed

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 7/18 proteorhodopsin from Modern Stromatolites

Fig 3. Phylogenetic tree inferred using the Maximum likelihood method of 21 amino acid sequences of Exiguobacterium spp. closely related with E17R. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1,000 replicates) are shown next to the branches. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. doi:10.1371/journal.pone.0154962.g003

the retinal oxime (inset in Fig 4)[57]. The ratio of protein (280 nm) to chromophore (524 nm) absorbance can be taken as a purity factor and it amounts to ~1.3 (taking into account increased scattering around 280 nm). The value is slightly higher than the theoretical limit of ~1 as calculated from the extinction coefficients for the protein moiety and the chromophore (Fig 4). We followed pH-induced absorbance changes in E17R between pH 2 and pH 10 in order to identify de- or reprotation of amino acids in close proximity to the chromophore (Fig 5). As reported for the closely related ESR [58], pH-titration reveals an ongoing red-shift of the absorption maximum upon acidification. Usually, this colour change is attributed to the protonation of the counterion of the protonated Schiff-base (D85 in ESR). pH-Titration experiments with E17R showed a similar absorbance change suggesting that D86 is the Schiff base counterion in E17R. pH-values lower than 2.0 could not be reached, as the protein lost stability irreversibly. In the accessible pH range, we observed a complex titration behaviour, similar to that described for ESR. For this protein, the data was interpreted as a strong interaction between the counter-ion D85 (D86 for E17R) and a conserved histidine residue (H57 in ESR, H58 in E17R) located in the second transmembrane helix (TMH B) [55]. For E17R, the fit of

the absorbance changes shifting from 525 nm to 535 nm or higher gives a pK2~3 (Fig 5). Note that this is an approximate value, as the final state at acidic pH (<2) was not reached. The

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 8/18 proteorhodopsin from Modern Stromatolites

Fig 4. Absorption spectrum of E17R in 20 mM Hepes, 100 mM NaCl, 0.03% DDM. The inset shows the difference spectrum after bleaching with hydroxylamine. The extinction coefficient was calculated using the absorbance of the oxime product ΔA(oxime) in comparison to the bleached rhodopsin absorbance ΔA(Rh) as reference. doi:10.1371/journal.pone.0154962.g004

second transition is observed with a pK1 ~10 accompanied by smaller spectral changes (from 517 nm to 525 nm). The fluorescence data is in line with the argument given in Balashov et al. (2012) [58] for the assignment of the titratable groups. The major increase in fluorescence intensity observed at wavelengths longer than 620 nm can be attributed to the protonation of the proton acceptor D86 and takes place only at pH values <3 (data not shown). On the other hand, there is almost no change in fluorescence intensity at pH values >8. Therefore, the second transition might be assigned to deprotonation of the His interacting strongly with D86. It was proposed that this interaction keeps the counter-ion in its deprotonated state even at low pH values. Such behaviour is different from, e.g., the very first identified green-absorbing PR

that has a pKa ~7.5 for its counter-ion D100 [58].

Flash photolysis Flash photolysis as a time-resolved detection method for absorption changes allows following —after a short laser pulse (routinely 10 ns)—absorbance changes at a wide wavelength range on a time domain from few microseconds into the seconds range. Detection of absorbance

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 9/18 proteorhodopsin from Modern Stromatolites

Fig 5. A. Absorbance spectra of E17R at three selected pH values. B. Titration curve as determined from the absorbance maxima. doi:10.1371/journal.pone.0154962.g005

changes, followed by fit procedures provides details on the lifetime of intermediate states formation and decay. In most cases, the global fit approach is used where all detected absorbance changes over time and wavelengths are subjected to one single fit model. This data evaluation yields lifetimes and also spectral shape of the intermediates through which the protein travels. Here, we followed the built-up and decay of photo-intermediates from time-resolved absorbance changes in the microsecond to second time range after excitation of E17R’s initial dark state by a nanosecond laser flash (Fig 6A). The results allow us to deduce a photocycle with a

Fig 6. Transient absorbance changes of E17R. A. The kinetics of the absorbance changes are shown for selected wavelengths (399 nm, 517 nm and 598 nm). B. The decay-associated spectra are depicted as obtained from the global fit. The corresponding time constants are given in the figure. doi:10.1371/journal.pone.0154962.g006

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 10 / 18 proteorhodopsin from Modern Stromatolites

lifetime of 82 ms for the recovery of the parental dark state. Global fit analysis of the kinetic traces clearly identifies formation and decay of different intermediates and the bleaching and recovery of the ground state (Fig 6B). The first intermediate is red-shifted compared to the ground state. Its formation is below the time resolution of our system. It decays partly in a fast τ μ process ( 1 = 3.5 s), followed by the formation of another red-shifted intermediate. This spe- τ μ τ cies is built in a two-step process with small amplitudes ( 2 =84 s and 3 = 11 ms). A lifetime τ of 4 = 82 ms corresponds to the final decay of all red-shifted intermediates and the concomitant formation of the parental 524 nm-form. These flash photolysis experiments did not reveal the presence of a deprotonated species (M-state). Such behaviour is reminiscent of the ESR that also has no noticeable transient concentration of an M-state [58–60]. It was argued for the ESR protein that the absence of this intermediate is solely due to kinetics parameters, i.e., caused by similarly fast formation and decay kinetics.

Light-induced proton transport currents The protein pumping activity for the closely related ESR protein was measured in E. coli cell suspensions [58]. For E17R, we used a more direct proof employing the black lipid membrane (BLM) method [61]. E17R was reconstituted in proteoliposomes that were attached to BLM yielding capacitive coupling. Light activation generated a current of the compound membrane (Fig 7A). Doping the membrane with the protonophore FCCP makes the membrane

Fig 7. Electrometric record of E17R (in 20 mM Hepes, pH 7.4) with the BLM system. A. Transient currents and B. after the addition of the protonophore FCCP. Black bars indicate illumination with a 75 W XBO long-pass filtered at >495 nm. The grey bar shows the additional excitation of the M-state (>380 nm). doi:10.1371/journal.pone.0154962.g007

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 11 / 18 proteorhodopsin from Modern Stromatolites

proton selective and turns it into a DC coupled system, now also allowing measurement of stationary transport currents as a result of the light-driven proton transport (Fig 7B). The data clearly illustrates the light-driven proton translocation by E17R, but, interestingly, this method displays a blue-light effect in the photocycle after additional irradiation with blue light (Fig 6B, illumination >380 nm), and thus gives evidence for the presence of an M-like intermediate. Its presence becomes visible by the concomitant excitation of the ground state and the M-state that absorbs in the blue-/near UV region. The excitation of the M-state provides a shunt in the photocycle abolishing proton pumping which leads to a current ampli- tude decrease.

Discussion In the course of a genomic project from extremophilic microbes from HAAL, a gene (E17R) predicted as belonging to the bacteriorhodopsin family was found in the genome of Exiguobac- terium sp. S17. This halotolerant, UV-resistant microorganism [48] was originally isolated from modern stromatolite at L. Socompa (3,570 m; Fig 1)[41] being able to grow at concentrations up to 100 mM As[V] and 5mM As[III] [45, 50]. Neighbor-Joining and maximum likelihood trees based on 16SrRNA gene sequences and MLSA of Exiguobacterium sp. S17 indicated that it is closely related and branched together with other mesophile strains (Exiguobacterium sp. 8-11-1, Exiguobacterium pavilonensis RW-2, Exiguobacterium sp. AT1b, Exiguobacterium sp. N139) clearly distant from psycrophilic strains clade: i.e. Exiguobacterium sibiricum 255–15 and Exiguobacterium antarticum B7 [45]. Sequence comparison and phylogenetic analysis of E17R indicated that it was a green-tuned proteorhodopsin, being the case also to all proteorhodopsins found in the so far reported Exi- guobacterium genomes. Accordingly, the original environments of most isolates are shallow or surficial such as microbial mats, soil, rhizosphere, glacier ice, or hot springs (S1 File). Neverthe- less, PRs were also present in Exiguobacterium strains isolated from ecological niches where light is less abundant such as marine sediments, industrial waste and even animal gastrointesti- nal tracts (S1 File). The ecological significance of PRs for survival under stressful conditions was evident for a Vibrio sp. strain in which deletion of the PR gene reduced survival during car- bon starvation [64]. Light-dependent promotion of starvation survival and enhanced PR expression was also established for the marine dinoflagellate Oxyrrhis marina [36]. Bohorquez et al., (2012) [16] found diverse PR-like genes in microorganisms from shallow, oligotrophic hot spring waters at the Colombian high-Andes indicating that rhodopsin photosystems could be advantageous in these acidic hot springs by contributing to survival in ecosystems that receive abundant sunlight and where alternative energy sources may vary or be scarce. For the Socompa stromatolite community, the use of light to counteract stressful conditions could also be an interesting strategy especially to those communities situated at the top layer where light is fully available but where UV stress is maximal [41]. In accordance, all four available genomes isolated from stromatolite indigenous strains (from the top 5 mm) contained microbial rhodopsin-likes genes even despite belonging to distant taxonomical groups (firmicutes, proteobacteria and bacterioidetes) [18, 48, 49]. Although isolated from completely different niches (Andean stromatolites in hypersaline lakes vs. Siberian Permafrost), E17R shows highest homology to ESR for which recently a three-dimensional structure was reported [21]. Both PRs share the essential pump amino acids and show clear differences to the blue light-absorbing BPR. Common between all three proteins compared here (Fig 2) is the proton acceptor (D86, D85 in ESR), and also the close interaction with His58 (His57 in ESR, and His76 in BPR) that keeps the proton acceptor in its deprotonated state during the photocycle. Different, however, is a lysine as putative proton

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 12 / 18 proteorhodopsin from Modern Stromatolites

donor (Lys97 in E17R, Lys96 in ESR). Here, BPR resembles bacteriorhodopsin as both proteins carry an acidic amino acid as proton donor (D96 in BR, and E109 in BPR). Apparently, the protonated side chain of lysine (better coined an ammonium group) undergoes transient deprotonation and therefore functionally complements the normally present acidic side chain of, e.g., D96 in BR. Another feature shared by both green-absorbing PRs (different from BPR and BR) is an extension of the loop between helices three and four. Seven additional, nearly fully conserved amino acids- NGGFTQL (N is exchanged for A in ESR)-are present in the green absorbing PRs. Inspection of the sequences of E17R and ESR identify positions at which amino acid are different (vertical arrows in Fig 2). In most cases aliphatic amino acids are exchanged against each other. Alternatively, polar (Thr, Ser) and acidic amino acids (Asp, Glu) were exchanged. Over- all, the exchanges found in ESR are of a more polar nature compared to the corresponding positions in E17R (numbering follows E17R, exchanges in ESR are given after the dash: T76-D, G147—S, N183—K, V224—T). A more drastic change is found in ESR that carries an ion pair (K172 and D178), located at the entrance into the sixth helix at a position, where E17R carries a serine (173) and a lysine (179). Highest conservation between all three aligned PRs is found in helix three with the proton donor and acceptor residues, except for the exchange of the acidic proton donor in BR and BPR for a lysine. In order to accomplish its function, the pK value of this lysine has to be low- ered significantly compared to aqueous solution (as that of aspartate or glutamate is remark- ably increased for maintaining their protonated state). Helix seven also shows a high degree of sequence similarity, except for the exchange V220 ! F (E17R to ESR) and V224 ! T in direct vicinity to the retinal-binding lysine. Flash photolysis of E17R showed a photocycle comprised only of long wavelengths absorbing intermediates. From the initial photoproduct on (detection of its formation being below the time resolution of the instrument), all further intermediates show continuously τ μ μ increasing life times ( i =3.5 s, 84 s, and 11 ms). The initial PR524 state is then formed with τ 4 = 82 ms. At first glance the apparent absence of the M-intermediate in the transient absorbance changes at the measuring pH value of 7.4 is puzzling, as a de- and reprotonation is instrumental for a vectorial proton transport. In fact, the BLM measurements clearly demon- strate the presence of a photocurrent and blue light irradiation opens a short circuit pathway leading to reduce current. Thus, the lack of an M-intermediate at the pH-value chosen for this experiment is due to purely kinetic reasons (see also the pH-dependent measurements in [58]).

Conclusions In this work, we have successfully cloned and heterologously overexpressed a proteorhodopsin gene from a halotolerant, UV- and arsenic resistant strain originally isolated from a high-altitude Andean Lake (ca. 3,600 m). To our knowledge, this is the first report of a functional proteorhodopsin isolated from modern stromatolites. Its presence in a stromatolite and its functional, light-driven proton transport activity identifies Exiguobacterium sp. S17 as an important member of this complex microbial community. Further studies including knock-out strategies on E17R are in progress to link light-harvesting functionality to microbial physiol- ogy. High-throughput sequencing metagenomic projects on Socompa stromatolites are also being conducted to ultimately reveal the ecological significance of these versatile novel ion pumps in the structure and function of the microbial extremophilic communities. The studies offer a high window to study ancient phototrophic systems in homologous stromatolites of their Precambrian counterparts.

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 13 / 18 proteorhodopsin from Modern Stromatolites

Supporting Information S1 File. Isolation source and other relevant information on Exiguobacterium spp. strains in which PR –like genes were found. (DOCX)

Acknowledgments The authors acknowledge the generous financial support by the PICT 2010 2991 Project (FON- CyT, Argentina), PIP-CONICET 2013–0519 Project, the Proalar Agreement (MINCYT- DAAD) DA/13/01, and the Max-Planck-Society (MPI for Biophysics and MPI for Chemical Energy Conversion, Germany). A Marie Curie FP7-People-2010-IIF EXTREMOPHIL (273831) in Germany and its return phase in Argentina (PIIFR-GA-2010-910831-EXTREMO- PHIL) supported research work of VHA. MEF and VHA are researchers from the National Research Council (CONICET) in Argen- tina. IK is a postdoctoral fellow from CONICET. Electron micrographs used in this paper were taken at the Center for Electron Microscopy (CIME) belonging to UNT and CCT, CONICET, Tucumán, Argentina.

Author Contributions Conceived and designed the experiments: CB PGW VHA WG. Performed the experiments: CB PGW IK. Analyzed the data: CB WG VHA. Wrote the paper: VHA CB WG. Obtained funding for the original project idea: MEF EB WG PGW VHA. Conceived the original HAAL project: MEF. Revised the manuscript: EB PGW. Performed HAAL sampling expeditions and provided the stromatolites and strains and the S17 genome sequence: MEF. Provide reagents/materials and analysis tools: EB PGW WG. All authors read and approved this manuscript: VHA IK CB PGW EB MEF WG.

References 1. Spudich JL. The multitalented microbial sensory rhodopsins. Trends in Microbiology. 14(11):480–7. doi: 10.1016/j.tim.2006.09.005 PMID: 17005405 2. Oesterhelt D, Stoeckenius W. Rhodopsin-like protein from the purple membrane of Halobacterium halobium. Nature: New biology. 1971; 233(39):149–52. Epub 1971/09/29. PMID: 4940442. 3. Bamann C, Bamberg E, Wachtveitl J, Glaubitz C. Proteorhodopsin. Biochimica et biophysica acta. 2014; 1837(5):614–25. Epub 2013/09/26. doi: 10.1016/j.bbabio.2013.09.010 PMID: 24060527. 4. Beja O, Aravind L, Koonin EV, Suzuki MT, Hadd A, Nguyen LP, et al. Bacterial rhodopsin: evidence for a new type of phototrophy in the sea. Science. 2000; 289(5486):1902–6. Epub 2000/09/16. PMID: 10988064. 5. Beja O, Spudich EN, Spudich JL, Leclerc M, DeLong EF. Proteorhodopsin phototrophy in the ocean. Nature. 2001; 411(6839):786–9. Epub 2001/07/19. doi: 10.1038/35081051 PMID: 11459054. 6. Giovannoni SJ, Bibbs L, Cho JC, Stapels MD, Desiderio R, Vergin KL, et al. Proteorhodopsin in the ubiquitous marine bacterium SAR11. Nature. 2005; 438(7064):82–5. Epub 2005/11/04. doi: 10.1038/ nature04032 PMID: 16267553. 7. Sabehi G, Loy A, Jung KH, Partha R, Spudich JL, Isaacson T, et al. New insights into metabolic properties of marine bacteria encoding proteorhodopsins. PLoS biology. 2005; 3(8):e273. Epub 2005/07/13. doi: 10.1371/journal.pbio.0030273 PMID: 16008504; PubMed Central PMCID: PMC1175822. 8. Sabehi G, Massana R, Bielawski JP, Rosenberg M, Delong EF, Beja O. Novel Proteorhodopsin variants from the Mediterranean and Red Seas. Environmental microbiology. 2003; 5(10):842–9. Epub 2003/09/27. PMID: 14510837. 9. Frigaard NU, Martinez A, Mincer TJ, DeLong EF. Proteorhodopsin lateral gene transfer between marine planktonic Bacteria and Archaea. Nature. 2006; 439(7078):847–50. Epub 2006/02/17. doi: 10.1038/ nature04435 PMID: 16482157.

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 14 / 18 proteorhodopsin from Modern Stromatolites

10. Stingl U, Desiderio RA, Cho JC, Vergin KL, Giovannoni SJ. The SAR92 clade: an abundant coastal clade of culturable marine bacteria possessing proteorhodopsin. Applied and environmental microbiology. 2007; 73(7):2290–6. Epub 2007/02/13. doi: 10.1128/AEM.02559-06 PMID: 17293499; PubMed Central PMCID: PMC1855678. 11. Campbell BJ, Waidner LA, Cottrell MT, Kirchman DL. Abundant proteorhodopsin genes in the North Atlantic Ocean. Environmental microbiology. 2008; 10(1):99–109. Epub 2008/01/24. doi: 10.1111/j. 1462-2920.2007.01436.x PMID: 18211270. 12. Nguyen D, Maranger R, Balague V, Coll-Llado M, Lovejoy C, Pedros-Alio C. Winter diversity and expression of proteorhodopsin genes in a polar ocean. The ISME journal. 2015; 9(8):1835–45. Epub 2015/02/24. doi: 10.1038/ismej.2015.1 PMID: 25700336; PubMed Central PMCID: PMC4511940. 13. Atamna-Ismaeel N, Sabehi G, Sharon I, Witzel KP, Labrenz M, Jurgens K, et al. Widespread distribution of proteorhodopsins in freshwater and brackish ecosystems. The ISME journal. 2008; 2(6):656–62. Epub 2008/03/29. doi: 10.1038/ismej.2008.27 PMID: 18369329. 14. Sharma AK, Zhaxybayeva O, Papke RT, Doolittle WF. Actinorhodopsins: proteorhodopsin-like gene sequences found predominantly in non-marine environments. Environmental microbiology. 2008; 10 (4):1039–56. Epub 2008/01/26. doi: 10.1111/j.1462-2920.2007.01525.x PMID: 18218036. 15. Sharma AK, Sommerfeld K, Bullerjahn GS, Matteson AR, Wilhelm SW, Jezbera J, et al. Actinorhodop- sin genes discovered in diverse freshwater habitats and among cultivated freshwater Actinobacteria. The ISME journal. 2009; 3(6):726–37. Epub 2009/02/27. doi: 10.1038/ismej.2009.13 PMID: 19242530. 16. Bohorquez LC, Ruiz-Perez CA, Zambrano MM. Proteorhodopsin-like genes present in thermoacidophi- lic high-mountain microbial communities. Applied and environmental microbiology. 2012; 78(21):7813– 7. Epub 2012/09/04. doi: 10.1128/AEM.01683-12 PMID: 22941077; PubMed Central PMCID: PMC3485705. 17. Zhu W, Lan Y, Lou X, Han N, Ran T, Xu L, et al. Isolation of proteorhodopsin-bearing bacterium JL-3 from fresh water and characterization of the proteorhodopsin. FEMS microbiology letters. 2013; 344 (1):10–7. Epub 2013/04/05. doi: 10.1111/1574-6968.12144 PMID: 23551202. 18. Gorriti MF, Dias GM, Chimetto LA, Trindade-Silva AE, Silva BS, Mesquita MM, et al. Genomic and phe- notypic attributes of novel salinivibrios from stromatolites, sediment and water from a high altitude lake. BMC genomics. 2014; 15:473. Epub 2014/06/15. doi: 10.1186/1471-2164-15-473 PMID: 24927949; PubMed Central PMCID: PMC4094778. 19. Koh EY, Atamna-Ismaeel N, Martin A, Cowie RO, Beja O, Davy SK, et al. Proteorhodopsin-bearing bacteria in Antarctic sea ice. Applied and environmental microbiology. 2010; 76(17):5918–25. Epub 2010/07/06. doi: 10.1128/AEM.00562-10 PMID: 20601510; PubMed Central PMCID: PMC2935036. 20. Petrovskaya LE, Lukashev EP, Chupin VV, Sychev SV, Lyukmanova EN, Kryukova EA, et al. Predicted bacteriorhodopsin from Exiguobacterium sibiricum is a functional proton pump. FEBS letters. 2010; 584(19):4193–6. Epub 2010/09/14. doi: 10.1016/j.febslet.2010.09.005 PMID: 20831870. 21. Gushchin I, Chervakov P, Kuzmichev P, Popov AN, Round E, Borshchevskiy V, et al. Structural insights into the proton pumping by unusual proteorhodopsin from nonmarine bacteria. Proceedings of the National Academy of Sciences of the United States of America. 2013; 110(31):12631–6. Epub 2013/07/ 23. doi: 10.1073/pnas.1221629110 PMID: 23872846; PubMed Central PMCID: PMC3732926. 22. Atamna-Ismaeel N, Finkel OM, Glaser F, Sharon I, Schneider R, Post AF, et al. Microbial rhodopsins on leaf surfaces of terrestrial plants. Environmental microbiology. 2012; 14(1):140–6. Epub 2011/09/ 03. doi: 10.1111/j.1462-2920.2011.02554.x PMID: 21883799; PubMed Central PMCID: PMC3608849. 23. Yoshizawa S, Kawanabe A, Ito H, Kandori H, Kogure K. Diversity and functional analysis of proteorhodopsin in marine Flavobacteria. Environmental microbiology. 2012; 14(5):1240–8. Epub 2012/02/15. doi: 10.1111/j.1462-2920.2012.02702.x PMID: 22329552. 24. Kwon SK, Kim BK, Song JY, Kwak MJ, Lee CH, Yoon JH, et al. Genomic makeup of the marine flavo- bacterium Nonlabens (Donghaeana) dokdonensis and identification of a novel class of rhodopsins. Genome biology and evolution. 2013; 5(1):187–99. Epub 2013/01/08. doi: 10.1093/gbe/evs134 PMID: 23292138; PubMed Central PMCID: PMC3595038. 25. Lami R, Cottrell MT, Campbell BJ, Kirchman DL. Light-dependent growth and proteorhodopsin expression by Flavobacteria and SAR11 in experiments with Delaware coastal waters. Environmental microbiology. 2009; 11(12):3201–9. Epub 2009/08/21. doi: 10.1111/j.1462-2920.2009.02028.x PMID: 19691506. 26. Gonzalez JM, Fernandez-Gomez B, Fernandez-Guerra A, Gomez-Consarnau L, Sanchez O, Coll- Llado M, et al. Genome analysis of the proteorhodopsin-containing marine bacterium Polaribacter sp. MED152 (Flavobacteria). Proceedings of the National Academy of Sciences of the United States of America. 2008; 105(25):8724–9. Epub 2008/06/17. doi: 10.1073/pnas.0712027105 PMID: 18552178; PubMed Central PMCID: PMC2438413.

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 15 / 18 proteorhodopsin from Modern Stromatolites

27. de la Torre JR, Christianson LM, Beja O, Suzuki MT, Karl DM, Heidelberg J, et al. Proteorhodopsin genes are distributed among divergent marine bacterial taxa. Proceedings of the National Academy of Sciences of the United States of America. 2003; 100(22):12830–5. Epub 2003/10/21. doi: 10.1073/ pnas.2133554100 PMID: 14566056; PubMed Central PMCID: PMC240704. 28. McCarren J, DeLong EF. Proteorhodopsin photosystem gene clusters exhibit co-evolutionary trends and shared ancestry among diverse marine microbial phyla. Environmental microbiology. 2007; 9 (4):846–58. Epub 2007/03/16. doi: 10.1111/j.1462-2920.2006.01203.x PMID: 17359257. 29. Ran T, Ozorowski G, Gao Y, Sineshchekov OA, Wang W, Spudich JL, et al. Cross-protomer interaction with the photoactive site in oligomeric proteorhodopsin complexes. Acta crystallographica Section D, Biological crystallography. 2013; 69(Pt 10):1965–80. Epub 2013/10/09. doi: 10.1107/ S0907444913017575 PMID: 24100316. 30. Luecke H, Richter HT, Lanyi JK. Proton transfer pathways in bacteriorhodopsin at 2.3 angstrom resolution. Science. 1998; 280(5371):1934–7. Epub 1998/06/25. PMID: 9632391. 31. Spudich JL, Yang CS, Jung KH, Spudich EN. Retinylidene proteins: structures and functions from archaea to humans. Annual review of cell and developmental biology. 2000; 16:365–92. Epub 2000/10/ 14. doi: 10.1146/annurev.cellbio.16.1.365 PMID: 11031241. 32. Man D, Wang W, Sabehi G, Aravind L, Post AF, Massana R, et al. Diversification and spectral tuning in marine proteorhodopsins. The EMBO journal. 2003; 22(8):1725–31. Epub 2003/04/19. doi: 10.1093/ emboj/cdg183 PMID: 12682005; PubMed Central PMCID: PMC154475. 33. Wang WW, Sineshchekov OA, Spudich EN, Spudich JL. Spectroscopic and photochemical characterization of a deep ocean proteorhodopsin. The Journal of biological chemistry. 2003; 278(36):33985–91. Epub 2003/06/25. doi: 10.1074/jbc.M305716200 PMID: 12821661. 34. Sabehi G, Kirkup BC, Rozenberg M, Stambler N, Polz MF, Beja O. Adaptation and spectral tuning in divergent marine proteorhodopsins from the eastern Mediterranean and the Sargasso Seas. The ISME journal. 2007; 1(1):48–55. Epub 2007/11/29. doi: 10.1038/ismej.2007.10 PMID: 18043613. 35. Jung JY, Choi AR, Lee YK, Lee HK, Jung KH. Spectroscopic and photochemical analysis of proteorhodopsin variants from the surface of the Arctic Ocean. FEBS letters. 2008; 582(12):1679–84. Epub 2008/04/26. doi: 10.1016/j.febslet.2008.04.025 PMID: 18435930. 36. Guo Z, Zhang H, Lin S. Light-promoted rhodopsin expression and starvation survival in the marine dinoflagellate Oxyrrhis marina. PloS one. 2014; 9(12):e114941. Epub 2014/12/17. doi: 10.1371/journal. pone.0114941 PMID: 25506945; PubMed Central PMCID: PMC4266641. 37. DeLong EF, Beja O. The light-driven proton pump proteorhodopsin enhances bacterial survival during tough times. PLoS biology. 2010; 8(4):e1000359. Epub 2010/05/04. doi: 10.1371/journal.pbio.1000359 PMID: 20436957; PubMed Central PMCID: PMC2860490. 38. Albarracin VH, Gartner W, Farias ME. Forged Under the Sun: Life and Art of Extremophiles from Andean Lakes. Photochemistry and photobiology. 2016; 92(1):14–28. Epub 2015/12/10. doi: 10.1111/ php.12555 PMID: 26647770. 39. Albarracin VH, Kurth D, Ordonez OF, Belfiore C, Luccini E, Salum GM, et al. High-Up: A Remote Reser- voir of Microbial Extremophiles in Central Andean Wetlands. Frontiers in microbiology. 2015; 6:1404. Epub 2016/01/07. doi: 10.3389/fmicb.2015.01404 PMID: 26733008; PubMed Central PMCID: PMC4679917. 40. Kurth D, Belfiore C, Gorriti MF, Cortez N, Farias ME, Albarracin VH. Genomic and proteomic evidences unravel the UV-resistome of the poly-extremophile Acinetobacter sp. Ver3. Frontiers in microbiology. 2015; 6:328. Epub 2015/05/09. doi: 10.3389/fmicb.2015.00328 PMID: 25954258; PubMed Central PMCID: PMC4406064. 41. Farias ME, Rascovan N, Toneatti DM, Albarracin VH, Flores MR, Poire DG, et al. The discovery of stromatolites developing at 3570 m above sea level in a high-altitude volcanic lake Socompa, Argentinean Andes. PloS one. 2013; 8(1):e53497. Epub 2013/01/12. doi: 10.1371/journal.pone.0053497 PMID: 23308236; PubMed Central PMCID: PMC3538587. 42. Bequer Urbano S, Albarracin VH, Ordonez OF, Farias ME, Alvarez HM. Lipid storage in high-altitude Andean Lakes extremophiles and its mobilization under stress conditions in Rhodococcus sp. A5, a UV-resistant actinobacterium. Extremophiles: life under extreme conditions. 2013; 17(2):217–27. Epub 2013/01/04. doi: 10.1007/s00792-012-0508-2 PMID: 23283521. 43. Farias ME, Contreras M, Rasuk MC, Kurth D, Flores MR, Poire DG, 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: life under extreme conditions. 2014; 18(2):311–29. Epub 2014/01/21. doi: 10.1007/s00792-013-0617-6 PMID: 24442191. 44. Rasuk MC, Kurth D, Flores MR, Contreras M, Novoa F, Poire D, et al. Microbial characterization of microbial ecosystems associated to evaporites domes of gypsum in Salar de Llamara in Atacama

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 16 / 18 proteorhodopsin from Modern Stromatolites

desert. Microbial ecology. 2014; 68(3):483–94. Epub 2014/05/27. doi: 10.1007/s00248-014-0431-4 PMID: 24859438. 45. Ordoñez OF, Lanzarotti EO, Kurth DG, Cortez N, Farias ME, Turjanski AG. Genome comparison of two Exiguobacterium strains from high altitude andean lakes with different arsenic resistance: Identification and 3D modeling of the Acr3 efflux pump. Frontiers in Environmental Science. 2015; 3:50. 46. Rasuk MC, Fernández AB, Kurth D, Contreras M, Novoa F, Poiré D, et al. Bacterial diversity in microbial mats and sediments from the Atacama Desert. Microbial ecology. 2016; 71(1):44–56. doi: 10.1007/ s00248-015-0649-9 PMID: 26224164 47. Rascovan N, Maldonado J, Vazquez MP, Eugenia Farias M. Metagenomic study of red biofilms from Diamante Lake reveals ancient arsenic bioenergetics in haloarchaea. The ISME journal. 2015. doi: 10. 1038/ismej.2015.109 48. Ordonez OF, Lanzarotti E, Kurth D, Gorriti MF, Revale S, Cortez N, et al. Draft Genome Sequence of the Polyextremophilic Exiguobacterium sp. Strain S17, Isolated from Hyperarsenic Lakes in the Argen- tinian Puna. Genome announcements. 2013; 1(4). Epub 2013/07/28. doi: 10.1128/genomeA.00480-13 PMID: 23887911; PubMed Central PMCID: PMC3735063. 49. Farias ME, Revale S, Mancini E, Ordonez O, Turjanski A, Cortez N, et al. Genome sequence of Sphin- gomonas sp. S17, isolated from an alkaline, hyperarsenic, and hypersaline volcano-associated lake at high altitude in the Argentinean Puna. Journal of bacteriology. 2011; 193(14):3686–7. Epub 2011/05/ 24. doi: 10.1128/JB.05225-11 PMID: 21602338; PubMed Central PMCID: PMC3133308. 50. Belfiore C, Ordonez OF, Farias ME. Proteomic approach of adaptive response to arsenic stress in Exi- guobacterium sp. S17, an extremophile strain isolated from a high-altitude Andean Lake stromatolite. Extremophiles: life under extreme conditions. 2013; 17(3):421–31. Epub 2013/03/26. doi: 10.1007/ s00792-013-0523-y PMID: 23525943. 51. Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC bioinformatics. 2004; 5:113. Epub 2004/08/21. doi: 10.1186/1471-2105-5-113 PMID: 15318951; PubMed Central PMCID: PMC517706. 52. Dereeper A, Audic S, Claverie JM, Blanc G. BLAST-EXPLORER helps you building datasets for phylogenetic analysis. BMC evolutionary biology. 2010; 10:8. Epub 2010/01/14. doi: 10.1186/1471-2148-10- 8 PMID: 20067610; PubMed Central PMCID: PMC2821324. 53. Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic acids research. 2008; 36(Web Server issue):W465–9. Epub 2008/04/22. doi: 10.1093/nar/gkn180 PMID: 18424797; PubMed Central PMCID: PMC2447785. 54. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic biology. 2003; 52(5):696–704. Epub 2003/10/08. PMID: 14530136. 55. Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annota- tions using subsystems technology. BMC genomics. 2008; 9:75. Epub 2008/02/12. doi: 10.1186/1471- 2164-9-75 PMID: 18261238; PubMed Central PMCID: PMC2265698. 56. Markowitz VM, Chen IM, Palaniappan K, Chu K, Szeto E, Pillay M, et al. IMG 4 version of the integrated microbial genomes comparative analysis system. Nucleic acids research. 2014; 42(Database issue): D560–7. Epub 2013/10/30. doi: 10.1093/nar/gkt963 PMID: 24165883; PubMed Central PMCID: PMC3965111. 57. Scharf B, Pevec B, Hess B, Engelhard M. Biochemical and photochemical properties of the photopho- bic receptors from Halobacterium halobium and Natronobacterium pharaonis. European journal of biochemistry / FEBS. 1992; 206(2):359–66. Epub 1992/06/01. PMID: 1597180. 58. Balashov SP, Petrovskaya LE, Lukashev EP, Imasheva ES, Dioumaev AK, Wang JM, et al. Aspartate- histidine interaction in the retinal schiff base counterion of the light-driven proton pump of Exiguobacter- ium sibiricum. Biochemistry. 2012; 51(29):5748–62. Epub 2012/06/29. doi: 10.1021/bi300409m PMID: 22738070; PubMed Central PMCID: PMC3415241. 59. Balashov SP, Petrovskaya LE, Imasheva ES, Lukashev EP, Dioumaev AK, Wang JM, et al. Breaking the carboxyl rule: lysine 96 facilitates reprotonation of the Schiff base in the photocycle of a retinal protein from Exiguobacterium sibiricum. The Journal of biological chemistry. 2013; 288(29):21254–65. Epub 2013/05/23. doi: 10.1074/jbc.M113.465138 PMID: 23696649; PubMed Central PMCID: PMC3774394. 60. Dioumaev AK, Petrovskaya LE, Wang JM, Balashov SP, Dolgikh DA, Kirpichnikov MP, et al. Photo- cycle of Exiguobacterium sibiricum rhodopsin characterized by low-temperature trapping in the IR and time-resolved studies in the visible. The journal of physical chemistry B. 2013; 117(24):7235–53. Epub 2013/05/31. doi: 10.1021/jp402430w PMID: 23718558; PubMed Central PMCID: PMC3813598. 61. Bamberg E, Apell H-J, Dencher N, Sperling W, Stieve H, Läuger P. Photocurrents generated by bacteriorhodopsin on planar bilayer membranes. Biophysics of structure and mechanism. 1979; 5(4):277–92.

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 17 / 18 proteorhodopsin from Modern Stromatolites

62. Albarracin VH, Pathak GP, Douki T, Cadet J, Borsarelli CD, Gartner W, et al. Extremophilic Acinetobac- ter strains from high-altitude lakes in Argentinean Puna: remarkable UV-B resistance and efficient DNA damage repair. Origins of life and evolution of the biosphere: the journal of the International Society for the Study of the Origin of Life. 2012; 42(2–3):201–21. Epub 2012/05/31. doi: 10.1007/s11084-012- 9276-3 PMID: 22644565. 63. Albarracin VH, Simon J, Pathak GP, Valle L, Douki T, Cadet J, et al. First characterisation of a CPD- class I photolyase from a UV-resistant extremophile isolated from High-Altitude Andean Lakes. Photo- chemical & photobiological sciences: Official journal of the European Photochemistry Association and the European Society for Photobiology. 2014; 13(5):739–50. Epub 2014/03/19. doi: 10.1039/ c3pp50399b PMID: 24637630. 64. Gomez-Consarnau L, Gonzalez JM, Coll-Llado M, Gourdon P, Pascher T, Neutze R, et al. Light stimu- lates growth of proteorhodopsin-containing marine Flavobacteria. Nature. 2007; 445(7124):210–3. Epub 2007/01/12. doi: 10.1038/nature05381 PMID: 17215843.

PLOS ONE | DOI:10.1371/journal.pone.0154962 May 17, 2016 18 / 18