Benthic protists and fungi of Mediterranean deep hypsersaline anoxic basin redoxcline sediments

Some of the most extreme marine habitats known are the Mediterranean deep hypersaline anoxic basins (DHABs; water depth ∼3500 m). Brines of DHABs are nearly saturated with salt, leading many to suspect they are uninhabitable for eukaryotes. While diverse bacterial and protistan communities are reported from some DHAB water-column haloclines and brines, the existence and activity of benthic DHAB protists have rarely been explored. Here, we report findings regarding protists and fungi recovered from sediments of three DHAB (Discovery, Urania, L’ Atalante) haloclines, and compare these to communities from sediments underlying normoxic waters of typical Mediterranean salinity. Halocline sediments, where the redoxcline impinges the seafloor, were studied from all three DHABs. Microscopic cell counts suggested that halocline sediments supported denser protist populations than those in adjacent control sediments. Pyrosequencing analysis based on ribosomal RNA detected eukaryotic ribotypes in the halocline sediments from each of the three DHABs, most of which were fungi. Sequences affiliated with Ustilaginomycotina Basidiomycota were the most abundant eukaryotic signatures detected. Benthic communities in these DHABs appeared to differ, as expected, due to differing brine chemistries. Microscopy indicated that only a low proportion of protists appeared to bear associated putative symbionts. In a considerable number of cases, when prokaryotes were associated with a protist, DAPI staining did not reveal presence of any nuclei, suggesting that at least some protists were carcasses inhabited by prokaryotic scavengers.


INTRODUCTION
Redox boundaries in marine sediments can have significant geochemical gradients, transitioning from fully aerated to anoxic (lack of detectable dissolved oxygen) conditions within short vertical distances (e.g., Cai and Sayles, 1996;Yucel, 2013). These chemoclines are zones of intense biogeochemical cycling, involving all major elements including carbon, oxygen, nitrogen, sulfur, and hydrogen as well as iron and manganese. While considerable effort has been dedicated to studying the biogeochemistry of redox boundaries of marine sediments in neritic zones (e.g., Yucel, 2013), silled basins (e.g., Reimers et al., 1996;Bernhard et al., 2003), and hydrocarbon seeps (e.g., Orcutt et al., 2005), less is known regarding chemoclines in the deep bathyal to hadal zones.
The deep bathyal Mediterranean has numerous deep hypersaline anoxic basins (DHABs; Figure 1 and also see Figure 1 in Stock et al., 2013b), which are brine-filled bathymetric depressions formed from the dissolution of subterranean Miocene salt deposits exposed to seawater after tectonic activity (Camerlenghi, 1990). Due to the sequence of different chemical ions precipitating from seawater as it evaporates, different layers of those salt deposits are characterized by differing chemistries. It follows that DHABs differ in brine chemistry. In the three DHABs studied here, Urania brine is the highest in free sulfide and methane, Discovery brine is highest in chlorine and magnesium, and L' Atalante brine is highest in sodium, potassium, and sulfate ( Table 1). Further details on the brine chemistries are available in van der Wielen et al. (2005). These DHABs have been the topic of attention over the last decade. In particular, the water columns overlying these bathymetric features have received considerable study, for both prokaryote communities (van der Wielen and Heijs, 2007;Ferrer et al., 2012) and eukaryote communities (e.g., Alexander et al., 2009;Edgcomb et al., 2009Edgcomb et al., , 2011dStock et al., 2013a). The sediments overlain by the brines have also been investigated for their prokaryotic (e.g., van der Wielen and Heijs, 2007;Yakimov et al., 2007;Ferrer et al., 2012) as well as metazoan constituents (Danovaro et al., 2010). The benthic redoxclines of DHABs have received less study (e.g., the study of Red Sea brine pool benthic eukaryotes by Wang et al., 2014), probably due to their limited spatial extent. Haloclines of DHABs are typically ∼2 m in vertical thickness (van der Wielen et al., 2005;Hallsworth et al., 2007;Yakimov et al., 2007). Thus, when the halocline impinges the seafloor, the pitch of the bathymetric slope directly impacts the lateral manifestation of the halocline (redoxcline) on the seafloor (Figure 2).
The goal of this study was to investigate the microbial eukaryotes in the redoxcline of Mediterranean DHAB sediments. Here we report results describing the eukaryotic community composition and abundances. We include results from classic morphologic analyses as well as from RNA-based sequencing efforts and a fluorescence-labeling approach to identify active eukaryotes. For comparison, we include data on microbial eukaryote communities obtained from adjacent normoxic sediments.

