Original Research ARTICLE
Fuerstia marisgermanicae gen. nov., sp. nov., an Unusual Member of the Phylum Planctomycetes from the German Wadden Sea
- 1Leibniz Institut Deutsche Sammlung Von Mikroorganismen und Zellkulturen, Braunschweig, Germany
- 2School of Biology, Newcastle University, Newcastle, UK
- 3Helmholtz Centre for Infectious Disease, Braunschweig, Germany
Members of the phylum Planctomycetes are ubiquitous bacteria that dwell in aquatic and terrestrial habitats. While planctomycetal species are important players in the global carbon and nitrogen cycle, this phylum is still undersampled and only few genome sequences are available. Here we describe strain NH11T, a novel planctomycete obtained from a crustacean shell (Wadden Sea, Germany). The phylogenetically closest related cultivated species is Gimesia maris, sharing only 87% 16S rRNA sequence identity. Previous isolation attempts have mostly yielded members of the genus Rhodopirellula from water of the German North Sea. On the other hand, only one axenic culture of the genus Pirellula was obtained from a crustacean thus far. However, the 16S rRNA gene sequence of strain NH11T shares only 80% sequence identity with the closest relative of both genera, Rhodopirellula and Pirellula. Thus, strain NH11T is unique in terms of origin and phylogeny. While the pear to ovoid shaped cells of strain NH11T are typical planctomycetal, light-, and electron microscopic observations point toward an unusual variation of cell division through budding: during the division process daughter- and mother cells are connected by an unseen thin tubular-like structure. Furthermore, the periplasmic space of strain NH11T was unusually enlarged and differed from previously known planctomycetes. The complete genome of strain NH11T, with almost 9 Mb in size, is among the largest planctomycetal genomes sequenced thus far, but harbors only 6645 protein-coding genes. The acquisition of genomic components by horizontal gene transfer is indicated by the presence of numerous putative genomic islands. Strikingly, 45 “giant genes” were found within the genome of NH11T. Subsequent analysis of all available planctomycetal genomes revealed that Planctomycetes as such are especially rich in “giant genes”. Furthermore, Multilocus Sequence Analysis (MLSA) tree reconstruction support the phylogenetic distance of strain NH11T from other cultivated Planctomycetes of the same phylogenetic cluster. Thus, based on our findings, we propose to classify strain NH11T as Fuerstia marisgermanicae gen. nov., sp. nov., with the type strain NH11T, within the phylum Planctomycetes.
Planctomycetes are a phylum of ubiquitous, environmentally important bacteria that play key roles in global carbon and nitrogen cycles (Fuerst and Sagulenko, 2011; Kartal et al., 2013). Together with the Verrucomicrobia and the Chlamydia, planctomycetes belong to the PVC superphylum (Wagner and Horn, 2006). While Chlamydia are intracellular pathogens (Stephens et al., 1998), Verrucomicrobia (Petroni et al., 2000), and Planctomycetes (Fuerst et al., 1997; Lage and Bondoso, 2014) can be closely associated with eukaryotes as well. The Planctomycetes are divided into four distinct orders. While the orders Phycispherales, Tepidisphaerales, and Planctomycetales are based on axenic cultures, the order Brocadiales is formed by well described enrichment cultures, the so-called anammox-Planctomycetes. These organisms are capable of anaerobic ammonium oxidation, a trait extremely useful for wastewater treatment (Kartal et al., 2013). In particular, members of the order Planctomycetales were found to encode numerous secondary metabolite-related genes and gene clusters (Jeske et al., 2013), that were recently found to be active under conditions that chemically mimicked the interaction with eukaryotes (Jeske et al., 2016). Consequently, Planctomycetes were postulated to be “talented producers” of small molecules and they might represent a yet untapped resource of novel bioactive molecules (Jeske et al., 2016). Furthermore, the biotechnological application of planctomycetal enzymes such as sulfatases as biocatalysts was demonstrated (Wallner et al., 2005). Thus, Planctomycetes are environmentally important and of general biotechnological interest.
Planctomycetes were further proposed to comprise conspicuous cell biological features such as endocytosis-like uptake of proteins (Fuerst and Sagulenko, 2010) and compartmentalization of their cytosol (Fuerst, 2005). In addition, they were believed to lack peptidoglycan (PG) in their cell walls (König et al., 1984). Recently some of these unique features were questioned. For example, the presence of PG was demonstrated (Jeske et al., 2015; van Teeseling et al., 2015). However, in particular the unusual cell division of Planctomycetales through polar budding makes them unique among bacteria, as they lack the otherwise universal bacterial cell division protein FtsZ (van Niftrik et al., 2009; Jogler et al., 2012).
Despite their importance for environmental microbiology, biotechnology and cell biology, only 30 planctomycetal species were obtained as axenic cultures and only 10 completely closed genome sequences are available through NCBI GenBank. Thus, from a phylogenetic point of view, the phylum Planctomycetes is heavily undersampled and only few representatives of this phylum are taxonomically characterized in detail (Ward, 2010; Fuerst and Sagulenko, 2011).
In this study, we selectively target eukaryote-associated species from a marine habitat. We revisit the German North Sea, which is a well-known resource for the cultivation of Planctomycetes (Glöckner et al., 2003; Winkelmann and Harder, 2009). To focus on potentially eukaryote-associated species, we collected crab shells as the successful isolation of novel Planctomycetes from Crustacea was previously reported (Fuerst et al., 1997). As planctomycetes are known to be resistant to multiple antibiotics such as carbenicillin (Cayrou et al., 2010; Aghnatios et al., 2015), we employed a 96 well-based high-throughput cultivation procedure and subsequently screened hundreds of carbenicillin-resistant bacterial cultures with a PCR targeting the 16S rRNA gene. The obtained strain NH11T was selected for further analysis and full-length 16S rRNA gene sequencing revealed its affiliation with the phylum Planctomycetes. Here we show how strain NH11T differs from other planctomycetal species and therefore propose the new genus Fuerstia gen. nov., with the type species Fuerstia marisgermanicae sp. nov.
Materials and Methods
The strain was isolated from a crustacean shell (crab) which was collected on the 22nd May 2012 at low tide (8:30 a.m.) at 16°C ambient temperature in some residue water on the tidal mud flat of the German Wadden Sea, Neuharlingersiel (53°42′15.5″ N, 7°42′15.9″ E).
Isolation and Maintenance
In a first step the surface of a crustacean shell was scraped-off and the shell material was transferred to sterile artificial sea water (ASW) supplemented with 10-fold diluted HD medium, to serve as biofilm suspension. ASW was prepared modified after Levring (1946) consisting of (per liter distilled water): 23.6 g NaCl; 0.64 g KCl; 4.53 g MgCl2 · 6 H2O; 5.94 g MgSO4 · 7 H2O; 1.3 g CaCl2 · 2 H2O; 10 mg Na2PO4 · 2 H2O; 2.1 mg NH4NO3. To avoid precipitation, the CaCl2 solution was sterilized separately (Bruns et al., 2003). HD medium was composed of 0.25 g/l yeast extract, 0.1 g/l glucose, 0.5 g/l peptone and 2.38 g/l HEPES; the pH was adjusted to 7.3. The suspension of the bacterial biofilm was used to inoculate fifty 96 well plates employing a multidrop device as previously described (Jogler et al., 2013). Cultures were subsequently transferred to fresh 96 well plates with the same medium but supplemented with 2 mg/ml carbenicillin. Cultures that survived this treatment were screened employing a 16S rRNA gene targeting PCR, using the primer set 8f (5′–AGA GTT TGA TCM TGG CTC AG–3′) and 1492r (5′–GGY TAC CTT GTT ACG ACT T–3′) modified from Lane (1991). PCR amplifications were performed employing a Veriti 96-Well Thermal Cycler (Applied Biosystems) applying the following conditions: initial denaturation at 94°C for 5 min, followed by 10 cycles of denaturation at 94°C for 30 s, annealing at 59°C for 30 s and elongation at 72°C for 60 s. This first 10 cycles were followed by 20 cycles of denaturation at 94°C for 30 s, annealing at 54°C for 30 s, elongation at 72°C for 60 s and a final elongation at 72°C for 5 min. Amplification products were subject to 16S rRNA gene sequencing and NCBI database comparison to identify novel strains. Based on this analysis, strain NH11T was selected for subsequent experiments. Further cultivation was performed employing a modified M2 culture broth (M2mod) previously used for Rhodopirellula baltica (Jeske et al., 2013) containing 0.5 g/l peptone, 0.25 g/l glucose, 250 ml/l artificial sea water (ASW), 5 ml/l vitamin solution (double concentrated) and 20 ml/l mineral salt solution. The medium was buffered with 2.38 g/l HEPES at pH 7.0. Artificial sea water was composed of 46.94 g/l NaCl, 7.84 g/l Na2SO4, 21.28 g/l MgCl2 · 6 H2O, 2.86 g/l CaCl2 · 2 H2O, 0.384 g/l NaHCO3, 1.384 g/l KCl, 0.192 g/l KBr, 0.052 g/l H3BO3, 0.08 g/l SrCl2 · 6 H2O, and 0.006 g/l NaF. The vitamin solution was composed of 4 mg/l biotin, 4 mg/l folic acid, 20 mg/l pyridoxine-HCl, 10 mg/l riboflavin, 10 mg/l thiamine-HCl · 2 H2O, 10 mg/l nicotinamide, 10 mg/l D-Ca-pentothenate, 0.2 mg/l vitamin B12, and 10 mg/l p-aminobenzoic acid. Mineral salt solution was composed of 10 g/l nitrilotriacetic acid (NTA), 29.70 g/l MgSO4 · 7 H2O, 3.34 g/l CaCl2 · 2 H2O, 12.67 mg/l Na2MoO4 · 2 H2O, 99 mg/l FeSO4 · 7 H2O, and 50 ml/l metal salt sol. “44”. Metal salt sol. “44” was composed of 250 mg/l Na-EDTA, 1.095 g/l ZnSO4 · 7 H2O, 0.5 g/l FeSO4 · 7 H2O, 154 mg/l MnSO4 · H2O, 39.20 mg/l CuSO4 · 5 H2O, 24.80 mg/l Co(NO3)2 · 6 H2O, and 17.70 mg/l Na2B4O7 · 10 H2O. Solid medium was prepared by adding 15 g/l of three times prewashed agar (Becton, Dickinson and Company) to the medium.
Phase contrast (Phaco) and differential interference contrast (DIC) analysis were performed employing a Nikon Eclipse Ti inverted microscope with a Nikon N Plan Apochromat λ 100x/1.45 Oil objective and a Nikon DS-Ri2 camera. Specimens were immobilized in MatTek glass bottom dishes (35 mm, No. 1.5) employing a 1% agarose cushion. Images were analyzed using the Nikon NIS-Elements software (Version V4.3).
For field emission scanning electron microscopy (FESEM) bacteria were fixed in 1% formaldehyde in HEPES buffer (3 mM HEPES, 0.3 mM CaCl2, 0.3 mM MgCl2, 2.7 mM sucrose, pH 6.9) for 1 h on ice and washed one time employing the same buffer. Cover slips with a diameter of 12 mm were coated with a poly–L–lysine solution (Sigma–Aldrich) for 10 min, washed in distilled water and air–dried. Fifty microliter of the fixed bacteria solution was placed on a cover slip and allowed to settle for 10 min. Cover slips were then fixed in 1% glutaraldehyde in TE buffer (20 mM TRIS, 1 mM EDTA, pH 6.9) for 5 min at room temperature and subsequently washed twice with TE–buffer before dehydrating in a graded series of acetone (10, 30, 50, 70, 90, and 100%) on ice for 10 min at each concentration. Samples from the 100% acetone step were brought to room temperature before placing them in fresh 100% acetone. Samples were then subjected to critical–point drying with liquid CO2 (CPD 300, Leica). Dried samples were covered with a gold/palladium (80/20) film by sputter coating (SCD 500, Bal–Tec) before examination in a field emission scanning electron microscope (Zeiss Merlin) using the Everhart Thornley HESE2–detector and the inlens SE–detector in a 25:75 ratio at an acceleration voltage of 5 kV. TEM micrographs of NH11T cells were taken after negative staining with aqueous 0.1–2% uranyl acetate, employing a Zeiss transmission electron microscope EM 910 at an acceleration voltage of 80 kV at calibrated magnification as previously described (Wittmann et al., 2014). Prior TEM analysis cells were fixated for 2 h with 3% glutaraldehyde in 3 mM EM-HEPES buffer.
Thin sections of strain NH11T were prepared by high pressure freezing and freeze substitution as previously described (Jogler et al., 2011). Sections were subsequently analyzed employing a JEOL 1200EX–80kV TEM microscope.
Physiological tests such as salinity, pH, and temperature tolerance were performed in liquid medium M2mod. To test the NaCl tolerance, medium M2mod was prepared employing ASW devoid of NaCl. The required NaCl concentrations were adjusted prior inoculation, using a 30% NaCl (w/v) solution. The ASW tolerance was tested using M2mod without ASW. Again, the required concentration of ASW was adjusted prior inoculation. Growth was detected by monitoring the optical density at 600 nm using a Photometer Ultrospec II (LKB Biochrom). Carbon source utilization was tested using the GN2 MicroPlate™ system (Biolog), while enzymatic activities were tested using the API® ZYM method (bioMérieux). The physical features of the cell wall were analyzed by Gram staining, KOH test and Bactident® Aminopeptidase test strips. Growth under anoxic conditions was investigated using the api® 20 NE system (bioMérieux) according to the manufacturer's instructions.
Phylogenetic Analysis and Tree Reconstruction
For phylogenetic analysis version 6.0.2 of the ARB software package (Ludwig et al., 2004) was used together with the SILVA database SSURef_NR99 (Version 119, released on 14.07.2014; Pruesse et al., 2007). The full length 16S rRNA gene sequence was imported into ARB and then aligned using the fast aligner tool of the ARB software package. The resulting alignment was further edited manually to improve alignment quality. During the phylogenetic tree reconstruction, different type strains of the phylum Planctomycetes were used as reference sequences, while type strains of the phylum Verrucomicrobia served as out-group. All sequences that were included are listed in Table S1. Tree reconstruction was performed with the ARB software package (Ludwig et al., 2004). The Maximum Likelihood RAxML module was used with the rate distribution model GTR GAMMA running the rapid bootstrap analysis algorithm. The Neighbor Joining tool was employed with Felsenstein correction, while the Maximum Parsimony analysis was achieved with the Phylip DNAPARS module. Bootstrap values for all three methods were calculated with 1000 resamplings including the E. coli 16S rRNA gene positions 101–1371.
To facilitate the taxonomic classification of strain NH11T, a cluster analysis was performed employing version 6.0.6 of the ARB software package (Ludwig et al., 2004) together with the database SSURef_NR99 (Version 123.1, released on 03.03.16) using a 87.65% sequence identity cutoff (threshold for the minimal sequence identity within a taxonomic family after Yarza et al., 2014) and E. coli 16S rRNA gene positions 112–1393. This analysis included all sequences of the family Planctomycetaceae that were available in the database as well as sequences of uncultivated planctomycetes, while the three clusters obtained for Rubinisphaera brasiliensis DSM 5305T, Gimesia maris 534-30T and strain NH11T were selected for further phylogenetic tree reconstruction. For the NH11T containing cluster, the highest, lowest, average and median sequence identities were determined by calculating a distance matrix based on the E. coli 16S rRNA position 112–1393.
Genomic DNA of strain NH11T was extracted using the Genomic-tip 20/G kit (Qiagen, Germany) and a >10 kb SMRTbell™ template library was prepared according to the manufacturer instructions (PacificBiosciences, USA). In brief, ~10 μg of genomic DNA was end-repaired and ligated to hairpin adapters overnight (DNA/Polymerase Binding Kit 2.0, Pacific Biosciences, USA). The SMRTbell™ template was exonuclease treated for removal of incompletely formed reaction products. Conditions for annealing of sequencing primers and binding of polymerase to purified SMRTbell™ template were assessed with the calculator in RS Remote (PacificBiosciences, USA). Five SMRT cells were sequenced on the PacBio RS/RSII (PacificBiosciences, USA), each covered with a 90-min movie. Eight additional SMRT cells were used applying the DNA/Polymerase Binding Kit P4 (PacificBiosciences, USA) in order to collect larger read lengths. For those cells, 180-min movies were taken and a sub read lengths up to 20 kb was observed.
In addition, 2 × 9,143,110 paired-end Illumina reads were obtained employing an Illumina HiSeq 2500 for 101 cycles in both directions using the TruSeq DNA Sample Prep Kit v2 according to the manufacturer's instructions (Illumina Inc., San Diego, CA, USA). The genome assembly was performed using the RS_HGAP_Assembly.2 protocol included in SMRT Portal version 2.2.0 utilizing 520,016 reads from all 13 SMRT cells applying standard parameters with exception of the genome size, which was set to 9,100,000. Thus, one final contig could be obtained, which afterwards was error-corrected using a subset of 2 × 4,000,000 Illumina reads using BWA (Li and Durbin, 2009) with subsequent variant and consensus calling employing the CLC Genomics Workbench 7.03 (http://www.clcbio.com). Visual inspection of 41 called variants has been performed using IGV (Thorvaldsdóttir et al., 2013) in addition to prediction of rRNAs using RNAmmer (Lagesen et al., 2007). Hereby, the error-corrected consensus was trimmed, circularized and adjusted to dnaA as the first gene. For the circular genome of Fuerstia marisgermanicae NH11T the finishing quality is estimated to be at 99.9999% (QV 60) confirmed by final PacBio resequencing. Genome annotation has been performed using PROKKA 1.8 (Seemann, 2014).
Genomic islands and CRISPR regions were identified using IslandViewer3 (Dhillon et al., 2015) and CRT (Bland et al., 2007), respectively. The planctomycetal genomes for the gene content analysis were derived from NCBI and IMG (Markowitz et al., 2012) in April 2016 and had to match the following criteria upon CheckM analysis (Parks et al., 2015): completeness >90, contamination <5 and strain heterogeneity <20. Orthologs were detected with Proteinortho5 (Lechner et al., 2011), a tool that identifies the reciprocal best hits from the given protein sequences. The genome plot was then generated with BRIG (Alikhan et al., 2011). The G+C content of the DNA was estimated by analyzing the genome data using the software Artemis (Rutherford et al., 2000).
The gene size was calculated from annotations and analyzed with R (R Core Team, 2015) using the packages ggplot2 (Wickham, 2009), reshape (Wickham, 2007), xlsx (Dragulescu, 2014). Threshold for large genes were set to 5 kb according to Reva and Tümmler (2008). The biggest giant genes of NH11T with a size over 20 kb were further analyzed using the InterProScan web service (Jones et al., 2014; Mitchell et al., 2015). The genome of strain NH11T was screened for putative secondary metabolite clusters using the tool antiSMASH 3.0 (Weber et al., 2015).
Multilocus Sequence Analysis (MLSA)
Orthologs were identified using Proteinortho5 (Lechner et al., 2011) with the “-selfblast” option (which enables paralog-detection) enabled. Only genes present exclusively in single copy in all compared genomes were selected for MLSA. Alignments of the respective gene product amino acid sequences were generated individually for each ortholog group using MUSCLE v3.8.31 (Edgar, 2004a,b) and subsequently concatenated. Unalignable regions, caused by e.g., unique N-terminal or C-terminal sequence overhangs, were filtered from the concatenated alignment using Gblocks v.0.91b (Castresana, 2000). Phylogenetic relationships were inferred from the remaining unambiguous alignment positions by Neighbor Joining clustering with 1000 bootstrap iterations, using ARB 6.0.5 (Ludwig et al., 2004) and by Maximum Likelihood calculation using RAxML v. 8.0.26 (Stamatakis, 2014).
Fatty Acid Analysis
For fatty acid analysis, biomass of the three compared strains (NH11T, R. brasiliensis DSM 5305T and G. maris 534-30T) was obtained from liquid cultures, grown in M2mod culture broth, at 20°C and slight agitation in baffled flasks. The obtained biomasses were processed according to the standards of the Identification Service of the German Collection of Microorganisms and Cell Cultures (DSMZ) (Miller, 1982; Kuykendall et al., 1988; Kämpfer and Kroppenstedt, 1996).
Nucleotide Sequence Accession Numbers
The accession numbers of the nucleotide sequences used for giant gene analysis or phylogenetic reconstruction are shown in Table S1. The NCBI genome accession number is CP017641.
Based on 16S rRNA gene phylogeny (Figure 1), strain NH11T clusters together with G. maris within the family Planctomycetaceae (bootstrap support values: 69% for Maximum Likelihood, 43% for Maximum Parsimony, 46% for Neighbor Joining). However, it shares only 87.1% 16S rRNA gene sequence identity with G. maris based on BLAST analysis. Its next closest relative besides G. maris is R. brasiliensis (84.5% 16S rRNA gene sequence identity). An overview of the identities, generated with the manually curated SILVA alignment of all species of the phylum Planctomycetes, based on the 1270 base pairs used for phylogenetic reconstruction, is given in Table S2. In this analysis the Felsenstein correction was used, to take evolutionary events that occurred during speciation into account. According to this analysis strain NH11T shares 85.4% identity with G. maris and 84.1% with R. brasiliensis, differing from BLAST search results described above. The performed cluster analysis revealed, that NH11T, G. maris and R. brasiliensis form three distinct clusters within the family Planctomycetaceae with an identity of at least 87.65% within each cluster, which is identical with the minimum sequence identity threshold for taxonomic families (Yarza et al., 2014; Figures S1–S3). In particular, within the cluster of strain NH11T the minimal sequence identity is 88.3% and the maximal identity is 99.9%. The average identity in this cluster is 94.48% with a median of 94.4%, matching the threshold of 92.25% for median sequence identity for taxonomic families.
Figure 1. Phylogeny of strain NH11T. A maximum-likelihood 16S rRNA gene based tree of the phylogenetic position of strain NH11T. Strain NH11T is highlighted in bold letters and in comparison planctomycetal type strains are shown, while selected Verrucomicrobia served as outgroup. The bootstrap percentages of 1000 resamplings of three different tree building methods (ML, Maximum Likelihood; NJ, Neighbor Joining; MP, Maximum Parsimony) are incorporated in this tree, indicated by different shaded dots. Black dots indicate support values above 70% for all three methods while white dots refer to a bootstrap support below 50% for at least one method. Branches that were not supported by all three methods show no dot.
Cells of strain NH11T are pear to ovoid shaped (Figure 2) and form cream colored colonies on solid medium and motile swimmer cells in liquid culture. Cells are 1.2–2.5 × 0.9–1.7 μm in size. Its major phenotypic characteristics compared to those of G. maris and R. brasiliensis are listed in Table 1. No rosettes or stalks are formed in culture in contrast to the closest cultivated relatives G. maris and R. brasiliensis (Figures 2, 3). The surface of strain NH11T is smooth and seems to lack crateriform structures except of polar regions were fiber-like structures seem to emerge from crateriform pits (Figure 3). Furthermore, the cell architecture of strain NH11T and its mode of division seem to differ from all other planctomycetes described thus far (Figure 2). High pressure frozen and freeze substitution sections of strain NH11T display a condensed nucleoid that was previously described in other planctomycetes, too (Figures 4A–C, white arrow). However, some cells display exceptional patterns of cytoplasmic invaginations (Figures 4A,B), while others, except for the condensed nucleoid, comprise a cell envelope comparable to E. coli (Figure 4C). Despite of the degree of invaginations with characteristics different from other planctomycetes, ribosomes seem to localize frequently in close proximity to the cytoplasmic membrane in a rope-of-pearls-like pattern (Figures 4A–C, white arrowheads). Cells of strain NH11T divide through budding and unlike all other planctomycetes described thus far, the mother- and daughter-cells seem to be connected via a tubular structure during cell division (Figures 4D–F). In particular, a dark structure at the division plane seems to parallel the cell division ring previously described for anammox Planctomycetes (van Niftrik et al., 2009; Figure 4D, black arrow).
Figure 2. Morphology of strain NH11T. Representative differential interference contrast (DIC), Phase-contrast (Phaco) and scanning electron microscopic (SEM) micrographs of strain NH1. Shown are individual representative cells during division through polar budding (A–C) and -aggregate formation (D–F). In addition, an overview of multiple cells from a representative liquid culture is provided (G–I). Bar, 1 μm.
Figure 3. Crateriform structures of strain NH11T. TEM analysis of negative stained NH11T cells revealed crateriform structures solely on the cell pole (A, white arrow) and fiber like structures seem to arise from such pits (B, black arrow). Bar, 0.5 μm.
Figure 4. Subcellular organization and cell division of strain NH11T. High pressure frozen and freeze substituted cells of strain NH11T were subject to thin-sectioning and TEM analysis. Three representative sections of individual cells correspond to differences in subcellular organization of strain NH11T (A–C). Cells always comprise a condensed nucleoid which is visible as a black structure (white arrow), while their cytoplasmic membrane can show different degrees of invaginations (A,B, black arrowheads). Some cells even lack any invagination of the cytoplasmic membrane (C). Some ribosomes (white arrowheads) are always located at the cytoplasmic membrane despite its degree of invagination (A–C). Dividing cells are interconnected by a tubular like structure (D–F). Between mother and daughter cell, sometimes a black structure, potentially correlating with a division ring, became visible (D, black arrow). This structure could only be observed in some sections, while it was absent in others (E–F). Bar, 0.1 μm.
Strain NH11T is able to grow between pH 6 and 10, with an optimal growth at pH 7 (Figure S4A). NH11T tolerates up to 5% NaCl supplementation (Figure S4B) and requires at least 27.5% ASW for detectable growth (Figure S4C). In contrast, up to 230% ASW can be tolerated, while optimal growth conditions are in the range of 50–117.5% ASW. Temperature-wise, strain NH11T can grow between 20 and 30°C, with optimal growth at 28°C (Figure S4D). Thus, strain NH11T is a mesophilic organism. It is capable of utilizing a variety of carbon sources, listed in Table 2. In particular N-acetyl-D-galactosamine, N-acetyl-D-glucosamine, L-arabinose, D-cellobiose, D-galactose, gentiobiose, α-D-glucose, α-D-lactose, lactulose, maltose, D-mannose, D-melibiose, β-methyl-D-glucoside, sucrose, D-trehalose, turanose, succinic acid mono-methyl-ester, acetic acid, γ-hydroxybutyric acid, itaconic acid, propionic acid and glycerol were utilized. In contrast to the closest related G. maris and R. brasiliensis, strain NH11T was not able to use L-rhamnose as carbon source. The enzymatic repertoire of strain NH11T comprises alkaline phosphatase, esterase, esterase lipase, leucinearylamidase, valinearylamidase, cysteine arylamidase, trypsin, α-chymotrypsin, acid phosphatase, naphthol-AS-BI-phosphohydrolase, α-galactosidase, and α-glucosidase. All enzymatic features are listed in Table 3. No growth under anoxic conditions was observed.
The fatty acids of strain NH11T and its closest relatives G. maris and R. brasiliensis were analyzed and compared. The major fatty acids of strain NH11T consist of 59.16% C16:1 ω6c/16:1 ω7c (summed feature 3), 19.83% C18:1 ω6c/18:1 ω7c (summed feature 8) and 15.12% C16:0. In comparison, G. maris comprised 26.63% of C16:1ω6c/16:1 ω7c (summed feature 3), 23.50% of C16:0, and 12.98% of C16:0 10-methyl/Iso-C17:1 ω6c (summed feature 9). The fatty acid profile of R. brasiliensis is composed of 47.10% C16:0 and 45.77% C16:1ω6c/16:1 ω7c (summed feature 3). The complete fatty acid profiles of strain NH11T, G. maris and R. brasiliensis are listed in Table 4.
Taken together, the general fatty acid profile of strain NH11T is similar, yet distinct compared to its closest relatives G. maris and R. brasiliensis and differs mainly in the components proportions.
The complete chromosome of strain NH11T comprises 8,920,478 bp, with a GC content of 55.9%. Within the genome, 6732 genes were annotated, of which 6645 were identified as protein-coding genes. As known for other planctomycetes (Fuerst and Sagulenko, 2011), only for 54.3% of these genes a function could be predicted, while the remaining genes were annotated as hypothetical proteins or proteins with unknown function. Notable islands of unique gene content were found to especially comprise regions of predicted genomic islands and regions containing so-called “giant genes” (Reva and Tümmler, 2008). The genomic features of strain NH11T are summarized in Figure 5. In total, the genome of strain NH11T comprises 45 giant genes, with a size greater than 5 kb (Figure S5). By the time of the first description of giant genes in the genome of R. baltica (Reva and Tümmler, 2008) no other planctomycetal genome sequence was available. Thus, we analyzed the general distribution of giant genes amongst Planctomycetes and found all 29 available planctomycetal genomes to encode such genes (Figure 6). Fifteen strains encoded genes with a size >20 kb, two of which Rhodopirellula sp. K833 and NH11T, even encode five >20 kb genes (Figure S6). Thus, genes >20 kb could be found in less than 50% of the available planctomycetal genomes. Furthermore, genes >30 kb were exclusively found in strain NH11T and G. maris. Thus strain NH11T is somewhat unique if gene-length is compared to other planctomycetes. Consequently, we analyzed the potential function of five NH11T protein-coding genes >20 kb. The results, obtained from the InterProScan web service, points toward potential functions in cell adhesion, carbohydrate binding and partial location on the outer membrane of NH11T.
Figure 5. Genomic features of strain NH11T. The circular plot of strain NH11T's 8.9 Mb chromosome summarizes multiple genomic features: the outer circles display protein (light and dark blue), tRNA (turquoise), and rRNA (pink) encoding genes as well as predicted genomic islands (gray) and giant genes (purple). The innermost circle shows the GC plot (gray). In between, ortholog genes from planctomycetal strains, for which high quality genomes were available, were identified by reciprocal BLAST and are depicted in green, yellow, and orange, according to the color scheme of Figure 7. Notable islands of unique gene content in NH11T are visualized as gaps and they are found to in regions of predicted genomic islands and/or giant genes.
Figure 6. Planctomycetal giant genes. The gene size distribution of different Planctomycetes with available genomes is ordered by declining phylogenetic relationship to NH11T. Median size of genes, 25 and 75% percentiles are shown and whiskers represents the 1.5* interquartile range and dots outliers. Fuerstia marisgermanicae comprises 45 giant genes (>5 kb), the biggest gene is 36,239 bp in size.
Because TEM analysis pointed toward an anammox Planctomycetes-like cell division ring (van Niftrik et al., 2009; Figure 4D, black arrow), the genome was analyzed toward the presence of the putative cell division ring protein of Candidatus “K. stuttgartiensis” (locus tag kustd1438, NCBI Accession Number CAJ72183). The best match was a hypothetical protein with 36% sequence identity but only 19% query coverage (Fuma_00096), suggesting that this is presumably not a homolog of the putative protein found in Candidatus “K. stuttgartiensis.” However, with a size of 8699 bp, the gene that encodes for this hypothetical protein of strain NH11T belongs to the identified giant genes. Analysis of this hypothetical protein, using the InterProScan web service, revealed the presence of Lamin Tail- (IPR001322), Ig-like- (IPR032812), and PapD-like domains (IPR008962), indicating putative involvement of this protein in cell shape maintenance or sub cellular organization.
Besides the giant genes, further analysis of the genome revealed the presence of a complete CRISPR type II system (cas9: Fuma_04885, cas1: Fuma_4887, cas2: Fuma_4888; Taylor et al., 2015). Only 2 other of the 35 analyzed planctomycetal genomes, corresponding to organisms Blastopirellula marina SH 106 and Candidatus “Brocadia sinica JPN1,” comprise such a system while the closest relative of strain NH11T, G. maris lacks equivalent genes. This indicates that bacteriophages can infect strain NH11T as the CRISPR type II system is primarily used for phage defense. Furthermore, absence in G. maris indicates differences in terms of phage susceptibility and/ or defense strategy in comparison to strain NH11T.
As Planctomycetes are a known source of potential bioactive compounds (Jeske et al., 2013), the genome of strain NH11T was analyzed using antiSMASH 3.0 (Weber et al., 2015). In total, four secondary metabolite-associated gene clusters were found. In contrast, G. maris, as closest relative to NH11T, encodes 9 putative secondary metabolite related genes and gene clusters. Thus, again, strain NH11T differes from its peers.
However, the four identified genes and clusters might be related to bacteriocins (2), terpenes (1), and ectoines (1) based on the antiSMASH prediction. Ectoins are compatible solutes, playing a role in the salt stress response, e.g., in halophilic eubacteria (Peters et al., 1990; Bernard et al., 1993; Pastor et al., 2010). However, physiological tests gave no evidence for an increased salt tolerance [>5% NaCl (w/v)] of strain NH11T (see Table 1).
The genome sequence of strain NH11T was further used to verify the 16S rRNA gene sequence based phylogeny via Multi Locus Sequence Analysis (MLSA). This method enables to increase the phylogenetic resolution and to avoid the bias caused by single marker gene-based phylogenies. This was achieved by considering sequences of multiple independent single copy genes shared by the comparison organisms (Glaeser and Kämpfer, 2015). The MLSA-based phylogeny is shown in Figure 7 and further supports our 16S rRNA gene based placement.
Figure 7. Multilocus sequence analysis (MLSA) of strain NH11T. The phylogenetic tree is based on Neighbor Joining analysis, with bootstrap values (presented by black, gray and white dots) based on 1000 permutations. Additionally, Maximum Likelihood (ML) analysis was performed and is also indicated by dots. The performed calculations where based on 68,695 unambiguous amino acid sequence alignment positions of 143 concatenated orthologous single-copy protein gene products shared by all 37 included genomes. The ortholog selection was based on bidirectional BLAST analysis as implemented in Proteinortho5 (Lechner et al., 2011). Sequences were aligned individually for each ortholog group and subsequently concatenated. Unalignable regions were filtered from the alignments using Gblocks (Castresana, 2000). The genomes of Opitutus terrae PB90-1 and Verrucomicrobium spinosum DSM 4136 served as outgroup. To enhance clarity, several monophyletic groups of reference organisms are shown collapsed as others are indicated by orange and green boxes, according to the color scheme of Figure 5.
Strain NH11T is in many ways an unusual planctomycete. Based on 16S rRNA gene sequence comparison, it branches within the former genus Planctomyces which was recently found to be polyphyletic. Consequently, it was further subdivided into the separate genera Gimesia, Planctopirus, Rubinisphaera, and Planctomyces, which all belong to the order Planctomycetales and the family Planctomycetaceae (Scheuner et al., 2014). Similar results were obtained by a MLSA phylogenetic analysis, while both approaches provided rather low bootstrap support. However, since both methods gave similar results, strain NH11T is likely to cluster within the family Planctomycetaceae.
To allow a detailed taxonomic classification, we performed further 16S rRNA gene sequence comparison (Table 1). We found, if evolutionary events were considered, that strain NH11T belongs to a novel family based on current thresholds (Yarza et al., 2014). In contrast, direct comparison, without acknowledging evolutionary events, led to 87.1% 16S rRNA gene sequence identity to its closest relative G. maris, while the current taxonomic thresholds of the bacterial family border is 86.5%. However, within a family, a minimum sequence identity of 87.65% should be given between the separate species (Yarza et al., 2014), while strain NH11T and G. maris comprise only 87.1% sequence identity. Thus without taking evolutionary events into account, strain NH11T falls right in between both thresholds (Yarza et al., 2014). To solve this conflict, a 16S rRNA sequence alignment based cluster analysis was performed. The results (Figures S1–S3) clearly demonstrate that all available 16S rRNA gene sequences from cultivated and uncultivated members of the family Planctomycetaceae form multiple clusters, while strain NH11T and the genera Gimesia and Rubinisphaera belong to three distinct clusters. Thus from a phylogenetic perspective strain NH11T belongs to a novel family.
This conclusion is consistent with other types of evidence obtained in this study: the morphological features of strain NH11T differ from those of its closest relatives. For example, G. maris is described to stain Gram negative (Bauld and Staley, 1976) whereas employing the same Gram staining method to strain NH11T delivered no clear result. However, others and we recently demonstrated the presence of peptidoglycan in a phylogenetically representative set of planctomycetal model species (Jeske et al., 2015; van Teeseling et al., 2015). Consequently, we suggest to exclude this classical Gram staining test as a valid tool for future characterization of novel planctomycetal strains, as it is of limited explanatory power for Planctomycetes. Nevertheless, other aspects of strain NH11T differ from its peers as well. Another aspect that differs between strain NH11T and members of the family Planctomycetaceae such as P. limnophila is the localization of the crateriform structures. They are hardly visible (Figure 3) and seem to be associated with fiber-like structures. However, such fiber-like structures might be artifacts from EM fixation methods and thus require further investigation. If such structures are associated with cell attachment to surfaces remains as well enigmatic. While G. maris and R. brasiliensis form stalks and rosettes in culture, strain NH11T lacks both features. In addition, strain NH11T forms a tubular connection between daughter and mother cell (Figures 4D–F). Within this tubes in some sections a dark structure became visible that was previously identified as cell division ring in anammox Planctomycetes (van Niftrik et al., 2009). However, anammox bacteria divide, in contrast to budding members of the order Planctomycteales, through binary fission. While both types of organisms lack the otherwise universal cell division protein FtsZ, only for anammox Planctomycetes a potential substitute is known (locus tag kustd1438, NCBI Accession Number CAJ72183). While no homolog protein could be determined in the genome of strain NH11T, the molecular mechanism of its unusual cell division remains enigmatic. One might speculate that Fuma_00096, a protein similar, yet distinct from kustd1438, might be involved in the formation of the unique tubular connection between mother and daughter cell and might correlate with the anecdotally observed structure (Figure 4D, black arrow). This idea is based on the Lamin Tail domain of Fuma_00096, which corresponds to an intermediate filament (IF). Such IFs are frequently found in eukaryotic cells to provide mechanical strength and support for fragile tubulin structures (Goldman et al., 1996). However, this is only a hypothesis at this stage, that requires experimental verification.
Besides morphological differences, the genome of strain NH11T shows distinct features if compared to other Planctomycetes as well. First, it contains areas of unique genes, if compared against the available high quality planctomycetal genomes. Second, these areas were found to comprise both, genomic islands - likely acquired through horizontal gene transfer - and giant genes. In total, 19 giant genes with a size >10 kb were detected, which is the highest number amongst all analyzed planctomycetal genomes (Figure S7). Representatives of the family Planctomycetaceae in contrast possess a maximum of 13 giant genes of this dimension. The protein products of giant genes, if synthetized, would represent a huge metabolic burden, suggesting that they have an important cellular function. Thus, 6–19 additional giant genes >10 kb compared to its peers suggests a huge difference in strain NH11T's cell surface or secondary metabolism. This is because more than 90% of giant genes encode either surface proteins or polyketide/nonribosomal peptide synthetases (PKS/NRPS) (Reva and Tümmler, 2008). Surprisingly, our antiSMASH analysis revealed only four secondary metabolite-associated gene clusters. One of these gene clusters is predicted to encode the compatible solute ectoine, giving a hint toward an increased salt tolerance of the strain. However, since ASW is tolerated up to 230% (v/v) and NaCl up to 5% (w/v), it is rather adapted to moderate salt concentration. In contrast to the four clusters found in strain NH11T's genome, its closest relatives R. brasiliensis and G. maris comprise 8 and 9 of such clusters (Jeske et al., 2013). Furthermore, a significant linear correlation between genome size and number of secondary metabolite related genes has been described (Jeske et al., 2013). With a genome size of almost 9 MB strain NH11T is an exception to this observation. Two potential scenarios could be envisioned: either the giant genes of strain NH11T fulfill other, maybe structural functions, or they might be involved in the formation of yet unknown secondary metabolites. While the latter hypothesis fits to the recent observation of novel antibiotic small molecules from Planctomycetes (Jeske et al., 2016), at least some giant genes are predicted to encode proteins with structural function such as the discussed Fuma_00096. While resolving this issue is beyond the scope of this study, we conclude that the genome of strain NH11T differs in important aspects from its peers.
Taken together, the 16S rRNA gene sequence identity, the unusual mode of cell division, its unique cell plan and its unusual planctomycetal genome require to place this novel and exceptional species in a new family. However, the family Planctomycetaceae would not be monophyletic anymore and based on our cluster analysis requires division into three distinct families. To prevent back- and forth renaming of species, genera and families, we describe strain NH11T as novel genus and species for now. We soon will present multiple novel planctomycetal strains and only in the light of such isolates rewriting the planctomycetal taxonomy would make sense. Thus, rearrangements within the family Planctomycetaceae will be revisited in the future (Jogler C, personal communication).
Description of Fuerstia gen. nov.
Fuerstia (named in honor of John Fuerst, an Australian microbiologist from University of Queensland, who played a key role in planctomycetal research). The pear to ovoid shaped cells form aggregates in liquid culture, but no rosettes. Daughter cells are motile, while mother cells are non-motile and no stalk formation was observed. The surface is smooth, crateriform stuctures are limited to one pole and cells reproduce by polar budding while mother- and daughter cells are connected by a thin tubular-like structure. The lifestyle is heterotrophic, obligatory aerobic and mesophilic. The major fatty acids are C16:1 ω6c/16:1 ω7c (Summed feature), C18:1ω6c/18:1 ω7c (Summed feature), and C16:0. Member of the phylum Planctomycetes, class Planctomycae, order Planctomycetales, family Planctomycetaceae. The type species is Fuerstia marisgermanicae.
Description of Fuerstia marisgermanicae sp. nov.
Fuerstia marisgermanicae (ma′ris L. n. maris of the sea and ger.ma'ni.cae L. adj. German, pertaining to the German North Sea from which the type strain was isolated). In addition to the features described for the genus, the species exhibits the following properties. Colonies on solid medium are cream colored. Cells are 1.2–2.5 × 0.9–1.7 μm in size. Non motile mother cells spawn motile, swimming, daughter cells. Gram staining delivers no clear result. KOH test and aminopeptidase test are negative, while oxidase and catalase tests positive. The organism is able to degrade a wide range of carbon sources, in particular N-acetyl-D-galactosamine, N-acetyl-D-glucosamine, L-arabinose, D-cellobiose, D-galactose, gentiobiose, α-D-glucose, α-D-lactose, lactulose, maltose, D-mannose, D-melibiose, β-methyl-D-glucoside, sucrose, D-trehalose, turanose, succinic acid mono-methyl-ester, acetic acid, γ-hydroxybutyric acid, itaconic acid, propionic acid, and glycerol were utilized. The enzymatic repertoire of the species, tested with API ZYM, is listed in Table 3. Growth occurs between pH 6 and 10 with an optimum at pH 7. At least 27.5% ASW is needed for growth, The optimal temperature for growth is 28°C (range between 20 and 30°C). The type strain is NH11T (= DSM 27554 = LMG 27831) isolated from a crustacean shell.
TK analyzed the data and wrote the manuscript. AH performed the physiological experiments and proofread the manuscript. MJ was involved in planning the project, aided in cultivation and helped writing and proofread the manuscript. JV, AK developed the pipeline for the genome analysis. CB prepared thin sections of the strain and performed TEM imaging, was involved in writing and proofreading the manuscript. JO, BB did the genome sequencing of the strain and were involved writing and proofreading the manuscript. PR accomplished phylogenetic analyses, helped writing and proofreading the manuscript. DB performed light-microscopic observations. IG performed the initial cultivation and isolated the stain. HF analyzed the genome for giant genes and proofread the manuscript. HK conducted the fatty acid analysis and proofread the manuscript. MR did the scanning electron microscopy of the strain, TEM imaging of negatively stained cells and proofread the manuscript. SW was involved in the genome analysis and did the multilocus sequence analysis and proofread the manuscript. CJ designed the project and wrote the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fmicb.2016.02079/full#supplementary-material
Figure S1. Maximum likelihood tree of the detected clusters of NH11T, Gimesia maris, and Rubinispharea brasiliensis using a minimum sequence identity threshold of 87.65% within a cluster corresponding to phylogenetic families. All three cultivated species form distinct clusters with sequences of uncultivated planctomycetes with at least 87.65% similarity within their cluster. The corresponding clusters are marked with brackets. The bootstrap percentages of 1000 resamplings are given for each node. Anammox planctomycetes served as outgroup.
Figure S2. Maximum parsimony tree of the detected clusters of NH11T, Gimesia maris, and Rubinispharea brasiliensis using a minimum sequence identity threshold of 87.65% within a cluster corresponding to phylogenetic families. All three cultivated species form distinct clusters with sequences of uncultivated planctomycetes with at least 87.65% similarity within their cluster. The corresponding clusters are marked with brackets. The bootstrap percentages of 1000 resamplings are given for each node. Anammox planctomycetes served as outgroup.
Figure S3. Neighbor joining tree of the detected clusters of NH11T, Gimesia maris, and Rubinispharea brasiliensis using a minimum sequence identity threshold of 87.65% within a cluster corresponding to phylogenetic families. All three cultivated species form distinct clusters with sequences of uncultivated planctomycetes with at least 87.65% similarity within their cluster. The corresponding clusters are marked with brackets. The bootstrap percentages of 1000 resamplings are given for each node. Anammox planctomycetes served as outgroup.
Figure S4. Temperature, pH, ASW, and NaCl optimum of strain NH11T. To determine the pH (A), NaCl (B), and ASW (C) tolerance as well as the optimal growth temperature (D), the optical density was measured at 600 nm (OD 600 nm) and slope values, corresponding to change of OD 600 nm over time during exponential growth phase, were plotted against the corresponding pH, NaCl, ASW, or temperature value. Growth is best at pH 7, NaCl concentrations up to 5% are tolerated and ASW is tolerated from 27.5 to 230%. Growth is best at 28°C. Each dot represents the mean of triplicate measurements.
Figure S5. Number of giant genes >5 kb amongst the analyzed strains ordered by declining phylogenetic relationship. Strain NH11T comprises 45 giant genes >5 kb, which is the second highest number after Zavarzinella Formosa with 60 giant genes above this size.
Figure S6. Number of giant genes >20 kb amongst the analyzed strains ordered by declining phylogenetic relationship. Strain NH11T comprises five giant genes >20 kb, which is the second highest number after Rhodopirellula sp. K833 with seven giant genes above this threshold.
Figure S7. Number of giant genes >10 kb amongst the analyzed strains ordered by declining phylogenetic relationship. With 19 giant genes >10 kb, Strain NH11T comprises the highest number above this threshold.
This work was funded by the Deutsche Forschungsgemeinschaft to CJ (JO 893/3-1) and to JO (TRR 51/2 TA07). We thank Dr. Cathrin Spröer for support and Simone Severitt, Nicole Heyer, and Gabriele Pötter for skilled technical assistance. We thank Elizabeth Benecci from the Harvard Medical School's EM Facility for technical assistance.
Aghnatios, R., Cayrou, C., Garibal, M., Robert, C., Azza, S., Raoult, D., et al. (2015). Draft genome of Gemmata massiliana sp. nov, a water-borne Planctomycetes species exhibiting two variants. Stand. Genomic Sci. 10, 120. doi: 10.1186/s40793-015-0103-0
Bernard, T., Jebbar, M., Rassouli, Y., Himdi-Kabbab, S., Hamelin, J., and Blanco, C. (1993). Ectoine accumulation and osmotic regulation in Brevibacterium linens. Microbiology 139, 129–136. doi: 10.1099/00221287-139-1-129
Bland, C., Ramsey, T. L., Sabree, F., Lowe, M., Brown, K., Kyrpides, N. C., et al. (2007). CRISPR recognition tool (CRT): a tool for automatic detection of clustered regularly interspaced palindromic repeats. BMC Bioinformatics 8:209. doi: 10.1186/1471-2105-8-209
Bruns, A., Philipp, H., Cypionka, H., and Brinkhoff, T. (2003). Aeromicrobium marinum sp. nov., an abundant pelagic bacterium isolated from the German Wadden Sea. Int. J. Syst. Evol. Microbiol. 53, 1917–1923. doi: 10.1099/ijs.0.02735-0
Cayrou, C., Raoult, D., and Drancourt, M. (2010). Broad-spectrum antibiotic resistance of Planctomycetes organisms determined by Etest. J. Antimicrob. Chemother. 65, 2119–2122. doi: 10.1093/jac/dkq290
Dhillon, B. K., Laird, M. R., Shay, J. A., Winsor, G. L., Lo, R., Nizam, F., et al. (2015). IslandViewer 3: more flexible, interactive genomic island discovery, visualization and analysis. Nucleic Acids Res. 43, W104–W108. doi: 10.1093/nar/gkv401
Dragulescu, A. A. (2014). Read, Write, Format Excel 2007 and Excel 97/2000/XP/2003 Files. R package version 0.5.7 ed. Available online at: https://cran.r-project.org/web/packages/xlsx/
Fuerst, J. A., Gwilliam, H. G., Lindsay, M., Lichanska, A., Belcher, C., Vickers, J. E., et al. (1997). Isolation and molecular identification of planctomycete bacteria from postlarvae of the giant tiger prawn, Penaeus monodon. Appl. Environ. Microbiol. 63, 254–262.
Fuerst, J. A., and Sagulenko, E. (2010). Protein uptake by bacteria: an endocytosis-like process in the planctomycete Gemmata obscuriglobus. Commun. Integr. Biol. 3, 572–575. doi: 10.4161/cib.3.6.13061
Glöckner, F. O., Kube, M., Bauer, M., Teeling, H., Lombardot, T., Ludwig, W., et al. (2003). Complete genome sequence of the marine planctomycete Pirellula sp. strain 1. Proc. Natl. Acad. Sci. U.S.A. 100, 8298–8303. doi: 10.1073/pnas.1431443100
Jeske, O., Jogler, M., Petersen, J., Sikorski, J., and Jogler, C. (2013). From genome mining to phenotypic microarrays: Planctomycetes as source for novel bioactive molecules. Antonie Van Leeuwenhoek 104, 551–567. doi: 10.1007/s10482-013-0007-1
Jeske, O., Surup, F., Ketteniß, M., Rast, P., Förster, B., Jogler, M., et al. (2016). Developing techniques for the utilization of Planctomycetes as producers of bioactive molecules. Front. Microbiol. 7:1242. doi: 10.3389/fmicb.2016.01242
Jogler, C., Glöckner, F. O., and Kolter, R. (2011). Characterization of Planctomyces limnophilus and development of genetic tools for its manipulation establish it as a model species for the phylum Planctomycetes. Appl. Environ. Microbiol. 77, 5826–5829. doi: 10.1128/AEM.05132-11
Jogler, C., Waldmann, J., Huang, X., Jogler, M., Glöckner, F. O., Mascher, T., et al. (2012). Identification of proteins likely to be involved in morphogenesis, cell division, and signal transduction in Planctomycetes by comparative genomics. J. Bacteriol. 194, 6419–6430. doi: 10.1128/JB.01325-12
Jogler, M., Chen, H., Simon, J., Rohde, M., Busse, H. J., Klenk, H. P., et al. (2013). Description of Sphingorhabdus planktonica gen. nov., sp. nov. and reclassification of three related members of the genus Sphingopyxis in the genus Sphingorhabdus gen. nov. Int. J. Syst. Evol. Microbiol. 63, 1342–1349. doi: 10.1099/ijs.0.043133-0
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., Mcanulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kartal, B., de Almeida, N. M., Maalcke, W. J., Op den Camp, H. J., Jetten, M. S., and Keltjens, J. T. (2013). How to make a living from anaerobic ammonium oxidation. FEMS Microbiol. Rev. 37, 428–461. doi: 10.1111/1574-6976.12014
König, E., Schlesner, H., and Hirsch, P. (1984). Cell-wall studies on budding bacteria of the Planctomyces/ Pasteuria group and on a Prosthecomicrobium sp. Arch. Microbiol. 138, 200–205. doi: 10.1007/BF00402120
Kuykendall, L. D., Roy, M. A., Apos Neill, J. J., and Devine, T. E. (1988). Fatty acids, antibiotic resistance, and deoxyribonucleic acid homology groups of Bradyrhizobium japonicum. Int. J. Syst. Evol. Microbiol. 38, 358–361. doi: 10.1099/00207713-38-4-358
Lagesen, K., Hallin, P., Rødland, E. A., Staerfeldt, H. H., Rognes, T., and Ussery, D. W. (2007). RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 35, 3100–3108. doi: 10.1093/nar/gkm160
Lechner, M., Findeiss, S., Steiner, L., Marz, M., Stadler, P. F., and Prohaska, S. J. (2011). Proteinortho: detection of (co-)orthologs in large-scale analysis. BMC Bioinformatics 12:124. doi: 10.1186/1471-2105-12-124
Markowitz, V. M., Chen, I. M., Palaniappan, K., Chu, K., Szeto, E., Grechkin, Y., et al. (2012). IMG: the Integrated Microbial Genomes database and comparative analysis system. Nucleic Acids Res. 40, D115–D122. doi: 10.1093/nar/gkr1044
Mitchell, A., Chang, H. Y., Daugherty, L., Fraser, M., Hunter, S., Lopez, R., et al. (2015). The InterPro protein families database: the classification resource after 15 years. Nucleic Acids Res. 43, D213–D221. doi: 10.1093/nar/gku1243
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114
Pastor, J. M., Salvador, M., Argandoña, M., Bernal, V., Reina-Bueno, M., Csonka, L. N., et al. (2010). Ectoines in cell stress protection: uses and biotechnological production. Biotechnol. Adv. 28, 782–801. doi: 10.1016/j.biotechadv.2010.06.005
Petroni, G., Spring, S., Schleifer, K. H., Verni, F., and Rosati, G. (2000). Defensive extrusive ectosymbionts of Euplotidium (Ciliophora) that contain microtubule-like structures are bacteria related to Verrucomicrobia. Proc. Natl. Acad. Sci. U.S.A. 97, 1813–1817. doi: 10.1073/pnas.030438197
Pruesse, E., Quast, C., Knittel, K., Fuchs, B. M., Ludwig, W., Peplies, J., et al. (2007). SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 35, 7188–7196. doi: 10.1093/nar/gkm864
R Core Team (2015). R: A Language Environment for Statistical Computing. Available online at: https://www.r-project.org/
Rutherford, K., Parkhill, J., Crook, J., Horsnell, T., Rice, P., Rajandream, M. A., et al. (2000). Artemis: sequence visualization and annotation. Bioinformatics 16, 944–945. doi: 10.1093/bioinformatics/16.10.944
Scheuner, C., Tindall, B. J., Lu, M., Nolan, M., Lapidus, A., Cheng, J. F., et al. (2014). Complete genome sequence of Planctomyces brasiliensis type strain (DSM 5305(T)), phylogenomic analysis and reclassification of Planctomycetes including the descriptions of Gimesia gen. nov., Planctopirus gen. nov. and Rubinisphaera gen. nov. and emended descriptions of the order Planctomycetales and the family Planctomycetaceae. Stand. Genomic. Sci. 9:10. doi: 10.1186/1944-3277-9-10
Stephens, R. S., Kalman, S., Lammel, C., Fan, J., Marathe, R., Aravind, L., et al. (1998). Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis. Science 282, 754–759.
Taylor, D. W., Zhu, Y., Staals, R. H., Kornfeld, J. E., Shinkai, A., van der Oost, J., et al. (2015). Structures of the CRISPR-Cmr complex reveal mode of RNA target positioning. Science 348, 581–585. doi: 10.1126/science.aaa4535
Thorvaldsdóttir, H., Robinson, J. T., and Mesirov, J. P. (2013). Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 14, 178–192. doi: 10.1093/bib/bbs017
van Niftrik, L., Geerts, W. J., van Donselaar, E. G., Humbel, B. M., Webb, R. I., Harhangi, H. R., et al. (2009). Cell division ring, a new cell division protein and vertical inheritance of a bacterial organelle in anammox planctomycetes. Mol. Microbiol. 73, 1009–1019. doi: 10.1111/j.1365-2958.2009.06841.x
van Teeseling, M. C., Mesman, R. J., Kuru, E., Espaillat, A., Cava, F., Brun, Y. V., et al. (2015). Anammox Planctomycetes have a peptidoglycan cell wall. Nat. Commun. 6, 6878. doi: 10.1038/ncomms7878
Wagner, M., and Horn, M. (2006). The Planctomycetes, Verrucomicrobia, Chlamydiae and sister phyla comprise a superphylum with biotechnological and medical relevance. Curr. Opin. Biotechnol. 17, 241–249. doi: 10.1016/j.copbio.2006.05.005
Wallner, S. R., Bauer, M., Würdemann, C., Wecker, P., Glöckner, F. O., and Faber, K. (2005). Highly enantioselective sec-alkyl sulfatase activity of the marine planctomycete Rhodopirellula baltica shows retention of configuration. Angew. Chem. Int. Ed. Engl. 44, 6381–6384. doi: 10.1002/anie.200501955
Ward, N. L. (2010). “Phylum XXV. Planctomycetes Garrity and Holt 2001, 137 emend,” Ward (this Volume) in Bergey's Manual® of Systematic Bacteriology: Volume Four The Bacteroidetes, Spirochaetes, Tenericutes (Mollicutes), Acidobacteria, Fibrobacteres, Fusobacteria, Dictyoglomi, Gemmatimonadetes, Lentisphaerae, Verrucomicrobia, Chlamydiae, and Planctomycetes, eds N. R. Krieg, J. T. Staley, D. R. Brown, B. P. Hedlund, B. J. Paster, N. L. Ward, W. Ludwig, and W. B. Whitman (New York, NY: Springer), 879–925.
Weber, T., Blin, K., Duddela, S., Krug, D., Kim, H. U., Bruccoleri, R., et al. (2015). antiSMASH 3.0-a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res. 43, W237–W243. doi: 10.1093/nar/gkv437
Wittmann, J., Dreiseikelmann, B., Rohde, C., Rohde, M., and Sikorski, J. (2014). Isolation and characterization of numerous novel phages targeting diverse strains of the ubiquitous and opportunistic pathogen Achromobacter xylosoxidans. PLoS ONE 9:e86935. doi: 10.1371/journal.pone.0086935
Yarza, P., Yilmaz, P., Pruesse, E., Glöckner, F. O., Ludwig, W., Schleifer, K. H., et al. (2014). Uniting the classification of cultured and uncultured bacteria and archaea using 16S rRNA gene sequences. Nat. Rev. Microbiol. 12, 635–645. doi: 10.1038/nrmicro3330
Keywords: planctomycetes, Fuerstia marisgermanicae, giant genes, cell division, animal associated
Citation: Kohn T, Heuer A, Jogler M, Vollmers J, Boedeker C, Bunk B, Rast P, Borchert D, Glöckner I, Freese HM, Klenk H-P, Overmann J, Kaster A-K, Rohde M, Wiegand S and Jogler C (2016) Fuerstia marisgermanicae gen. nov., sp. nov., an Unusual Member of the Phylum Planctomycetes from the German Wadden Sea. Front. Microbiol. 7:2079. doi: 10.3389/fmicb.2016.02079
Received: 30 September 2016; Accepted: 08 December 2016;
Published: 22 December 2016.
Edited by:Damien Paul Devos, Pablo de Olavide University, Spain
Reviewed by:Svetlana N. Dedysh, Winogradsky Institute of Microbiology, Russian Academy of Science, Russia
Johannes Gijsbrecht Kuenen, Delft University of Technology, Netherlands
Copyright © 2016 Kohn, Heuer, Jogler, Vollmers, Boedeker, Bunk, Rast, Borchert, Glöckner, Freese, Klenk, Overmann, Kaster, Rohde, Wiegand and Jogler. 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) or licensor 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: Christian Jogler, email@example.com