In situ production of branched glycerol dialkyl glycerol tetraethers in a great basin hot spring (USA).

Branched glycerol dialkyl glycerol tetraethers (bGDGTs) are predominantly found in soils and peat bogs. In this study, we analyzed core (C)-bGDGTs after hydrolysis of polar fractions using liquid chromatography-atmospheric pressure chemical ionization-mass spectrometry and analyzed intact P-bGDGTs using total lipid extract (TLE) without hydrolysis by liquid chromatography-electrospray ionization-multiple stage mass spectrometry. Our results show multiple lines of evidence for the production of bGDGTs in sediments and cellulolytic enrichments in a hot spring (62-86°C) in the Great Basin (USA). First, in situ cellulolytic enrichment led to an increase in the relative abundance of hydrolysis-derived P-bGDGTs over their C-bGDGT counterparts. Second, the hydrolysis-derived P- and C-bGDGT profiles in the hot spring were different from those of the surrounding soil samples; in particular, a monoglycosidic bGDGT Ib containing 13,16-dimethyloctacosane and one cyclopentane moiety was detected in the TLE but it was undetectable in surrounding soil samples even after sample enrichments. Third, previously published 16S rRNA gene pyrotag analysis from the same lignocellulose samples demonstrated the enrichment of thermophiles, rather than mesophiles, and total bGDGT abundance in cellulolytic enrichments correlated with the relative abundance of 16S rRNA gene pyrotags from thermophilic bacteria in the phyla Bacteroidetes, Dictyoglomi, EM3, and OP9 ("Atribacteria"). These observations conclusively demonstrate the production of bGDGTs in this hot spring; however, the identity of organisms that produce bGDGTs in the geothermal environment remains unclear.

Introduction

Recent advances in liquid chromatography-mass spectrometry (LC-MS) have significantly expanded our view of the occurrence of unique lipid biomarkers in the natural environment (Hopmans et al., 2000; Schouten et al., 2000, 2013; Sturt et al., 2004; Liu et al., 2012). Branched glycerol dialkyl glycerol tetraethers (bGDGTs; Figure A1) are unusual lipids that have been detected in a variety of natural settings using LC-MS, including soil and peat bogs (Hopmans et al., 2004; Weijers et al., 2006, 2007a; Peterse et al., 2009, 2010, 2012; Liu et al., 2010), lakes (Sinninghe Damsté et al., 2009; Tierney and Russell, 2009; Tierney et al., 2011; Sun et al., 2011; Wang et al., 2012), rivers and estuaries (Kim et al., 2010, 2012; Zhu et al., 2011; Zhang et al., 2012; Yang et al., 2013), and continental margin sediments (Weijers et al., 2007b,c; Schouten et al., 2008; Rueda et al., 2009; Bendle et al., 2010). While bGDGTs show clear structural similarities to known archaeal membrane lipids (Schouten et al., 2000), the source of bGDGTs was uncertain until recent studies demonstrated production of a single bGDGT, bGDGT I (see Figure A1), by two members of the phylum Acidobacteria (Sinninghe Damsté et al., 2011). Although other bGDGTs may possibly be produced by mesophilic Acidobacteria growing optimally in acidic and anoxic environments (Weijers et al., 2006; Sinninghe Damsté et al., 2011), this has yet to be determined unequivocally.

Figure A1

Schematic structures of core branched glycerol dialkyl glycerol tetraethers (bGDGTs).

Despite their largely uncertain origin, bGDGT distributions tend to correlate with annual mean air temperature or soil pH (Weijers et al., 2007a; Peterse et al., 2009), thus having the potential to record paleocontinental temperatures (Weijers et al., 2007b,c; Ballantyne et al., 2010; Peterse et al., 2011, 2012; Weijers et al., 2011; Zhou et al., 2011) or soil pH (Weijers et al., 2007b; Tyler et al., 2010; Fawcett et al., 2011; Zhou et al., 2011). In combination with the archaeal biomarker crenarchaeol, bGDGTs are also used to estimate soil organic contribution in the marine environment (see review by Schouten et al., 2013).