SAMPLE COLLECTION AND SALINITY
Samples were collected using the ROV Jason from 24 November to 6 December 2011. Most sediments for this study were collected with the typical Deep Submergence Lab (www.whoi.edu/groups/DSL/) Alvin-type pushcores (6.35-cm diameter; hereafter referred to as "cores") configured with a seal to prevent contamination during ascent. Cores were collected from nearby (within 100 m) aerated sediments, just upslope from each DHAB. Thus, control samples at Urania and Discovery were consistently from 1 to 3 m shallower water depths than the halocline samples. Because ROV Jason is neutrally buoyant, it was not possible to core deep in the DHAB brine. Thus, cores of deep halocline sediments were obtained by Jason reaching toward the brine with the ROV manipulator. The depth differential between halocline samples did not exceed 2 m. Given the waters depths for all samples were >3400 m, these depth differences are negligible. Coordinates and water depths reflect our approximate sampling area for each DHAB. Upon ROV Jason recovery, all cores were taken as quickly as possible (within 5 min) to the RV Atlantis environmental room set at 9.5 • C (±0.5 • C), until further processing. Cores for microelectrode profiling were placed in an N 2 -flushed glove bag and profiled as described below. Cores for quantification (counts) and imaging were first sampled for their bottom-water salinity, measured via hand-held refractometer. Many samples needed to be diluted up to 10-fold to register on the refractometer. The handheld refractometer was not expected to provide precise salinity data, especially for the very high MgCl 2 brine of Discovery, but was intended as a relative indicator between samples from different halocline horizons in any given DHAB. After salinity was recorded, each of these cores was sectioned into 1-cm intervals to a depth of 3 cm. Each sediment slice was preserved in either 3.8% formaldehyde or 3% TEM-grade glutaraldehyde in a 0.1 cacodylic acid sodium salt buffer. In some cases, there was no overlying bottom water due to full penetration or over-penetration of the corer. In these cases, salinity data are not reported. In instances where no salinity data was obtained from any given core, we estimated salinity using published literature ( Table 1) and confirmed by our own water column high-range CTD data (collected for a separate water-column study, data not shown), extrapolated for specific core location relative to others with known salinity and the halocline boundary. We designate halocline samples by their relative location (depth, refractometer data) in the following descending order: upper halocline, mid-halocline, lower halocline.

Frontiers in Microbiology | Extreme Microbiology
To avoid the possibility that pressure and temperature changes during ascent may cause death of some fraction of the eukaryotic community, a subset of cores were incubated in situ using the viability indicator CellTracker™ Green CMFDA (chlorofluoromethyl fluorescein diacetate; Life Technologies). CellTracker Green relies on active esterases, which are known to function in at least one DHAB brine (Ferrer et al., 2005). For these samples, push corers were emplaced on the seafloor and left for ≥24 h prior to collection. These pushcores were modified to allow injection of the CellTracker Green using the ROV manipulator. Similar "injector" cores have been employed before for seafloor incubations using the same fluorescent probe (Bernhard et al., 2009). During the next Jason dive (2-3 days later), these cores were removed from the seafloor per standard procedure and brought to the sea surface. Because only one dive was executed in L' Atalante, those incubations were begun early on the dive and terminated after recovery later that day, resulting in incubation times of approximately 12.5 h. The injector cores were sectioned and preserved in glutaraldehyde as noted above.

