Evaluating Production of Cyclopentyl Tetraethers by Marine Group II Euryarchaeota in the Pearl River Estuary and Coastal South China Sea: Potential Impact on the TEX86 Paleothermometer

TEX86 [TetraEther indeX of glycerol dialkyl glycerol tetraethers (GDGTs) with 86 carbon atoms] has been widely applied to reconstruct (paleo-) sea surface temperature. Marine Group I (MG-I) Thaumarchaeota were thought to be the primary source of GDGTs constituting the TEX86 formula; however, recent research has suggested that Marine Group II (MG-II) Euryarchaeota may also contribute significantly to the GDGT pool in the ocean. Little is known regarding the potential impact of MG-II Euryarchaeota-derived GDGTs on TEX86 values recorded in marine sediments. In this study, we assessed the relationship between distributions of GDGTs and MG-II Euryarchaeota and evaluated its potential effect on the TEX86 proxy. Lipid and DNA analyses were performed on suspended particulate matter and surface sediments collected along a salinity gradient from the lower Pearl River (river water) and its estuary (mixing water) to the coastal South China Sea (SCS, seawater). TEX86-derived temperatures from the water column and surface sediments were significantly correlated and both were lower than satellite-based temperatures. The ring index (RI) values in these environments were higher than predicted from the calculated TEX86-RI correlation, indicating that the GDGT pool in the water column of the PR estuary and coastal SCS comprises relatively more cyclopentane rings, which thereby altered TEX86 values. Furthermore, the abundance of MG-II Euryarchaeota 16S rRNA gene in the mixing water was two to three orders of magnitude higher than those observed in the river or seawater. Significant linear correlations were observed between the gene abundance ratio of MG-II Euryarchaeota to total archaea and the fractional abundance of GDGTs with cyclopentane rings. Collectively, these results suggest that MG-II Euryarchaeota likely produce a large proportion of GDGTs with 1–4 cyclopentane moieties, which may bias TEX86 values in the water column and sediments. As such, valid interpretation of TEX86 values in the sediment record, particularly in coastal oceans, should consider the contribution from MG-II Euryarchaeota.