bGDGTs have been reported from terrestrial hot springs in Yellowstone National Park, where the major source of bGDGTs was suggested to be soil runoff (Schouten et al., 2007). On the other hand, the presence of bGDGTs in mesophilic bacteria is recognized to be possibly a relict feature from thermophilic ancestors, as ether bonds as well as membrane-spanning core lipids have previously been reported in some thermophilic bacteria (Langworthy et al., 1983; DeRosa et al., 1988; Huber et al., 1992, 1996).

The Great Basin in western United States is an endorheic region with widely distributed geothermal activity. The hot springs of the Great Basin are characterized by low-inorganic energy yielding species such as ammonia, hydrogen sulfide, or hydrogen (Zhang et al., 2008), which are in contrast with more inorganic energy-rich geothermal systems fueled by subsurface volcanism (e.g., Yellowstone, Kamchatka, and Italy). The biological research in Great Basin hot springs has recently made important findings in lipid biomarker biogeochemistry and microbial carbon and nitrogen cycling processes (Pearson et al., 2004, 2008; Zhang et al., 2006, 2007; Huang et al., 2007; Costa et al., 2009; Dodsworth et al., 2011, 2013; Cole et al., 2013). Here we show multiple lines of evidence that bGDGTs are produced in situ in Great Boiling Spring (GBS) in the Great Basin.

Materials and methods

Sampling

GBS is a large geothermal spring located in the US Great Basin near the town of Gerlach, Nevada [N40°39.689′ W119°21.968′; 9.15 m deep, 7.6 m diameter; described in Costa et al. (2009) and Dodsworth et al. (2011)]. GBS has a relatively well-mixed, oxic water column and a relatively uniform clay bottom composed primarily of smectite, illite, kaolinite, quartz, and zeolite (Costa et al., 2009). Sediment samples (top ~1 cm of sediment/water interface) were collected at five locations (Sites A, B, C, D, and E; Figure 1) in February 2010; at each location sediment was homogenized on site in a sterile pie tin. Subsamples of the sediment homogenate were separated into a 50-mL polypropylene tube for lipid analysis. Other subsamples were collected for a variety of other analyses, including identification of predominant minerals and 16 S rRNA gene pyrosequencing. Temperature and pH were measured at the precise location of sampling prior to sample collection using a LaMotte pH 5 meter (LaMotte, Chestertown, MD). The current paper focuses on analysis of GDGTs. The details of the mineralogy, 16S rRNA gene pyrosequencing, and field chemistry were reported in detail previously (Cole et al., 2013).

Figure 1

Sites in GBS where sediments were collected and . The soil transect started at the edge of the hot spring at the lower right corner of the photo.

In addition, eight enrichments designed to stimulate growth of cellulolytic organisms were incubated in situ. Nylon bags (100 micron pore size, 10 × 10 cm) were filled with 20 g of either aspen shavings (AS) or ammonia fiber explosion (AFEX)-treated corn stover (CS). Bags were loaded into 20 × 12 × 5 cm polypropylene boxes punctured with ~100 0.5 cm holes to allow water exchange and incubated either suspended in spring water or buried ~1 cm deep in the sediment. The polypropylene boxes were anchored into the sediment or to structures adjacent to the spring by using stainless steel wire. Lignocellulose materials were incubated for 64 days (Figure 1, Site C; ~77°C) or 92 days (Figure 1, Site A; ~85°C). The difference in the incubation times was based on the time required to observe visual changes to the lignocellulose substrates that were consistent with cellulolysis. After incubation, the bags were removed, homogenized, and distributed into a sterile 50 mL polypropylene tube as described above for sediments. Details of lignocellulose degradation and 16S rRNA gene pyrosequencing were reported previously (Peacock et al., 2013).

The sample code for cellulolytic enrichments consists of three parameters, temperature (77°C or 85°C), cellulose substrate (A, aspen shavings; C, corn stover), and incubation environment (W, suspended in water; S, buried in sediment) as described earlier (Peacock et al., 2013). As described in Peacock et al. (2013), the temperature at each incubation site (within 0.5 meters of each lignocellulose enrichment) was tracked for the majority of the duration of the incubation using high-temperature iButtons (Maxim Integrated, San Jose, CA). During the time the temperature was tracked, temperatures ranged from 68°C to 82°C at site C (mean 77°C) and 74–88°C at site A (mean 85°C).