SEDIMENT PROFILING AND SAMPLING FOR RNA EXTRACTION
In the environmental room on RV Atlantis, designated cores containing overlying site water were profiled for dissolved oxygen with a vertical resolution of 0.25 mm in a N 2 -filled glove bag using polarographic microelectrodes (Visscher et al., 1991). These cores were then subsampled from each core center using a sterile 20-ml syringe (∼1.4 cm inside diameter) with the Luer end removed. These syringe cores were frozen at −80 • C and transported to WHOI, where they were slightly thawed to enable slicing of the subcore. The surface 2 cm of the subcore was www.frontiersin.org placed in conical centrifuge tubes and stored at −80 • C until RNA extraction.

CHARACTERIZATION OF SEDIMENT EUKARYOTES USING SMALL SUBUNIT RIBOSOMAL RNA (SSU rRNA)
RNA from approximately 8 g of the top 2 cm from each frozen syringe subcore was extracted using an optimized protocol with the RNA Power Soil kit (MoBio, USA). Major modifications included introduction of three cycles of freeze-thaw (−80 • C, 5 min., 65 • C, 5 min.), bead beating with 2 × 5 min. intervals on a horizontal vortexer, an overnight nucleic acid precipitation, and a 1-hour centrifugation during the precipitation step. In addition, we introduced two DNAase treatments using TurboDNAase (Ambion, USA). Removal of DNA was confirmed by PCR using general eukaryotic primers. RNA was purified using the MEGAclear kit (Ambion, USA). Reverse transcription of the purified RNA samples was performed using the QuantiTect kit (Qiagen, USA). Tagpyrosequencing of the eukaryotic small subunit rRNA (18S rRNA) gene was performed using PCR amplification of the V4 region of the 18S rRNA gene and the primer pair TAReuk454FWD1 (5 -CCAGCA(G/C)C(C/T)GCGGTAATTCC-3 ) and TAReukREV3 (5 -ACTTTCGTTCTTGAT(C/T)(A/G)A-3 ). In brief, a one-step, 30-cycle PCR reaction was performed using GoTaq polymerase (Promega, USA). PCR conditions included: 94 • C for 3 min, followed by 30 cycles of 94 • C for 30 s; 55 • C for 40 s; 72 • C for 1 min; and a final elongation step at 72 • C for 5 min. Following PCR, all amplicon products (ca. 450 bp) from different samples were purified using the MinElute Reaction Cleanup kit (Qiagen, USA). In parallel, blank extractions and amplification reactions were performed to check for contamination. Distinct tags (multiplex identifiers) were used for each of the samples. Amplicon libraries were sequenced utilizing the Roche 454 FLX Titanium platform and reagents following manufacturer's guidelines at the MR DNA (Molecular Research LP, Shallowater, TX, USA) sequencing facility.
Denoising of the flowgrams was performed using Acacia (Bragg et al., 2012). Processing of the resulting sequences, i.e., trimming and quality control was performed using QIIME (Caporaso et al., 2010). Sequences with ≥380 bp and no ambiguous base calls and no homopolymers ≥6 bp were included in further analysis. All sequences were binned into operational taxonomic units (OTUs) and were clustered (average neighbor algorithm) at 97% sequence similarity identity. Taxonomic assignments were made using BLAST within QIIME. Sequences have been deposited in the GenBank SRA archive under the accession number SRP049010.

