Impact Factor 4.235 | CiteScore 6.4
More on impact ›

ORIGINAL RESEARCH article

Front. Microbiol., 03 July 2020 | https://doi.org/10.3389/fmicb.2020.01506

Quantification of Cable Bacteria in Marine Sediments via qPCR

Jeanine S. Geelhoed1*, Sebastiaan J. van de Velde2 and Filip J. R. Meysman1,3*
  • 1Department of Biology, University of Antwerp, Antwerp, Belgium
  • 2Department of Earth and Planetary Sciences, University of California, Riverside, Riverside, CA, United States
  • 3Department of Biotechnology, Delft University of Technology, Delft, Netherlands

Cable bacteria (Deltaproteobacteria, Desulfobulbaceae) are long filamentous sulfur-oxidizing bacteria that generate long-distance electric currents running through the bacterial filaments. This way, they couple the oxidation of sulfide in deeper sediment layers to the reduction of oxygen or nitrate near the sediment-water interface. Cable bacteria are found in a wide range of aquatic sediments, but an accurate procedure to assess their abundance is lacking. We developed a qPCR approach that quantifies cable bacteria in relation to other bacteria within the family Desulfobulbaceae. Primer sets targeting cable bacteria, Desulfobulbaceae and the total bacterial community were applied in qPCR with DNA extracted from marine sediment incubations. Amplicon sequencing of the 16S rRNA gene V4 region confirmed that cable bacteria were accurately enumerated by qPCR, and suggested novel diversity of cable bacteria. The conjoint quantification of current densities and cell densities revealed that individual filaments carry a mean current of ∼110 pA and have a cell specific oxygen consumption rate of 69 fmol O2 cell–1 day–1. Overall, the qPCR method enables a better quantitative assessment of cable bacteria abundance, providing new metabolic insights at filament and cell level, and improving our understanding of the microbial ecology of electrogenic sediments.

Introduction

Cable bacteria are long, multicellular, filamentous bacteria that transport electrons from cell-to-cell along the filament over distances of up to several centimeters (Pfeffer et al., 2012; Bjerg et al., 2018; Meysman, 2018; Meysman et al., 2019). This long-distance electron transport enables an electrical coupling of sulfide oxidation in deeper anoxic sediment layers to oxygen reduction near the sediment surface (Nielsen et al., 2010). Electrogenic sulfur oxidation (e-SOx) by cable bacteria makes an important contribution to sulfur cycling in aquatic sediments (Risgaard-Petersen et al., 2012; Malkin et al., 2014; Burdorf et al., 2017; Sandfeld et al., 2020), and has a strong imprint on the local geochemistry of the seafloor (Risgaard-Petersen et al., 2012; Meysman et al., 2015). The sulfide oxidizing half-reaction results in an acidification of the pore water, leading to the dissolution of minerals, such as carbonates and sulfides, thereby greatly affecting the cycling of elements like Ca, Fe, Mn, P and trace elements (Risgaard-Petersen et al., 2012; Rao et al., 2016; Sulu-Gambari et al., 2016; van de Velde et al., 2016, 2017). The dissolution of iron sulfide supplies additional sulfide to the cable bacteria and results in a strong increase of the ferrous iron (Fe2+) concentration in the pore water. Part of this aqueous Fe2+ will diffuse to the oxic sediment surface, where it will reoxidize and reprecipitate as iron(III) (hydr)oxides (Seitaj et al., 2015; Rao et al., 2016). This precipitated iron(III) (hydr)oxide plays a critical role in coastal ecosystems experiencing seasonal hypoxia, because it intercepts sulfide that is released from the sediment to the bottom water, and hence reduces the risk of sulfide toxicity to fauna in seasonally stratified coastal systems (Seitaj et al., 2015). Accordingly, given this strong environmental impact of cable bacteria, it is adamant to accurately quantify both their abundance and activity in aquatic sediments.

Phylogenetic studies of the 16S rRNA gene of cable bacteria show affiliation to the family Desulfobulbaceae, in which cable bacteria form a sister clade of the genus Desulfobulbus (Pfeffer et al., 2012; Trojan et al., 2016; Kjeldsen et al., 2019). Two genera of cable bacteria have been identified in aquatic surface sediments. The genus Candidatus Electrothrix contains sequences originating from marine and salt marsh sites and comprises four species, Ca. Electrothrix aarhusiensis, E. communis, E. japonicus and E. marina. The genus Candidatus Electronema encompasses sequences of cable bacteria from freshwater sediments and comprises two species, Ca. Electronema nielsenii and En. palustris (Trojan et al., 2016). Cable bacteria have also been reported from a groundwater aquifer. These subsurface sequences were also affiliated to the Desulfobulbaceae, but clustered more closely with Desulfurivibrio alkaliphilus with 91% sequence identity (Müller et al., 2016, 2019). Currently, no cable bacteria isolates are available in pure culture.

The natural occurrence of cable bacteria has been documented in sediments from a range of marine and freshwater environments (Malkin et al., 2014; Risgaard-Petersen et al., 2015; Burdorf et al., 2017). However, the actual abundance of cable bacteria has only been quantified in a few locations. Quantification of cable bacteria abundance has thus far been performed using fluorescence in situ hybridization (FISH), and the DSB706 probe is most commonly used to quantify cable bacteria densities (Pfeffer et al., 2012; Malkin et al., 2014, 2017; Schauer et al., 2014; Risgaard-Petersen et al., 2015; Vasquez-Cardenas et al., 2015; Burdorf et al., 2016; van de Velde et al., 2016; Hermans et al., 2019). The DSB706 probe targets most members of the family Desulfobulbaceae and the genus Thermodesulforhabdus (Loy et al., 2002), and therefore, only the hybridized filaments and not the single cells are counted as cable bacteria. The ELF645 probe has been developed to more narrowly target the monophyletic clade of surface sediment cable bacteria (Pfeffer et al., 2012; Marzocchi et al., 2014), while additional species-level probes are targeting specific species within the Ca. Electrothrix and Ca. Electronema genera (Trojan et al., 2016; Marzocchi et al., 2018).

The quantification of cable bacteria by FISH has been shown to suffer from substantial variation, leading to standard errors of 50–80% when using three replicates (Risgaard-Petersen et al., 2015). This variation is likely caused by the small volume of sediment that is used in sample preparation to avoid high background of sediment particles and non-hybridized cells. Typically, only ∼1 μl of wet sediment is deposited onto a filter of 25 mm diameter, and this small amount of sample allows for strong variability in the FISH analysis between samples. An alternative approach for quantification is the use of qPCR using primers targeting the 16S rRNA gene of specific phyla or targeting specific functional genes (Agogue et al., 2008; Kato et al., 2009; Miyashita et al., 2009; Blazejak and Schippers, 2011; Pohlner et al., 2017).

Here we explore the use of qPCR as a method to quantify the abundance of cable bacteria in marine and salt marsh sediments. We used primers that are selective either for cable bacteria of the genus Ca. Electrothrix or for the larger Desulfobulbaceae family. We enumerated cable bacteria, Desulfobulbaceae and total bacteria in sediments from laboratory sediment incubations with and without the presence of cable bacteria. Also, we compared the qPCR results with abundance estimates obtained using 16S rRNA amplicon sequencing. Filament and cell densities of cable bacteria are subsequently combined with electrical current densities derived from a detailed assessment of the sediment geochemistry. This allows to estimate the electron flow and metabolic activity at the single filament and single cell level.

Materials and Methods

Sediment Collection and Laboratory Sediment Incubations

Sediment for laboratory incubations was collected from Lake Grevelingen (Netherlands), a coastal basin in which the presence of cable bacteria has been extensively documented (e.g., Malkin et al., 2014; Seitaj et al., 2015). The sampling site (Den Osse basin, site S1, 51.747°N, 3.890°E, salinity 33, water temperature 11.8°C) is at 34 m water depth and sediment was taken using gravity cores (60 mm inner diameter, 60 cm long; Uwitec, Austria).