Finally, soil samples were collected at 10-, 20-, 30-, 50-, 150-, 200-, 300- and 500 cm distance from the edge (zero cm) of the spring (Figure 1) in order to provide a contrast in bGDGT profiles between the soil and the hot spring, which serves to evaluate possible soil contamination into the hot spring. This is a potential concern because of previous report of bGDGTs in soil next to hot springs in the Great Basin (Peterse et al., 2009). The soil temperature was determined by inserting a stainless steel temperature probe ~3 cm into the soil where the sample was collected. The soil pH was determined in the lab following a previously described procedure (Zhang et al., 2012).

All samples were frozen on dry ice in the field and stored at −80°C before analysis. DNA from each sample was extracted by using the FastDNA Spin Kit for Soil (MP Biomedicals, Solon, OH); raw data were presented in Cole et al. (2013) and Peacock et al. (2013).

LC-MS analysis of polar bGDGTs (hydrolysis method)

Lipid extraction, fractionation, and separation of core (C)- and hydrolysis-derived polar (P)- bGDGTs followed a sonication method described in Zhang et al. (2012), in which the P-bGDGTs were calculated as the difference between the hydrolyzed and non-hydrolyzed polar fractions. The GDGTs were analyzed on an Agilent 1200 liquid chromatography equipped with an automatic injector coupled to QQQ 6460 MS and Mass Hunter LC-MS manager software using a procedure modified from Hopmans et al. (2004). Detection was performed using the Agilent 6460 triple-quadrupole spectrometer MS with an atmospheric pressure chemical ionization (APCI) ion source (Zhang et al., 2012). Separation of peaks was achieved using a Prevail Cyano column (2.1 mm×150 mm, 3 μm; Alltech Deerfield, IL, USA) maintained at a temperature of 40°C (Zhang et al., 2012). The detection limit of the LC-MS was 0.8 pg (Zhang et al., 2012).

LC-MS analysis of intact polar bGDGTs (non-hydrolysis method)

While the hydrolysis method gives total abundance of all polar bGDGTs, it does not identify the types of polar bGDGTs. To identify specific head groups of the intact polar lipid bGDGTs, total lipid extracts (TLEs) were also analyzed by a reverse phase liquid chromatography-electrospray ionization-multiple stage mass spectrometry (RP-ESI-MSn) at University of Bremen, Germany (Zhu et al., in review). In brief, the analysis of TLEs was performed on a Dionex Ultimate 3000 ultra-high pressure liquid chromatograph (UHPLC) coupled to a Bruker maXis Ultra High Resolution orthogonal accelerated quadrupole—time-of-flight (qTOF) tandem MS/MS, equipped with an electrospray ionization source (ESI) in positive ionization mode (Bruker Daltonik, Bremen, Germany). Ether lipids were eluted through an ACE3 C18 column (3 μm, 2.1 × 150 mm; Advanced Chromatography Technologies Ltd., Aberdeen, Scotland), starting with 100% eluent A isocratically for 10 min, followed by a gradient to 24% eluent B in 5 min, and then to 65% eluent B in 55 min at a flow rate of 0.2 mL/min, where the eluent A was 100:0.04:0.10 of methanol/formic acid/14.8 M NH3aq and B was 100:0.04:0.10 of 2-propanol/formic acid/14.8 M NH3(aq). The column was washed with 90% eluent B for 10 min and subsequently re-equilibrated with 100% A for another 10 min. Ether lipids were scanned from m/z 100 to 2000 in a positive mode at a scan rate of 1 Hz with automated data-dependent fragmentation of the three most abundant ions. To ensure mass accuracy, an internal lock mass (m/z 922.0077) and tuning mixture solution (m/z 322.0481, 622.0290, 922.0098, 121.9906, 1521.9715, and 1821.9523) were infused directly into the ion source throughout a complete run and at the near end of the run, respectively. Lipids were detected as protonated [M+H]+, ammoniated [M+NH4]+, and sodiated [M+Na]+ molecular ions and identified by retention time, accurate masses (better than 1 ppm), and diagnostic fragments (Weijers et al., 2006; Liu et al., 2010).