QUANTIFICATION AND IMAGING OF PROTISTS USING MICROSCOPY
At WHOI, pushcore samples preserved in either glutaraldehyde or formalin on the ship were used for morphologic documentation and/or quantification. Percoll density gradients were used to isolate the eukaryotes (Starink et al., 1994), subjecting to each established gradient 0.5 ml of sample from the surface 1 cm. The eukaryote-laden supernatant of each density gradient was filtered onto one 2-μm blackened polycarbonate filters and mounted with Citifluor AF1: Vectashield: 1× PBS (Phosphate buffered saline; 11:2:1) media containing 1 mg/ml DAPI (4 ,6-diamidino-2-phenylindole) on a microscope slide under a coverslip. Eukaryotes were counted from each filter [up to 60 fields were counted at 400× magnification on a Zeiss Axio Imager M2 upright microscope with differential interference contrast (DIC) and phase contrast optics, until 350 cells per filter were counted]. For each sample, eukaryotes were isolated from three separate Percoll gradients. Thus, three replicate filters were counted for each sample.
Because the filter impedes imaging, the specimen-laden side of a subset of filters from selected samples was blotted onto a drop of mounting media atop a microscope slide, covered with a coverslip and then sealed. The mounted slide was scanned for eukaryotes appropriate for imaging (clean, larger, mostly unobstructed specimens). Images were taken either with an Olympus BX51 upright epifluorescence microscope with DIC and phase contrast optics and an Olympus DP70 color digital camera or the Zeiss Axio Imager M2 and a Canon EOS Rebel color digital camera.
Non-quantitative aliquots of samples incubated on the seafloor with CellTracker Green were similarly subjected to Percoll extraction, DAPI staining and/or imaging, although specimens were also examined under 480-nm excitation and 520-nm emission to determine if they had active esterases at the time of incubation.

PORE-WATER DISSOLVED OXYGEN
Dissolved oxygen in control cores remained detectable to depths deeper than 1 cm ( Figure 3A). In all but one of the halocline cores, oxygen was detectable in overlying waters (Figures 3B,C). At the sediment-water interface of halocline cores, [O 2 ] was typically below 21 μM. Oxygen remained detectable to at least 2.0 cm depth in 7 of the 12 halocline cores. In some halocline cores (i.e., L' Atalante dive 611 core 1 and core 6; Discovery dive 609 core 10), oxygen remained detectable at the deepest measured horizon (2.5 cm). Oxygen was near the detection limit (<3 μM) throughout the top 2.5 cm in core 608 c 11 (Urania); this core was not analyzed for eukaryotes.

BOTTOM-WATER SALINITY
Salinity in bottom waters within 10 cm of the sediment-water interface of cores from Discovery and L' Atalante varied widely. Salinity was not measured from cores collected in and near Urania DHAB. Control cores had bottom-water salinities of 41-50. The highest salinity measured was 220 (L' Atalante; Table 2); the highest salinity measured in Discovery was 146. We reiterate that these salinities do not represent brine characteristics per se but allow relative comparisons between pushcores collected within each halocline.

PROTIST MORPHOTYPES
Protist morphotypes were detected from every habitat (normoxic control, upper halocline, mid-halocline, lower halocline) in each sample examined for this purpose ( Table 2). A number of morphotypes were observed (Figures 4-8). For example, what appeared to be a thecate foraminifera was observed in a control sample (Figure 4A), but none of these were noted Frontiers in Microbiology | Extreme Microbiology  in halocline samples. Small (3-5 micron) eukaryotes that were seen under fluorescence microscopy appeared to be consistent with flagellates, dominated microscopic counts in terms of abundance, but the most photogenic protists were ciliates (Figures 5-7). Not all protists had discernable nuclei. When no nuclei or few prokaryotes via DAPI were observed, the specimen was deemed dead (Figures 4C,D). In some cases, nuclei were visible in both DIC and DAPI images (e.g., Figures 5A,B,D,E and 7A,C,D,E but DAPI staining of nucleus not shown). In other cases, even though prokaryotic associates (see below) were discernable by DAPI staining, nuclei that should also be discernable via DAPI staining were not obvious (e.g., www.frontiersin.org     Figures 6D,E). Some samples (i.e., 610 c14; 610 cL; 611 c10; 611 cC) had numerous filamentous bodies that were likely fungal hyphae or, perhaps, filamentous bacteria similar to Beggiatoa, although well-developed Beggiatoa mats were not observed. Moreover, Beggiatoa sequences were not detected in our samples (Kormas et al., unpublished data). Metazoa that were found (Bernhard et al., unpublished data) are not considered in this contribution.