INTRODUCTION
TEX 86 is a popular temperature proxy applied in paleoclimatological studies, which is based on the relative distribution of cyclopental rings among isoprenoid glycerol dialkyl glycerol tetraether (GDGT; Figure S1) lipids produced by archaea in marine and terrestrial environments (see review by Schouten et al., 2013). Global core-top calibrations of TEX 86 values were empirically correlated with the annual mean sea surface temperature (SST; Schouten et al., 2002;Kim et al., 2008Kim et al., , 2010. However, mounting evidence indicates anomalies of TEX 86 -derived SST in coastal seas and the open ocean, which have been attributed to multiple inputs of GDGTs from terrestrial (e.g., Weijers et al., 2006) or bathypelagic sources (e.g., Lee et al., 2008), as well as production in marine sediments (e.g., Liu X. L. et al., 2011).
A great deal of effort has been made to assess TEX 86 accuracy in marine and lake sediments. For example, application of the TEX 86 proxy is cautioned under the following circumstances: a branched and isoprenoid tetraether (BIT) index > 0.2 (Zhu et al., 2011), a ratio of GDGT-2/crenarchaeol > 0.4 (Weijers et al., 2011a), a Methane Index > 0.5 , a ratio of GDGT-0/crenarchaeol > 2 ( Blaga et al., 2009), or when %GDGT-2 > 45 (Sinninghe Damsté et al., 2012). Recently, based on an assessment of the relationship between the weighted average number of cyclopentane rings in GDGTs (ring index, RI) and published TEX 86 data from core-top sediments, Zhang et al. (2016) established a significant correlation between TEX 86 and RI [RI = 3.32 × (TEX 86 ) 2 − 0.77 × TEX 86 + 1.59; ±2σ ∼ 0.3]. This relationship was expected and reflects the physiological response of marine archaea to synthesize GDGTs with more rings (higher RI values) at higher temperatures (higher TEX 86 values). Deviations from this relationship suggest that temperature is not a dominant factor governing GDGT ring distribution, or, alternatively, the relationship between GDGT ring distribution and temperature is different from the modern analog as defined by the global core-top dataset.
The TEX 86 -related GDGTs (GDGTs-1, -2, -3 and the crenarchaeol regioisomer) in the water column are thought to primarily derive from the Marine Group I (MG-I) Thaumarchaeota (e.g., Schouten et al., 2008;Pitcher et al., 2011a,b), as it is one of the dominant groups of planktonic archaea in the ocean (Karner et al., 2001). In particular, crenarchaeol, containing one cyclohexane and four cyclopentane moieties, is accepted as a specific biomarker for MG-I Thaumarchaeota Schouten et al., 2008). Marine Group II (MG-II) Euryarchaeota are another group of planktonic archaea that predominantly inhabit coastal water and (near) surface waters of the open ocean (e.g., DeLong, 1992;Galand et al., 2010;Hugoni et al., 2013). Recently, this cosmopolitan group was also invoked as another major source of GDGTs (including crenarchaeol) in the ocean (Lincoln et al., 2014a), which supports an earlier hypothesis about GDGT-producing MG-II Euryarchaeota (Turich et al., 2007). However, concrete evidence of crenarchaeol production by MG-II is lacking due to the inability to obtain a pure culture or isolate, and the exact composition of GDGTs produced by these organisms remains uncertain (Lincoln et al., 2014b;Schouten et al., 2014).
To further evaluate the relationship between MG-II Euryarchaeota and the distribution of GDGTs, we quantified the abundance of MG-II 16S rRNA gene, determined the distribution of GDGT core and intact polar lipids (CL and IPL, respectively), and assessed TEX 86 and Ring Index (RI) values from suspended particulate matter and surface sediments collected from the lower Pearl River (PR) and its estuary to the coastal South China Sea (SCS). Our results provide a mechanistic explanation for deviations of the TEX 86 paleothermometer and have important implications for the sources of GDGTs in marine environments and probing past changes in global climate.

Sample Collection
Sampling locations and other information for suspended particulate matter (SPM) and surface sediments are shown in Figure 1 and Table 1. SPM samples (n = 18) and surface sediments (n = 8) were collected along a salinity gradient from the lower Pearl River and its estuary to the coastal South China Sea in the summer of 2011. SPM samples were collected from the surface (stations R1-R6) and the bottom (stations R1 and R2) water in the lower Pearl River, from three water depths (surface, middle, and bottom) and during three tidal periods (high tide, slack tide, and low tide) at station M located in the PR estuary, and at four water depths (surface, subsurface, middle, and bottom) at station S in the seawater of the coastal SCS (Figure 1). The depth of the sampling layers in the water column is given in Table 1. About 4-103 liters of water were filtered onto combusted (450 • C, overnight) glass-fiber filters (Whatman GF/F, 0.7 µm, 142 mm diameter) using an in situ submersible pump. The pH, temperature, salinity, and depth were determined in situ by a Horiba instrument (W-20XD, Kyoto, Japan; Table 1). Surface sediments (top ca. 10 cm) were collected at all stations ( Figure 1; Table 1) using a grab sampler. All samples were frozen immediately in liquid nitrogen and kept at −80 • C in the laboratory before analysis.