Statistical analyses

Mann-Whitney U tests and Wilcoxon signed rank tests were used as non-parametric alternatives to independent- and paired samples t-tests to explore relationships between bGDGT and experimental conditions. Linear regressions were calculated to quantify relationships between temperature and bGDGT fractions from sediment samples. These analyses were all calculated at the 0.05 level of significance.

Cluster analysis was performed on C-bGDGTs and hydrolysis-derived P-bGDGTs from soil and the hot spring samples using the base program in R 2.12.1. The relative abundances of C-bGDGTs and P-bGDGTs from all samples were imported into R and the Euclidean method was used to compute the distance matrix and generate a hierarchical clustering tree.

Spearman's rho, non-parametric correlation coefficients, were calculated to identify positive relationships between total bGDGTs (normalized to ng DNA) and relative abundance of phyla based on quality-filtered pyrotag sequence reads from cellulose enrichments. Analyses were completed for phyla that occurred at ≥1% relative abundance in one or more of the in situ cellulose enrichments. Subsequently, the same statistical framework was applied to individual Operational Taxonomic Units (OTUs) defined at 97% within phyla (Peacock et al., 2013) that were positively correlated with bGDGT abundance. Results are reported for 1-tailed significance.

Results and discussion

Abundance of C- and hydrolysis-derived P-bGDGTs

Sediment samples from the hot spring had C-bGDGTs ranging from 12 ng/g dry sediment to 280 ng/g dry sediment (Table 1), while P-bGDGTs were about 4-5-fold less abundant than C-bGDGTs (Table 1). Linear regression analyses indicated statistically significant, negative relationships between temperature and all bGDGT fractions (when normalized to gram dry sediment) from the hot spring sediments, with r2 values ranging from 0.80 to 0.82 (p < 0.05, Figure 2). The soil samples were also dominated by C-bGDGTs (up to 2–318 ng/g) with P-bGDGTs being less than 10 ng/g in five out of six samples (Table 1).

Table 1

Core- and hydrolysis-derived polar bGDGTs from sediments, .

			bGDGTs (ng/g dry wt)
Sample	Temp (°C)	pH	Core	Polar	% Polarb
SEDIMENTS
Site A	82.0	6.6	12	3	21
Site B	79.2	6.8	50	13	21
Site C	72.1	7.1	38	8	17
Site D	82.0	7.2	25	6	19
Site E	61.9	7.2	282	53	16
ENRICHMENTSa
77AS	78.8	6.8	22	57	72
77AW	78.8	6.8	6	9	59
77CS	78.8	6.8	63	88	58
77CW	78.8	6.8	4	11	72
85AS	85.5	6.9	3	6	64
85AW	85.5	6.9	3	2	45
85CS	85.5	6.9	28	12	29
85CW	85.5	6.9	7	3	34
SOILS
Soil-0 cm	NAd	8.2	NA	NA	NA
Soil-10 cm	33.6	7.8	149	8	5
Soil-20 cm	34.0	7.8	231	0	0
Soil-30 cm	30.9	7.4	319	7	2
Soil-50 cm	NA	NA	72	76	105
Soil-100 cm	31.7	8.2	NA	NA	NA
Soil-150 cm	NA	NA	2	0	6
Soil-200 cm	30.0	6.4	NA	NA	NA
Soil-300 cm	NA	NA	6	1	13
Soil-500 cm	30.5	9.2	NA	NA	NA
MAATc	10.7

Enrichment sample codes consist of three parameters: average temperature (77°C or 85°C), cellulose substrate (A, aspen shavings; C, core stover), and incubation location (W, suspended in water; S, buried in sediment); temperature and pH listed are those measured at the time of sample collection.

Hydrolysis-derived polar bGDGTs as a percentage of total (core + polar) bGDGTs.