Sediment was sieved (<500 μm) and homogenized, left to settle overnight and repacked into plexiglass core liners of 40 mm diameter with the sediment surface level at the top of the core liner. The sediment cores were incubated at a temperature of 16°C in containers containing artificial seawater at salinity 30 (Instant Ocean Sea Salt). The containers were kept dark to prevent photosynthesis and the overlying water was continuously aerated. Four cores were incubated in a container that only contained freshly prepared artificial seawater (denoted “unamended incubation”). In a second treatment (denoted the “stimulated incubation”), four cores were incubated in a container that already contained “inoculum cores” with active cable bacteria enrichments that were started 32 days earlier (Supplementary Figure S1).

Microsensor Profiling

Pore water concentration depth profiles of free sulfide (H2S), oxygen (O2) and pH were recorded using commercially available micro-electrodes (Unisense, Denmark) after 22 days of incubation. A single depth profile of H2S and pH, and duplicate profiles of O2, were taken in 2 cores from each treatment following the procedure outlined in van de Velde et al. (2016). Sensors were calibrated following the standard calibration procedures as described previously (Malkin et al., 2014) (H2S: 5 point standard curve using Na2S standards; O2: 2 point calibration using 100% in air bubbled seawater and the anoxic zone of the sediment; pH: 3 NBS standards (pH = 4, 7, 10) and TRIS buffer). The pH data are reported on the total pH scale and ΣH2S was calculated from H2S based on the pH values measured at the same depth. pH values are reported on the total scale.

Flux Calculations

The sedimentary oxygen consumption was calculated from the oxygen depth profiles using Fick’s 1st law: Jo2=ϕDO21-2ln(ϕ)d[O2]dx. Here, the diffusive oxygen flux (JO_2) across the sediment-water interface is calculated using the slope of the oxygen depth profile (d[O2]dx) and the molecular diffusion coefficient (DO_2) corrected for tortuosity (1−2ln⁡(ϕ)) and the mean porosity of the oxic zone (ϕ = 0.88). The molecular diffusion coefficient was calculated in R using the marelac package (Soetaert et al., 2010) taking into account the temperature (16°C) and salinity (30) of the incubations. Similarly the upward flux of ΣH2S was calculated using the slope of the ΣH2S depth profile using Fick’s first law at the base of the suboxic zone using the diffusion coefficient for H2S. Finally, upward fluxes were estimated from the iron pore water profiles using a similar approach and the molecular diffusion coefficient for ferrous iron (Fe+2). Flux data are presented as the mean ± SD with n = 4 for O2 and n = 2 for sulfide and ferrous iron fluxes.

Sediment Sampling and DNA Extraction

On day 24, two cores per treatment were transferred into an anaerobic glove box (N2 atmosphere, Coy Laboratory Products, United States) and the sediment was sectioned at 0.5 cm resolution. Each slice of sediment was centrifuged at 4000 × g for 7.5 min, and the extracted pore water was collected and filtered through 0.22 μm cellulose acetate filters. Pore water samples were immediately stabilized with 50 μl ml–1 bidistilled HNO3 (65%, Suprapure, Merck) and stored for ∼3 weeks at 4°C until further analysis of Ca, Fe and Mn using ICP-OES (precision was <2% for all analytes).

The two remaining cores per treatment were sectioned in 0.5 cm thick sediment slices in ambient air. Per slice, the sediment was homogenized and subsamples were flash frozen in liquid N2 and stored at −80°C until DNA extraction. Another subsample was fixed with an equal volume of ethanol (Ethanol for molecular biology, Merck) and stored at −20°C for analysis by fluorescence in situ hybridization (FISH).

DNA was extracted from a weighed amount (∼0.5 g) of wet sediment following an amended protocol of Zhou et al. (1996) that is effective for cable bacteria lysis and isolating DNA from sediment. The sediment was incubated with lysozyme (final concentration 2.5 mg ml–1) and RNase (100 μg ml–1) in DNA extraction buffer (0.1 M Tris–HCl, 0.1 M sodium EDTA, 0.1 M sodium phosphate, 1.5 M NaCl, 1% hexadecyltrimethylammonium bromide and 1% poly[vinylpolypyrrolidone] (pH 8)) at 37°C for 30 min, followed by addition of SDS (to a final concentration of 0.5%) and proteinase K (final concentration 350 μg ml–1) and 30 min incubation at 37°C. Subsequently, the SDS concentration in the extraction liquid was increased to 1% and the incubation continued at 60°C for 2 h. Samples were mixed by inversion every 10–15 min during the whole incubation period. DNA contained in the aqueous phase was purified by two extractions with chloroform/isoamyl alcohol (24:1 v/v) and precipitated with 1 volume of 10% polyethylene glycol (molecular mass ∼ 8000 g/mol) overnight at 4°C. The precipitated DNA was collected by centrifugation, washed twice with ice-cold 70% ethanol and dissolved in 10 mM Tris–HCl (pH 8) buffer.

qPCR Primers and Cycling Conditions

The abundance of (1) total bacteria, (2) bacteria belonging to the Desulfobulbaceae family, and (3) marine cable bacteria was estimated from the quantity of 16S rRNA gene copies determined by qPCR with target-specific primers (Table 1). Universal primers for total bacteria were Eub338 (Lane, 1991) and Eub518 (Muyzer et al., 1993). For quantification of Desulfobulbaceae, the forward primer DSBB280wF designed for the family Desulfobulbaceae was combined with the reverse primer SRB385R targeting sulfate-reducing bacteria (Amann et al., 1990; Sass et al., 1998). Primer DSBB280wF is an extension of primer DSBB280F targeting the family Desulfobulbaceae (Kjeldsen et al., 2007) allowing for a variation in the sequence (A instead of G at position 10) present in Ca. Electrothrix marina and Ca. E. aarhusiensis. For marine cable bacteria, forward primer ELF645wF was combined with reverse primer CB836wR. The ELF645wF primer is an extension of primer ELF645F targeting cable bacteria (Pfeffer et al., 2012; Marzocchi et al., 2014) to accommodate C instead of T at position 2 as present in the 16S rRNA gene sequence of Ca. E. marina. Reverse primer CB836wR is a newly designed primer that targets the genus Ca. Electrothrix, but also few Chloroflexi (see section “Results and Discussion”). Since this reverse primer is intended to be combined with forward primer ELF645wF it should have a similar melting temperature, in addition to the general requirements: GC content 40–60%, preferably start with G or C, no GC-rich end, no hairpin or self-dimer formation. The primer pair should be specific for cable bacteria and produce a relatively short amplicon (<200 bp) suitable for qPCR, and yield no hetero-dimer formation. A multiple alignment of the 16S rRNA gene sequence of cable bacteria (Trojan et al., 2016) and other Desulfobulbaceae (cf. Figure 4) was made with Muscle (Edgar, 2004) and manually inspected. Potential primer sequences specific for cable bacteria were identified and further assessed for specificity of the primer with Test Probe v 3.0, and specificity of the primer pair with TestPrime v 1.0, both against the Silva SSU reference database v 138 (Klindworth et al., 2013).

TABLE 1
www.frontiersin.org

Table 1. Primers used for qPCR and for PCR-amplification of qPCR standards.

qPCR analysis was performed with a realtime PCR analyzer (Bio-Rad CFX, Hercules, CA, United States and Rotor-Gene 6000, Corbett Research, Sydney, Australia) using ABsolute QPCR SYBR Green mix (Thermo Fisher scientific). The mix contained 3 mM MgCl2 and primer concentrations were 0.2 μM. Following hot start polymerase activation for 15 min at 95°C, the cycling conditions were: denaturation at 95°C for 15 s, annealing at 60°C (63°C for total bacteria) for 30 s, and extension at 72°C for 20 s. A series of 40 cycles was followed by 95°C for 30 s, then 60°C for 30 s and subsequently increasing the temperature to 95°C in steps of 0.5°C for 10 s/step, thus creating a melting curve for the obtained amplicons. Calibration curves were made using standards with at least 4 consecutive dilution steps (each step comprising a 10-time dilution). Sample DNA was analyzed in duplicate, and was diluted 104 times for quantification of total bacteria and 102 times for quantification of Desulfobulbaceae and cable bacteria. The detection limit of qPCR for Desulfobulbaceae and cable bacteria calculated using 102 dilution of sample DNA and 10 copies μl–1 detected in a qPCR run was approximately 100 16S copies g–1 wet sediment.