GDGT Extraction and Separation
The SPM samples (n = 18) and surface sediments (n = 8) were freeze-dried and extracted using a modified Bligh and Dyer method (Bligh and Dyer, 1959); the separation of core lipids and intact polar lipids followed the procedure described in Weijers et al. (2011b). Briefly, the total lipid extract (TLE) was obtained by ultrasonic extraction (10 min each, 6 times) of SPM (1 filter) or sediments (5 g) with a single-phase solvent mixture of methanol, dichloromethane (DCM), and phosphate buffer (2:1:0.8, v/v/v; pH 7.4). The TLE was separated over an activated silica gel column eluted with n-hexane/ethyl acetate (1:1, v/v) and methanol for CL and IPL, respectively. For GDGT quantification, a known amount of an internal C46 GDGT standard was added into the CL fraction or IPL fraction (Huguet et al., 2006). The CL fraction was directly measured. The IPL fraction was determined by measuring the CL after cleavage of the polar head groups via acid or base hydrolysis (Pitcher FIGURE 1 | Map showing locations of sites (dark circles) in the lower Pearl River (PR), the PR estuary, and coastal South China Sea (SCS). R, River water, which is followed by the station #1-6; M, mixing water; S, seawater. Station S includes four water layers (surface, subsurface, middle, and bottom). Station M includes three water layers (surface, middle, and bottom). The surface water at station M was also collected during the high tide-, slack tide-, and low tide-periods. Stations R1 and R2 include two water layers (surface and bottom). Surface sediments were collected at each sampling site. Weijers et al., 2011b). Briefly, 1/3 IPL fraction (nonhydrolyzed IPL fraction) was directly condensed; another 1/3 IPL fraction was hydrolyzed (2 h) in 1.5 N HCl in methanol, which was called the acid-hydrolyzed IPL fraction (total IPL). DCM and MilliQ water were added, and the DCM fraction was collected (repeated 4 times). The DCM fraction was rinsed (6 times) with MilliQ water in order to remove acid and dried under N 2 gas. The last 1/3 IPL fraction was subjected to base hydrolysis (2 h) in a 1N KOH in methanol/H 2 O mixture (95:5, v/v), which was called the base-hydrolyzed IPL fraction (phospho IPL). Together with the condensed CL fraction, the CL fraction, two fractions of IPL-derived core lipids and non-cleaved IPL fraction were dissolved in n-hexane/isopropanol (99:1, v/v), and filtered using PTFE filters (pore diameter of 0.45 µm). The acidhydrolyzed IPL reflected total IPL, which, prior to hydrolysis, were attached to both phosphatidyl and glycosidic head groups; The base-hydrolyzed IPL represented IPLs with phosphatidyl head groups only (phospho IPL; Weijers et al., 2011b). Analysis of the non-hydrolyzed IPL fraction was performed to determine any carryover of CL into the IPL fraction.

GDGT-Based Indices
Indices based on the fractional abundance of GDGTs were calculated as follows:  (1) ( 3) with the GDGT numbers corresponding to the GDGT structures in Figure S1. Note that RI 2 was originally developed for the current study and modified from RI 1 (Zhang et al., 2016), in which the fractional abundance of crenarchaeol was replaced by GDGT-4 in order to eliminate the influence of crenarchaeol on weighted average number of cyclopentane rings in GDGTs.

DNA Extraction and the Quantitative Polymerase Chain Reaction (qPCR)
The SPM samples (n = 12) from station R1 (river water), station M (mixing water), and station S (seawater) were selected for the DNA analysis. The frozen filters were washed 3 times by phosphate buffered saline (pH 7.4). The supernatants were centrifuged under 11,000 g for 10 min. The DNA was extracted following the protocol of FastDNA SPIN Kit. The DNA samples were dissolved with a final dilution in 100-µL deionized water and preserved at −80 • C until further processing. The DNA concentrations were quantified in duplicate with a Nano-Drop spectrophotometer (Thermo Fisher Scientific Inc., Wilmington, DE, USA). The quantitative PCR primers were Arch_334F (5 ′ ACGGGGCGCAGCAGGCGCGA3 ′ )/Arch_ 518R (5 ′ ATTACCGCGGCTGCTGG3 ′ ) for total archaeal 16S gene quantification (Bano et al., 2004), GII-554f (5 ′ GTCGMTTTTATTGGGCCTAA3 ′ ), and Eury806-r (5 ′ CACAGCGTTTACACCTAG3 ′ ) for MG-II Euryarchaeota 16S gene quantification (Massana et al., 1997;Teira et al., 2004). The qPCR analysis was performed at 95 • C for 30 s and 40 cycles at 94 • C for 30 s, 55 • C for total archaea and 53 • C for MG-II Euryarchaeota for 30 s and 68 • C for 1 min. Triplicate measurements were run for each sample and standard.
PCR bands of 16S rRNA gene and MG-II gene were amplified from SPM samples in Station M. They were recovered by a Gel Extraction Kit (omega) and sequenced on the 3730sequencing platform. The sequences were annotated as the corresponding target genes, which demonstrated the specificity of the chosen qPCR primers. A dilution series of purified DNA from those positive clones were used as standards. A melting curve analysis was performed to demonstrate that the fluorescent signal obtained in a given reaction was consistent with the expected profile for specific PCR products. The R 2 values of standard curves were >0.99. The efficiency of each qPCR was between 87 and 99%.
Amplicon Sequencing of the Archaeal 16S rRNA Gene SPM samples collected in river water, mixing water and seawater were selected to conduct MiSeq pyrosequencing targeting the archaeal 16S rRNA gene. In contrast to the qPCR primers, these primers targeted longer sequences to increase the precision of phylogenetic analysis. The primers were Arch_344F (5 ′ ACGGGGCGCAGCAGGCGCGA3 ′ ) and Arch_915R (5 ′ GTGCTCCCCCGCCAATTCCT3 ′ ; Gantner et al., 2011). The MiSeq sequencing was conducted on the MiSeq platform (2 × 250 PE, Illumina) at the Shanghai Personalbio Biotechnology (Shanghai, China). Mothur (version 1.29.2; Schloss et al., 2009) was applied to filter the raw pyrosequencing data. The selected sequences were analyzed using the QIIME standard pipeline (Caporaso et al., 2010). Taxonomy was assigned according to the Ribosomal Database Project (RDP) classifier 2.2 (minimum confidence of 80%; Cole et al., 2009). The GenBank accession numbers are PRJNA38421 for those archaeal 16S rRNA genes.

Satellite-Derived Surface Water Temperature (SWT)
The satellite-derived SWT was determined with a spatial resolution of 4 km from the NOAA advanced very-highresolution radiometer (AVHRR; version 5.2; http://www.nodc. noaa.gov/SatelliteData/pathfinder4km/). The June mean SWT was obtained from the daily averaged values of 30 days in June 2011 (sampling month). The annual mean SWT and winter mean SWT represented 8-year mean values of annual mean temperature (2004-2011) and monthly mean temperature (December-February), respectively, as the surface sediment (top ca. 10 cm) collected in this study might represent a deposition of 6-10 years, based on an estimation from Strong et al. (2012).

TEX 86 -Derived Temperature and Ring Index
The TEX 86 -derived temperature was calculated based on the calibration of Kim et al. [SST = 68.4 × log (TEX 86 ) + 38.6; (Kim et al., 2010)]. The CL-TEX 86 temperatures derived from either the SPM or the surface sediments were close to the satellitebased annual mean SWT in the river water and mixing water, station R and M, respectively; whereas CL-TEX 86 temperatures for the seawater station S were lower than the winter mean SWT (Figure 2). The correspondence between the CL-TEX 86 temperatures derived from SPM and sediment samples (R 2 = 0.70, P < 0.01) along the salinity gradient indicated that the TEX 86 signal in the sediment predominantly reflected archaea from the water column, which is consistent with previous studies in the PR estuary (Wang et al., 2015) and other coastal settings (Herfort et al., 2006;Zell et al., 2014).
Total IPL-TEX 86 temperatures in the water column were consistently lower than both the June mean SWT and in situ measurements, although the sampling season was summer (Figure 2). In the mixing water and seawater, the phospho IPL-TEX 86 was lower than the total IPL-TEX 86 whereas it was very close to the CL-TEX 86 (Figure 2; Table 1). Since phosphatidyl head groups can be degraded faster than the glycosidic head groups, the phospho IPL is considered to be a better reflection of the living microorganisms (Harvey et al., 1986;Schouten et al., 2010), which may explain the deviation between these IPL pools. Relatively rapid conversion of phospho IPL to CL may result in more similar ring distributions and thus TEX 86 values between phosphor IPL and CL. Furthermore, in the same study area, the variability in TEX 86 was suggested to be due to changes in archaeal community composition in the water column, in which the unusually low TEX 86 -derived temperature in the coastal SCS was speculated to be linked to MG-II Euryarchaeota (Wang et al., 2015).
Since TEX 86 can be influenced by factors other than temperature, the ring index was proposed to evaluate the accuracy of TEX 86 in marine sediments (Zhang et al., 2016). Here, CL-, phosphor IPL-, and total IPL-TEX 86 values were plotted against RI 1 (Equation 2; Figure 3) using data derived from the SPM and surface sediment samples collected during the current and previous studies (Wei et al., 2011;Ge et al., 2013;Zhang et al., 2013;Wang et al., 2015). Most SPM and surface sediment samples from the open South China Sea plotted within the RI 1 -TEX 86 -confined zone [RI 1 = 3.32 × (TEX 86 ) 2 − 0.77 × TEX 86 + 1.59, ±2σ ∼ 0.3; (Zhang et al., 2016)], whereas the majority of samples from coastal SCS and the PR estuary fell above the calibration zone (Figure 3), exhibiting higher ring index values than those from the open SCS. This implies that the GDGT pool in the water column of the PR estuary and coastal SCS comprised relatively more cyclopentane rings than predicted from the measured SWT. If it is assumed that GDGTs produced only by Thaumarchaeota underpin the relationship between SWT and TEX 86 (e.g., Schouten et al., 2002), then it follows that either Thaumarchaeota in the PR estuary respond differently to temperature than marine strains, or there is another source of cycloalkyl-containing GDGTs in the estuary.
To further assess the distribution of RI values in the SPM and to explore other contributor(s) to the cyclopentyl GDGT pool in the study area, mean values of CL-, total IPL-, and phospho IPLring index were examined (Figure 4). Note that crenarchaeol, a biomarker for MG-I Thaumarchaeota, was excluded from the ring index calculation (RI 2 , Equation 3) in order to limit its overwhelming influence on the index. The re-defined RI 2 equation is more sensitive to the variation of cyclopentanecontaining GDGTs that might be contributed by other archaea. Compared with the river water and seawater, the highest RI 2 value for either CL (avg. 0.39 ± 0.08) or IPL (avg. 0.48 ± 0.07 for total IPL; avg. 0.47 ± 0.10 for phospho IPL) occurred at station M in the mixing water (Figure 4; Table 1), suggesting that the PR estuary (station M) appeared to be a hot spot of production of GDGTs with cyclopentane moieties. Further confirmation came from the comparison of %GDGT 1-4 in different water FIGURE 4 | Distribution of the mean values of Ring Index (RI 2 ) compared with the abundance of MG-II Euryarchaeota 16S rRNA gene and the gene abundance ratio of MG-II Euryarchaeota to total archaea along the salinity gradient from the river water to seawater. RI (Equation 3) was calculated from CL (red bars), total-IPL (yellow bars), and phosphor-IPL (blue bars). Details are shown in Table 1. settings (Table S1), which showed that the sum of the fractional abundances of the GDGTs with 1-4 cyclopentane moieties in the mixing water was significantly higher than that in the seawater or river water. In the mixing water station, the mean values of total IPL-RI and phospho IPL-RI were not significantly different; both were higher than the CL-RI (Figure 4). In the seawater station, however, the total IPL-RI (0.34 ± 0.07) was more elevated than the phospho IPL-RI (0.22 ± 0.01) and the CL-RI (0.24 ± 0.04; Figure 4). This ring index distribution pattern at station S corresponded to the TEX 86 -temperature distribution in the SPM and sediment samples (Figure 2), suggesting that cyclopentanecontaining GDGTs altered the TEX 86 record in the water column, and the vertical transportation of GDGTs from the water column appeared to be a predominant source to the sediment.

Relationship between MG-II Euryarchaeota and Cyclopentane-Containing GDGTs
Previous studies reported that a significant proportion of MG-II Euryarchaeota was diversely present in the estuarine and coastal regions, including the Pearl River estuary (Liu et al., 2014;Wang et al., 2015), the Yangtze River estuary (Liu M. et al., 2011), and the Jiulong River estuary (Hu et al., 2015). In this study, the MG-II Euryarchaeota 16S rRNA gene averaged 5.4 ± 5.9 × 10 8 copies L −1 (n = 6) at the mixing water station, which was two to three orders of magnitude higher than that in river water station (avg. 1.5 ± 2.1 × 10 5 copies L −1 , n = 2) and seawater station (avg. 4.9 ± 9.4 × 10 6 copies L −1 , n = 4; Figure 4). Considering the heterotrophic life style of MG-II, which have been demonstrated by the former genetic analysis (Iverson et al., 2012;Li et al., 2015;Martin-Cuadrado et al., 2015) and cultivation experiment (Orsi et al., 2015(Orsi et al., , 2016, the high abundances of MG-II in the mixing water station seem to be due to the high phototrophs that enhanced by nutrients input from upper river (Gan et al., 2014).
The ratio of MG-II Euryarchaeota 16S rRNA gene abundance to total archaeal 16S rRNA gene abundance ([MG-II 16S]/[Archaea 16S]) in mixing water (avg. 0.25 ± 0.08, n = 6) was significantly higher than in the seawater (avg. 0.10 ± 0.004, n = 4; P < 0.01), whereas it was negligible (avg. <0.0001) in the river water (Figure 4). This observation was also supported by pyrosequencing analysis, which exhibited a linear correlation with the qPCR-based ratio of [MG-II]/[Archaea 16S] (Figure S2; Table S2). These results further confirmed that the PR estuary (mixing zone, salinity avg. 16.6) provided a habitat to sustain a natural enrichment of planktonic MG-II Euryarchaeota.
The presence of (more labile) phospho IPL-GDGTs implied in situ production of isoprenoidal GDGTs in the water column along the entire salinity gradient from the Pearl River to coastal SCS. The elevated phospho IPL-RI in the mixing water FIGURE 5 | RI 2 (A), fractional abundance of GDGTs (B-E, G) and Crenarchaeol and isomer (H,F, respectively), and the ratio of GDGT-2 to GDGT-3 (I) vs. the ratio of the MG-II Euryarchaeota 16S rRNA genes to the total archaeal 16S rRNA genes. The black points represent SPM samples collected from the lower Pearl River, the PR estuary, and the coastal SCS.
Frontiers in Microbiology | www.frontiersin.org additionally suggests that higher relative proportions of GDGTs with 1-4 cyclopentane moieties were produced in the PR estuary. A linear regression analysis confirmed the positive relationship between phospho IPL-derived RI and the [MG-II 16S]/[Archaea 16S] ratio in the water column along the salinity gradient (R 2 = 0.72, P < 0.01; Figure 5A). Therefore, it is reasonable to hypothesize that MG-II Euryarchaeota preferentially synthesized GDGTs with 1-4 cyclopentane moieties in this region. Moreover, production of GDGTs-1 to -4 by MG-II Euryarchaeota could represent the missing source needed to explain the elevated value of RI and deviation from the RI-TEX 86 relationship driven by the preservation of Thaumarchaeota GDGTs in global core-top sediments.
To further constrain the relationship between MG-II Euryarchaeota and cyclopentane-containing GDGTs, linear regression analysis was conducted between the fractional abundance of GDGTs and the [MG-II 16S]/[Archaea 16S] ratio in the SPM along the salinity gradient. Results exhibited a significantly positive linear correlation between the [MG-II 16S]/[Archaea 16S] ratio and the fractional abundance of phospho IPL-based GDGT with cyclopentane moieties (Figures 5B-F). In contrast, we observed no correlation between the [MG-II 16S]/[Archaea 16S] ratio and the fractional abundances of phospho IPL-based GDGT-0 or crenarchaeol (Figures 5G,H; Table 2). Similar trends of the linear correlations were also observed between the [MG-II 16S]/[Archaea 16S] ratio and the CL-and total IPL-GDGTs with 1-4 cyclopentane moieties, with a less significant correlation between the [MG-II 16S]/[Archaea 16S] ratio and the total IPL-based crenarchaeol regioisomer ( Table 2). Although MG-II Euryarchaeota were suggested to be an alternative source of crenarchaeol in the ocean (Lincoln et al., 2014a), our study showed an absence of a significant correlation between the distribution of MG-II Euryarchaeota and crenarchaeol ( Figure 5H). However, members of MG-II Euryarchaeota living in the North Pacific Subtropical Gyre (i.e., those targeted by Lincoln et al., 2014a) may be different from those living in the coastal zone of the PR estuary. This is consistent with the phylogenetic distribution of MG-II reported by Wang et al. (2015), which showed diverse groups of MG-II living in this region.
By comparison of the R 2 values and slopes of the regression equations (Figures 5B-F), GDGT-1 exhibited not only the strongest correlation with MG-II Euryarchaeota in the study area, but also largest relative enrichment (i.e., slope of 14.06). If MG-II Euryarchaeota preferentially synthesized GDGT-1, additional contribution of GDGT-1 to the water column of the PR estuary and the coastal SCS would be reflected by an increase in RI values and a substantial decrease in TEX 86 values. Offsets in TEX 86 values have been similarly proposed to result from the decreased ratio of GDGT-2 to GDGT-3 in the surface sediments of this area (Wang et al., 2015), whereas the increased ratio of GDGT-2 over GDGT-3 in the deep-water column seems to be responsible for a warm bias of TEX 86derived temperature in other marine environments (Taylor et al., 2013;Hernandez-Sanchez et al., 2014). In this study, however, no correlation is exhibited between GDGT-2/3 ratio and [MG-II 16S]/[Archaea 16S] ratio ( Figure 5I). A recent study by Kim et al. (2015) suggested that coincident increases in GDGT-2 and the crenarchaeol regioisomer and decreases in GDGT-1 and GDGT-3 shifted TEX 86 -derived temperatures toward higher values in the deep-water surface sediments of the Mediterranean Sea. In combination with our data (Figures 4, 5), these observations are consistent with the interpretation that planktonic Euryarchaeota have the potential to bias TEX 86 by changing the distribution of TEX 86 -related GDGTs (especially those with 1-3 cyclopentane rings) in different marine environments.

CONCLUSIONS
This study assessed the relationship between TEX 86 -related GDGTs and MG-II Euryarchaeota along a salinity gradient from river water to seawater. The fractional abundance of MG-II Euryarchaeota was correlated with %GDGTs with cyclopentane moieties as well as ring index values, implying that MG-II Euryarchaeota may contribute ringed GDGTs to the total GDGT pool. This source would thus increase the ring index value and potentially bias the TEX 86 proxy. However, MG-II Euryarchaeota living in the estuary and coastal region did not seem to be a significant source of crenarchaeol. These and other findings based on environmental distributions provide indirect evidence of the lipid profile of MG-II Euryarchaeota, which cannot be validated until a pure culture is available.

AUTHOR CONTRIBUTIONS
JW and CZ designed this study. JW extracted and analyzed lipids. WX extracted and analyzed DNA. JW, WX, YZ, TM, and CZ wrote the paper.
on an earlier version of the manuscript. An early version of this manuscript was discussed online at the Biogeosciences Discussion forum (doi: 10.5194/bgd-12-12455-2015); however, the final paper was declined for publication by Biogeosciences. This research was supported by the National Key Basic Research Program of China grant #2013CB955703 and 2016YFA0601101 (CZ), and the National Science Foundation of China Grant # 41530105 and 41673073 (CZ). This study is also a contribution to the international IMBER project. Comments from the two reviewers and the editor significantly improved the quality of the paper and are greatly appreciated.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb. 2017.02077/full#supplementary-material Figure S1 | Structures of archaeal core GDGTs described in the text. Figure S2 | Relationship of fractional abundances of MG-II Euryarchaeota derived from quantitative polymerase chain reaction (qPCR) and pyrosequencing. % MG-II (qPCR) is calculated based on the ratio of MG-II 16S rRNA gene to Archaea 16S rRNA gene; % MG-II (sequencing) refers to the fractional abundance of the MG-II OTUs to the total archaeal OTUs.