Aureliella helgolandensis gen. nov., sp. nov., a novel Planctomycete isolated from a jellyfish at the shore of the island Helgoland

A novel planctomycetal strain, designated Q31aT, was isolated from a jellyfish at the shore of the island Helgoland in the North Sea. The strain forms lucid white colonies on solid medium and displays typical characteristics of planctomycetal strains, such as division by budding, formation of rosettes, presence of crateriform structures, extracellular matrix or fibre and a holdfast structure. Q31aT is mesophilic (temperature optimum 27 °C), neutrophilic (pH optimum 7.5), aerobic and heterotrophic. A maximal growth rate of 0.017 h− 1 (generation time of 41 h) was observed. Q31aT has a genome size of 8.44 Mb and a G + C content of 55.3%. Phylogenetically, the strain represents a novel genus and species in the recently introduced family Pirellulaceae, order Pirellulales, class Planctomycetia. We propose the name Aureliella helgolandensis gen. nov., sp. nov. for the novel species, represented by Q31aT (= DSM 103537T = LMG 29700T) as the type strain.


Introduction
Planctomycetes is a bacterial phylum displaying exceptional physiological and morphological features (Fuerst 1995;Lage et al. 2019;Staley et al. 1992;Ward 2010;Wiegand et al. 2020). Members of this phylum can be found in a variety of different habitats on earth, while the majority of species characterised so far have been isolated from aquatic environments (Wiegand et al. 2018). Phylogenetically, the phylum Planctomycetes, along with Chlamydiae, Verrucomicrobia and others, forms the PVC superphylum, which is of medical and biotechnological relevance Wagner and Horn 2006). Planctomycetes have large genomes of up to 12.4 Mb and usually 40-50% of the annotated genes are of unknown function (Jogler et al. 2012;Ravin et al. 2018;Wiegand et al. 2020). The taxonomy of the phylum Planctomycetes was recently revised (Dedysh et al. 2019). No changes were made to the original division of the phylum into the classes Candidatus Brocadiae, Phycisphaerae and Planctomycetia, while the latter is now further subdivided into the orders Gemmatales, Isosphaerales, Pirellulales and Planctomycetales.
Species belonging to the class Planctomycetia have been isolated from various marine biotic and abiotic surfaces (Boersma et al. 2019;Bondoso et al. 2014Bondoso et al. , 2017Kallscheuer et al. 2020;Peeters et al. 2020;Vollmers et al. 2017), on which they can be highly abundant (Bengtsson and Øvreås 2010). Due to the oligotrophic nature of marine environments, such species are suggested to digest complex carbon substrates, e.g. from biotic surfaces to which they frequently attach (Jeske et al. 2013;Lachnit et al. 2013). The observed dominance of planctomycetal species e.g. on algal surfaces is astonishing given their slow growth compared to other natural competitors in this ecological niche, e.g. members of the 'Roseobacter group' (Frank et al. 2014;Wiegand et al. 2018). The underlying mechanisms allowing species of the class Planctomycetia to compensate for slower growth are not understood, but might include the capability to produce bioactive small molecules (Graca et al. 2016;Jeske et al. 2016;Kallscheuer et al. 2019c), their resistance against several antibiotics (Cayrou et al. 2010;Godinho et al. 2019) and a specialised machinery for the uptake and intracellular digestion of complex polysaccharides. The latter is suspected to be facilitated by unique pili-forming crateriform structures and an extremely enlarged periplasmic space .
In the last decade, novel microscopic techniques and genetic tools for Planctomycetes (Jogler et al. 2011;Jogler and Jogler 2013;Rivas-Marin et al. 2016) allowed for a more detailed analysis of the cell envelope architecture of these bacteria. Planctomycetes were shown to possess peptidoglycan (Jeske et al. 2015;van Teeseling et al. 2015), supporting the assumption that all free-living bacteria have a peptidoglycan cell wall. The cell envelope architecture of Planctomycetes is therefore similar to that of Gramnegative bacteria Devos 2014). However, the phylum Planctomycetes is still exceptional. Characterised members were found to lack canonical divisome proteins including the otherwise universal FtsZ (Jogler et al. 2012;Pilhofer et al. 2008). Members of the class Phycisphaerae divide by binary fission, while budding is performed by species in the class Planctomycetia .
To extend the current collection of axenic cultures and as a basis to further study the interesting cell biology and metabolism of Planctomycetes here we describe strain Q31a T isolated from a jellyfish close to the island Helgoland in the North Sea. According to the results of our analyses, the strain represents a novel species and genus in the recently proposed family Pirellulaceae, order Pirellulales in the class Planctomycetia (Dedysh et al. 2019).