Standards for qPCR

Plasmid DNA with an inserted PCR-amplified partial 16S rRNA gene product was used to produce a calibration curve of 16S rRNA copy concentration vs. cycle quantification (Cq) value with the Desulfobulbaceae and cable bacteria primer sets. DNA from the incubated sediment was PCR amplified with primers ELF645F-CB836R (for the cable bacteria standard, CB) at an annealing temperature (TA) of 60°C or with DSBB280F-SRB385R (for the Desulfobulbaceae standard, DSB) at TA = 62°C using Taq DNA polymerase (New England Biolabs). Each PCR product was purified (EZNA Cycle Pure kit, Omega Bio-tek, Norcross, GA, United States) and ligated with the pCRTM 4-TOPO TA vector and cloned into One Shot TOP10 Competent E. coli cells (Invitrogen) according to the manufacturer’s instructions. Insert-containing E. coli colonies were transferred to liquid medium, incubated overnight at 37°C and the cells harvested for plasmid isolation. For 20 clones each with DSB and CB inserts, the insert was PCR amplified using primers M13F and M13R and the PCR product sequenced (BaseClear, Leiden, Netherlands). This showed incorporation of an insert with the correct qPCR primer sites. Plasmid DNA was linearized by digestion with PstI, purified and quantified by fluorescence (Qubit 3.0 and Qubit dsDNA HS assay kit; Life Technologies, Thermo Fisher Scientific).

For qPCR quantification of total bacteria, a PCR-amplified 16S rRNA gene product was used as the standard. DNA from incubated sediment (Lake Grevelingen) was amplified with general primers 27F (5′-AGA GTT TGA TCM TGG CTC AG-3′) and 1492R (5′-GGY TAC CTT GTT ACG ACT T-3′) using TA = 55°C. The product was purified (EZNA Cycle Pure kit) and the DNA concentration quantified by absorbance at 260 nm (NanoDrop, Thermo Fisher scientific).

Fluorescence in situ Hybridization

FISH analysis of cable bacteria using sediment from 0 to 0.5 cm depth was performed as described in Seitaj et al. (2015). One filter was prepared per sediment sample and the length of filaments hybridized to probe DSB706 was measured in 200 fields (105 × 141 μm).

16S rRNA Amplicon Sequence Analysis

The V4-V5 region of the 16S rRNA gene was amplified using DNA extracted from sediment with primers 515F-Y (5′-GTGYCAGCMGCCGCGGTAA-3′) and 926R (5′-CCGYCAATTYMTTTRAGTTT-3′) (Parada et al., 2016) with Illumina adapters added onto the target sequences. PCR was performed in triplicate using Q5 Hot start high-fidelity DNA polymerase (New England Biolabs). Cycling parameters were: initial denaturation at 98°C for 2 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 55°C for 20 s, extension at 72°C for 15 s, and final extension at 72°C for 2 min. Amplified DNA was checked for size on gel and triplicate products were pooled. Library preparation and sequencing (Illumina MiSeq, 2 × 300 bp) was performed at Eurofins Genomics, Ebersberg, Germany. The data for this study have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under accession number PRJEB37045.

Amplicon sequence variants (ASVs) were analyzed using the R-package dada2 (Callahan et al., 2016). Only the forward reads, covering the V4 region of the 16S rRNA gene, were included in the analysis. The reads were filtered and chimeras were removed following the default settings of the dada2 pipeline. Singletons were ignored because of the higher chance that these originate from sequencing errors. The taxonomy of ASVs was assigned using the Silva small subunit rRNA nr database v132 (Quast et al., 2013).

16S Phylogeny

16S rRNA gene sequences of cable bacteria, other Desulfobulbaceae and sequences of clones from sediments were aligned with the software Muscle (Edgar, 2004) and manually inspected. Phylogenetic trees were calculated using maximum likelihood implemented in RAxML (Stamatakis, 2014) with Gamma rate heterogeneity and 1000 bootstraps. The short V4 amplicon sequences were added using the evolutionary placement algorithm, also implemented in RAxML.

Filament Abundance of Cable Bacteria

The volumetric filament density NF, volume (i.e., the length of cable bacterium filaments per cm3 of wet sediment) is calculated from the 16S copy density CqPCR (i.e., the number of 16S copies per gram of wet sediment) determined by qPCR via the relation NF, volume = CqPCR × Lcell × ρw/ncopy, where Lcell is the mean length of an individual cell in a filament, ρw is the wet bulk density of the sediment (in g cm–3) and ncopy is the number of 16S rRNA gene copies present per cell. The aerial filament density NF,area (i.e., the length of cable bacteria filaments per cm2 of sediment surface) is subsequently calculated as NF,area = NF,volume × LCB, where LCB is the depth to which cable bacteria grow downward in the sediment.

Results and Discussion

qPCR Method Development: Primer Sets, Specificity and Calibration

Quantification of cable bacteria via qPCR requires an assessment of the performance of the primers that are used. Known cable bacteria sequences from surface sediments belong to two genera (Trojan et al., 2016): the genus Ca. Electrothrix encompassing the marine cable bacteria (4 species) and the genus Ca. Electronema, encompassing the freshwater cable bacteria (2 species). Ideal primer combinations for qPCR should allow annealing temperatures of 60°C or higher and generate a relatively small product (≤200 nucleotides). Based on these prerequisites, no suitable primer sequences were identified that could target both marine and freshwater cable bacteria with one single primer set. Therefore, we focused our quantification efforts on the genus Ca. Electrothrix that encompasses the marine cable bacteria.

The previously used primer ELF645F (Marzocchi et al., 2014) targets the genus Ca. Electrothrix, but for species Ca. E. marina, there is one mismatch at primer position 2. We allowed for this mismatch in the adapted primer ELF645wF. Analysis of primer ELF645wF using TestProbe showed that the primer targets only Ca. Electrothrix. When 1 mismatch is taken into account also other Desulfobulbaceae were targeted, with virtually all hits belonging to uncultured clones that were affiliated to the genus Desulfobulbus, comprising 9.2% of sequences in the Desulfobulbus lineage. We combined ELF645wF with a newly designed reverse primer CB836wR. Sequence CB836R matches to the species Ca. E. aarhusiensis, E. communis and E. japonica, and has 1 mismatch to Ca. E. marina that has been included in primer CB836wR. Analysis of CB836wR with TestProbe showed that it targets Ca. Electrothrix, and that there were also few perfect matches with uncultured clones affiliated to the order Ardenticatenales (Chloroflexi, 3 sequences) and to the genus Desulfobulbus (5 sequences). Primer CB836wR is a less specific primer compared to ELF645wF with in total 596 hits with 1 mismatch in the database, corresponding to sequences classified predominantly as Desulfobulbaceae, 55% of the lineage Desulfobulbus, but also Chloroflexi (7.4% of sequences affiliated to the order Ardenticatenales) and Gammaproteobacteria (10.2% of Cardiobacterales). The primer pair ELF645wF-CB836wR makes a more specific combination compared to the individual primers and is highly specific for Ca. Electrothrix. The primer combination does not target Chloroflexi and Gammaproteobacteria and, when allowing 1 mismatch for each primer, the only hits outside the Ca. Electrothrix lineage are 2 sequences affiliated to the lineage Desulfobulbus (out of 390 sequences affiliated to the genus Desulfobulbus in the database).

Members of the Desulfobulbaceae family were quantified using a specific primer targeting the family level [DSBB280F; (Kjeldsen et al., 2007)] that was extended to allow a mismatch with the 16S rRNA gene sequences of Ca. E. marina and E. aarhusiensis (DSBB280wF). The variant sequence had 31 hits in TestProbe, which were all affiliated to the Desulfobulbales. The forward primer DSBB280wF was combined with a more general primer for the class Deltaproteobacteria [SRB385R; (Sass et al., 1998)].

Quantitative PCR reactions using primer sets ELF645wF-CB936wR for cable bacteria (CB), DSB280wF-SRB385R for Desulfobulbaceae (DSB) and Eub338-Eub518 for total bacteria (TB) and the corresponding standards, all produced reproducible calibration curves with a calculated amplification efficiency of 92–104%, showing the suitability of the primer sets for qPCR-based enumeration. Hence, we used these primer sets for the quantification of CB, DSB, and TB in incubated marine sediment.