PROTIST VIABILITY SUGGESTED BY CELLTRACKER GREEN
In the aliquots examined for both CellTracker Green labeling and DAPI labeling, esterase-positive (i.e., CellTracker Green labeled) protists were observed that had easily distinguishable nuclei, either in addition to or exclusive from any associated ecto-and/or endobionts (Figure 8). Such specimens were observed in normoxic and some halocline samples (i.e., 607-608 cB and cE; 609-610 cL).

MOLECULAR SIGNATURES OF PROTISTS AND FUNGI
After quality trimming and removal of non-protistan sequences and fungal sequences, the number of eukaryotic reads obtained from each sample analyzed for eukaryote diversity varied Frontiers in Microbiology | Extreme Microbiology drastically, with the most reads (6546) in the Discovery upper halocline, followed by both L' Atalante halocline samples (539 reads in lower halocline, 390 reads in upper halocline), and Urania halocline (82 reads) and Urania normoxic control (75 reads). For two of the samples analyzed, the number of curated reads was too low and these were excluded from the graphical representation. Discovery mid-halocline had a total of 8 reads, where most were Ustilaginomycotina Basidiomycota (not shown). Discovery lower halocline had a total of 6 reads, where most were Chlorophyceae (not shown). A higher number of orders of eukaryotic taxa based on SSU rRNA signatures was detected in the Urania aerated control sample compared to the halocline samples, implying that alpha diversity may be lower inside these DHABs. The only exception was L' Atalante upper halocline, which also had relatively high diversity (Figure 10). In four of these five samples, the majority of eukaryotic reads were fungal, with most being affiliated to Ustilaginomycotina Basidiomycota. L' Atalante upper halocline was dominated by Chlorophyceae. Other groups with high diversity were Syndiniophyceae (L' Atalante upper halocline, Urania aerated) and Dinophyceae (L' Atalante upper halocline) Dinophyta. Rhizarians were also present, sometimes in considerable proportions (13%, L' Atalante upper halocline), mostly as Radiolaria.