Isolation of the novel strain and cultivation
For the isolation and cultivation of strain Q31a T , M1H NAG ASW medium was used (Kallscheuer et al. 2019a). Strain Q31a T was isolated from a dead common jellyfish (Aurelia aurita) found at the shore of Helgoland Island (exact location 54.188 N 7.875 E) on the 5th of June 2013. A piece of the tentacles was cut off and then swabbed over a M1H NAG ASW plate containing 8 g/L gellan gum, 1000 mg/L streptomycin, 200 mg/L ampicillin and 20 mg/L cycloheximide, which was subsequently incubated at 20°C for three weeks. The 16S rRNA gene of obtained colonies was amplified by PCR and sequenced following an established protocol ). This step was performed in order to check whether isolated strains represent members of the phylum Planctomycetes. DNA extraction and genome sequencing are described in a previously published study ).

Determination of pH and temperature optimum
Cultivation for determination of the pH optimum was performed in M1H NAG ASW medium and for ensuring a stable pH 100 mM HEPES was used for cultivations at pH 7, 7.5 and 8. For cultivation at pH 5 and 6 HEPES was replaced by 100 mM 2-(Nmorpholino)ethanesulfonic acid (MES), whereas 100 mM N-cyclohexyl-2-aminoethanesulfonic acid (CHES) served as a buffering agent at pH 9 and 10. Cultivations for determination of the pH optimum were performed at 28°C. Cultivations for determination of the temperature optimum were performed in standard M1H NAG ASW medium at pH 7.5.

Microscopy protocols
Phase contrast and field emission scanning electron microscopy were performed as previously described (Boersma et al. 2019).

Genome information
Genome and 16S rRNA gene sequence of strain Q31a T are available from GenBank under accession numbers CP036298 and MK559992, respectively. Numbers of carbohydrate-active enzymes were obtained from the CAZY database (Lombard et al. 2014). Gene clusters potentially involved in the production of secondary metabolites were determined using antiSMASH 4.0 (Blin et al. 2017).

Phylogenetic analysis
16S rRNA gene-based phylogeny was computed for Q31a T , the type strains of all described planctomycetal species (assessed in January 2020) and all isolates recently described (Boersma et al. 2019;Kallscheuer et al. 2019aKallscheuer et al. , b, d, 2020Kohn et al. 2019;Peeters et al. 2020;Rensink et al. 2020). The 16S rRNA gene sequences were aligned with SINA (Pruesse et al. 2012) and the phylogenetic inference was calculated with RAxML (Stamatakis 2014) with a maximum likelihood approach with 1000 bootstraps, nucleotide substitution model GTR, gamma distributed rate variation and estimation of proportion of invariable sites (GTRGAMMAI option). Three 16S rRNA genes of bacterial strains from the PVC superphylum but outside of the phylum Planctomycetes (GenBank accession numbers AJ229235, KC665948 and NR_027571) were used as outgroup. For the multilocus sequence analysis (MLSA) the unique singlecopy core genome of the analysed genomes (GenBank acc. no. CP036298) was determined with pro-teinortho5 (Lechner et al. 2011) with the 'selfblast' option enabled. The protein sequences of the resulting orthologous groups were aligned using MUSCLE v.3.8.31 (Edgar 2004). After clipping, partially aligned C-and N-terminal regions and poorly aligned internal regions were filtered using Gblocks (Castresana 2000). The final alignment was concatenated and clustered using the maximum likelihood method implemented by RaxML (Stamatakis 2014) with the 'rapid bootstrap' method and 500 bootstrap replicates.
Four planctomycetal genomes from different families were used as outgroup. The average nucleotide identity (ANI) was calculated using OrthoANI (Lee et al. 2016). The average amino acid identity (AAI) was calculated using the aai.rb script of the enveomics collection (Rodriguez-R and Konstantinidis 2016) and the percentage of conserved proteins (POCP) was calculated as described (Qin et al. 2014). The rpoB nucleotide sequences were taken from publicly available planctomycetal genome annotations and the sequence identities were determined as described (Bondoso et al. 2013). Upon extracting only those parts of the sequence that would have been sequenced with the described primer set, the alignment and matrix calculation was performed with Clustal Omega (Sievers et al. 2011).