Metabolic Activity of Cable Bacteria in Sediment Incubations

The metabolic activity of cable bacteria can be detected through their strong imprint on the pore water geochemistry, which can be verified by microsensor depth profiling (Meysman et al., 2015). Depth profiles of O2 and H2S obtained after 22 days of incubation showed that in the unamended incubation there was no development of a suboxic zone (i.e., a zone where both O2 and H2S are below detection limit) (Figure 1A), and no subsurface pH maximum occurred, thus indicating no detectable metabolic activity of cable bacteria. In the stimulated incubation, H2S was only detectable in deeper sediment layers, resulting in a suboxic zone that was ∼10 mm wide (Figure 1E). The oxygen penetration depth (1.01 ± 0.16 mm, n = 4) was also shallower compared to the unamended incubation (1.60 ± 0.25 mm, n = 4), indicating higher O2 consumption. Microsensor depth profiles of pH showed a subsurface peak coinciding with the O2 penetration depth (pH values 0.16–0.31 units above the pH of the overlying water, n = 2) and a decreased pH at depth (minimum pH values 6.0–6.4, at approximately 1.0–1.4 cm depth in the sediment, n = 2). This particular combination of O2, H2S, and pH depth profiles functions as the unique fingerprint of cable bacteria activity (Pfeffer et al., 2012; Meysman et al., 2015), and results from the spatial segregation of the anodic sulfide oxidation and cathodic oxygen reduction half-reactions in electrogenic sulfur oxidation.

FIGURE 1
www.frontiersin.org

Figure 1. Pore water geochemistry of sediment in unamended incubations (upper row, open symbols) and stimulated incubations (bottom row, filled symbols). (A,E) Microsensor profiles of O2, ΣH2S and pH (t = 22 days). In the stimulated incubations a suboxic zone of ∼10 mm wide was present where both oxygen and sulfide were below detection limit, indicative of the activity of cable bacteria. (B–D,F–H) Pore water concentrations (t = 24 days) of Ca, Fe, and Mn were much higher in the stimulated compared to the unamended incubations, indicating dissolution of minerals stimulated by acidification of the suboxic zone through anodic sulfide oxidation by cable bacteria.

In each treatment, pore water Ca, Fe and Mn was measured in the same cores as the microsensor depth profiles. Pore water concentrations of Ca, Fe and Mn showed a strong increase in the stimulated incubation compared to the unamended incubation (Figures 1B–D,F–H). Previous studies examining the effect of cable bacteria on sediment geochemistry have shown that the subsurface accumulation of Ca, Mn, and Fe is linked to the dissolution of Ca and Mn-carbonates and Fe sulfide. This mineral dissolution is a direct response to the acidification of the pore water (Figure 1E), resulting from proton production by the sulfide oxidation half-reaction performed by cable bacteria (Risgaard-Petersen et al., 2012; Rao et al., 2016). Together the microsensor depth profiles and the pore water data revealed a substantial metabolic activity of cable bacteria in the stimulated incubation and similar sets of depth profiles have been recorded in sediments with high cable bacteria activity (van de Velde et al., 2016; Malkin et al., 2017). No sign of cable bacteria activity was apparent in the unamended incubation.

Together, these data suggest that the cable bacteria activity developed far more rapidly in the four sediment cores of the stimulated incubation compared to those of the unamended incubation (Figures 1, 2). In the stimulated incubation, “inoculum cores” were present that had already developed active cable bacteria, thus suggesting that a transfer occurred through the oxygenated overlying water. Such a transfer of cable bacteria via the water phase has not been documented previously, and the actual mechanism of transfer remains unclear. Yet, it suggests that cable bacteria can disperse through the water column and that they survive at least a short period in oxygenated water.

FIGURE 2
www.frontiersin.org

Figure 2. qPCR quantification of 16S rRNA gene copies of (A) total bacteria (TB), (B) Desulfobulbaceae (DSB) and (C) cable bacteria (CB). (D) CB 16S copies relative to DSB 16S copies and (E) CB 16S copies relative to TB 16S copies. Data is for two replicate cores each from unamended incubations (open symbols) and stimulated incubations (filled symbols). Data in panels (A–C) shows the mean with error bars denoting the standard deviation (n = 2). 16S rRNA gene copies are expressed per gram of wet sediment.

qPCR Quantification of Bacterial Abundances

Quantification of cable bacteria, Desulfobulbaceae and total bacteria via qPCR was performed for each incubation treatment. The two replicate cores showed very similar abundances for all clades and treatments investigated, thus demonstrating good reproducibility of qPCR results (Figure 2). Total bacteria 16S rRNA gene copy numbers (“TB 16S copies”) were very similar in the unamended and stimulated incubations and were largest in the top 0–0.5 cm, 1.3 ± 0.2 × 1010 (g wet sediment)–1, compared to 5.4 ± 0.7 × 109 (g wet sediment)–1 at >0.5 cm depth (Figure 2A). The increase of total bacteria in the top layer likely resulted from the growth of aerobic microorganisms at the sediment surface.

Quantification of Desulfobulbaceae in the unamended incubation showed a presence of approximately 1.5 ± 0.2 × 108 16S copies (g wet sediment)–1 below 1 cm depth (Figure 2B, open symbols). The DSB 16S copy number did not increase in the anoxic layers of the unamended incubation, indicating there was no enrichment of sulfate-reducing Desulfobulbaceae during the time of the incubation. In the stimulated incubation, the number of DSB 16S copies in the deepest layers (2–3 cm depth) was similar to the unamended incubation, but increased to 4.6 ± 0.8 × 108 (g wet sediment)–1 in the 0–1 cm depth layer (Figure 2B, filled symbols). This increase in Desulfobulbaceae coincides with the suboxic zone and can be attributed to the presence of cable bacteria (see below).

The deeper layers (>1 cm) of the unamended incubation contained 9.1 ± 2.0 × 104 cable bacteria 16S copies (g wet sediment)–1 (Figure 2C, open symbols), which we interpret to be the low background concentration in the homogenized sediments without an enriched and active population of cable bacteria. In one replicate of the unamended incubation, the CB 16S copy number remained constant throughout the depth profile, indicating cable bacteria had not developed. In the other replicate (core N1), CB 16S copies were 78 times higher in the upper 0–0.5 cm of the sediment [7.0 ± 0.3 × 106 (g wet sediment)–1], an increase that corresponds to ∼6 doubling times. Assuming the growth rate is the same for all cells in a cable bacterium filament and using the doubling time of ∼20 h as reported in Schauer et al. (2014), this indicates that cable bacteria development likely started ∼5 days before core sectioning.

In the stimulated incubation, the number of detected CB 16S copies was 40–2400 times higher compared to the unamended incubation, amounting to 2.4–3.2 × 108 16S copies (g wet sediment)–1 in the upper 0.5 cm of sediment. The CB 16S copy numbers decreased to 0.4–2.3 × 106 in the layer of 2.5–3 cm depth (Figure 2C, filled symbols). This confirms the accelerated development of cable bacteria in the stimulated incubation compared to the unamended incubation. The 16S CB copy number was constant within the suboxic zone (0–1.5 cm depth horizon), and rapidly declined below the depth at which free sulfide accumulated in the pore water (Figures 1E, 2C). A similar depth distribution has previously been observed in depth profiles of cable bacteria abundances quantified by FISH (Schauer et al., 2014; Vasquez-Cardenas et al., 2015; van de Velde et al., 2016).

Based on qPCR quantifications, CB 16S copies comprised up to 64% of DSB 16S copies in the stimulated incubation (Figure 2D). The difference between the DSB 16S copy numbers in the unamended and stimulated incubations matches the difference in CB 16S copy numbers between these treatments, demonstrating that the increase in the Desulfobulbaceae population is entirely due to the development of cable bacteria. Cable bacteria contributed up to 5.2% of the 16S copies of the total bacterial population within the suboxic zone (depth 0.5–1.5 cm), while a lower contribution of 1.6–2.7% is seen in the oxic zone (0–0.5 cm layer) due to the larger number of total bacteria in the oxic surface layer (Figures 2A,E). Studies of the microbial community composition of marine environments using 16S amplicon sequencing have found cable bacteria abundances up to 0.6% in Baltic Sea sediments (Klier et al., 2018), and up to 4.5% in sediment from a shallow lagoon in Aarhus Bay (Klier et al., 2018; Otte et al., 2018). The largest relative abundance at the Aarhus Bay site was detected in sediment sampled at 3–10 mm depth, consistent with the results here.