Mean annual air temperature at Gerlach, Nevada (http://www.wunderground.com/weather-forecast/US/NV/Gerlach.html).

NA, not available.

Figure 2

Regression analyses indicated negative, linear relationships between the absolute abundance of C- and hydrolysis-derived P- bGDGTs in hot spring sediments and temperature (sig. < 0.05). bGDGT abundance was normalized to g dry mass of sample.

bGDGTs were not detected in cellulosic substrates before incubation (data not shown) and in situ cellulolytic enrichments had C- and P-bGDGTs ranging from 2.0 ng/g solids to 88 ng/g solids. The cellulolytic enrichments had significantly higher percentages of P-bGDGTs than the hot spring sediments (Table 1; Mann-Whitney U test, p = 0.003), indicating enrichment of bGDGT-producing bacteria among cellulolytic consortia. Within the cellulolytic enrichments lower temperature sites (~77°C) had higher total bGDGT concentrations (C- plus P- bGDGTs) as compared to their corresponding higher temperature sites (Wilcoxon signed rank test, p = 0.068), which is consistent with the temperature relationships observed in the hot spring sediments (Figure 2). These observations suggest that organisms producing the bGDGTs tend to have higher biomass at lower temperatures (e.g., 60°C) in GBS.

Comparisons of paired cellulolytic enrichments showed that incubations within the sediment contained a higher absolute abundance of bGDGTs than their corresponding water column incubations (Wilcoxon test, p = 0.068), suggesting that anaerobic conditions favored bGDGT-producing organisms in this hot spring. Genetic data suggest the enrichments buried in the sediments were anaerobic. For example, the majority of 16S rRNA gene pyrotags described in Peacock et al. (2013) were strict anaerobes (e.g., dominant groups are Archaeoglobales, Thermotogales, Thermofilaceae, etc.). Paired comparisons of cellulosic enrichment substrates revealed that the absolute abundance of bGDGTs was elevated in corn stover as compared to aspen shavings (Wilcoxon test, p = 0.068), indicating a potential preference of bGDGT organisms for corn stover.

Schouten et al. (2007) first reported C-bGDGTs in hot springs, which accounted for up to 64% of total GDGTs. However, absolute concentrations of bGDGTs were not reported, precluding any comparisons with our data. Peterse et al. (2009), on the other hand, reported C-bGDGTs from geothermally heated soil near two hot springs in the Great Basin, which had the majority of C-bGDGTs (0.9–760 ng/g dry wt) in the same range as what is reported here (Table 1). More recently, two studies have reported the presence of bGDGTs associated with hydrothermal deposits in the mid-ocean ridges (Hu et al., 2012; Lincoln et al., 2013). Collectively, these observations suggest that bGDGTs can be produced at elevated temperatures, which may indicate a possible origin of bGDGTs in thermophilic bacteria.

Composition of bGDGTs in sediments and enrichments

Hot spring sediment samples contained bGDGT I, Ib, Ic, II, IIb, IIc, and III (Table 2). In most sediment samples, bGDGT I was the dominant bGDGT, but bGDGT Ib, Ic, II, and III were also present in significant amounts (>15% relative abundance) in one or more sediment samples. Cellulose enrichment generally led to simplification of bGDGT profiles, with bGDGT I comprising up to 100% of total bGDGTs in some enrichments with aspen shavings, and bGDGT II as the only other lipid detected in enrichments with aspen shavings. Enrichments with corn stover were more complex. bGDGT I was the dominant lipid in most corn stover enrichments, but bGDGT Ib, Ic, II, IIc, and III were also present in significant amounts (>15%) in one or more corn stover enrichments. The dominance of bGDGT I is consistent with the high relative abundance of bGDGT I in other hot spring environments (Schouten et al., 2007), in soils of warmer climate (Weijers et al., 2007a), and in some geothermally heated soils, although bGDGT I and II, with various degrees of cyclization, were present in roughly equal amounts in others (Peterse et al., 2009). The different composition of bGDGTs in the cellulolytic enrichments compared with those of the hot spring sediment samples (Sites A–E, Figures 1, 3), along with the increase in P-bGDGTs over C-bGDGTs, suggests that the enrichments stimulated growth of a distinct population of b-GDGT-producing thermophiles. This is also supported by cluster analysis based on relative abundance of C- and P-bGDGTs in soil samples as the majority of them were distinct from those collected from the hot spring environment (Figure 3).

Table 2

Relative abundances of bGDGTs and bGDGTs-derived proxies for cellulolytic enrichments and natural samples from GBS.

Relative abundance of bGDGTs, normalized to 100% for core bGDGTs and for P-bGDGTs															bGDGT proxiesa
	I		Ib		Ic		II		IIb		IIc		III		MBT		CBT		TMBT/CBT (°C)		pHCBT
Sample	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar	Core	Polar
SEDIMENTS
Site A	56.5	21.5	17.9	24.6	9.0	14.3	9.4	16.4	2.6	2.7			4.6	20.5	0.83	0.60	0.51	0.14	30.9	22.8	7.4	8.4
Site B	50.5	34.3	16.3	15.7	4.3	13.3	16.8	23.3	6.3	0.3			5.7	13.1	0.71	0.63	0.47	0.56	25.1	20.4	7.5	7.3
Site C	72.8	39.1	10.4	12.0	2.6	10.7	8.6	18.0	3.0	8.6			2.7	11.5	0.86	0.62	0.78	0.44	29.4	20.7	6.7	7.6
Site D	46.6	47.4	18.4	17.9	5.8	15.6	17.2	2.2	5.3	11.0			6.8	5.8	0.71	0.81	0.43	0.24	25.3	32.2	7.6	8.1
Site E	72.3	68.1	14.2	9.8	10.5	7.6	1.7	8.6	0.2	2.2			1.2	3.7	0.97	0.86	0.71	0.81	35.7	29.1	6.9	6.6
ENRICHMENT SAMPLES
77AS	90.1	84.3					9.9	15.7
77AW	100.0	100.0
77CS	75.0	56.2	4.6				5.4	25.0	14.9					18.9	0.80	0.56	0.61		28.0	22.0	7.1	8.8
77CW	63.2	52.6					36.8	18.1						29.3
85AS	63.9	42.6					36.1	57.4
85AW	100.0	100.0
85CS	17.7	15.3	28.8	17.0	22.7	15.3	14.1	28.1	11.0	5.4	5.6	19.0			0.69	0.48	0.10	0.29	29.5	15.0	9.0	8.0
85CW	66.2	57.4					33.8	42.6

MBT = (I+Ib+Ic)/(I+Ib+ Ic+II + IIb+ IIc + III), in which IIIb and IIIc are all zero values and not included; CBT = −log[(Ib+IIb)/(I+II)]; TMBT/CBT(°C) = (MBT−0.122–0.187*CBT)/0.02; pHCBT = (3.33-CBT)/0.38.

All equations are derived from Weijers et al. (2007a).

Figure 3

Cluster analysis based on the relative abundance of C-bGDGTs (A) and hydrolysis-derived P-bGDGTs (B) from hot spring and soil samples, which distinguish the lipid profiles in hot spring sediments and cellulolytic enrichments from the majority of surrounding soil samples.

We also calculated the methylation index of branched tetraethers (MBT) and cyclization ratio of branched tetraethers (CBT), and the temperature and pH estimates derived from them, according to Weijers et al. (2007a) (Table 2). In the hot spring sediment samples, the calculated pH values averaged 7.2 ± 0.4 (n = 5) for C-bGDGTs and 7.6 ± 0.4 (n = 5) for P-bGDGTs; the former was significantly (P < 0.05) lower than the average pH values (8.0 ± 0.6, n = 6) calculated from soil C-bGDGTs (Table A1), which is close to the average measured pH value of the soil samples (7.8 ± 0.9, n = 6) (Table 1).

Table A1

Relative abundances of Core (C)- bGDGTs and C-bGDGTs-derived proxies for soil samples collected along the transect from GBS.

Sample	Relative abundance of C-bGDGTs (%)									C-bGDGT proxiesa
	I	Ib	Ic	II	IIb	IIc	III	IIIb	IIIc	MBT	CBT	pHCBT	TMBT/CBT (°C)
GBS-10cm	57.7	12.0	2.2	22.5	2.6	0.3	2.6	0.2	0.0	0.72	0.74	6.8	22.9
GBS-20cm	29.8	16.3	1.9	24.9	16.4	1.5	7.2	1.6	0.3	0.48	0.22	8.2	15.8
GBS-30cm	12.1	7.9	1.8	29.0	22.3	1.7	20.2	4.2	0.9	0.22	0.14	8.4	3.5
GBS-50cm	10.7	7.6	1.8	27.8	15.0	1.4	31.7	3.7	0.5	0.2	0.23	8.1	1.7
GBS-150cm	20.9	15.6	8.0	18.8	9.4	0.0	27.3	0.0	0.0	0.44	0.2	8.2	14.3
GBS-300cm	13.7	11.4	6.4	17.1	13.7	0.0	31.5	6.1	0.0	0.32	0.09	8.5	8.9

MBT = (I+Ib+ Ic)/(I+Ib+ Ic+II + IIb+ IIc + III + IIIb + IIIc).

CBT = −log[(Ib + IIb)/(I + II)]; TMBT/CBT (°C) = (MBT−0.122 − 0.187*CBT)/0.02; pHCBT = (3.33–CBT)/0.38. All equations are derived from Weijers et al. (2007a).

The calculated temperatures in sediment samples were statistically indistinguishable (p > 0.2) between the C-bGDGTs (29.3 ± 4.4°C, n = 5) and P-bGDGTs (25.0 ± 5.3°C, n = 5). The average C-bGDGTs-derived temperature was close to the average measured soil temperature (31.8 ± 1.7°C, n = 6) at the time of sampling but was significantly higher than the soil C-bGDGTs-derived temperatures (1.7–22.9°C, average = 11.2 ± 8.0°C, n = 6; Table A1), which were on the other hand closer to the mean annual air temperature (MAAT) for the region (10.7°C; Table 1). Schouten et al. (2007) showed that one of the Yellowstone hot springs had MBT and CBT temperatures very close to the local MAAT and suggested that the bGDGTs derived from the surrounding soil area. While we cannot unequivocally exclude the possibility of soil contamination in the hot spring sediment samples in GBS, the MBT-CBT temperatures calculated from the enrichments conducted at high temperature are also in the range of 15–32°C (Table 2), which cannot be explained by soil contamination. Overall, these results suggest that the MBT-CBT proxies derived from soil environments may not be applicable in geothermal environments and that factors other than pH or temperature may control the relative abundance of bGDGTs in geothermal systems.

Intact polar lipid bGDGTs analyzed by RP-ESI-MSn

Intact glycosidic bGDGTs without hydrolysis have been detected in a peat bog where bGDGT-producing bacteria thrive (Liu et al., 2010). In the GBS, activity of bGDGT-producing bacteria was evidenced by the occurrence of monoglycosidic bGDGT Ib detected using total lipid extract (Figure A2), which accounted for ca. 6.5% of total C-bGDGTs (in terms of total peak areas of [M+Na]+, [M+H]+, [M+NH4]+). In contrast, intact polar bGDGTs were undetectable in the adjacent soil samples by RP-ESI-MSn, even after purification of soil extracts by preparative HPLC (data not shown), which is consistent with the low P-bGDGTs by the hydrolysis method (above). However, the reverse phase LC-MS method failed to detect other polar bGDGTs obtained by the hydrolysis method (Table 2).

Figure A2

Occurrence of monoglycosidic bGDGT Ib (1G-bGDGT Ib) at site C in Figure 1, which was determined by RP-ESI-MS The measured accurate masses of [M+H]+, [M+NH4]+, [M+Na]+ adduct ions. (B) Extracted ion chromatogram of 1199.0739 (the theoretical mass of m/z [M+NH4]+) ± 0.05. (C) MS2 spectra of [M+NH4]+ adduct ion (Note the ring location is based on Weijers et al., 2006).

16S rRNA gene pyrosequences

Previously published pyrosequencing data (Peacock et al., 2013) further supported the in situ production of bGDGTs in the hot spring environment. Peacock et al. (2013) demonstrated that each of the cellulolytic substrates led to enrichment of thermophilic, cellulolytic consortia, as documented by (1) significant increase in DNA yield over un-incubated cellulose substrates, (2) changes in lignocelluloses substrate composition consistent with cellulolysis, and (3) dramatic changes in microbial community composition, with dominant community members closely related to known cellulolytic and hemicellulolytic Thermotoga and Dictyoglomus, cellulolytic and sugar-fermenting Desulfurococcales, and sugar-fermenting and hydrogenotrophic Archaeoglobales.

The relative abundance of several bacterial phyla in the pyrotag datasets was significantly correlated with bGDGT abundance: Bacteroidetes, Dictyoglomi, and candidate phyla EM3 and OP9 (“Atribacteria” Dodsworth et al., 2013) (Figure 4; Table A2). Within the Bacteroidetes, only a single species-level OTU was positively correlated with bGDGT abundance (Table A3); however, that pyrotag could not be assigned to a class or any lower taxonomic level. Similarly, EM3 and OP9 were represented by one and two species-level OTUs that correlated with bGDGT abundance, respectively, but no cultures are available for those candidate phyla (Table A3). Although these groups are plausible sources of bGDGTs, definitive experiments to test whether they are sources of bGDGTs may await isolation and study of representative strains. The phylum Dictyoglomi was represented by two OTUs that correlated with bGDGT abundance that could be ascribed to the genus Dictyoglomus (Table A3). That genus is currently comprised of only two closely related species, Dictyoglomus thermophilum and D. turdigum, both of which are obligate fermenters that are known to decompose components of hemicellulose and cellulose (Saiki et al., 1985; Patel et al., 1987). Dictyoglomus spp. Rt46-B1 is known to produce predominantly C16:0 phospholipid derived fatty acids and traces (<10%) of C14:0, C15:0, C17:0, C18:0, C20:0, aC17:0, C16:1ω11, and C18:1ω9 (Patel et al., 1991). Our own analysis of lipid extracts from pure cultures of Dictyoglomus thermophilum DSM 3960T and D. turdigum DSM 6724T, both grown on DSMZ 388 medium with gentle agitation in the dark at 70°C for 2 days, failed to reveal bGDGTs (data not shown). This argues against Dictyoglomi as a source of bGDGTs, although it is possible that such organisms may produce bGDGTs in the natural environment but not under the laboratory conditions. Lastly, Acidobacteria that are commonly believed to be sources of soil/peat bog bGDGTs (Weijers et al., 2006) constituted less than 0.01% of total pyrotags (Peacock et al., 2013) and did not show a significant correlation with bGDGTs in the enrichment, even though the primers used for PCR each match >95% of Acidobacteria 16S rRNA gene sequences.

Figure 4

Spearman's rho correlation coefficients indicated statistically significant, positive relationships between bGDGT absolute abundance and the bacterial Bacteroidetes, Dictyoglomi, EM3, and OP9 (“Atribacteria;” Dodsworth et al., . Total bGDGT abundance, including core- and polar lipids, were normalized to ng DNA.

Table A2

Correlations between normalized bGDGT concentrations (pg lipids/ng DNA) and relative abundance of bacterial phyla.

aStatistically significant, positive relationships highlighted in pink (1-tailed sig. ≤ 0.05). No archaeal phyla showed significant relationships with bGDGTs.

bRelative abundance of phyla-based pyrotag sequence reads.

Table A3

Correlations between normalized bGDGT concentrations (pg lipids/ng DNA) and relative abundance of dominant bacterial OTUs.

aStatistically significant, positive relationships highlighted in pink (1-tailed sig. ≤ 0.05). No archaeal phyla showed significant relationships with bGDGTs.

bRelative abundance of OTUs-based pyrotag sequence reads (See Peacock et al., in review for more detail on OTU assignments).

In summary, we show three lines of evidence supporting the in situ production of bGDGTs at high temperature in a Great Basin hot spring: (1) the greater abundance of hydrolysis-derived P-bGDGTs over C-bGDGTs in enrichments vs. the opposite in the soil samples and the distinct composition of bGDGTs in the enrichments, (2) the presence of the intact monoglycosidic-bGDGT Ib identified by the RP-ESI-MSn in the hot spring and its absence in the surrounding soil, and (3) the significant correlations between bGDGTs and certain groups of thermophilic bacteria combined with the low abundance of Acidobacteria, as inferred from pyrotag data, in all samples.

Conflict of interest statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.