Phylogenetic analysis
Based on 16S rRNA gene phylogeny and whole genome-based MLSA, strain Q31a T groups within the planctomycetal family Pirellulaceae (Fig. 1). Within this family, its current closest neighbours on 16S rRNA gene level are Mariniblastus fucicola and Pirellula staleyi, further close neighbours are Blastopirellula sp. and Bremerella sp. However, supporting bootstrap values in this clade are sometimes rather low and no definitive closest neighbourhood could be determined by MLSA analysis. Therefore, strain Q31a T was compared to all described genera within the family Pirellulaceae (Fig. 2).
16S rRNA gene sequence identity analysis (Fig. 2) shows that all minimal identities between the novel strain Q31a T and the nine most related genera are notably below the genus threshold of 94.5% that might place Q31a T in any of these taxa (Yarza et al. 2014). The similarity values for strain Q31a T and its relatives are also below the genus threshold values used with rpoB nucleotide sequences identities (75.5-78.0%) (Kallscheuer et al. 2019d), AAI (60-80%) (Konstantinidis and Tiedje 2005) and POCP (50%) (Qin et al. 2014) (Fig. 2). With all used methods suggesting the placement of strain Q31a T in a novel genus, we conclude the strain represents a novel genus and species within the family Pirellulaceae, for which we propose the name Aureliella helgolandensis gen. nov., sp. nov.

Morphological and physiological analyses
Basic features of strain Q31a T regarding cell morphology and mechanism of cell division are summarised in Table 1. As we could not identify a current closest relative of Q31a T , the morphological and genomic features were compared to all strains identified as potential candidates during the phylogenetic analysis. Morphological features of Q31a T cells harvested during the exponential growth phase were analysed using phase contrast and scanning electron microscopy (Fig. 3). Cells of strain Q31a T are 1.9 ± 0.2 9 1.0 ± 0.2 lm in size and acorn-shaped ( Fig. 3a-c). Cells form aggregates of typically 8-25 cells, which in turn are often loosely connected to each other (Fig. 3d). Cells divide by polar budding (Fig. 3a). Extracellular matrix or fibre originates from one pole. At this pole crateriform structures can also be observed, which cover around 10-20% of the cell surface. Daughter cells of Q31a T have the same shape as the mother cell. The strain follows a dimorphic lifecycle involving sessile mother cells and flagellated daughter cells. Colonies of strain Q31a T lack pigmentation and have a lucid white colour.
Despite aggregate formation, measurement of optical densities (OD 600 ) in liquid cultures was possible. In M1H NAG ASW medium, Q31a T was found to be able to grow in a temperature range of 10-33°C and in a pH range of 6.0-8.0 (Fig. 4). Optimal growth was observed at 27°C and pH 7.5. The maximal growth rate observed in M1H NAG ASW medium was 0.017 h -1 (Fig. 4), corresponding to a generation time of 41 h. Q31a T is an aerobic heterotroph. During comparison of preferred temperature and pH of Q31a T with the close relatives M. fucicola FC18 T , Blastopirellula marina DSM 3645 T , Roseimaritima ulvae UC8 T , P. staleyi DSM 6068 T and Crateriforma conspicua Mal65 T considerable differences were observed. The temperature optimum of Q31a T is between the optima of P. staleyi (24°C)/M. fucicola (25°C) and the other three strains (30-36°C) ( Table 1). The pH range for growth of Q31a T is narrow compared to ranges of 6.0-10.0 observed for R. ulvae and C. conspicua, but comparable to M. fucicola (for the other two strains no data was available). The lucid white colony colour of Q31a T indicates the lack of carotenoid formation of the strain, which is a common feature of B. marina and P. staleyi, but separates Q31a T from the pink-pigmented R. ulvae and C. conspicua (Table 1).