DISCUSSION
Based on microscopic and sequencing analyses, protists and fungi were present in all our sediment samples, including lower halocline (higher salinity) samples. Eukaryote densities in the control (normoxic) sediments were roughly similar to densities of other protist communities in some oxygenated bathyal sediments (i.e., off Southern California, CA, USA; Bernhard et al., 2000), but higher than those in other normoxic bathyal sediments (i.e., Coral Sea; Alongi, 1987). Densities of protists in the DHAB halocline sediments were approximately equal to those in bathyal chemocline sediments of a silled basin off Baja California (i.e., Soledad Basin; Bernhard and Buck, 2004) but were only about one third those of the eukaryote community of bacterial mat chemocline sediments of Santa Barbara Basin, CA, USA (Bernhard et al., 2000). The detected densities might be considered surprisingly high given the oligotrophic nature of the Mediterranean (e.g., Kress et al., 2003;Pujo-Pay et al., 2011;Huertas et al., 2012;Tanhua et al., 2013). Even though none of the haloclines supported well-developed bacterial mats, our quantitative data suggests FIGURE 10 | Relative abundance (%) of rRNA signatures of unicellular eukaryotic groups of the DHAB sediments based on V4 SSU rRNA pyrotags clustered at 97% similarity. Taxonomic assignments were performed using BLAST implemented in QIIME against SILVA 111 database.
www.frontiersin.org that the haloclines of these DHABs are inhabited by diverse eukaryotes, seemingly supporting higher numbers of eukaryotic microbes than the surrounding normoxic deposits. In deep oxic and oligotrophic Mediterranean waters above DHABs, protist counts are very low, ∼15-223 cells mL −1 , but within the water-column haloclines (e.g., Urania halocline), protist numbers can rival those in near-surface waters, ∼1071 ± 117cells mL −1 (Pachiadaki et al., in press). Our quantitative observations from the three DHAB benthic halocline habitats had two orders of magnitude more protists than in documented DHAB watercolumn haloclines (Edgcomb et al., 2011d;Pachiadaki et al., in press).
To determine whether protist density showed a relationship with salinity or chemical constituent, we plotted our protist density data (Figure 9) with published data on the seawater chemistry of each DHAB (Table 1). In general, higher protist densities occurred with higher concentrations of HS − , Mg 2+ , and SO 4 2 to a threshold (Figure 11). Protist densities decreased with the highest concentrations of sulfide and magnesium. No relationship was apparent between salinity and protist density. Refractometer readings of salinity for athalassohaline brines such as Discovery are problematic, and therefore comparisons to thalassohaline brines such as Urania and L' Atalante are difficult. Protist density may relate to the differing cation chemistries in each of the DHABs (Table 1) and whether encountered specimens were living or dead (see below). Statistical analyses to test for trends were not performed as protist FIGURE 11 | Abundances of protists per cubic centimeter, calculated from microscopic counts, presented with respect to inferred chemical parameter in overlying waters. Chemical constituent data from van der Wielen et al. (2005). data from replicate cores at each habitat and halocline site were not available.
Although a minority, a considerable number of protists supported associated prokaryotes. These putative symbionts were more commonly endobionts than ectobionts. In one case, what appeared to be clumped endobionts were noted (Figures 7J,K). This pattern is somewhat similar to observations of another ciliate with multiple endobionts in a structured configuration, inhabiting chemocline sediments (Edgcomb et al., 2011c). Alternatively, these clumped endobionts may be in food vacuoles. It is possible that some of the DAPI signals detected in some eukaryotic specimens were actually of other smaller protists. For example Figures 6D,E may show numerous parasitic protists inside a remnant cell of another eukaryote of unknown affinity. We recovered some sequences of Syndiniales, which are parasitic dinoflagellates known to sometimes occur in deep-sea hydrothermal vent settings (e.g., Sauvadet et al., 2010).
Observed protists may not have been living individuals but merely locales of prokaryote infestation. Specifically, anaerobic bacteria and/or archaea may have inhabited protists as suggested by the lack of a nucleus/nuclei in some protists (Figures 7D,E and 8B,C) yet presence of copious endo-or ectobionts evidenced by DAPI staining. Prokaryote presence within degrading foraminifera cytoplasm has been documented in hydrocarbon seep habitats (Bernhard et al., 2010b). Ribosomal RNA can be preserved in the high salinity environments of DHAB haloclines and brines (Hallsworth et al., 2007), and therefore, an unknown fraction of taxa whose rRNA signatures were detected here may represent inactive or dead and preserved cells. The detection of rRNA signatures of the photosynthetic Chlorophyceae in halocline samples in this study supports this suggestion, although the same signatures were detected in oxic normal saline sediments, which argues against preservation, and for the possibility that they may still be active. Ribosomal RNA signatures of additional photosynthetic organisms have previously been reported from deep marine sediments (Orsi et al., 2013). The potential for preservation of rRNA, particularly under anoxia and hypersaline conditions means that we must interpret rRNA-based data from these habitats cautiously, and find support for activity from microscopical observations and other types of molecular data. Furthermore, although here we report positively labeled eukaryotes with CellTracker Green (i.e., esterase activity), bacteria also react to this fluorogenic probe  could cause the eukaryotes to appear positively labeled. In sum, the protists present in the halocline samples may have been dead and only appeared to be inhabitants of these extreme habitats. Indeed, it may be that eukaryote carcasses, or "pickled protists" have been in the DHAB haloclines and brines since their formation. Estimates of Eastern Mediterranean DHAB formation are, at a minimum, 2000 years (Wallmann et al., 1997).
The presence of pickled protists may explain some of the observed trends of protist abundance (densities) with respect to water chemistry (Figure 11), at least for some constituents. More specifically, for Mg 2+ and HS − , it is plausible that the observed protist peak abundances coincide with the approximate maximum tolerance of protists for the constituent being considered (Figures 11A,C). Regarding Mg 2+ , although an upper limit for prokaryotic life may be ∼2 M (Hallsworth et al., 2007), concentrations exceeding 0.2 M MgCl 2 may be inhospitable to protists (Figure 11A). Similarly, DHAB environments exceeding ∼1.5 mM HS − ( Figure 10C) and ∼200 mM SO 4 2− (Figure 11B) may not support living protist populations. Thresholds probably also depend on concentrations of additional constituents and their possible synergistic interactions in these brines, so generalizations involving habitats with typical marine salinities are likely invalid.
At least some protists were living in normoxic and halocline DHAB sediments, as evidenced by positive CellTracker Green labeling along with DAPI analysis indicating the lack of associated microbes but the presence of a nucleus or nuclei (e.g., Figure 8). Importantly, we do not assert positive DAPI staining indicates viability, as we are aware that DNA can be preserved in DHAB brines (e.g., Borin et al., 2008). Regardless, at least some protists lacking prokaryotic associates appear to inhabit upper, mid-, and lower DHAB halocline sediments.
Our assertion that filamentous items are most likely Fungi is supported by the recent microscopic observations documenting filamentous fungi in the oxycline and deep anoxic zones of the Black Sea (Sergeeva and Kopytina, 2014). The rRNA sequence data suggest that Fungi dominate the eukaryotic community in the MgCl 2 -dominated Discovery basin upper halocline, and in the L'Atlante basin lower halocline, confirming previous findings of high abundance of Fungi in the lower haloclines of DHAB water columns (Stock et al., 2012;Pachiadaki et al., 2014). Fungi were also reported to dominate the DNA libraries of the hypersaline sediments from a Red Sea brine pool (Wang et al., 2014). In that study, cultured representatives from those brine sediments were affiliated to Ascomycota, while in the present study the majority of retrieved ribotypes of Fungi from the halocline sediments were identified as Basidiomycota, suggesting that different halotolerant (or halophilic) fungal groups might have adapted to different local environments. The libraries from the lower haloclines of both Discovery and L' Atalante share the same most abundant OTU. This OTU showed high nucleotide identity (99% with 100% sequence coverage) with an uncultured Malassezia isolate [KC487833] from the hypersaline Lake Tyrrell, VIC, Australia (Heidelberg et al., 2013). Recently and Amend (2014) provided evidence that the genus Malassezia is cosmopolitan and probably of great importance in deep-sea extreme environments. Interestingly, we have retrieved the same OTU (99% of sequence similarity) in the anoxic lower halocline water layer of another DHAB, Thetis (JF308281; Stock et al., 2012), as well as in the anoxic marine Cariaco Basin (GU824665; Edgcomb et al., 2011a), in the anoxic fjord Saanich Inlet (HQ866112; Orsi et al., 2012) and in the deep subsurface marine sediment of Peru Margin (GU972524; Edgcomb et al., 2011b). The marine species of Malassezia remain uncultured, so we cannot yet infer much about the ecology and trophic status of this OTU found in our survey.
Our data also support the notion that increasing salinity (e.g., L' Atalante basin rRNA profiles) leads to a decrease in relative abundance of selected groups, e.g., a decrease in abundance within Dinophyceae, Cryptophyceae, Centroheliozoa, and Acantharea in the lower halocline sample from L' Atalante basin relative to the upper halocline. The dominance of chlorophyte signatures in the Urania basin upper halocline sample suggests decreased protist diversity possibly due to the relatively high sulfide concentrations there.
The application of metatranscriptomics to these communities may be the only feasible means at present to fully resolve eukaryote activity (i.e., gene expression) and their inhabitation of these DHAB haloclines (Pachiadaki et al., unpublished data). Future efforts may reveal novel taxa and heretofore unknown biochemical pathways in these spatially restricted benthic habitats.