The observed abundance data illustrate once more that cable bacteria can grow to dense filament networks in aquatic sediments (Schauer et al., 2014; Vasquez-Cardenas et al., 2015; Malkin et al., 2017). If we assume that individual cells are Lcell = 3 μm long (based on FISH imaging, cf. Figure 5), and there are two 16S rRNA gene copies present per cell [as for the genome sequence of Ca. Electrothrix aarhusiensis; (Kjeldsen et al., 2019)], and sediment wet bulk density is ρw = 1.32 g cm–3 sediment, then our observed cable bacteria abundance of CqPCR = 3 × 108 16S copies (g wet sediment)–1 corresponds to a volumetric filament density of 600 m per bulk cm3 of sediment. This filament length density is at the high end of densities obtained with FISH under field conditions, which are in the range of 2–290 m per cm3 bulk sediment (Malkin et al., 2014; Burdorf et al., 2016; van de Velde et al., 2016; Marzocchi et al., 2018). However, in laboratory sediment incubations, similar to the ones conducted here, filament lengths of cable bacteria can amount up to 2800 m per cm3 bulk sediment (Schauer et al., 2014; Vasquez-Cardenas et al., 2015; Malkin et al., 2017).

After integration with depth, assuming that cable bacteria filaments extend to a depth of 1.5 cm (as derived from the qPCR data in Figure 2C), we obtain an aerial filament density of 900 m per cm2 of sediment surface area. This is again larger than FISH based values under field conditions, which range up to 480 m per cm2 (Risgaard-Petersen et al., 2015; Seitaj et al., 2015; Burdorf et al., 2016; van de Velde et al., 2016; Marzocchi et al., 2018; Hermans et al., 2019). Assuming that all filaments are straight and vertically oriented in the sediment (i.e., filaments are 1.5 cm long), we find an areal density of 6 × 104 filaments cm–2, which provides an average horizontal distance of 45 μm between filaments.

Cable Bacteria and Desulfobulbaceae Abundance Assessed by Amplicon Sequencing

Amplicon sequencing of the V4 region of the 16S rRNA gene yielded ≥52,000 sequences per sample (Supplementary Figure S2 shows the microbial community composition). Using the qPCR TB 16S copy number per core and depth layer, the calculated detection limit for the amplicon sequencing variant (ASV) abundance ranged between 5 × 104 and 4 × 105 16S copies (g wet sediment)–1. The ASV detection limit was in the same range as the qPCR CB 16S copy number in sediment layers where cable bacteria were not enriched during incubation, and therefore, amplicon sequencing could not be used to accurately assess the abundance of cable bacteria in most depth layers of the unamended incubation.

Amplicon sequences related to Ca. Electrothrix were detected, but no sequences were related to Ca. Electronema, consistent with the known ecology of the two genera. Ca. Electrothrix is found in marine environments, similar to the incubations here, while Ca. Electronema occurs in freshwater environments (Trojan et al., 2016). The abundance of cable bacteria relative to total bacteria assessed by qPCR and amplicon sequencing corresponded well for the sediment layers in the stimulated incubation (slope = 0.77, r2 = 0.99; Figure 3A). However, there was one notable exception. In core N1 of the unamended incubation, qPCR detected only ∼20% of the amplicon sequences that were assigned to Ca. Electrothrix (Figure 3A inset), which was caused by primer mismatch (as discussed below). In core N2 of the unamended incubation, no Electrothrix related sequences were detected, in agreement with qPCR results.

FIGURE 3
www.frontiersin.org

Figure 3. Ratio of the relative abundance of (A) cable bacteria and (B) Desulfobulbaceae quantified by qPCR and by classification of amplicon sequences. For panel (A), relative abundance of cable bacteria by qPCR (y-axis) was calculated as CB 16S copies/TB 16 copies (%), and relative abundance by amplicon sequencing (x-axis) as Σreads classified as Ca. Electrothrix/Σreads classified as Bacteria (%). For panel (B), calculations were similar, using DSB 16S copies, and reads classified as Desulfobulbaceae. Each data point shows the ratio for one depth layer in one of the replicate cores of the unamended incubation (open symbols) and stimulated incubation (filled symbols). The inset in panel (A) shows a close-up for unamended incubation N1. Ratios for depth layers with cable bacteria present are shown in green, and for depth layers without cable bacteria in red.

The relative abundance of Desulfobulbaceae estimated by classification of amplicon sequences and qPCR showed discrepancies in sediment layers without cable bacteria (Figure 3B). In layers without cable bacteria, qPCR overestimated Desulfobulbaceae by approximately 70% (slope = 1.72, r2 = 0.14), while in sediment layers with cable bacteria, qPCR and amplicon sequencing data showed good agreement (slope = 1.10, r2 = 0.89). It is likely that PCR amplification of other sequences than Desulfobulbaceae caused this overestimation. Since the sequenced V4 region is outside the region targeted by the DSB primers, we cannot evaluate if sequences that were classified as Desulfobulbaceae have mismatches with the primers used for qPCR. A TestProbe query of the DSBB280wF primer shows that there are a few thousand sequences with 1 mismatch present in the database that are affiliated to other families in the Deltaproteobacteria, and which hence may anneal with the reverse primer SRB385R. The SRB385R primer itself has >100,000 hits in the database. The use of a reverse primer that is (more) specific for the Desulfobulbaceae could likely decrease the incidence of misidentification.

Metabolic Activity of Cable Bacteria at the Single Filament and Single Cell Level

The convergence of amplicon sequencing and qPCR estimates suggests we can attain reliable estimates of densities of cells and filaments of cable bacteria in our sediment incubations. If we combine these data with estimates of the rate of long-distance electron transport, we can obtain insight into the electron flow at the single filament and single cell level. The electrical current density due to cable bacteria activity can be derived from either the rate of cathodic oxygen reduction or the rate of anodic sulfide oxidation (Meysman et al., 2015; van de Velde et al., 2016). Both rates can be estimated from the geochemical flux data (Figures 1, 2), and in theory, the resulting value for the current density should be the same (Meysman et al., 2015).

The sedimentary oxygen consumption was 35.1 ± 7.7 mmol O2 m–2 d–1 in the stimulated incubations with cable bacteria and 21.6 ± 3.7 mmol O2 m–2 d–1 in the unamended incubations with a low density of cable bacteria. The difference between both, 13.6 ± 11.4 mmol O2 m–2 d–1 or 39% of the total O2 uptake, can be attributed to cathodic oxygen reduction by cable bacteria, matching previous field studies, which estimate that 5–50% of total O2 uptake can be attributed to e-SOx (Malkin et al., 2014; van de Velde et al., 2016). Knowing that 4 electrons are transferred to each O2 molecule during cathodic O2 reduction, this provides a current density of 61 ± 51 mA m–2 (1A = 1.036 × 10–5 mol e s–1). Current densities reported for different field studies range from 3.8 to 47 mA m–2 (Malkin et al., 2014, 2017; van de Velde et al., 2016, 2017) and go up to 74–96 mA m–2 for sediment incubations under laboratory conditions (Risgaard-Petersen et al., 2012; Schauer et al., 2014; Rao et al., 2016; van de Velde et al., 2017). Accordingly, our observed current density is at the high end of the field observations, and at the low end of the laboratory experiments.