Genome-based analysis of metabolic features
Numbers of carbohydrate-active enzymes and secondary metabolite-associated genes clusters were analysed based on the genome sequences of Q31a T and type species of related genera (Table 2). These numbers give a first impression on the potential of Q31a T for degradation of complex and highly decorated polysaccharides and for production of bioactive small molecules. In total, Q31a T harbours 159 carbohydrate-active enzyme as currently listed in the CAZY database. This number is comparable to R. ulvae UC8 T , which also has a similar genome size. Although having a genome 1.3 Mb smaller than Q31a T , the highest number of 217 carbohydrate-active enzymes was observed in C. conspicua Mal65 T . This difference can mainly be attributed to the glycoside hydrolase family, since 52 enzymes were found in Q31a T and 121 in C. conspicua Mal65 T . Q31a T contains the highest number of enzymes of the carbohydrate esterase family of the species used for comparison. During analysis of secondary metabolite-associated gene clusters a heterogeneous distribution for the investigated species was observed. While 2-3 terpenoid-related clusters were found in all species, other clusters putatively involved in the production of ectoine, resorcinol or non-ribosomal peptides seem to be restricted to individual genera. Similar results were obtained for comparison of type I and type III polyketide synthases (PKSs), while type II PKSs appear to be absent from the compared genomes. The total numbers of predicted clusters is between 5 and 10, while higher numbers are not reflected by larger genomes in this case. Q31a T has the largest genome, but is ranked 3rd with regard to the number of gene clusters. In contrast, B. marina DSM 3645 T has the highest number of predicted gene clusters, but is amongst the species with the smallest genomes.

Conclusions
The performed comparison of morphological, physiological and genomic features supports the results of the phylogenetic analysis that Q31a T does not belong to the genera Mariniblastus, Pirellula, Blastopirellula, Rhodopirellula, Novipirellula, Rubripirellula, Bremerella, Crateriforma or Roseimaritima, but instead represents a new species belonging to a novel genus.  The analysis is based on GenBank accession numbers CP036298 (Q31a T ), CP042912 (Mariniblastus fucicola FC18 T ), GCF_000153105 (Blastopirellula marina DSM 3645 T ), CP042914 (Roseimaritima ulvae UC8 T ), GCF_000025185 (Pirellula staleyi DSM 6068 T ) and CP036319 (Crateriforma conspicua Mal65 T ). The genome of B. marina was not listed in the CAZY database Thus, we propose the name Aureliella helgolandensis gen. nov., sp. nov. for Q31a T and propose this species as the type species of the genus and Q31a T as the type strain of the novel species.
Members of the genus are Gram-negative, aerobic, mesophilic, neutrophilic and heterotrophic. Cells are acorn-shaped, have crateriform structures at one pole and divide by polar budding. The genus belongs to the family Pirellulaceae, order Pirellulales, class Planctomycetia, phylum Planctomycetes. The type species is Aureliella helgolandensis.
Aureliella helgolandensis (hel.go.lan.den'sis N.L. fem. adj. helgolandensis of Helgoland; corresponding to the origin of the strain from the German island Helgoland).
Cells are 1.9 ± 0.2 lm 9 1.0 ± 0.2 lm in size and form aggregates. Matrix or fibre originates from the budding pole and a holdfast structure is present at the opposite pole. Grows at 10-33°C (optimum 27°C) and at pH 6.0-8.0 (optimum 7.5). Colonies are lucid white. The genome of the type strain has a G ? C content of 55.3%.
The type strain is Q31a T (DSM 103537 T = LMG 29700 T ), isolated from a jellyfish (Aurelia aurita) on the shore of Helgoland Island in June 2013. The type strain genome (8.44 Mb, acc. no. CP036298) and 16S rRNA gene sequence (acc. no. MK559992) are available from GenBank.