Compartmentalization in PVC super-phylum: evolution and impact
© The Author(s). 2016
Received: 13 May 2016
Accepted: 2 August 2016
Published: 9 August 2016
The PVC super-phylum gathers bacteria from seven phyla (Planctomycetes, Verrucomicrobiae, Chlamydiae, Lentisphaera, Poribacteria, OP3, WWE2) presenting different lifestyles, cell plans and environments. Planctomyces and several Verrucomicrobiae exhibit a complex cell plan, with an intracytoplasmic membrane inducing the compartmentalization of the cytoplasm into two regions (pirellulosome and paryphoplasm). The evolution and function of this cell plan is still subject to debate. In this work, we hypothesized that it could play a role in protection of the bacterial DNA, especially against Horizontal Genes Transfers (HGT). Therefore, 64 bacterial genomes belonging to seven different phyla (whose four PVC phyla) were studied. We reconstructed the evolution of the cell plan as precisely as possible, thanks to information obtained by bibliographic study and electronic microscopy. We used a strategy based on comparative phylogenomic in order to determine the part occupied by the horizontal transfers for each studied genomes.
Our results show that the bacteria Simkania negevensis (Chlamydiae) and Coraliomargarita akajimensis (Verrucomicrobiae), whose cell plan were unknown before, are compartmentalized, as we can see on the micrographies. This is one of the first indication of the presence of an intracytoplasmic membrane in a Chlamydiae. The proportion of HGT does not seems to be related to the cell plan of bacteria, suggesting that compartmentalization does not induce a protection of bacterial DNA against HGT. Conversely, lifestyle of bacteria seems to impact the ability of bacteria to exchange genes.
Our study allows a best reconstruction of the evolution of intracytoplasmic membrane, but this structure seems to have no impact on HGT occurrences.
This article was reviewed by Mircea Podar and Olivier Tenaillon.
The PVC super-phylum gathers seven bacterial phyla (Planctomycetes, Verrucomicrobiae, Chlamydiae, Lentisphaerae, Poribacteria, OP3, WWE2) [1–4], and comports 37 species of bacteria whose genome was entirely sequenced. The monophyly of this super phylum have been discussed a lot in the last years, due to the difficulty to obtain a consensual phylogeny [5–14]. Recently, it seems that a global consensus was reached, that include all these bacteria in a same super-phylum [2, 5]. This idea is confirmed by the discovery in 2012 and 2014 of a molecular signature conserved in all PVC bacteria [1, 15]. The phylogenetic relations between PVC and other bacteria are still a controversial subject [5, 13]. PVC bacteria are present in an important variety of environments: water, soils, vertebrates, amoeba, insects… [5, 16–21].
The bacteria of this super-phylum show different interesting characteristics [11, 22]. Some of these features can be found in other bacteria: some Planctomycetes are implicated in carbon and nitrogen cycles [23, 24] or synthesize special sterols , many Chlamydiae are pathogens of mammals [26, 27], some Chlamydiae and Verrucomicrobiae are symbionts . But PVC present also genetic and cellular features unusual among bacteria: compartmentalization in Verrucomicrobiae and Planctomycetes [28, 29], absence of tubulin like protein FtZ in Planctomycetes and Chlamydiae [30, 31], crateiform surface and budding reproduction  in Planctomycetes. Among these characteristics, we were specifically interested in the compartmentalization of bacteria. This feature is characterized by a specific cell plan, and concerns all the Planctomycetes [29, 33, 34], some Verrucomicrobiae , one Lentisphaera and one Poribacteria . An intracytoplasmic membrane (ICM) separates the cytoplasm of bacteria into two distinct compartments: the pirellulosome inside (containing the DNA ), and the paryphoplasm outside (the size and shape of these two compartments varied a lot among PVC bacteria). This ICM is a lipid bilayer in contact with proteins [29, 32, 33] presenting structural similarities with some proteins from eukaryotic membranes, such as clathrins [37, 38]. Some compartmentalized bacteria also present a more complex cell plan, with the presence of another compartment, like the anamoxosome in Candidatus Kuenenia Stuttgartiensis [32, 39, 40], or a double internal membrane with ribosomes  in the same compartment that DNA  in Gemmata obscuriglobus . Compartmentalization in PVC bacteria is still debatable and three opinions are defended actually : 1- As some features of PVC bacteria are current in Eukaryota , this observation allows some people to assume that the compartmentalization of PVC is the precursor of Eukaryotic nuclei, but this idea is not very popular [34, 43–48]. 2- Another proposition is that the compartment is structurally and functionally similar to a nucleus but that these two structures appeared independently . 3- Some people identify the membrane as an invagination of a Gram negative external membrane .
Considering the reality of the existence of a compartmentalization in PVC bacteria, it would be interesting to determine the function of this specific cell plan; indeed, no previous studies had been able to determine it. Here we assumed that this membrane could be a protection against Horizontal Gene Transfers (HGT). This hypothesis was based on observations of the impact of different membranes already studied : the role of the internal membrane in Eukaryotes, or the influence of host membrane on genomes from intracellular bacteria (which limits the contact between foreign elements and DNA). We used a phylogenomic strategy of HGT detection to reconstruct the history of HGT during the evolution of PVC bacteria, before and after appearance and disappearance of compartmentalization.
Lifestyle of species and evolution of cell plan
PVC bacteria present two different cell plan, compartmentalized by an intra cytoplasmic membrane (ICM) or non compartmentalized. They also present two different lifestyles: allopatric lifestyle, for bacteria isolated from other microorganisms (obligate intracellular bacteria living in non amoeba cells) and sympatric, for bacteria in contact with other microorganisms. The studies already published allowed determining cell plan of 26 PVC bacteria: PVC super phylum presented 18 compartmentalized species distributed in three phyla (Verrucomicrobiae, Planctomycetes and Lentisphaerae) and eight non compartmentalized bacteria in one phyla (Chlamydiae). These bacteria were used in our study, but we also selected seven bacteria with an unknown cell plan, distributed in two phyla (Verrucomicrobiae and Chlamydiae). These data allowed the reconstruction of cell plan evolution with a good reliability, however it was still difficult to conclude about ancestral state at different nodes of the tree, especially in Chlamydiae phylum (their common ancestor presented a probability of 80.2 % to be compartmentalized). It was interesting to add new information in the dataset, in order to improve this reconstruction.
Relation between cell plan and genomes contents
The compartmentalization of bacteria could disturb only some particular mechanisms of HGT processes, leading to a low impact on HGT global proportions. The study of mobilome revealed that no element present proportions varying significantly according cell plan (Additional file 1). The lifestyle of bacteria seems to clearly influence the variations of some mobilome elements; indeed, the intracellular bacteria present more abundant proportions of transposases than extracellular bacteria (Kruskal test : p-value = 2.0*10-2). Conjugation genes are significantly over represented in the genomes of Intracellular allopatric bacteria, and underrepresented in extracellular compartmentalized bacteria (Kruskal test : p-value = 3.0*10-2).
Lifestyle of species and evolution of cell plan
The results of electronic microscopy allowed to defined S.negevensis and C.akajimensis as compartmentalized bacteria. The Chlamydiae is characterized by an important presence of internal membranes. This observation is probably related to the folds of the intracytoplasmic membrane, as it is visible for the G.obscuriglobus in the micrographies and the three dimensional reconstruction published by Santarella-Mellwig et al in 2010  and 2013  respectively. Conversely, the cell plan observed is different from that identified in the cryotomographies realized by Pilhofer et al in 2014 . The differences could be related to the variations in methodology used and in the growth stage observed. The possible compartmentalization of S.negevensis is very interesting because the cell plan is unknown in a large majority of sympatric Chlamydiae. C.akajimensis presents a more classical compartmentalized cell plan, as it was already observed in some PVC bacteria by Santarella-Mellwig in 2013  and Kuo-Chang Lee and al in 2009 . These micrographies allowed an improving in reconstruction of cell plan evolution among PVC bacteria, especially in Chlamydiae. They permitted to analyze the impact of ICM in the phyla of intracellular sympatric Chlamydiae, with two bacteria presenting a known cell plan and the probable presence of compartmentalization in their last common ancestor (with a probability of 71 %, against 57 % without the information about S.negevensis). The micrography of C.akajimensis completes the studies already realized on the Verrucomicrobiae and reinforces the probability of compartmentalization in the last common ancestor of Verrucomicrobiae (probability of compartmentalization = 83 %, against 79 % without the cell plan of C.akajimensis). These observations allowed to conclude with fewer uncertainties our study concerning the HGT in compartmentalized bacteria.
Cell plan evolution reconstruction indicated that compartmentalization appeared at the root of PVC super-phylum and disappeared once, in the group of allopatric intracellular bacteria. This disappearance seems to occurred simultaneously with bacteria conversion to intracellular allopatric lifestyle. Intracellular allopatric bacteria are isolated from other bacteria, so their genomes are protected against HGT. Conversely, extracellular and intracellular sympatric bacteria are more exposed to exchanges because they are in contact with many microorganisms. An hypothesis could be that the intracytoplasmic membrane plays a role in protection of genomes against HGT, as host membrane of intracellular bacteria or eukaryotic membrane . This intracytoplasmic membrane became useless in allopatric bacteria, due to the absence of HGT, leading to its disappearance. However, this hypothesis needs to be confirmed by data concerning genomes evolution in compartmentalized and non compartmentalized bacteria.
Relation between cell plan and genomes contents
The analysis of HGT in the different groups of bacteria highlighted the role of lifestyle on HGT process and the absence of impact of compartmentalization. The important decrease of HGT proportions observed in Chlamydiae after their conversion to the intracellular allopatric lifestyle is probably due to the physical isolation from other microorganisms, which prevents the opportunities of HGT . This agrees with previous studies showing that the predominant evolutionary process in intracellular bacteria is genome reduction leading to small genome sizes [53, 54] with the bacteria from amoeba as the exception [55, 56]. The similarity of HGT proportions between compartmentalized and non compartmentalized bacteria suggests that the presence on an intracytoplasmic membrane has no impact on HGT process. We could imagine that this absence of difference is due to a problem of our strategy but three observations demonstrate the efficiency of our method: Chlamydiae present an important rate of ancient HGT with eukaryotes, significantly higher than other bacteria (ANOVA test p-value = 3.3 × 10-6). This observation is consistent with the existing literature dedicated to transfers between plants and Chlamydiae [57, 58] that supports a role of Chlamydiae in chloroplast endosymbiosis; The HGT proportions detected in two Planctomycetes (Planctomyces maris and Candidatus Kuenenia stuttgartiensis) correspond to these observed by Kamneva et al in 2012 ; Spirochaetes seem to exchanges more with Firmicutes as it was already detected in previous studies [59, 60].
A similar HGT proportion in compartmentalized and non-compartmentalized bacteria does not necessary contradict the hypothesis of protection against HGT by compartmentalization. Indeed, the intracytoplasmic membrane of PVC bacteria can represents a barrier to HGT by disturbing only some transfers mechanisms mediated by mobilome elements. This perturbation could lead or not to a decrease of HGT level, depending if the deficiency of some mechanisms are offset by the others. The absence of difference in mobilome elements quantities and proportions between the different groups of bacteria indicated that the compartmentalization has no impact, not only on HGT proportions, but also on elements implicated in HGT process. The lifestyle seems to have an influence on two mobilomes elements, proportions of transposases and conjugation genes. The high proportion of transposases in intracellular allopatric bacteria could be related to the importance of transposable elements in genomes of intracellular bacteria [61, 62]. However, conjugation genes present other uses in bacteria than conjugation [63, 64]. So we can assume that the abundance of conjugation genes in bacteria has probably no influence on the HGT.
The function of compartmentalization in PVC bacteria remains unknown and the causes of intracytoplasmic membrane disappearance in intracellular isolated bacteria is still unresolved.
The different pictures obtained for Verrucomicrobiae and Chlamydiae allow a best definition of ancestral and modern states of compartmentalization and so, permit to reconstruct intracytoplasmic membrane evolution in PVC bacteria. But this compartmentalization seems have no significant impact on HGT quantity, HGT proportion or partners of transfers.
Bacteria selection and genomes recovery
We selected bacteria from 4 phyla of PVC super-phylum: Planctomycetes, Verrucomicrobia, Lentisphaerae, and Chlamydiae. Bacteria from three phylogenetically close phyla were chosen as negative controls, based on a reference tree : Bacteroidetes, Spirochaetes and Chlorobi. We retrieved the proteomes of all bacteria selected in genomes NCBI database (Additional file 2).
Lifestyle and cell plan determination
The lifestyle of selected bacteria is determined by a bibliographic study of each bacterium. Two types of lifestyles are known for bacteria, according the possibility of exchanges with other microorganisms: allopatric or sympatric bacteria [65, 66]. Sympatric bacteria are bacteria living in community with other microorganisms, they exchange easily genes by lateral transfer through their interactions. Allopatric bacteria are living in cells and are not in contact with other microorganisms, more specialized, with a reduced genome size and less genetic exchanges. They are obligate intracellular bacteria, living in non amoeba cells.
Cell plans of bacteria are determined thanks to transmission electron microscopy pictures already available in bibliography [28, 29, 34, 35] and microscopic observations of the bacteria whose cell plan is unknown. First, the species Simkania negevensis (DSM27360T, type strain) was grown within a culture of its amoebal host A. castellanii. At H6, and H16 post-infection, cultures were centrifuged at 2000 rpm/min during 10 min, and the pellets were fixed for electron microscopy in a fixative solution (2.5 % glutaraldehyde in 0.1 M sodium cacodylate buffer). The second species, Coraliomargarita akajimensis (DSM 45221T, type strain) was cultivated directly into Bacto Marine Broth medium (Difco). Growth occurred after 6 days of cultivation, and when the bacterial suspension reached the exponential growth phase, it was centrifuged at 5500 rpm/min during 30 min, and fixed for electron microscopy in the same fixative solution. An aliquot of each culture was kept for DNA extraction, followed by a standard 16S rRNA PCR and sequencing, in order to confirm the bacterial identification. For both bacteria, Electron microscopic observations were realized in different steps: for embedding, cells were fixed for 1 h with glutaraldehyde 2.5 % in 0.1 M sodium cacodylate buffer and washed three times. Cells were post-fixed for 1 h with 1 % OsO4 diluted in 0.2 M Potassium hexa-cyanoferrate (III) / 0.1 M sodium cacodylate solution. After four 5 min washes with distilled water, cells were gradually dehydrated with ethanol by successive 10 min baths in 30, 50, 70, 96, 100 and 100 % ethanol. Substitution was achieved by successively placing the cells in 25, 50 and 75 % Epon solutions for 15 min. Cells were placed for 1 h in 100 % Epon solution and in fresh Epon 100 % over-night under vacuum at room-temperature. Polymerization took place with cells in fresh 100 % Epon for 24 h at 60°. Ultrathin 70 nm sections were cut with a UC7 ultramicrotome (Leica) and placed on HR25 300 Mesh Copper/Rhodium grids (TAAB, UK). Electron micrographs were obtained on a Tecnai G20 TEM operated at 200 keV equipped with a 4096 × 4096 pixels resolution Eagle camera (FEI)).
Phylogeny and cell plan evolution reconstruction
We reconstructed the species tree of PVC bacteria and Bacteroidetes-Chlorobi-Spirochaetes thanks to 12 markers common to the 65 species (three ribosomal proteins (16S, 23S, 30S), two elongation factors (Tu and Ts), two DNA polymerase subunits, CTP synthetase, and four tRNA ligases). The selected sequences were concatenated and Mega5  software is used to perform the phylogeny (alignment by Muscle algorithm, manual removal of non-conserved positions, tree building with NJ and ML methods with 150 bootstraps, comparison of phylogenies and selection of a consensus (Additional file 3)).
Thanks to information about compartmentalization in bacteria, we reconstructed the ancestral state of cell plan with parsimony method (Mesquite software  and Phytools package on R software . Phytools allowed the reconstruction of ancestral sates of discrete character by a method of Maximum likelihood). The tree obtained allowed us to date the different events occurred during evolution of the phyla studied (appearance and disappearance of compartmentalization, conversion to intracellular allopatric or sympatric lifestyle).
We studied 9 elements directly or indirectly implicated in the horizontal transfers: phages (complete or incomplete), conjugation genes and plasmids, these three elements are involved in entrance of foreign sequences in bacteria; Transposases, integrases and CRISP (candidate or confirmed), regulating positively or negatively the sequences integration in recipient genome; and tRNA used for sequences translation. Different databases were used as RepBase  and CRISPfinder  to detect the mobilome elements (the complete list of databases used is presented in the Additional file 4). For each element we counted the quantity in each species; if this quantity is related to the genome size (CRISPs, transposases, integrases and conjugation genes), we calculated the proportion of these elements in the genomes. We used the same statistical tests as those used for HGT analysis, in order to determine the elements overrepresented in the different classes of bacteria.
OrthoMCL  was used for construction of orthologous groups. Genes absent to all orthologous groups are either acquired by HGT or generated de novo (ORFans). Blast against nr database allowed the identification genes presenting no identifiable homologous genes (e-value < 10-4 and coverage > 50 %), considered as ORFans.
Phylopattern  pipeline allowed the detection of genetic events in four steps: comparison between species tree and each orthologous groups; detection of missing species in orthologous groups; reconstruction of ancestral state for each proteins, based on the pattern of presence/absence of sequences in modern species; identification of two types of genetic events: gains and losses. Gains could be HGT, de novo genes, or artifacts. We focused on gene gains and Blast  each sequences against nr database (NCBI). HGT candidates were identified thanks to filtering on Blast results: For each blast results, sequences derived from phylum where gene gain was detected were removed (for example, if gain was identified in Akkermansia muciniphila (Verrucomicrobiae), all sequences in Blast results derived from Verrucomicrobiae were removed). Sequences with e-value > 10-5, coverage < 60 % or identities < 30 % were also removed. Then, species of the first ten sequences remaining in Blast result were identified. If species identified do not belong to one of the sister phyla studied, gene gain is probably an HGT (for example if gain was identified in Akkermansia muciniphila, gain is considered as HGT only if the first ten hits do not belong to the phyla Planctomycetes or Chlamydiae). The position on gains in species tree allows to differentiate HGT occurred only in modern species (specific HGT) and HGT occurred in the ancestor (non specific HGT). We calculated the quantity of HGT and the percentage of proteins present in each species, due to transfers. For each HGT detected (recent or ancient), we determined with which organism the transfers were realized.
For each elements studied (HGT quantity, proportions, partners and mobilome), we performed a comparison of variances between the different groups of bacteria (based on cell plan). We tested the normality of data distribution by the shapiro-Wilk test and the homogeneity of variance by Levene test. These tests were followed by a comparison of variance between the different classes, in order to evaluate if there was a significant difference between them (Kruskal test for data without a normal distribution and ANOVA test for data with a normal distribution). Tests of Nemenuyi (data not normally distributed) or Tukey (data normally distributed) were performed to obtain a comparison of each pairs of classes. We realized also these statistical analysis with bacterial groups based on the phyla, to determine if classes based on phyla present the same results that classes based on cell plan. If it had been the case, it would have been impossible to determine if differences observed were related to cell plan or to phylogenetic relations.
Reviewer 1 (Mircea Podar)
“In their study, Pinos et al tested their hypothesis that cellular compartmentalization, observed in some representatives of the PVC super phylum, is correlated to the frequency of horizontal gene transfer. The authors performed ultrastructural analyses of a representative from the Chlamydia and one from Verrucomicrobiae, which were shown to be compartmentalized. That information was combined with results from literature on other PVC species and was correlated with phylogenomic analyses on HGT frequencies, based on completed genomes. The result indicates that compartmentalization is in fact not significantly impacting HGT, which has a stronger correlation with the organismal life-style.
The study is interesting and original, even though the conclusion did not support the original hypothesis. It does provide a foundation for additional tests and focused study of related organisms/genomes, which may increase the signal to link compartmentalization to more subtle evolutionary genomic effects.”
Reviewer 2 (Olivier Tenaillon)
“In the manuscript entitled “Compartmentalization in PVC super-phylum: evolution and impact”, Sandrine Pinos and colleagues use both microscopy and a phylogenetic approach to study Compartmentalization and its potential impact on horizontal gene transfer.
The idea is interesting, but the paper is not very clear in its present form. The usage of microscopy is not introduced properly, statistics and methods are over-simplified and impossible to reproduce. I would consider therefore that the paper needs very substantial rewriting and the phylogenetic analysis should include uncertainty in the phylogeny.”
Reviewer recommendations to authors
Reviewer 1 (Mircea Podar): The manuscript needs extensive work in correcting the many language errors (spelling, syntax, grammar). It also does not provide sufficient detail to evaluate the power of the phylogenomic analyses that were performed. While I like the concept of this study and the overall approach, the manuscript requires a complete overhaul.
Reviewer 2 (Olivier Tenaillon): How robust is the tree? Before being able to conclude on the ancestry of a trait, the robustness of the phylogeny has to be tested. In my experience a 16S phylogeny is not very robust and much more loci should be involved in the reconstruction, bootstrap values should be presented. Character reconstruction should also take into account the tree uncertainty. This is critical for both figures and both character mapping procedures.
Reviewer 2 (Olivier Tenaillon): line 219: “If species identified do not belong to one of the close phyla studied, gene gain is probably an HGT” how close is close, the threshold use should be clearly defined?”
Author’s response: The different comments and questions concerning our phylogenomic strategy indicate an important problem in our presentation of this method and its results. In order to correct these problems, we expanded the methodology part of phylogenomic strategy, by adding indications and some examples (we detailed some points, as what we want to say by “If species identified do not belong to one of the close phyla studied, gene gain is probably an HGT”, in this case, “close phyla” refers to species belonging to the sister phylum). We also added, in the Additional file 3, phylogenies of PVC bacteria and Bacteroidetes-Spirochaetes-Chlorobi, with the indication of bootstrap values. These phylogenies are more robust than those of the first version, due to the use of ten supplementary markers, for alignment and phylogenetic reconstruction. We want to underline that some results obtained by our phylogenomic method are convergent with results already observed in previous studies (line 155–161), this is an indication of the relative efficiency of our method.
The improvement of the phylogeny allowed a better reconstruction of cell plan evolution in PVC bacteria. However, In order to take into account the problem of phylogenetic uncertainties, we used another method of ancestral reconstruction that not only considers these uncertainties, but also provides the probabilities of each state of cell plan at each nodes of the tree. The results of this analysis were included in the Figure 2, they led to a better vision of PVC bacteria evolution and highlighted the importance of microscopic data.
Reviewer 2 (Olivier Tenaillon): Introduction - The link between horizontal gene transfer and compartmentalization is not clearly stated in the introduction. Two arguments are mentioned, the second one a physical limit is somehow simple to understand, but the first one remains enigmatic: “the presence or absence of compartmentalization in the different PVC bacteria and the reconstruction of ancestral states”. What is meant by this sentence is not clear to me. I think developing more on the subject would improve the introduction.
Author’s response: This sentence is difficult to explain without the knowledge about the compartmentalization and the lifestyle of PVC bacteria, however, it could be long to detail all these information in introduction. Consequently, we decided to remove this sentence from the introduction, we presented and explained it in the results part. This sentence resumed the idea that if Chlamydiae are the only non compartmentalized PVC bacteria sequenced and also the only intracellular bacteria of this super-phylum, if the compartmentalization protects genomes against HGT, its lost in the intracellular species could be related to the isolation of these bacteria that induces a low proportion of HGT, and leads to the uselessness of compartmentalization. So, the bibliographic knowledge about the presence or absence of compartmentalization in species suggests the hypothesis of a role of the ICM in protection against HGT.
Reviewer 2 (Olivier Tenaillon): Results - The result section starts with the description of microscopic structures. It is not introduced in anyway. The two first paragraphs should be inverted: first mention that the compartmentalization is unknown in some species, then present the evidence for the ones lacking information.
Discussion - again no clear explanation of the interest of the microscopy.
Author’s response: As the result of the comments of reviewers, we decided to reorganize the results of microscopic analyzes to improve our presentation. We insisted on the importance of these data in our study by adding explanations of their use in ancestral reconstruction (71–80) and their role in the study of the relation between HGT and compartmentalization (97–99). We realized a more detailed ancestral state reconstruction, with a more complex method (analysis performed thanks to R software, with phylotools package, by a maximum likelihood method) that provided the probabilities of each cell plan states at each nodes of the tree. These analyzes highlighted the importance to add our data to the study, indeed our supplement of information allowed to improve significantly the probability of compartmentalization for several nodes of tree.
Reviewer 2 (Olivier Tenaillon): Results - The term allopatric and sympatric should clearly be defined in the present context and not in the method. Being an evolutionary biologist rather than a microbiologist, I do not support the use of these terms that refer to speciation and not to accessibility of foreign DNA, but as the author have previously used them… It should nevertheless be explicitly defined, with the present definition it is not clear who has access to other microbes.
Author’s response: Indeed the terms allopatry and sympatry needed to be clearly defined, we added the complete definition in methodology part and we also put a short definition in the results part.
Reviewer 2 (Olivier Tenaillon): The mobilome analysis are not described enough to be reproducible.
Author’s response: We improved the paragraphs dedicated to the mobilome analysis by adding explanations about our strategy (elements selected and calculation of proportions). We completed these explanations by a table of the databases used to retrieve the mobilome elements, available in the Additional file 4.
Reviewer 2 (Olivier Tenaillon): “line 206 “Genes absent to all orthologous groups are either acquired by HT or generated de novo”: this sentence is not clear. Sometimes HT is used sometimes HGT.
Author’s response: The term HT indicates an horizontal transfer of sequence, this sequence could contain one or more genes (average of genes implicated by transfer is 1.5). The term HGT refers to a transfer of one gene. As the use of these two terms complicated our manuscript, we decided to use only the term of HGT, and we adapted our data in consequence.
Reviewer 2 (Olivier Tenaillon): line 206 NR database: the author should say what it is and use either nr or NR.
Author’s response: Indeed we used the nr database, we corrected this mistakes in manuscript.
Reviewer 2 (Olivier Tenaillon): Figure - species names should be written in full and in italic with no underscore. Line 70–75 Statistics should be given to compare “sympatric” groups.
Author’s response: The figures were corrected and we added the statistics for sympatric bacteria comparison in the paragraph.
We thank Olivier Chabrol for assistance in computer programming during the elaboration of HGT detection strategy, Manuela Royer Carenzi for the assistance in statistical analyzes and Pauline Fournerie for her implication in electronic microscopy. We also thank the Xegen company for the assistance in HGT detection by using of Phylopattern software. This work was supported by the Assistance Publique - Hopitaux de Marseille (Marseille Public University Hospital System).
This work was supported by the Assistance Publique - Hopitaux de Marseille (Marseille Public University Hospital System). The funders had no role in study design, data collection and interpretation, or the decision to submit the work.
Availability of data and materials
PS carried out the design of the study, the strategy elaboration and the collect of data, performed the statistical analysis and interpretation of results and drafted the manuscript. PP participated in strategy elaboration, data interpretation and revised the manuscript. RD conceived the study, participated in its design and coordination and revised the manuscript. BJP coordinated and realized the microscopic observations and helped in data interpretation. PI coordinated and participated in microscopic observations, helped in data interpretation and revised the manuscript. All authors read and approved the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Gupta RS, Bhandari V, Naushad HS. Molecular signatures for the PVC clade (planctomycetes, verrucomicrobia, chlamydiae, and lentisphaerae) of bacteria provide insights into their evolutionary relationships. Front Microbiol. 2012;3:327.PubMedPubMed CentralGoogle Scholar
- Wagner M, Horn M. The planctomycetes, verrucomicrobia, chlamydiae and sister phyla comprise a superphylum with biotechnological and medical relevance. Curr Opin Biotechnol. 2006;17(3):241–9.View ArticlePubMedGoogle Scholar
- Cho JC, Vergin KL, Morris RM, Giovannoni SJ. Lentisphaera araneosa gen. nov., sp. nov, a transparent exopolymer producing marine bacterium, and the description of a novel bacterial phylum, lentisphaerae. Environ Microbiol. 2004;6(6):611–21.View ArticlePubMedGoogle Scholar
- Siegl A, Kamke J, Hochmuth T, Piel J, Richter M, Liang C, … Hentschel U. Single-cell genomics reveals the lifestyle of Poribacteria, a candidate phylum symbiotically associated with marine sponges. ISME J. 2011;5(1):61–70.Google Scholar
- Kamneva OK, Knight SJ, Liberles DA, Ward NL. Analysis of genome content evolution in PVC bacterial super-phylum: assessment of candidate genes associated with cellular organization and lifestyle. Genome Biol Evol. 2012;4(12):1375–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Yarza P, Richter M, Peplies J, Euzeby J, Amann R, Schleifer KH, … Rosselló-Móra R. The All-Species Living Tree project: A 16S rRNA-based phylogenetic tree of all sequenced type strains. Syst. Appl. Microbiol. 2008;31(4):241–250.Google Scholar
- Griffiths E, Gupta RS. Phylogeny and shared conserved inserts in proteins provide evidence that verrucomicrobia are the closest known free-living relatives of chlamydiae. Microbiology. 2007;153(8):2648–54.View ArticlePubMedGoogle Scholar
- Brochier C, Philippe H. Phylogeny: a non-hyperthermophilic ancestor for bacteria. Nature. 2002;417(6886):244.View ArticlePubMedGoogle Scholar
- Wu D, Hugenholtz P, Mavromatis K, Pukall R, Dalin E, Ivanova NN, … Eisen JA. A phylogeny-driven genomic encyclopaedia of Bacteria and Archaea. Nat. 2009;462(7276):1056–1060.Google Scholar
- Jenkins C, Fuerst JA. Phylogenetic analysis of evolutionary relationships of the planctomycete division of the domain bacteria based on amino acid sequences of elongation factor Tu. J Mol Evol. 2001;52(5):405–18.PubMedGoogle Scholar
- Pilhofer M, Rappl K, Eckl C, Bauer AP, Ludwig W, Schleifer KH, Petroni G. Characterization and evolution of cell division and cell wall synthesis genes in the bacterial phyla verrucomicrobia, lentisphaerae, chlamydiae, and planctomycetes and phylogenetic comparison with rRNA genes. J Bacteriol. 2008;190(9):3192–202.View ArticlePubMedPubMed CentralGoogle Scholar
- Gupta RS, Griffiths E. Critical issues in bacterial phylogeny. Theor Popul Biol. 2002;61(4):423–34.View ArticlePubMedGoogle Scholar
- Jun SR, Sims GE, Wu GA, Kim SH. Whole-proteome phylogeny of prokaryotes by feature frequency profiles: an alignment-free method with optimal feature resolution. Proc Natl Acad Sci. 2010;107(1):133–8.View ArticlePubMedGoogle Scholar
- Ciccarelli FD, Doerks T, Von Mering C, Creevey CJ, Snel B, Bork P. Toward automatic reconstruction of a highly resolved tree of life. Science. 2006;311(5765):1283–7.View ArticlePubMedGoogle Scholar
- Lagkouvardos I, Jehl MA, Rattei T, Horn M. Signature protein of the PVC superphylum. Appl Environ Microbiol. 2014;80(2):440–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Buckley DH, Huangyutitham V, Nelson TA, Rumberger A, Thies JE. Diversity of planctomycetes in soil in relation to soil history and environmental heterogeneity. Appl Environ Microbiol. 2006;72(7):4522–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang J, Jenkins C, Webb RI, Fuerst JA. Isolation of gemmata-like and Isosphaera-like planctomycete bacteria from soil and freshwater. Appl Environ Microbiol. 2002;68(1):417–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Sangwan P, Kovac S, Davis KE, Sait M, Janssen PH. Detection and cultivation of soil verrucomicrobia. Appl Environ Microbiol. 2005;71(12):8402–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Horn M. Chlamydiae as symbionts in eukaryotes. Annu Rev Microbiol. 2008;62:113–31.View ArticlePubMedGoogle Scholar
- Zwart G, Crump BC, Kamst-van Agterveld MP, Hagen F, Han SK. Typical freshwater bacteria: an analysis of available 16S rRNA gene sequences from plankton of lakes and rivers. Aquat Microb Ecol. 2002;28(2):141–55.View ArticleGoogle Scholar
- Bergmann GT, Bates ST, Eilers KG, Lauber CL, Caporaso JG, Walters WA, … Fierer N. The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol. Biochem. 2011;43(7):1450–1455.Google Scholar
- Ward N, Staley JT, Fuerst JA, Giovannoni S, Schlesner H, Stackebrandt E. The order planctomycetales, including the genera planctomyces, pirellula, gemmata and Isosphaera and the candidatus genera brocadia, kuenenia and scalindua. New York: Springer; 2006. p. 757–93.Google Scholar
- Shu Q, Jiao N. Different planctomycetes diversity patterns in latitudinal surface seawater of the open sea and in sediment. J Microbiol. 2008;46(2):154–9.View ArticlePubMedGoogle Scholar
- Strous M, Pelletier E, Mangenot S, Rattei T, Lehner A, Taylor MW, … Barbe V. Deciphering the evolution and metabolism of an anammox bacterium from a community genome. Nature. 2006;440(7085):790–794.Google Scholar
- Pearson A, Budin M, Brocks JJ. Phylogenetic and biochemical evidence for sterol synthesis in the bacterium gemmata obscuriglobus. Proc Natl Acad Sci. 2003;100(26):15352–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Bébéar C, De Barbeyrac B. Genital chlamydia trachomatis infections. Clin Microbiol Infect. 2009;15(1):4–10.View ArticlePubMedGoogle Scholar
- Corsaro D, Greub G. Pathogenic potential of novel chlamydiae and diagnostic approaches to infections due to these obligate intracellular bacteria. Clin Microbiol Rev. 2006;19(2):283–97.View ArticlePubMedPubMed CentralGoogle Scholar
- Lee KC, Webb RI, Janssen PH, Sangwan P, Romeo T, Staley JT, Fuerst JA. Phylum verrucomicrobia representatives share a compartmentalized cell plan with members of bacterial phylum planctomycetes. BMC Microbiol. 2009;9(1):5.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindsay MR, Webb RI, Strous M, Jetten MS, Butler MK, Forde RJ, Fuerst JA. Cell compartmentalisation in planctomycetes: novel types of structural organisation for the bacterial cell. Arch Microbiol. 2001;175(6):413–29.View ArticlePubMedGoogle Scholar
- Chopra I, Storey C, Falla TJ, Pearce JH. Antibiotics, peptidoglycan synthesis and genomics: the chlamydial anomaly revisited. Microbiology. 1998;144:2673–8.View ArticlePubMedGoogle Scholar
- McCoy AJ, Maurelli AT. Building the invisible wall: updating the chlamydial peptidoglycan anomaly. Trends Microbiol. 2006;14:70–7.View ArticlePubMedGoogle Scholar
- Lage OM, Bondoso J, Lobo-da-Cunha A. Insights into the Ultrastructural morphology of novel planctomycetes. Antonie Van Leeuwenhoek. 2013;104(4):467–76.View ArticlePubMedGoogle Scholar
- Lindsay MR, Webb RI, Fuerst JA. Pirellulosomes: a new type of membrane-bounded cell compartment in planctomycete bacteria of the genus pirellula. Microbiology. 1997;143(3):739–48.View ArticleGoogle Scholar
- Fuerst JA, Sagulenko E. Beyond the bacterium: planctomycetes challenge our concepts of microbial structure and function. Nat Rev Microbiol. 2011;9:403–13.View ArticlePubMedGoogle Scholar
- Fieseler L, Horn M, Wagner M, Hentschel U. Discovery of the novel candidate phylum “poribacteria” in marine sponges. Appl Environ Microbiol. 2004;70(6):3724–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Fuerst JA, Webb RI. Membrane-bounded nucleoid in the eubacterium gemmatata obscuriglobus. Proc Natl Acad Sci. 1991;88(18):8184–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Santarella-Mellwig R, Pruggnaller S, Roos N, Mattaj IW, Devos DP. Three-dimensional reconstruction of bacteria with a complex endomembrane system. PLoS Biol. 2013;11(5):e1001565.View ArticlePubMedPubMed CentralGoogle Scholar
- Santarella-Mellwig R, Franke J, Jaedicke A, Gorjanacz M, Bauer U, Budd A, Devos DP. The compartmentalized bacteria of the planctomycetes-verrucomicrobia-chlamydiae superphylum have membrane coat-like proteins. PLoS Biol. 2010;8(1):e1000281.View ArticlePubMedPubMed CentralGoogle Scholar
- Van de Graaf AA, Mulder A, de Bruijn P, Jetten MS, Robertson LA, Kuenen JG. Anaerobic oxidation of ammonium is a biologically mediated process. Appl Environ Microbiol. 1995;61(4):1246–51.PubMedPubMed CentralGoogle Scholar
- Strous M, Fuerst JA, Kramer EH, Logemann S, Muyzer G, van de Pas-Schoonen KT, … Jetten MS. Missing lithotroph identified as new planctomycete. Nat. 1999;400(6743):446–449.Google Scholar
- Fuerst JA. Intracellular compartmentation in planctomycetes. Annu Rev Microbiol. 2005;59:299–328.View ArticlePubMedGoogle Scholar
- Lonhienne TG, Sagulenko E, Webb RI, Lee KC, Franke J, Devos DP, … Fuerst JA. Endocytosis-like protein uptake in the bacterium Gemmata obscuriglobus. Proc. Natl. Acad. Sci. 2010;107(29):12883–12888.Google Scholar
- Devos DP, Reynaud EG. Evolution. Intermediate steps. Science. 2010;330:1187–8.View ArticlePubMedGoogle Scholar
- McInerney JO, Martin WF, Koonin EV, Allen JF, Galperin MY, Lane N, … Embley TM. Planctomycetes and eukaryotes: a case of analogy not homology. Bioessays. 2011;33(11):810–817.Google Scholar
- Fuchsman CA, Rocap G. Whole-genome reciprocal BLAST analysis reveals that planctomycetes do not share an unusually large number of genes with eukarya and archaea. Appl Environ Microbiol. 2006;72(10):6841–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Budd A, Devos DP. Evaluating the evolutionary origins of unexpected character distributions within the bacterial planctomycetes-verrucomicrobia-chlamydiae superphylum. Front Microbiol. 2012;3:401.View ArticlePubMedPubMed CentralGoogle Scholar
- Forterre P. A new fusion hypothesis for the origin of eukarya: better than previous ones, but probably also wrong. Res Microbiol. 2011;162(1):77–91.View ArticlePubMedGoogle Scholar
- Fuerst JA, Sagulenko E. Keys to eukaryality: planctomycetes and ancestral evolution of cellular complexity. Front Microbiol. 2012;3(167):b94.Google Scholar
- Fuerst JA. The PVC superphylum: exceptions to the bacterial definition? Antonie Van Leeuwenhoek. 2013;104(4):451–66.View ArticlePubMedGoogle Scholar
- Devos DP. PVC bacteria: variation of, but not exception to, the gram-negative cell plan. Trends Microbiol. 2014;22(1):14–20.View ArticlePubMedGoogle Scholar
- Pilhofer M, Aistleitner K, Ladinsky MS, König L, Horn M, Jensen GJ. Architecture and host interface of environmental chlamydiae revealed by electron cryotomography. Environ Microbiol. 2014;16(2):417–29.View ArticlePubMedGoogle Scholar
- Hirt RP, Alsmark C, Embley TM. Lateral gene transfers and the origins of the eukaryote proteome: a view from microbial parasites. Curr Opin Microbiol. 2015;23:155–62.View ArticlePubMedPubMed CentralGoogle Scholar
- Merhej V, Royer-Carenzi M, Pontarotti P, Raoult D. Massive comparative genomic analysis reveals convergent evolution of specialized bacteria. Biol Direct. 2009;4(1):13.View ArticlePubMedPubMed CentralGoogle Scholar
- Wolf YI, Koonin EV. Genome reduction as the dominant mode of evolution. Bioessays. 2013;35(9):829–37.View ArticlePubMedPubMed CentralGoogle Scholar
- Molmeret M, Horn M, Wagner M, Santic M, Kwaik YA. Amoebae as training grounds for intracellular bacterial pathogens. Appl Environ Microbiol. 2005;71(1):20–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Moliner C, Fournier PE, Raoult D. Genome analysis of microorganisms living in amoebae reveals a melting pot of evolution. FEMS Microbiol Rev. 2010;34(3):281–94.View ArticlePubMedGoogle Scholar
- Subtil A, Collingro A, Horn M. Tracing the primordial chlamydiae: extinct parasites of plants? Trends Plant Sci. 2014;19(1):36–43.View ArticlePubMedGoogle Scholar
- Facchinelli F, Colleoni C, Ball SG, Weber AP. Chlamydia, cyanobiont, or host: who was on top in the ménage à trois? Trends Plant Sci. 2013;18(12):673–9.View ArticlePubMedGoogle Scholar
- Viswanathan VK. Spirochaetes and their twisted ways. Gut Microbes. 2012;3(5):399–400.View ArticlePubMedGoogle Scholar
- Bellgard MI, Wanchanthuek P, La T, Ryan K, Moolhuijzen P, Albertyn Z. Genome sequence of the pathogenic intestinal spirochete brachyspira hyodysenteriae reveals adaptations to its lifestyle in the porcine large intestine. PLoS One. 2009;4:e4641.View ArticlePubMedPubMed CentralGoogle Scholar
- Plague GR, Dunbar HE, Tran PL, Moran NA. Extensive proliferation of transposable elements in heritable bacterial symbionts. J Bacteriol. 2008;190(2):777–9.View ArticlePubMedGoogle Scholar
- Bordenstein SR, Reznikoff WS. Mobile DNA in obligate intracellular bacteria. Nat Rev Microbiol. 2005;3(9):688–99.View ArticlePubMedGoogle Scholar
- Cascales E, Christie PJ. The versatile bacterial type IV secretion systems. Nat Rev Microbiol. 2003;1(2):137–49.View ArticlePubMedGoogle Scholar
- Christie PJ, Atmakuri K, Krishnamoorthy V, Jakubowski S, Cascales E. Biogenesis, architecture, and function of bacterial type IV secretion systems. Annu Rev Microbiol. 2005;59:451–85.View ArticlePubMedGoogle Scholar
- Georgiades K. Genomics of epidemic pathogens. Clin Microbiol Infect. 2012;18(3):213–7.View ArticlePubMedGoogle Scholar
- Merhej V, Notredame C, Royer-Carenzi M, Pontarotti P, Raoult D. The rhizome of life: the sympatric rickettsia felis paradigm demonstrates the random transfer of DNA sequences. Mol Biol Evol. 2011;28(11):3213.View ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- MESQUITE : http://mesquiteproject.org/ consulted in June and July 2016. Accessed June 2015.
- Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3(2):217–23.View ArticleGoogle Scholar
- Repbase: http://www.girinst.org/repbase/ consulted in march 2015. Accessed March 2015.
- CRISPfinder: http://crispr.i2bc.paris-saclay.fr/ consulted in February 2015. Accessed February 2015.
- Fischer S, Brunk BP, Chen F, Gao X, Harb OS, Iodice JB, Shanmugam D, Roos DS, Stoeckert CJ. Using OrthoMCL to assign proteins to OrthoMCL-DB groups or to cluster proteomes into New ortholog groups. Current Protocols in Bioinformatics. 2011;35(6):12. 1–6.12.19.Google Scholar
- Gouret P, Thompson JD, Pontarotti P. PhyloPattern: regular expressions to identify complex patterns in phylogenetic trees. BMC Bioinf. 2009;10(1):298.View ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.View ArticlePubMedGoogle Scholar