- Research article
- Open Access
The CtrA phosphorelay integrates differentiation and communication in the marine alphaproteobacterium Dinoroseobacter shibae
© Wang et al.; licensee BioMed Central Ltd. 2014
- Received: 26 July 2013
- Accepted: 6 February 2014
- Published: 13 February 2014
Dinoroseobacter shibae, a member of the Roseobacter clade abundant in marine environments, maintains morphological heterogeneity throughout growth, with small cells dividing by binary fission and large cells dividing by budding from one or both cell poles. This morphological heterogeneity is lost if the quorum sensing (QS) system is silenced, concurrent with a decreased expression of the CtrA phosphorelay, a regulatory system conserved in Alphaproteobacteria and the master regulator of the Caulobacter crescentus cell cycle. It consists of the sensor histidine kinase CckA, the phosphotransferase ChpT and the transcriptional regulator CtrA. Here we tested if the QS induced differentiation of D. shibae is mediated by the CtrA phosphorelay.
Mutants for ctrA, chpT and cckA showed almost homogeneous cell morphology and divided by binary fission. For ctrA and chpT, expression in trans on a plasmid caused the fraction of cells containing more than two chromosome equivalents to increase above wild-type level, indicating that gene copy number directly controls chromosome number. Transcriptome analysis revealed that CtrA is a master regulator for flagellar biosynthesis and has a great influence on the transition to stationary phase. Interestingly, the expression of the autoinducer synthase genes luxI 2 and luxI 3 was strongly reduced in all three mutants, resulting in loss of biosynthesis of acylated homoserine-lactones with C14 side-chain, but could be restored by expressing these genes in trans. Several phylogenetic clusters of Alphaproteobacteria revealed a CtrA binding site in the promoters of QS genes, including Roseobacters and Rhizobia.
The CtrA phosphorelay induces differentiation of a marine Roseobacter strain that is strikingly different from that of C. crescentus. Instead of a tightly regulated cell cycle and a switch between two morphotypes, the morphology and cell division of Dinoroseobacter shibae are highly heterogeneous. We discovered for the first time that the CtrA phosphorelay controls the biosynthesis of signaling molecules. Thus cell-cell communication and differentiation are interlinked in this organism. This may be a common strategy, since we found a similar genetic set-up in other species in the ecologically relevant group of Alphaproteobacteria. D. shibae will be a valuable model organism to study bacterial differentiation into pleomorphic cells.
- Mutant Strain
- Quorum Sense
- Asymmetric Cell Division
- Position Weight Matrix
- Binary Fission
Bacteria have in the past been regarded as very simple organisms, dividing by binary fission into clones of identical daughter cells which can perform an unlimited number of cell divisions and thus were considered immortal in the absence of external killing events. This view has since been abandoned and replaced by the concept of asymmetric cell division in most if not all bacterial species, resulting in progeny of different cellular composition, different history, and different fate . Asymmetric cell division has first been recognized in those bacterial species that undergo complex differentiation processes, e.g. formation of filaments, buds, spores, swarmer cells and stalked cells. However, even binary fission is asymmetric in many bacteria, including E. coli which produces morphologically identical daughter cells. Each daughter cell contains an old and a young cell pole, and this orientation is maintained in subsequent divisions, resulting in an increasing probability of death with increasing number of cell divisions [2, 3]. Like in eukaryotic cells, the number of cell divisions that an individual cell is able to perform appears to be limited in prokaryotic cells, too . Asymmetric cell division through growth from one cell pole has frequently been observed in Alphaproteobacteria and recognized as a possibility to create progeny with dissimilar age and cell fates .
The freshwater bacterium Caulobacter crescentus became the model organism to study asymmetric cell division due to its strictly dimorphic lifestyle. A sessile stalked cell gives birth to a flagellated motile cell that cannot divide until it loses its flagellum and develops a stalk again. The complex molecular mechanisms controlling the cell cycle of C. crescentus have been thoroughly analysed [5, 6]. In brief, asymmetry of daughter cells is ensured through a phosphorylation gradient of the essential response regulator CtrA  between the two cell poles . This gradient is maintained by the trans-membrane histidine kinase CckA that is localized at the poles in the dividing cell . Different enzymatic complexes maintain its activities as a kinase at the swarmer pole and as a phosphatase at the stalked pole that are transmitted to CtrA via the diffusible phosphotransferase ChpT . This phosphorelay system is conserved in most Alphaproteobacteria and has been found to control divergent traits in addition to cell division. The underlying functional principles however are only poorly understood outside of the specialized Caulobacter model system .
Polar growth is conserved in the order Rhizobiales as are most of the cell-cycle regulating genes of C. crescentus. CtrA has been demonstrated to be essential in the plant symbiont Sinorhizobium meliloti and might play a role during differentiation from free-living, dividing cells into non-dividing, nitrogen fixing and endosymbiotic bacteroids . There is strong evidence that CtrA is also essential for replication and cell division in the plant pathogen Agrobacterium tumefaciens[15, 16] as well as in Brucella abortus, the causative agent for brucellosis in mammals including humans .
In contrast to the aforementioned examples, CtrA has not been found to be essential in the Rhodospirillales and Rhodobacterales strains examined so far. In both orders only a subset of C. crescentus cell-cycle genes is present . The differentiation of the phototroph Rhodospirillum centeum from flagellated swarmer cells into non-motile aggregating cysts is controlled by a ctrA homologue . Based on evidence from the literature and their finding that CtrA controls flagellar biosynthesis but not cell division in Magnetospirillum magneticum, Greene et al. suggested that regulation of motility is the primordial role of the CtrA phosphorelay and other functions like control of the cell cycle might have been acquired later during evolution . Indeed, flagella and chemotaxis gene expression but not growth and cell division are impaired in a ctrA knock-out mutant of Rhodobacter capsulatus[20, 21]. In this organism however, the CtrA phosphorelay also controls expression of a gene cluster coding for the gene transfer agent (GTA) mediating a virus-like exchange of DNA between organisms . As only a small fraction of the population expresses the GTA gene cluster, regulation by CtrA might be heterogeneous or bistable in this case . Another exciting new finding was that expression of ctrA depends on the quorum sensing (QS) system of R. capsulatus. QS refers to a form of cell-to-cell communication that involves the production, excretion and detection of small diffusible signaling molecules called autoinducers (AI). Thus, the CtrA phosphorelay is integrated into the communication system of this organism.
Dinoroseobacter shibae is a representative of the Roseobacter clade, an ecologically important phylogenetic cluster of marine Rhodobacterales. It was isolated from the dinoflagellate Prorocentrum lima and lives in symbiosis with marine algae . D. shibae relies on acylated homoserine-lactones (AHL) for communication like many Proteobacteria. Genome analysis revealed the presence of three LuxI type AHL synthase genes (termed luxI 1-3 ). luxI 1 and luxI 2 are located on the chromosome downstream of a gene encoding a LuxR-type transcriptional regulator, whereas luxI 3 is on the 86-kb plasmid downstream of a gene encoding an autoinducer but no DNA binding domain . Homologs to the genes ctrA, chpT and cckA have been identified, too .
Recently we found that QS induces differentiation of D. shibae into pleomorphic cells utilizing distinct types of cell division . A luxI 1 knock-out led to a complete loss of AHL biosynthesis, homogeneous cell morphology and faster rate of cell division. Expression of the genes ctrA, chpT and cckA was strongly reduced in this strain. Here we test the hypothesis that the QS induced differentiation of D. shibae is mediated by the CckA-ChpT-CtrA phosphorelay system.
Homologs of cell cycle related genes from Caulobacter crescentus NA1000 in D. shibae DFL-12
Locus Tag D. shibae
Locus Tag C.crescentus
Integral membrane sensor hybrid histidine kinase
Two component transcriptional regulator
ATP-dependent Clp protease, ATP-binding subunit
0E + 00
Cell cycle regulator
DNA methylase N-4/N-6 domain protein
Chromosomal replication initiator protein
Single gene knockouts of ctrA, cckA and chpT reduce the morphological heterogeneity of D. shibae
Growth rate and doubling time of D. shibae wild-type and mutant strains
Growth rate ± s.d. (h-1)
Doubling time ± s.d. (h)
0.16 ± 0.00
4.33 ± 0.07
0.27 ± 0.02
2.61 ± 0.14
0.25 ± 0.01
2.72 ± 0.08
0.25 ± 0.01
2.80 ± 0.09
0.27 ± 0.01
2.57 ± 0.09
Replicative diversification is controlled by ctrA, cckA and chpT
Comparative transcriptome analysis revealed substantial overlaps of gene expression in ctrA, cckA and chpT-deficient strains
To understand the global contribution of the CtrA phosphorelay to gene regulation in D. shibae, we performed microarray analysis to examine gene expression profiles of the three mutant strains compared to the wild-type in both mid-exponential and stationary phase of growth. Differential expression was assumed for genes with an absolute log2 fold change higher or equal to one and a Benjamini-Hochberg-adjusted p-value smaller than 0.05. The log2 fold-changes and p-values for all genes present on the microarray can be found in Additional file 1: Table S1, also including the microarray data for ∆luxI 1 versus wild-type  for comparison. The differential expression of four genes was also confirmed by quantitative reverse transcription PCR (qRT-PCR). This data can be found in supplementary Additional file 2: Figure S1.
In the stationary phase, the number of genes that were differentially expressed increased remarkably (Figure 3A). 526, 417 and 517 genes were differentially expressed in the ∆ctrA, ∆cckA and ∆chpT mutant, respectively. In all three cases around half of the genes were up-regulated compared to the wild-type. 349 genes were differentially expressed in all three mutants. The large overlap in gene expression of all three mutant strains suggests that CckA, ChpT and CtrA indeed act in a common pathway in D. shibae.
Differentially regulated genes were categorized by Clusters of Orthologous Groups (COG) designations, and the numbers of genes differentially regulated both in exponential and stationary phase in each category are shown in Figure 3B. Categories that showed a notable over-representation of down-regulated genes in both growth phases included those related to cell motility, intracellular trafficking and secretion, and transcription and signal transduction. Specific groups of the CtrA regulon will be discussed in later sections. The CtrA binding site was well conserved in D. shibae and could be found in 101 promoter regions (Figure 3C and Additional file 3: Table S2).
The ctrA phosphorelay influences the adaptation of D. shibae to stationary phase
The enormous increase in the number of differentially expressed genes in the stationary phase was mainly due to the higher expression of genes falling in the categories primary metabolite biosynthesis, transport and catabolism as well as genes with a function in transcriptional control and signal transduction (Figure 3B). In particular many genes involved in amino acid biosynthesis and transport showed a stronger expression in all three mutants. Remarkably genes encoding the cytochrome c complex, a cytochrome bd ubiquinol oxidase as well as two other cytochromes were higher expressed indicating that the redox state of mutant cells may be altered and they may adapt to stationary phase in a different way. The activation of various genes predicted to encode heat- and cold-shock proteins as well as chaperones indicates a stress response in the mutant strains (Additional file 1: Table S1). We asked next which regulatory genes might be involved in the different adaptation strategies. Two sigma factors (rpoH 2 and a putative sigma-24 factor) and one sigma/anti sigma factor pair (rpoE and chrR) showed up-regulation in all three mutants exclusively in stationary phase (Figure 3D). These bona fide regulators of the oxidative stress response  therefore represent likely candidates for the observed induction of gene expression in the stationary phase. None of them had a CtrA binding site in their promoter region. The only sigma factor with a respective binding site was rpoH 1 (Figure 3E) which showed reduced expression in all three mutants during both growth phases (Figure 3D). Thus, the differences between wild-type and mutants in the stationary phase might not be a result of direct regulation by the CtrA phosphorelay.
Could phenotypic observations help to interpret these data? We found that the wild-type cells were still heterogeneous in the stationary phase, although cell size was on average lower than in the exponential growth phase. In contrast, mutant strain cells – small rods in the exponential phase – were more cocci-shaped in the stationary phase (Additional file 4: Figure S2A). Flow cytometric analysis confirmed the overall size reduction of all strains (Additional file 4: Figure S2B). Determination of the number of chromosome equivalents per cell revealed a shift of the peaks to 20% lower intensity levels (Additional file 4: Figure S2C). As this shift was not stoichiometric we assume that it was caused not by a reduction of chromosome equivalents but rather by a reduction of the RNA content in the cells. RNA has not been removed prior to staining and is also stained by SybrGreen. Furthermore the fraction of cells with only one chromosome equivalent was substantially higher for all strains. It was highest for the ∆luxI 1 mutant, followed by ∆chpT, ∆cckA, ∆ctrA and the wild-type. The latter was also the only strain for which cells with more than two chromosome equivalents could still be found. Interestingly, some cells of the ∆ctrA and ∆cckA mutant strains formed long chains indicating that they could not divide properly (Additional file 4: Figure S2A). Both, transcriptome and phenotypic data suggest that the QS induced and CtrA phosphorelay mediated pleomorphism substantially alters the adaptation of D. shibae to nutrient limitation.
The ctrA phosphorelay regulates the expression of a gene transfer agent gene cluster
Gene transfer agents (GTAs) are phage-like particles that have first been discovered in Rhodobacter capsulatus. In D. shibae, we have found an approximately 15 kb GTA gene cluster (Dshi_2162-2178) with identical organization like in R. capsulatus. Our microarray results showed a slightly but consistently reduced expression of the complete GTA gene cluster in all three mutant strains. Two genes of this cluster (Dshi_2174, Dshi_2175) were significantly down-regulated in all mutants (Additional file 1: Table S1). In R. capsulatus, CckA and CtrA have been shown to positively regulate the expression of the GTA structural genes . However, only 0.15% of all cells actively expressed this cluster . Assuming a similar frequency in D. shibae, it is likely that down-regulation of the GTA cluster might not result in a strong signal for transcriptome analysis of the whole culture.
CtrA is a master regulator of flagellar biosynthesis
Plasmid gene expression in the ∆ctrA mutant strain
Genes with a significant change in expression in all three mutants were located on the chromosome, with four exceptions, among them the luxI 3 gene. This picture changed when only those genes were taken into consideration that were regulated exclusively in the ∆ctrA strain. Two of the five plasmids of D. shibae, pDshi01 (NC_009955) and pDshi03 (NC_009957), are in large parts syntenic (103 kb shared with 90% identity) and therefore regarded as sister plasmids of common origin . Despite this strong similarity, they differ with respect to their replication operons which are from different compatibility groups . We found 43 genes on pDshi01 being down-regulated only in D. shibae ∆ctrA. In particular, the replication and toxin/antidote system showed a strong loss in the microarray signal (Additional file 5: Figure S3A). Plotting the expression changes of all genes against a theoretical normal distribution showed that the log2 fold changes were biased towards lower expression of pDshi01 (Additional file 5: Figure S3B). The observed bias in the expression changes of this plasmid leaves room for two different explanations. It could mean that CtrA acts exclusively as an activator of pDshi01 gene expression which is highly unlikely given the high degree of identity between pDshi01 and pDshi03. More likely, this bias could also explained by a modulation of the copy number or loss of this plasmid in this mutant. PCR using pDshi01 specific primers and genomic DNA as a template demonstrated that indeed the plasmid has been lost in the ∆ctrA knock-out strain (Additional file 5: Figure S3C). Given the strong overlap in gene expression changes between and consistent phenotypic alterations of all three mutants we assumed that plasmid loss has no great impact on the experimental results regarding the ∆ctrA mutant. However, as plasmid loss can also occur randomly, it should not be ascribed to the loss of CtrA activity.
Regulation of C. crescentus cell-cycle homologs by the CtrA phosphorelay
The CtrA phosphorelay induces autoinducer biosynthesis
Recently we found that knock-out of the autoinducer synthase luxI 1 led to a reduced expression of ctrA, cckA, chpT, divL as well as the QS genes luxR 2 I 2 and the AHL synthase luxI 3 during exponential growth and stationary phase, resulting in a complete loss of autoinducer biosynthesis . Here, we found that the latter QS genes were also down-regulated in the CtrA phosphorelay mutants in both growth phases (Figure 5B). Expression of luxR 2 I 2 was most strongly reduced in the ∆ctrA strain and for all three strains expression was lower in the exponential phase. By contrast, the down-regulation of the AHL synthase luxI 3 and the autoinducer binding domain gene Dshi_4067 which are located on the 86 kb plasmid pDshi04 (NC_009958) was more pronounced in the stationary phase. The gene expression of luxI 1 (Dshi_0312) and its cognate luxR 1 regulator (Dshi_0311) was not affected by deletion of any of the three genes in the CtrA regulon. The same was true for the three orphan luxR-type regulators (Figure 5B) whose function is currently unknown. We found two perfect CtrA binding sites in the promoter of luxR 2 and one conserved binding site in front of Dshi_4067 suggesting that CtrA is tightly integrated into the control of AHL biosynthesis (Figure 5C). Therefore we investigated the AHL production using the biosensor strain P. putida F117 pKR-C12 which responds mainly to long-chain AHLs . Surprisingly, the extracts of all three mutant strains failed to induce fluorescence in the reporter strain (Figure 5D). In contrast, we were able to detect the production of AHLs in the genetic complementation strains of all the knock-out mutants, although with a reduced amount compared to the wild-type.
AHLs in different strains of D. shibae identified by GC-MS
Integration of the CtrA phosphorelay and quorum sensing might be present in other Alphaproteobacteria
Knock-out of ctrA, cckA and chpT resulted in a similar reduction of morphological heterogeneity and cell size as in the QS mutant confirming that they act in a common pathway controlling cell growth. We further demonstrated that all three genes are involved in the control of replication as their loss led to reduction of cells with a larger number of chromosomes in the population. This was not unexpected as cell size and replication are coordinated in bacteria . In contrast to ctrA and chpT genetic complementation of ∆cckA failed to shift the population towards cells with higher chromosome content. As the gene was transcribed and autoinducer synthesis could be restored in trans, it is unlikely that the introduced plasmid was not functional at all. The sensor histidine kinase is a trans-membrane protein. In C. crescentus it is specifically localized at the cell poles and creates a gradient of phosphorylated CtrA . In D. shibae, the six fold over-expression of cckA might lead to a mislocation of the protein and therefore a gradient that might be necessary for differentiation into cells with multiple chromosomes could not be formed.
The phenotype of the ∆chpT mutant – in contrast to ∆ctrA and ∆cckA – is almost identical to that of the ∆luxI 1 mutant. However, comparison of the transcriptome data from ∆chpT and ∆luxI 1  reveals only few genes that are differentially expressed in both strains but not in ∆ctrA and ∆cckA (Additional file 1: Table S1). Among them are genes encoding parts of an ATP synthase (Dshi_3028-3031) that are up-regulated and T4SS (Dshi_3975/3976) that are down-regulated in both mutant strains as well as few hypothetical genes. However, these data offer no explanation for the physiological differences between the aforementioned strains. In C. crescentus, CtrA is not the only target of phosphorylation by ChpT  and the differences between ∆chpT and ∆luxI 1 (in the latter the whole phosphorelay is down-regulated) on one side and ∆ctrA and ∆cckA on the other might be a result of altered phosphorylation states rather than gene expression.
The remarkable variety of cell shapes and sizes as well as chromosome copy numbers of D. shibae naturally raises three questions: How is this heterogeneity induced, are there other traits that are heterogeneously regulated as well and finally, how does the population benefit from such heterogeneity? At the moment we can only speculate on the mechanisms underlying the development of heterogeneous cells. It might be stochastic gene expression or a bistable regulatory feedback loop (reviewed in  and ). Single cell techniques have to be developed to answer this question. We showed that the control of flagellar biosynthesis by CtrA is also conserved in D. shibae. Recently it has been demonstrated that under the same growth conditions that have been used for experiments in this manuscript, only a fraction of D. shibae cells are motile. Thus flagellar biosynthesis seems to be heterogeneously distributed, too . Furthermore it seems reasonable that the GTA cluster in D. shibae is heterogeneously activated like in R. capsulatus.
Differentiation clearly represents a burden for population growth of D. shibae, as the doubling time is longer and the maximum number of cells is reduced . In this context one might ask if it also represents a burden during adaptation to stationary phase. Our observations suggest that D. shibae reduces its cell size in response to nutrient depletion like marine Vibrios . This cell-size reduction is of course easier to achieve if the cells are already small on entry to stationary phase, and should reduce the need for costly cellular reprogramming. However, our data also suggest that the fast dividing mutant cells might be subjected to higher oxidative stress than the slower growing wild-type population. Thus, fast growth in the exponential phase could also represent a burden during adaption to stationary phase. The role of the CtrA phosphorelay for adaptation to nutrient depletion clearly needs more in depth analysis.
Our recent paper on the role of the LuxI1 autoinducer synthase  together with the present work suggest a strong integration of QS and the CtrA phosphorelay. Our hypothesis was that LuxR1 and LuxI1 are at the top of a hierarchical regulatory system that activates expression of CtrA phosphorelay genes, the luxR 2 I 2 operon and the plasmid encoded luxI 3 . The present findings suggest that CtrA actually activates luxR 2 I 2 and luxI 3 expression; consequently the phosphorelay is not only regulated by but actually integrated into the QS system. Thus, the CtrA mediated differentiation leads to enhanced autoinducer production. This has not been reported for any other organism so far, but our promoter-analysis in Alphaproteobacteria predicts that this kind of integration of these two regulatory systems might indeed be widespread. The question how the cell or the whole population benefits from this regulatory constellation will be addressed in future research.
It is now recognized that QS is not only employed to coordinate the behavior of a bacterial population but often induces heterogeneous responses towards signaling molecules. This has for example also been demonstrated for the marine pathogen Vibrio harveyi[44, 45] and the marine symbiont Vibrio fischeri. D. shibae is an outstanding organism as its heterogeneity is manifested as a pronounced morphological pleomorphism. The fact that loss of morphological heterogeneity results in faster growth will be beneficial for screening libraries of mutants in order to identify genes that are involved in this differentiation process. Studying this organism in more detail will also enhance our understanding of the evolution of CtrA phosphorelay mediated control of cellular polarity, which in the case of C. crescentus leads to a tightly controlled cell cycle and in the case of D. shibae leads to high heterogeneity of replication and cell division.
Bacterial strains and growth conditions
The strains and plasmids used for this study are listed in Additional file 7: Table S3. D. shibae DFL-12 strains were grown at 30°C and 160 rpm in half-concentrated Marine Broth or on the same medium solidified with 1.5% agar. Alternatively, cells were grown in a chemically defined sea water medium (SWM) supplemented with 5 mM succinate, prepared as described previously . When appropriate, 150 μg/ml of gentamicin or 500 μg/ml of kanamycin was added. Escherichia coli strains were grown at 37°C and 180 rpm in Luria-Bertani (LB) broth or on LB agar supplemented with ampicillin (100 μg/ml), gentamicin (20 μg/ml), kanamycin (50 μg/ml), and/or aminolevulinic acid (50 μg/ml) if necessary.
Construction of D. shibae deletion mutants
For deletion of ctrA (Dshi_1508), cckA (Dshi_1644), and chpT (Dshi_1470), each full-length ORF with flanking regions (approximately 800 bp on either side) was amplified by PCR using D. shibae DFL12 genomic DNA and appropriate primers (Additional file 8: Table S4). These fragments were cloned into the multiple cloning site of pEX18Ap, a suicide plasmid carrying an ampicillin resistance cassette. From the resulting plasmids, the target ORFs were replaced with a gentamicin resistance cassette from pBBR1MCS-5  using the CloneEZ PCR Cloning Kit (GenScript, Piscataway, USA) according to the manufacturer’s instructions. The resulting constructs were transferred to D. shibae DFL12 using conjugation with E. coli ST18 as a helper strain  and integrated into the wild-type strain by double-homologous recombination. The transconjugants were selected on half-concentrated Marine agar plates containing gentamicin. All mutant strains constructed in this study were confirmed by PCR and sequencing.
Complementation of mutants
Complementation of the deleted genes under the control of their native upstream sequences was performed with the broad-host-range vector pBBR1MCS-2 . The complementing fragments were amplified by PCR using appropriate primers (Table 2) and cloned into pBBR1MCS-2, generating plasmids pHW1, pHW2 and pHW3. These plasmids were transferred into the respective mutant strain via conjugation using E. coli ST18 . All complementations were also confirmed by PCR and sequencing.
Cell morphology was observed with an Olympus BX60 microscope with a Plan 100×/1.25 oil immersion objective. All images were acquired using the Olympus OM-4 Ti camera and the cell B image acquisition software (Olympus, Japan).
For growth curves, cell material from half-concentrated MB plates was inoculated into 20 ml of SWM + 5 mM succinate and incubated overnight at 30°C and 160 rpm to obtain a pre-culture. Optical density at 600 nm (OD600) was measured with spectrophotometer Ultrospec 3100 pro (Biochrom Ltd, Cambridge CB4 0FJ England). Cultures were diluted to 0.01 OD600 and then 200 μl were placed into each well of a Honeycomb 2 plate (100 wells each plate, Oy Growth Curves Ab Ltd, Helsinki, Finland). OD600 was monitored every 30 min for 36 hours in the automated microbiology growth analysis system Bioscreen C (Oy Growth Curves Ab Ltd, Helsinki, Finland). Growth curves were plotted in ln scale using the R program.
Cultures for the flow cytometry and microarray analysis were inoculated into 100 ml medium in 300 ml Erlenmeyer flasks with an initial OD600 of 0.01 from overnight cultures and incubated for 36 hours. Samples of the wild-type and all mutant strains were taken in the middle of the exponential phase at an OD600 of approximately 0.4 and in the stationary phase 6 h after OD600 had reached a maximum.
Samples for flow cytometry analyses were fixed at a final concentration of 2% glutaraldehyde and stained with SYBR Green I (Molecular Probes, Leiden, The Netherlands) by diluting the stock reagent 1: 10,000 into the samples. For each experiment, the DNA content was measured in a population of 50, 000 cells using the BD FACS Canto analyzer. The data were analysed using the ‘flowCore’ package  of the R BioConductor project.
Cells were collected by centrifugation at 12,000 × g for 1 min at 4°C, covered with 1 ml Trizol reagent (Ambion, Germany), immediately frozen in liquid N2 and stored at -70°C until processing. For RNA extraction, cells were homogenized with ~ 0.3 g of glass beads in the FastPrep-24 instrument (MP Biomedicals, California, USA) at 6.0 m/s for 3 min and then incubated for 5 min at room temperature. Samples were centrifuged at 12,000 × g for 10 min at 4°C and the supernatants were transferred to fresh tubes, followed by the addition of 100 μl of 1-bromo-3-chloropropane (BCP, Sigma, Germany) and incubation for 10 min at room temperature. Samples were centrifuged at 12,000 × g for 10 min at 4°C, after which the aqueous phase was transferred to new tubes and mixed with 500 μl of absolute ethanol. Extracts were applied to RNeasy spin columns (RNeasy mini kit, Qiagen, Hilden, Germany) and processed according to the manufacturer’s instructions. In addition, samples were treated with DNAse I (Qiagen, Hilden, Germany). Removal of genomic DNA was verified via PCR. The concentration of the RNA was quantified using a NanoDrop spectrophotometer (Peqlab, Erlangen, Germany) and the RNA integrity was assessed using Bioanalyzer 2100 (Agilent, Santa clara, USA).
2 μg of total RNA was labeled with Cy3 or Cy5 using the ULS-system (Kreatech, Amsterdam, The Netherlands) according to the manufacturer's protocol. 600 ng of labeled RNA was fragmented and hybridized to the microarray according to Agilent's two-color microarray protocol. Two biological replicates were applied in this experiment.
Microarray data analysis
Microarray slides were scanned using the Agilent DNA Microarray Scanner. Median foreground and background signals of the Cy3 and Cy5 channel were loaded into the R environment and processed using the LIMMA package . Spots were down-weighted if 2 or more quality flags were set by the scanner software. Background signals were subtracted using the normexp method , Cy3 and Cy5 signals were Loess normalized and finally quantile normalisation was performed on all microarrays from one dataset. Signals from replicate probes for single genes were averaged. A linear model was fitted for each comparison of interest. The obtained p-values were adjusted for false discovery rate (fdr) using the method by Benjamini and Hochberg . For further analysis, only those genes with a false-discovery-rate-adjusted p-value < 0.05 and an absolute log2 fold change > 1 under at least one condition were taken into account. Raw and processed microarray data have been deposited at the gene expression omnibus database under the accession number GSE47451.
Validation of microarray data by qRT-PCR
Reverse transcription of the same RNA samples as those used for microarray analysis was performed with the QuantiTect Reverse Transcription Kit (Qiagen, Hilden, Germany), according to the manufacturer’s protocol. qPCR analysis was performed using the LightCycler 480 (Roche, Mannheim, Germany) with the QuantiTect SYBR Green PCR Kit (Qiagen, Hilden, Germany). Primer sequences are listed in Additional file 8: Table S4. Relative expression was calculated by the ∆∆Ct method and normalized to the dapB gene coding for dihydrodipicolinate reductase. dapB was selected as reference gene based on the microarray result, as it shows similar expression in all strains. The experiment was performed in duplicate.
Identification of transcription factor binding sites (TFBS)
For the identification of binding sites for the transcription factor CtrA, promoter-regions ranging from 400 bp upstream to 50 bp downstream of the translation start were searched for matches with a position weight matrix (PWM) obtained from a comparative phylogenetic analysis of cell cycle regulation in Alphaproteobacteria. A TFBS was considered as identified if its score reached at least 85% of the maximum weight of the PWM.
Extraction and detection of AHLs
The production of AHLs was detected using the biosensor strain P. putida F117 pKR-C12 as previously described  with the following modifications. All D. shibae strains were grown in 100 ml of SWM with 5 mM succinate and 4% of adsorber resins (Amberlite XAD-16, Rohm & Haas) at 30°C and 160 rpm for 48 hours. Adsorber resins were harvested after 2 days by filtration and transferred to a separating funnel together with 50 mL distilled water, and then those resins were extracted with 25 ml of dichloromethane (DCM). The mixture was shaken vigorously for 1 min and the phases were allowed to separate. The DCM fraction was removed and another 25 ml of DCM was added. The whole extraction process was repeated three times. The combined DCM fractions were concentrated to 2 ml using a rotary evaporator (Heidolph VV2001, Schwabach, Germany) and stored at -20°C until the bioassay was conducted. For bioassays, 20 μl of the concentrated extract was applied to each well of a 96 well polypropylene micro-titer plate (PlateOne, Starlab, Hamburg, Germany). After air-drying, the micro-titer plate was overlaid with 100 μl of medium and 100 μl of the sensor strain and then incubated at 30°C with gentle agitation for 24 hours. Green fluorescence and optical density at 620 nm were determined in a Victor3 1420 Multilabel counter (PerkinElmer, Waltham, USA). The presence of AHLs in the strains was calculated by dividing the specific fluorescence (gfp535/OD620) of extracts by that of the negative control (DCM).
A volume of 1 μL of the extract was injected into an Agilent GC 7890A gas chromatograph connected to a 5975C mass-selective detector (Agilent) fitted with a HP-5 MS fused silica capillary column (30 m × 0.25 mm i.d., 0.22 μm film; Hewlett-Packard, Wilmington, USA). Conditions were as follows: inlet pressure: 67.5 kPa, He 24.2 ml min-1; injection volume: 1 μl; injector: 250°C; transfer line: 300°C; electron energy: 70 eV. The GC was programmed as follows: 50°C (5 min isothermic), increasing at 5°C min-1 to 320°C, and operated in splitless mode; carrier gas (He): 1.2 ml min-1. Retention indicies I were determined from a homologous series of n-alkanes (C8-C32).
Distribution of CtrA-binding sites in QS gene promoters in Alphaproteobacteria
Finished genomes of Alphaproteobacteria with at least one gene with an autoinducer synthesis domain (pfam00765) and a ctrA homolog present were selected from IMG (http://img.jgi.doe.gov/). Promoter regions ranging from 400 bp up to 50 bp downstream of the start codon were downloaded for all genes with an autoinducer synthesis (pfam00765) or an autoinducer binding (pfam03475) domain (Additional file 9: Dataset S1). These regions were searched for putative ctrA binding sites the same way as described for the D. shibae genes.
CtrA phylogeny in Alphaproteobacteria
The evolutionary history of CtrA protein sequences found in Additional file 10: Dataset S2 was inferred using the Neighbor-Joining method . The optimal tree with the sum of branch length = 3.83257521 is shown. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Poisson correction method  and are in the units of the number of amino acid substitutions per site. The analysis involved 93 amino acid sequences. All positions containing gaps and missing data were eliminated. There were a total of 179 positions in the final dataset. Evolutionary analyses were conducted in MEGA5 .
This work was supported by the DFG in Transregio-SFB TR51. We thank our reviewers whose comments helped to improve this manuscript.
- Ksiazek K: Bacterial aging: from mechanistic basis to evolutionary perspective. Cell Mol Life Sci. 2010, 67: 3131-3137. 10.1007/s00018-010-0417-4.PubMedView ArticleGoogle Scholar
- Stewart EJ, Madden R, Paul G, Taddei F: Aging and death in an organism that reproduces by morphologically symmetric division. PLoS Biol. 2005, 3: e45-10.1371/journal.pbio.0030045.PubMed CentralPubMedView ArticleGoogle Scholar
- Wang P, Robert L, Pelletier J, Dang WL, Taddei F, Wright A, Jun S: Robust growth of Escherichia coli. Curr Biol. 2010, 20: 1099-1103. 10.1016/j.cub.2010.04.045.PubMed CentralPubMedView ArticleGoogle Scholar
- Hallez R, Bellefontaine AF, Letesson JJ, De BX: Morphological and functional asymmetry in alpha-proteobacteria. Trends Microbiol. 2004, 12: 361-365. 10.1016/j.tim.2004.06.002.PubMedView ArticleGoogle Scholar
- Tsokos CG, Laub MT: Polarity and cell fate asymmetry in Caulobacter crescentus. Curr Opin Microbiol. 2012, 15: 744-750. 10.1016/j.mib.2012.10.011.PubMed CentralPubMedView ArticleGoogle Scholar
- McAdams HH, Shapiro L: A bacterial cell-cycle regulatory network operating in time and space. Science. 2003, 301: 1874-1877. 10.1126/science.1087694.PubMedView ArticleGoogle Scholar
- Quon KC, Marczynski GT, Shapiro L: Cell cycle control by an essential bacterial two-component signal transduction protein. Cell. 1996, 84: 83-93. 10.1016/S0092-8674(00)80995-2.PubMedView ArticleGoogle Scholar
- Chen YE, Tropini C, Jonas K, Tsokos CG, Huang KC, Laub MT: Spatial gradient of protein phosphorylation underlies replicative asymmetry in a bacterium. Proc Natl Acad Sci U S A. 2011, 108: 1052-1057. 10.1073/pnas.1015397108.PubMed CentralPubMedView ArticleGoogle Scholar
- Jacobs C, Domian IJ, Maddock JR, Shapiro L: Cell cycle-dependent polar localization of an essential bacterial histidine kinase that controls DNA replication and cell division. Cell. 1999, 97: 111-120. 10.1016/S0092-8674(00)80719-9.PubMedView ArticleGoogle Scholar
- Biondi EG, Reisinger SJ, Skerker JM, Arif M, Perchuk BS, Ryan KR, Laub MT: Regulation of the bacterial cell cycle by an integrated genetic circuit. Nature. 2006, 444: 899-904. 10.1038/nature05321.PubMedView ArticleGoogle Scholar
- Brilli M, Fondi M, Fani R, Mengoni A, Ferri L, Bazzicalupo M, Biondi EG: The diversity and evolution of cell cycle regulation in alpha-proteobacteria: a comparative genomic analysis. BMC Syst Biol. 2010, 4: 52-10.1186/1752-0509-4-52.PubMed CentralPubMedView ArticleGoogle Scholar
- Brown PJ, de Pedro MA, Kysela DT, Van der HC, Kim J, De BX, Fuqua C, Brun YV: Polar growth in the Alphaproteobacterial order Rhizobiales. Proc Natl Acad Sci U S A. 2012, 109: 1697-1701. 10.1073/pnas.1114476109.PubMed CentralPubMedView ArticleGoogle Scholar
- Barnett MJ, Hung DY, Reisenauer A, Shapiro L, Long SR: A homolog of the CtrA cell cycle regulator is present and essential in Sinorhizobium meliloti. J Bacteriol. 2001, 183: 3204-3210. 10.1128/JB.183.10.3204-3210.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Becker A, Berges H, Krol E, Bruand C, Ruberg S, Capela D, Lauber E, Meilhoc E, Ampe F, de Bruijn FJ, Fourment J, Francez-Charlot A, Kahn D, Kuster H, Liebe C, Puhler A, Weidner S, Batut J: Global changes in gene expression in Sinorhizobium meliloti 1021 under microoxic and symbiotic conditions. Mol Plant Microbe Interact. 2004, 17: 292-303. 10.1094/MPMI.2004.17.3.292.PubMedView ArticleGoogle Scholar
- Kahng LS, Shapiro L: The CcrM DNA methyltransferase of Agrobacterium tumefaciens is essential, and its activity is cell cycle regulated. J Bacteriol. 2001, 183: 3065-3075. 10.1128/JB.183.10.3065-3075.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim J, Heindl JE, Fuqua C: Coordination of division and development influences complex multicellular behavior in Agrobacterium tumefaciens. PLoS One. 2013, 8: e56682-10.1371/journal.pone.0056682.PubMed CentralPubMedView ArticleGoogle Scholar
- Bellefontaine AF, Pierreux CE, Mertens P, Vandenhaute J, Letesson JJ, De BX: Plasticity of a transcriptional regulation network among alpha-proteobacteria is supported by the identification of CtrA targets in Brucella abortus. Mol Microbiol. 2002, 43: 945-960. 10.1046/j.1365-2958.2002.02777.x.PubMedView ArticleGoogle Scholar
- Bird TH, MacKrell A: A CtrA homolog affects swarming motility and encystment in Rhodospirillum centenum. Arch Microbiol. 2011, 193: 451-459. 10.1007/s00203-011-0676-y.PubMedView ArticleGoogle Scholar
- Greene SE, Brilli M, Biondi EG, Komeili A: Analysis of the CtrA pathway in Magnetospirillum reveals an ancestral role in motility in alphaproteobacteria. J Bacteriol. 2012, 194: 2973-2986. 10.1128/JB.00170-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Mercer RG, Callister SJ, Lipton MS, Pasa-Tolic L, Strnad H, Paces V, Beatty JT, Lang AS: Loss of the response regulator CtrA causes pleiotropic effects on gene expression but does not affect growth phase regulation in Rhodobacter capsulatus. J Bacteriol. 2010, 192: 2701-2710. 10.1128/JB.00160-10.PubMed CentralPubMedView ArticleGoogle Scholar
- Lang AS, Beatty JT: A bacterial signal transduction system controls genetic exchange and motility. J Bacteriol. 2002, 184: 913-918. 10.1128/jb.184.4.913-918.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Fogg PC, Westbye AB, Beatty JT: One for all or all for one: heterogeneous expression and host cell lysis are key to gene transfer agent activity in Rhodobacter capsulatus. PLoS One. 2012, 7: e43772-10.1371/journal.pone.0043772.PubMed CentralPubMedView ArticleGoogle Scholar
- Leung MM, Brimacombe CA, Beatty JT: Transcriptional regulation of the Rhodobacter capsulatus response regulator CtrA. Microbiology. 2013, 159: 96-106. 10.1099/mic.0.062349-0.PubMed CentralPubMedView ArticleGoogle Scholar
- Wagner-Döbler I, Thiel V, Eberl L, Allgaier M, Bodor A, Meyer S, Ebner S, Hennig A, Pukall R, Schulz S: Discovery of complex mixtures of novel long-chain quorum sensing signals in free-living and host-associated marine alphaproteobacteria. Chembiochem. 2005, 6: 2195-2206. 10.1002/cbic.200500189.PubMedView ArticleGoogle Scholar
- Biebl H, Allgaier M, Tindall BJ, Koblizek M, Lünsdorf H, Pukall R, Wagner-Döbler I: Dinoroseobacter shibae gen. nov., sp. nov., a new aerobic phototrophic bacterium isolated from dinoflagellates. Int J Syst Evol Microbiol. 2005, 55: 1089-1096. 10.1099/ijs.0.63511-0.PubMedView ArticleGoogle Scholar
- Wagner-Döbler I, Ballhausen B, Berger M, Brinkhoff T, Buchholz I, Bunk B, Cypionka H, Daniel R, Drepper T, Gerdts G, Hahnke S, Han C, Jahn D, Kalhöfer D, Kiss H, Klenk HP, Kyrpides N, Liebl W, Liesegang H, Meincke L, Pati A, Petersen J, Piekarski T, Pommerenke C, Pradella S, Pukall R, Rabus R, Stackebrandt E, Thole S, Thompson L, et al: The complete genome sequence of the algal symbiont Dinoroseobacter shibae: a hitchhiker's guide to life in the sea. ISME J. 2010, 4: 61-77. 10.1038/ismej.2009.94.PubMedView ArticleGoogle Scholar
- Patzelt D, Wang H, Buchholz I, Rohde M, Grobe L, Pradella S, Neumann A, Schulz S, Heyber S, Munch K, Munch R, Jahn D, Wagner-Dobler I, Tomasch J: You are what you talk: quorum sensing induces individual morphologies and cell division modes in Dinoroseobacter shibae. ISME J. 2013, 7: 2274-2286. 10.1038/ismej.2013.107.PubMed CentralPubMedView ArticleGoogle Scholar
- Tzur A, Moore JK, Jorgensen P, Shapiro HM, Kirschner MW: Optimizing optical flow cytometry for cell volume-based sorting and analysis. PLoS One. 2011, 6: e16053-10.1371/journal.pone.0016053.PubMed CentralPubMedView ArticleGoogle Scholar
- Tomasch J, Gohl R, Bunk B, Diez MS, Wagner-Dobler I: Transcriptional response of the photoheterotrophic marine bacterium Dinoroseobacter shibae to changing light regimes. ISME J. 2011, 5: 1957-1968. 10.1038/ismej.2011.68.PubMed CentralPubMedView ArticleGoogle Scholar
- Marrs B: Genetic recombination in Rhodopseudomonas capsulata. Proc Natl Acad Sci U S A. 1974, 71: 971-973. 10.1073/pnas.71.3.971.PubMed CentralPubMedView ArticleGoogle Scholar
- Lang AS, Beatty JT: Genetic analysis of a bacterial genetic exchange element: the gene transfer agent of Rhodobacter capsulatus. Proc Natl Acad Sci U S A. 2000, 97: 859-864. 10.1073/pnas.97.2.859.PubMed CentralPubMedView ArticleGoogle Scholar
- Smith TG, Hoover TR: Deciphering Bacterial Flagellar Gene Regulatory Networks in the Genomic Era. 2009, Elsevier Inc., 257-295.Google Scholar
- Laub MT, Chen SL, Shapiro L, McAdams HH: Genes directly controlled by CtrA, a master regulator of the Caulobacter cell cycle. Proc Natl Acad Sci U S A. 2002, 99: 4632-4637. 10.1073/pnas.062065699.PubMed CentralPubMedView ArticleGoogle Scholar
- Petersen J, Brinkmann H, Pradella S: Diversity and evolution of repABC type plasmids in Rhodobacterales. Environ Microbiol. 2009, 11: 2627-2638. 10.1111/j.1462-2920.2009.01987.x.PubMedView ArticleGoogle Scholar
- Riedel K, Hentzer M, Geisenberger O, Huber B, Steidle A, Wu H, Hoiby N, Givskov M, Molin S, Eberl L: N-acylhomoserine-lactone-mediated communication between Pseudomonas aeruginosa and Burkholderia cepacia in mixed biofilms. Microbiology. 2001, 147: 3249-3262.PubMedView ArticleGoogle Scholar
- Neumann A, Patzelt D, Wagner-Döbler I, Schulz S: Identification of new N-acylhomoserine lactone signalling compounds of Dinoroseobacter shibae DFL-12 by overexpression of luxI genes. Chembiochem. 2013, 14: 2355-2361. 10.1002/cbic.201300424.PubMedView ArticleGoogle Scholar
- Pappas KM, Winans SC: A LuxR-type regulator from Agrobacterium tumefaciens elevates Ti plasmid copy number by activating transcription of plasmid replication genes. Mol Microbiol. 2003, 48: 1059-1073. 10.1046/j.1365-2958.2003.03488.x.PubMedView ArticleGoogle Scholar
- Hill NS, Kadoya R, Chattoraj DK, Levin PA: Cell size and the initiation of DNA replication in bacteria. PLoS Genet. 2012, 8: e1002549-10.1371/journal.pgen.1002549.PubMed CentralPubMedView ArticleGoogle Scholar
- Angelastro PS, Sliusarenko O, Jacobs-Wagner C: Polar localization of the CckA histidine kinase and cell cycle periodicity of the essential master regulator CtrA in Caulobacter crescentus. J Bacteriol. 2010, 192: 539-552. 10.1128/JB.00985-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Avery SV: Microbial cell individuality and the underlying sources of heterogeneity. Nat Rev Microbiol. 2006, 4: 577-587. 10.1038/nrmicro1460.PubMedView ArticleGoogle Scholar
- Smits WK, Kuipers OP, Veening JW: Phenotypic variation in bacteria: the role of feedback regulation. Nat Rev Microbiol. 2006, 4: 259-271. 10.1038/nrmicro1381.PubMedView ArticleGoogle Scholar
- Rex R, Bill N, Schmidt-Hohagen K, Schomburg D: Swimming in light: a large-scale computational analysis of the metabolism of Dinoroseobacter shibae. PLoS Comput Biol. 2013, 9: e1003224-10.1371/journal.pcbi.1003224.PubMed CentralPubMedView ArticleGoogle Scholar
- Novitsky JA, Morita RY: Morphological characterization of small cells resulting from nutrient starvation of a psychrophilic marine vibrio. Appl Environ Microbiol. 1976, 32: 617-622.PubMed CentralPubMedGoogle Scholar
- Anetzberger C, Pirch T, Jung K: Heterogeneity in quorum sensing-regulated bioluminescence of Vibrio harveyi. Mol Microbiol. 2009, 73: 267-277. 10.1111/j.1365-2958.2009.06768.x.PubMedView ArticleGoogle Scholar
- Anetzberger C, Schell U, Jung K: Single cell analysis of Vibrio harveyi uncovers functional heterogeneity in response to quorum sensing signals. BMC Microbiol. 2012, 12: 209-10.1186/1471-2180-12-209.PubMed CentralPubMedView ArticleGoogle Scholar
- Perez PD, Hagen SJ: Heterogeneous response to a quorum-sensing signal in the luminescence of individual Vibrio fischeri. PLoS One. 2010, 5: e15473-10.1371/journal.pone.0015473.PubMed CentralPubMedView ArticleGoogle Scholar
- Kovach ME, Elzer PH, Hill DS, Robertson GT, Farris MA, Roop RM, Peterson KM: Four new derivatives of the broad-host-range cloning vector pBBR1MCS, carrying different antibiotic-resistance cassettes. Gene. 1995, 166: 175-176. 10.1016/0378-1119(95)00584-1.PubMedView ArticleGoogle Scholar
- Piekarski T, Buchholz I, Drepper T, Schobert M, Wagner-Doebler I, Tielen P, Jahn D: Genetic tools for the investigation of Roseobacter clade bacteria. BMC Microbiol. 2009, 9: 265-10.1186/1471-2180-9-265.PubMed CentralPubMedView ArticleGoogle Scholar
- Thoma S, Schobert M: An improved Escherichia coli donor strain for diparental mating. FEMS Microbiol Lett. 2009, 294: 127-132. 10.1111/j.1574-6968.2009.01556.x.PubMedView ArticleGoogle Scholar
- Hahne F, LeMeur N, Brinkman RR, Ellis B, Haaland P, Sarkar D, Spidlen J, Strain E, Gentleman R: FlowCore: a bioconductor package for high throughput flow cytometry. BMC Bioinforma. 2009, 10: 106-10.1186/1471-2105-10-106.View ArticleGoogle Scholar
- Smyth GK: Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Limma: Linear Models for Microarray Data. 2005, New York: Springer, 397-420.Google Scholar
- Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.PubMedView ArticleGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate - a practical and powerful approach to multiple testing. J R Stat Soc Ser B-Methodological. 1995, 57: 289-300.Google Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Zuckerkandl E, Pauling L: Molecules as documents of evolutionary history. J Theor Biol. 1965, 8: 357-366. 10.1016/0022-5193(65)90083-4.PubMedView ArticleGoogle 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-2739. 10.1093/molbev/msr121.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.