A sulfide balance for the stimulated incubation allows a second estimate of the current density. There is a sulfide flux of 4.8 ± 1.3 mmol H2S m–2 d–1 from the deeper sulfidic zone toward the suboxic zone (Figure 1E), which is generated by sulfate reduction in deeper sediment (sediment zone = 59 mm wide). If we assume that the volumetric rate of sulfate reduction is similar in the suboxic zone (sediment zone = 10 mm wide), we can estimate that sulfate reduction produces 0.8 ± 0.2 [ = (10/59) × 4.8] mmol H2S m–2 d–1 in the suboxic zone. Anodic sulfide oxidation is the sum of consumption of sulfide generated by sulfate reduction in the suboxic zone (0.8 mmol H2S m–2 d–1), sulfide generated by sulfate reduction in the deeper sediment (4.8 mmol H2S m–2 d–1) and sulfide generated by net FeS dissolution (−2.0 ± 1.3 mmol H2S m–2 d–1 derived from the upward flux of dissolved Fe; the downward Fe flux is not accounted for as we assume this precipitates as FeS). Anodic oxidation of sulfide (total rate 7.6 mmol H2S m–2 d–1) produces 8 electrons when sulfide is oxidized to sulfate, and hence, the current density is estimated at 68 ± 25 mA m–2, which closely matches our estimate based on cathodic O2 consumption.

Based on the above values for the current density (61–68 mA m–2) and the aerial cable bacteria density (6 × 108 filaments m–2), we can calculate the metabolic activity of an average individual cable bacterium filament as 100−115 pA per filament. When cable bacteria were connected directly with electrodes with an applied voltage bias of 100 mV, measured currents ranged from 6 to 1350 pA per cable bacterium filament (Meysman et al., 2019). It should be noted, however, that a direct comparison is difficult, as the distance over which currents were transmitted was different. Nevertheless, the magnitude of currents observed under in situ conditions matches that of previous in silico electrical measurements.

The pathway by which oxygen is reduced in cable bacteria is currently unknown. Recent 13C and 15N labeling studies have shown that cells in the oxic zone of the sediment show very little or no carbon or nitrogen assimilation, implying that oxygen reduction in cable bacteria is not coupled to energy conservation (Geerlings et al., 2020). Moreover, no terminal oxidases could be identified in cable bacteria genomes and it was hypothesized that cable bacteria reduce oxygen via periplasmic cytochromes (Kjeldsen et al., 2019).

The cell specific oxygen consumption rate of cable bacteria in the oxic zone calculated from the cathodic oxygen consumption, the oxygen penetration depth (∼1 mm; Figure 1E) and qPCR enumerated cable bacteria is 69 ± 58 fmol O2 cell–1 day–1, which is within the same order of magnitude as the value of 36 fmol O2 cell–1 day–1 reported by Schauer et al. (2014). Measurements of the oxygen reduction capacity of cable bacteria using cyclic voltammetry resulted in a rate of 6.1 pmol O2 cell–1 day–1 (at 50 μM O2; Geerlings et al., 2020), hence the potential capacity for oxygen reduction is approximately a hundred times larger than the actual rates in sediment incubations.

Cable Bacteria Diversity and qPCR Quantification

Up to present, four different candidate species (communis, aarhusiensis, japonica and marina) have been identified within the marine cable bacteria genus Ca. Electrothrix (Trojan et al., 2016), and field research shows that these different species of cable bacteria may co-exist within the same sediment environment (Marzocchi et al., 2018). In the unamended and stimulated incubations together, four amplicon sequence variants (ASV1-ASV4) were identified that classified as Ca. Electrothrix. ASV1 showed 100% sequence identity to both Ca. E. communis and E. aarhusiensis, ASV2 was most closely related to E. marina with 97.2% identity, ASV3 shared 100% sequence identity with E. japonica, and ASV4 shared 97.2% sequence identity with both E. communis and E. aarhusiensis. Placement of the ASVs in a phylogenetic tree of 16S rRNA sequences of cable bacteria, other Desulfobulbaceae and sediment clones (Figure 4) showed that ASV1 and ASV3 clustered with the currently described species of marine cable bacteria. Placement of ASV1 with Ca. E. communis and E. aarhusiensis is equally likely, since the sequences have identical V4 regions. ASV2 was placed with a group of clones that clustered more distant from the described Ca. Electrothrix species, but still in the Electrothrix lineage. ASV4 was placed with another group of sediment clones, which clustered outside of this lineage, however, with low bootstrap support. The sediment clone sequences included in the tree were obtained from sediment incubations with cable bacteria, using sediment from Aarhus Bay, Denmark (Pfeffer et al., 2012; Schauer et al., 2014). Sequence identity between ASV2 and ASV4 and the related clones is 99–100%.

FIGURE 4
www.frontiersin.org

Figure 4. Phylogenetic tree of 16S rRNA gene sequences of cable bacteria, cultured members of the Desulfobulbaceae and selected clones from sediment (the origin of the sediments is indicated in square brackets) and the amplicon sequence variants (ASVs) classified as Ca. Electrothrix detected in the incubations of marine sediment (in bold). The ASV sequences (indicated with * followed by the length of the sequence) were inserted using the evolutionary placement algorithm. Bootstrap values (%) were calculated for 1000 re-samplings; only values ≥60 are shown.

Quantitative PCR is a relatively rapid and low-cost method for quantification of microbial populations, but its success depends inherently on the conservation of primer sites in the targeted populations. The sequenced V4 region contains the ELF645wF primer site, which showed full sequence identity for ASV1 and ASV3, but two mismatches for ASV4 (positions 1 and 3) and three mismatches for ASV2 (positions 2, 3 and 5). The ASV1 contributed up to 7.0% of bacterial V4 sequences, and ASV3 up to 0.62%. The ASVs 2 and 4 made up small proportions of sequences classified as Electrothrix in sediments of the stimulated incubation. However, in core N1 from the unamended incubation, ASV2 was the only variant detected. The presence of sequence mismatches in the ELF645wF primer site of ASV2 explains the underestimation of cable bacteria by qPCR in incubation N1 (Figure 3A).

So the question arises whether the ASVs 2 and 4 are indeed cable bacteria? FISH analysis of sediment from the N1 incubation (which only has ASV2) showed hybridization of filamentous bacteria to probe DSB706, supporting that ASV2 belongs to cable bacteria (Figure 5). The placement of sequence ASV2 in the cable bacteria lineage in the 16S rRNA phylogenetic tree (Figure 4) also supports that this ASV is representative of a cable bacterium strain. For sequence ASV4 and also clones previously obtained from marine sediment incubations (Figure 4; Pfeffer et al., 2012; Schauer et al., 2014), it is yet unclear if they belong to cable bacteria.

FIGURE 5
www.frontiersin.org

Figure 5. FISH images of filamentous bacteria hybridized with DSB706 from sediment of core N1 of the unamended incubation. The images show the overlay of false colored images for the DSB706 (red) and DAPI (blue) signal. The scale bars denote 10 μm.

For sediment incubations with nitrate, Marzocchi et al. (2014) observed that of bacterial filaments that hybridized with probe DSB706, only 21 ± 13% also hybridized with probe ELF645. FISH probe ELF645 targets the same sequence as primer ELF645F. Hence, based on the observations of Marzocchi et al. (2014), the diversity of marine cable bacteria appears larger than currently identified, which could lead to an underestimation of qPCR-based abundance data of cable bacteria. Therefore, a combination of different primers targeting different subgroups of cable bacteria may be needed for qPCR quantification, or primers for other targets, e.g., the 23S rRNA sequence, or genes with a function in the metabolism of cable bacteria.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article.

Author Contributions

JG and FM designed the study. JG designed the qPCR primers, performed qPCR analysis, and analyzed the amplicon sequencing data. SV performed the sediment incubations and the pore water analysis. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) [ERC Grant 306933 to FM], the Research Foundation Flanders [FWO project grant G031416N], the Netherlands Organisation for Scientific Research [VICI grant 016.VICI.170.072 to FM], the NASA Postdoctoral Program [post-doctoral fellowship to SV], and the Flanders Marine Institute [Brilliant Marine Research Idea grant].

Conflict of Interest

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.

Acknowledgments

The authors thank Anton Tramper for assistance with the sediment incubation experiments, and Silvia Hidalgo-Martinez and Francis Criens for the FISH analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01506/full#supplementary-material

References

Agogue, H., Brink, M., Dinasquet, J., and Herndl, G. J. (2008). Major gradients in putatively nitrifying and non-nitrifying Archaea in the deep North Atlantic. Nature 456, 788–U772.

Google Scholar

Amann, R. I., Binder, B. J., Olson, R. J., Chisholm, S. W., Devereux, R., and Stahl, D. A. (1990). Combination of 16S ribosomal-RNA-targeted oligonucleotide probes with flow-cytometry for analyzing mixed microbial-populations. Appl. Environ. Microbiol. 56, 1919–1925. doi: 10.1128/aem.56.6.1919-1925.1990

CrossRef Full Text | Google Scholar

Bjerg, J. T., Boschker, H. T. S., Larsen, S., Berry, D., Schmid, M., Millo, D., et al. (2018). Long-distance electron transport in individual, living cable bacteria. Proc. Natl. Acad. Sci. U.S.A. 115, 5786–5791. doi: 10.1073/pnas.1800367115

PubMed Abstract | CrossRef Full Text | Google Scholar

Blazejak, A., and Schippers, A. (2011). Real-time PCR quantification and diversity analysis of the functional genes aprA and dsrA of sulfate-reducing prokaryotes in marine sediments of the Peru Continental Margin and the Black Sea. Front. Microbiol. 2:253. doi: 10.3389/fmicb.2011.00253

PubMed Abstract | CrossRef Full Text | Google Scholar

Burdorf, L. D. W., Hidalgo-Martinez, S., Cook, P. L. M., and Meysman, F. J. R. (2016). Long-distance electron transport by cable bacteria in mangrove sediments. Mar. Ecol. Progr. Ser. 545, 1–8. doi: 10.3354/meps11635

CrossRef Full Text | Google Scholar

Burdorf, L. D. W., Tramper, A., Seitaj, D., Meire, L., Hidalgo-Martinez, S., Zetsche, E.-M., et al. (2017). Long-distance electron transport occurs globally in marine sediments. Biogeosciences 14, 683–701. doi: 10.5194/bg-14-683-2017

CrossRef Full Text | Google Scholar

Callahan, B. J., Mcmurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13:581. doi: 10.1038/nmeth.3869

PubMed Abstract | CrossRef Full Text | Google Scholar

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32, 1792–1797. doi: 10.1093/nar/gkh340

PubMed Abstract | CrossRef Full Text | Google Scholar

Geerlings, N. M. J., Karman, C., Trashin, S., As, K. S., Kienhuis, M. V. M., Hidalgo-Martinez, S., et al. (2020). Division of labor and growth during electrical cooperation in multicellular cable bacteria. Proc. Natl. Acad. Sciences 117, 5478–5485. doi: 10.1073/pnas.1916244117

PubMed Abstract | CrossRef Full Text | Google Scholar

Hermans, M., Lenstra, W. K., Hidalgo-Martinez, S., Van Helmond, N., Witbaard, R., Meysman, F. J. R., et al. (2019). Abundance and biogeochemical impact of cable bacteria in Baltic Sea sediments. Environ. Sci. Technol. 53, 7494–7503. doi: 10.1021/acs.est.9b01665

PubMed Abstract | CrossRef Full Text | Google Scholar

Kato, S., Kobayashi, C., Kakegawa, T., and Yamagishi, A. (2009). Microbial communities in iron-silica-rich microbial mats at deep-sea hydrothermal fields of the Southern Mariana Trough. Environ. Microbiol. 11, 2094–2111. doi: 10.1111/j.1462-2920.2009.01930.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kjeldsen, K. U., Loy, A., Jakobsen, T. F., Thomsen, T. R., Wagner, M., and Ingvorsen, K. (2007). Diversity of sulfate-reducing bacteria from an extreme hypersaline sediment, Great Salt Lake (Utah). FEMS Microbiol. Ecol. 60, 287–298. doi: 10.1111/j.1574-6941.2007.00288.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Kjeldsen, K. U., Schreiber, L., Thorup, C. A., Boesen, T., Bjerg, J. T., Yang, T. T., et al. (2019). On the evolution and physiology of cable bacteria. Proc. Natl. Acad. Sci. U.S.A. 116, 19116–19125.

Google Scholar

Klier, J., Dellwig, O., Leipe, T., Jurgens, K., and Herlemann, D. P. R. (2018). Benthic bacterial community composition in the oligohaline-marine transition of surface sediments in the Baltic Sea based on rRNA analysis. Front. Microbiol. 9:236. doi: 10.3389/fmicb.2018.00236

PubMed Abstract | CrossRef Full Text | Google Scholar

Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2013). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41:e1. doi: 10.1093/nar/gks808

PubMed Abstract | CrossRef Full Text | Google Scholar

Lane, D. J. (1991). “16S/23S rRNA sequencing,” in Nucleic Acid Techniques in Bacterial Systematics, ed. E. G. Stackebrandt (Chichester: John Wiley & Sons), 115–175.

Google Scholar

Loy, A., Lehner, A., Lee, N., Adamczyk, J., Meier, H., Ernst, J., et al. (2002). Oligonucleotide microarray for 16S rRNA gene-based detection of all recognized lineages of sulfate-reducing prokaryotes in the environment. Appl. Environ. Microbiol. 68, 5064–5081. doi: 10.1128/aem.68.10.5064-5081.2002

PubMed Abstract | CrossRef Full Text | Google Scholar

Malkin, S. Y., Rao, A. M. F., Seitaj, D., Vasquez-Cardenas, D., Zetsche, E.-M., Hidalgo-Martinez, S., et al. (2014). Natural occurrence of microbial sulphur oxidation by long-range electron transport in the seafloor. ISME J. 8, 1843–1854. doi: 10.1038/ismej.2014.41

PubMed Abstract | CrossRef Full Text | Google Scholar

Malkin, S. Y., Seitaj, D., Burdorf, L. D. W., Nieuwhof, S., Hidalgo-Martinez, S., Tramper, A., et al. (2017). Electrogenic sulfur oxidation by cable bacteria in bivalve reef sediments. Front. Mar. Sci. 4:28. doi: 10.3389/fmars.2017.00028

CrossRef Full Text | Google Scholar

Marzocchi, U., Bonaglia, S., Van De Velde, S., Hall, P. O. J., Schramm, A., Risgaard-Petersen, N., et al. (2018). Transient bottom water oxygenation creates a niche for cable bacteria in long-term anoxic sediments of the Eastern Gotland Basin. Environ. Microbiol. 20, 3031–3041. doi: 10.1111/1462-2920.14349

PubMed Abstract | CrossRef Full Text | Google Scholar

Marzocchi, U., Trojan, D., Larsen, S., Meyer, R. L., Revsbech, N. P., Schramm, A., et al. (2014). Electric coupling between distant nitrate reduction and sulfide oxidation in marine sediment. ISME J. 8, 1682–1690. doi: 10.1038/ismej.2014.19

PubMed Abstract | CrossRef Full Text | Google Scholar

Meysman, F. J. R. (2018). Cable bacteria take a new breath using long-distance electricity. Trends Microbiol. 26, 411–422. doi: 10.1016/j.tim.2017.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Meysman, F. J. R., Cornelissen, R., Trashin, S., Bonne, R., Martinez, S. H., Van Der Veen, J., et al. (2019). A highly conductive fibre network enables centimetre-scale electron transport in multicellular cable bacteria. Nat. Commun. 10:8.

Google Scholar

Meysman, F. J. R., Risgaard-Petersen, N., Malkin, S. Y., and Nielsen, L. P. (2015). The geochemical fingerprint of microbial long-distance electron transport in the seafloor. Geochim. Cosmochim. Acta 152, 122–142. doi: 10.1016/j.gca.2014.12.014

CrossRef Full Text | Google Scholar

Miyashita, A., Mochimaru, H., Kazama, H., Ohashi, A., Yamaguchi, T., Nunoura, T., et al. (2009). Development of 16S rRNA gene-targeted primers for detection of archaeal anaerobic methanotrophs (ANMEs). FEMS Microbiol. Lett. 297, 31–37. doi: 10.1111/j.1574-6968.2009.01648.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, H., Bosch, J., Griebler, C., Damgaard, L. R., Nielsen, L. P., Lueders, T., et al. (2016). Long-distance electron transfer by cable bacteria in aquifer sediments. ISME J. 10, 2010–2019. doi: 10.1038/ismej.2015.250

PubMed Abstract | CrossRef Full Text | Google Scholar

Müller, H., Marozava, S., Probst, A. J., and Meckenstock, R. U. (2019). Groundwater cable bacteria conserve energy by sulfur disproportionation. ISME J. 14, 1–12.

Google Scholar

Muyzer, G., Dewaal, E. C., and Uitterlinden, A. G. (1993). Profiling of complex microbial-populations by denaturing gradient gel-electrophoresis analysis of polymerase chain reaction-amplified genes-coding for 16S ribosomal-RNA. Appl. Environ. Microbiol. 59, 695–700. doi: 10.1128/aem.59.3.695-700.1993

CrossRef Full Text | Google Scholar

Nielsen, L. P., Risgaard-Petersen, N., Fossing, H., Christensen, P. B., and Sayama, M. (2010). Electric currents couple spatially separated biogeochemical processes in marine sediment. Nature 463, 1071–1074. doi: 10.1038/nature08790

PubMed Abstract | CrossRef Full Text | Google Scholar

Otte, J. M., Harter, J., Laufer, K., Blackwell, N., Straub, D., Kappler, A., et al. (2018). The distribution of active iron-cycling bacteria in marine and freshwater sediments is decoupled from geochemical gradients. Environ. Microbiol. 20, 2483–2499. doi: 10.1111/1462-2920.14260

PubMed Abstract | CrossRef Full Text | Google Scholar

Parada, A. E., Needham, D. M., and Fuhrman, J. A. (2016). Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414. doi: 10.1111/1462-2920.13023

PubMed Abstract | CrossRef Full Text | Google Scholar

Pfeffer, C., Larsen, S., Song, J., Dong, M., Besenbacher, F., Meyer, R. L., et al. (2012). Filamentous bacteria transport electrons over centimetre distances. Nature 491, 218–221. doi: 10.1038/nature11586

PubMed Abstract | CrossRef Full Text | Google Scholar

Pohlner, M., Degenhardt, J., Von Hoyningen-Huene, A. J. E., Wemheuer, B., Erlmann, N., Schnetger, B., et al. (2017). The biogeographical distribution of benthic Roseobacter group members along a Pacific transect is structured by nutrient availability within the sediments and primary production in different oceanic provinces. Front. Microbiol. 8:2550. doi: 10.3389/fmicb.2017.02550

PubMed Abstract | CrossRef Full Text | Google Scholar

Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596.

Google Scholar

Rao, A. M. F., Malkin, S. Y., Hidalgo-Martinez, S., and Meysman, F. J. R. (2016). The impact of electrogenic sulfide oxidation on elemental cycling and solute fluxes in coastal sediment. Geochim. Cosmochim. Acta 172, 265–286. doi: 10.1016/j.gca.2015.09.014

CrossRef Full Text | Google Scholar

Risgaard-Petersen, N., Kristiansen, M., Frederiksen, R. B., Dittmer, A. L., Bjerg, J. T., Trojan, D., et al. (2015). Cable bacteria in freshwater sediments. Appl. Environ. Microbiol. 81, 6003–6011.

Google Scholar

Risgaard-Petersen, N., Revil, A., Meister, P., and Nielsen, L. P. (2012). Sulfur, iron-, and calcium cycling associated with natural electric currents running through marine sediment. Geochim. Cosmochim. Acta 92, 1–13. doi: 10.1016/j.gca.2012.05.036

CrossRef Full Text | Google Scholar

Sandfeld, T., Marzocchi, U., Petro, C., Schramm, A., and Risgaard-Petersen, N. (2020). Electrogenic sulfide oxidation mediated by cable bacteria stimulates sulfate reduction in freshwater sediments. ISME J. 14, 1233–1246. doi: 10.1038/s41396-020-0607-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Sass, H., Wieringa, E., Cypionka, H., Babenzien, H. D., and Overmann, J. (1998). High genetic and physiological diversity of sulfate-reducing bacteria isolated from an oligotrophic lake sediment. Arch. Microbiol. 170, 243–251. doi: 10.1007/s002030050639

PubMed Abstract | CrossRef Full Text | Google Scholar

Schauer, R., Risgaard-Petersen, N., Kjeldsen, K. U., Bjerg, J. J. T., Jorgensen, B. B., Schramm, A., et al. (2014). Succession of cable bacteria and electric currents in marine sediment. ISME J. 8, 1314–1322. doi: 10.1038/ismej.2013.239

PubMed Abstract | CrossRef Full Text | Google Scholar

Seitaj, D., Schauer, R., Sulu-Gambari, F., Hidalgo-Martinez, S., Malkin, S. Y., Burdorf, L. D. W., et al. (2015). Cable bacteria generate a firewall against euxinia in seasonally hypoxic basins. Proc. Natl. Acad. Sci. U.S.A. 112, 13278–13283. doi: 10.1073/pnas.1510152112

PubMed Abstract | CrossRef Full Text | Google Scholar

Soetaert, K., Petzoldt, T., and Meysman, F. (2010). Marelac: Tools for Aquatic Sciences, R Package Version.

Google Scholar

Stamatakis, A. (2014). RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 1312–1313. doi: 10.1093/bioinformatics/btu033

PubMed Abstract | CrossRef Full Text | Google Scholar

Sulu-Gambari, F., Seitaj, D., Meysman, F. J. R., Schauer, R., Polerecky, L., and Slomp, C. P. (2016). Cable bacteria control iron-phosphorus dynamics in sediments of a coastal hypoxic basin. Environ. Sci. Technol. 50, 1227–1233. doi: 10.1021/acs.est.5b04369

PubMed Abstract | CrossRef Full Text | Google Scholar

Trojan, D., Schreiber, L., Bjerg, J. T., Boggild, A., Yang, T. T., Kjeldsen, K. U., et al. (2016). A taxonomic framework for cable bacteria and proposal of the candidate genera Electrothrix and Electronema. Syst. Appl. Microbiol. 39, 297–306. doi: 10.1016/j.syapm.2016.05.006

PubMed Abstract | CrossRef Full Text | Google Scholar

van de Velde, S., Callebaut, I., Gao, Y., and Meysman, F. J. R. (2017). Impact of electrogenic sulfur oxidation on trace metal cycling in a coastal sediment. Chem. Geol. 452, 9–23. doi: 10.1016/j.chemgeo.2017.01.028

CrossRef Full Text | Google Scholar

van de Velde, S., Lesven, L., Burdorf, L. D. W., Hidalgo-Martinez, S., Geelhoed, J. S., Van Rijswijk, P., et al. (2016). The impact of electrogenic sulfur oxidation on the biogeochemistry of coastal sediments: a field study. Geochim. Cosmochim. Acta 194, 211–232. doi: 10.1016/j.gca.2016.08.038

CrossRef Full Text | Google Scholar

Vasquez-Cardenas, D., Van De Vossenberg, J., Polerecky, L., Malkin, S. Y., Schauer, R., Hidalgo-Martinez, S., et al. (2015). Microbial carbon metabolism associated with electrogenic sulphur oxidation in coastal sediments. ISME J. 9, 1966–1978. doi: 10.1038/ismej.2015.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J. Z., Bruns, M. A., and Tiedje, J. M. (1996). DNA recovery from soils of diverse composition. Appl. Environ. Microbiol. 62, 316–322. doi: 10.1128/aem.62.2.316-322.1996

CrossRef Full Text | Google Scholar

Keywords: cable bacteria, quantitative PCR, Desulfobulbaceae, marine sediment, amplicon sequencing, current density, oxygen consumption rate

Citation: Geelhoed JS, van de Velde SJ and Meysman FJR (2020) Quantification of Cable Bacteria in Marine Sediments via qPCR. Front. Microbiol. 11:1506. doi: 10.3389/fmicb.2020.01506

Received: 30 March 2020; Accepted: 10 June 2020;
Published: 03 July 2020.

Edited by:

Tony Gutierrez, Heriot-Watt University, United Kingdom

Reviewed by:

Luke McKay, Montana State University, United States
Tingting Yang, University of Southern California, Los Angeles, United States

Copyright © 2020 Geelhoed, van de Velde and Meysman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Jeanine S. Geelhoed, jeanine.geelhoed@uantwerpen.be; Filip J. R. Meysman, f.j.r.meysman@tudelft.nl