Montserrat Elías-Arnanz, Marcos Peñalver-Mellado, S. Padmanabhan, and Francisco J. Murillo*
Departamento de Genética y Microbiología, Universidad de Murcia, 30100 Murcia, Spain
Received 21 January 2003/ Accepted 25 March 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The only reported prokaryotic example of an HMGA-type protein is CarD in the bacterium Myxococcus xanthus (28, 29, 31). Like HMGA proteins, CarD has multifunctional roles in vivo and has been shown to be involved in at least two processes: light-induced carotenogenesis and fruiting body formation (29). The randomly structured HMGA-like C-terminal domain of CarD is made up of a basic region of multiple AT hooks and a flanking highly acidic segment. This domain is very similar to human HMGA proteins in its physical, structural and DNA-binding properties (28, 31). CarD is, however, considerably larger and contains an N-terminal stretch of around 180 amino acids that is absent in eukaryotic HMGA proteins. This N-terminal domain, whose function remains to be elucidated, has defined secondary and tertiary structure, in contrast to the C-terminal HMGA-like region (31).
HMGA proteins are ubiquitous in higher eukaryotes (3, 8). However, as pointed out above, M. xanthus CarD is the only such protein identified so far in a prokaryote. The objective of this study was to examine whether HMGA-type proteins occur in other bacteria. As a first step we have done so for the bacterium Stigmatella aurantiaca, which belongs to the same taxonomic subgroup as M. xanthus and exhibits similar behavioral and developmental characteristics (40). We have identified a gene in S. aurantiaca that codes for a protein highly similar to M. xanthus CarD. Like CarD, its S. aurantiaca counterpart (CarDSa) contains an HMGA-like C-terminal region, as well as the N-terminal stretch of approximately 180 residues that is absent in eukaryotic HMGA proteins. The N-terminal regions of the two proteins are almost identical in sequence, while the HMGA-like C-terminal regions are less so. In CarDSa the latter region has one fewer AT hook and only one CKII phosphorylation site in the acidic part. In vitro analysis using purified proteins showed that CarDSa exhibits the same DNA-binding specificity as CarD and that CKII phosphorylation of CarDSa occurs but to a considerably lower extent than in CarD. In vivo, CarDSa can replace CarD in carotenogenesis and fruiting body formation. Our findings have been exploited to gain additional insights into the molecular bases that govern CarD activity in M. xanthus. We present evidence that the existence of HMGA domains in bacteria appears to be restricted to some other myxobacteria. The N-terminal domain in CarDSa and CarD, on the other hand, exists in various other bacteria as an independent module.
| MATERIALS AND METHODS |
|---|
|
|
|---|
was used for plasmid constructions, and strain BL21(DE3) was used for protein overexpression (44); both were grown in Luria broth medium (38). E. coli strains Q358 (2) and Y1090 (51) were used for infection by the
DASH and
gt11 phage libraries, respectively.
|
Southern hybridization analysis. Genomic DNA was isolated by using the Promega Wizard genomic DNA purification kit. Southern analysis was done by standard procedures (38). A 1.13-kb DNA fragment containing the entire carD gene (948 bp) or a fragment spanning nucleotides 354 to 948 of this gene was used as a probe. Low-stringency hybridization was performed at 59°C in the following buffer: 6x SSC (0.9 M sodium chloride plus 0.09 mM sodium citrate), 0.1% SDS, 5x Denhardt's solution, and 5% dextran sulfate. Blots were washed twice at room temperature and once at 59°C in 2x SSC-0.1% SDS.
Cloning of the S. aurantiaca gene homologous to carD.
Two independent genomic DNA libraries constructed in phages
DASH and
gt11 (generous gifts from H. U. Schairer, University of Heidelberg) were used to screen for the gene homologous to carD in S. aurantiaca (which is referred to as carDSa). The
DASH and
gt11 libraries were constructed with genomic DNA digested with SalI and HpaII, respectively. About 16,000 plaques generated from E. coli strain Q358 (for the
DASH phage library) or Y1019 (for the
gt11 phage library) were screened under low-stringency conditions with the 1.13-kb DNA probe containing carD by standard procedures (38). The most intense positive signals were selected and subjected to a second round of hybridization. Five positive clones were chosen from the
DASH library, and three were chosen from the
gt11 library. Phage DNA from each of the selected clones was isolated, cut with SalI (for
DASH) or EcoRI (for
gt11), electrophoresed, and tested again with the carD probe. Each of the five
DASH clones gave a single positive hybridization band of
5 kb which was purified and cloned into the vector pUC9-2 to generate plasmid pMAR541. Single hybridization bands, but of different sizes, were obtained with the three
gt11 clones, and of these, one of
2 kb was cloned into pUC19 to generate plasmid pMAR542. A DNA fragment of about 1.2 kb from each clone was sequenced twice along both strands.
Construction of an M. xanthus strain with gene carD deleted. pMAR975 is a pBJ114 derivative (22) which lacks the EcoRI site and contains a Kmr gene for positive selection and a galactose sensitivity (Gals) gene for negative selection (46). A 3.3-kb DNA fragment containing carD and about 1.2 kb of flanking DNA on each side was cloned into pMAR975 to obtain plasmid pMR2598. To generate a complete in-frame deletion of carD, pMR2598 was used as template for inverse PCR with the Expand long-template PCR system (Roche Applied Science) and two oligonucleotide primers with EcoRI overhangs (underlined): carD-DEL1 (5'-AAAGGAATTCGTCCCCCTCACGGGTGAGGT-3'), and carD-DEL2 (5'-AAAGGAATTCTGACAGCCCCATGGACCGAC-3'). The PCR-amplified fragment was cut with EcoRI and self-ligated to generate plasmid pMR2603, in which the entire carD gene is precisely deleted (from the ATG start codon to the nucleotide immediately upstream of the stop codon) and replaced by an EcoRI site. This in-frame deletion of carD is referred to as the carD3 allele. pMR2603 was electroporated into M. xanthus (strain DK1622), where it can be maintained only after integration into the chromosome by homologous recombination. Stable kanamycin-resistant transformants are thus merodiploids carrying wild-type carD as well as the carD3 allele. Cells having lost one of the two copies and the vector DNA through intramolecular recombination events were selected for on CTT plates supplemented with 10 mg of galactose per ml. Several Galr Kms colonies were picked and tested by PCR for the presence or absence of the carD deletion to isolate the carD deletion strain MR1900.
Complementation analyses.
For complementation analyses the following plasmids were constructed. (i) pMR2698 contains the carDSa gene inserted into the EcoRI site of pMR2603. The insert (initiator ATG to stop codon) was generated by PCR with S. aurantiaca genomic DNA as the template and the following primers with EcoRI overhangs (underlined): SaCD-RI-N (5'-AAAGAATTCATGCCAGAAGGACTCCAGCTC-3') and SaCD-RI-C (5'-AAAGAATTCTCACTACTCGGTCTCACCCTC-3'). (ii) pMR2745 is derived from pMR2698 by site-directed mutagenesis of the codon for Ser198 in carDSa to that for Ala by using the PCR overlap extension method (18). Compared to the normal arrangement for carD as in pMR2598, both pMR2698 and pMR2745 have the initiator ATG codon separated from the carD ribosomal binding site by an additional 6 bp (the EcoRI site). Consequently, in the complementation analyses carried out with pMR2698 and pMR2745, the following positive control plasmid, pMR2696, was used. (iii) pMR2696 contains the carD' allele, that is, the carD gene from the initiator ATG to the stop codon, inserted into the EcoRI site of pMR2603 (and thus separated from the ribosomal binding site by an additional 6 bp, as in constructs pMR2698 and pMR2745). The insert was generated by PCR with pMR2598 as template DNA and the following primers with EcoRI overhangs (underlined): MxCD-RI-N (5'-AAAGAATTCATGCCTGAAGGGTCCGCGTCA-3') and MxCD-RI-C (5'-AAAGAATTCTCAGCTCTCACCCTCGGGCGG-3'). (iv) pMR2768 is a derivative of pMR2598 containing a carD allele (carD
N) in which the DNA region coding for the N-terminal residues 2 to 178 of the protein has been deleted by following the same protocol described above for pMR2603. The PCR primers used in this case were carD-NDEL1 (5'-AAAGAATTCCATGTCCCCCTCACGGGTGAG-3') and carD-NDEL2 (5'-ACCGCAGCCGAATTCCAGCGCG-3', where GAA in the EcoRI site shown underlined replaces the wild-type GAG, both triplets coding for Glu179). All PCR-derived constructs described in this and other sections were verified by DNA sequencing.
Complementation analyses were performed with the carD deletion strain MR1900 (see above) as the recipient for electroporation with each of the four plasmids described above, none of which replicate in M. xanthus but which can integrate into the bacterial genome by homologous recombination. The resulting merodiploids (Kmr) bear both the carD deletion copy and the incoming allele. When complementation was observed, the merodiploids were further processed by using the same procedure described above for strain MR1900 to generate a strain with a copy of the incoming allele alone (strains MR1902 and MR1903 and the control strain MR1901 [Table 1]). These Kms strains could then be electroporated with pDAH217 (Kmr), which carries a lacZ transcriptional probe fused to the carD-dependent light-inducible carQRS promoter (19, 29), in order to quantify the ability of the complementing gene to replace carD. Promoter activity with and without illumination was assessed qualitatively by monitoring lacZ reporter gene expression through colony color formation on plates containing 40 µg of X-Gal (5-bromo-4-chloro-3-indolyl-ß-D-galactoside) per ml and quantitatively by measurements of ß-galactosidase activity as described previously (4).
Overexpression and purification of CarDSa.
The carDSa gene was PCR amplified with S. aurantiaca genomic DNA as the template to generate a fragment with an NdeI site containing the initiator Met and a BamHI site immediately downstream of the stop codon. The PCR-amplified fragment was purified and then cloned into the NdeI-BamHI sites of the overexpression vector pET11b or into pET15b to produce His6-tagged protein (44). Overexpression of CarDSa and its purification by ion-exchange chromatography, first from phosphocellulose and then by MonoS high-performance liquid chromatography (HPLC) (AKTA-Amersham Biosciences), were as described for CarD (31). His6-tagged proteins were purified by using TALON metal affinity resin and the accompanying purification protocol for native conditions (Clontech, Palo Alto, Calif.) and by MonoS HPLC. For concentrations determined from the absorbance at 280 nm, the values used for
280 (in molar-1 centimeter-1) are 8,480 and 9,970 for CarD and CarDSa, respectively (30).
Analytical size exclusion chromatography.
The apparent molecular mass for CarDSa was estimated by analytical size exclusion chromatography at room temperature with a Superdex-200 HPLC column (Amersham Biosciences). One hundred microliters of protein (
10 µM) was injected into the column, which was previously equilibrated with buffer A containing 200 mM NaCl, and its elution was tracked by the absorbance at 280, 235, and 220 nm at a flow rate of 0.3 ml/min. The molecular mass was estimated from the elution volume Ve and the calibration curve log molecular mass = 7.91 - 0.23Ve, generated as described elsewhere (31). The identity of eluted CarDSa was confirmed by Coomassie blue staining after SDS-PAGE and Western blotting with anti-CarD antibodies.
In vitro DNA-binding and CKII phosphorylation assays.
DNA binding was examined by electrophoretic mobility shift assay with a radiolabeled 169-bp DNA probe that spans positions -117 to +52 relative to the transcription start site and contains the CarD-binding site in the carQRS promoter region (28, 31). The fragment was generated by PCR with pDAH231 as the template (26) and the synthetic oligonucleotide primers carQRS5 (5'-GGGCAGGACGGGATGCTGCTG-3') and carQRS6 (5'-CCGTCGCGAAACCGTTCCATGA-3'). The primer carQRS5 was labeled with [
-32P]ATP and T4 polynucleotide kinase prior to its addition to the PCR mix. The electrophoretic mobility shift assay was carried out in 20-µl-total reaction volumes containing 1 to 3 pM end-labeled probe (
13,000 cpm), 500 nM protein, and 1 µg of double-stranded poly(dA-dT), poly (dG-dC), or poly(dI-dC) as nonspecific competitor DNA in binding buffer (50 mM NaCl, 15 mM HEPES, 4 mM Tris [pH 7.9], 1 mM dithiothreitol, 10% glycerol, 1 mg of bovine serum albumin per ml, and 0.1% Nonidet P-40). After a 30-min equilibration period at 4°C, DNA binding was analyzed by electrophoresis for 1 to 1.5 h in nondenaturing 4% polyacrylamide gels (acrylamide-bisacrylamide, 37.5:1) that were prerun at 200 V and 4°C for 30 min in 0.5x TBE buffer (45 mM Tris base, 45 mM boric acid, 1 mM EDTA). Gels were dried and analyzed by autoradiography.
CKII phosphorylation of purified CarD and CarDSa was examined by using recombinant human CKII (New England BioLabs). The protein (0.3 to 1 µM) was treated with 0.5 U of CKII and 1µCi of [
32-P]ATP in DNA-binding buffer for 1 h at 30°C. After removal of unincorporated [
32-P]ATP through Sephadex G-50, the samples were subjected to SDS-PAGE and analyzed by autoradiography. To compare relative amounts of protein used in the assay, mock CKII phosphorylation reactions were carried out under identical conditions but without labeled ATP and examined by Coomassie blue staining after SDS-PAGE.
Sequence comparisons and analysis. Sequence database searches were performed with the BLAST suite of programs provided at the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/BLAST), Baylor College of Medicine search launcher (http://searchlauncher.bcm.tmc.edu) and/or European Bioinformatics Institute (http://www.ebi.ac.uk). Similarity searches of the nonredundant protein sequence database were done with the gapped BLASTP program (1). Bacterial proteins with the AT hook motif were also retrieved from the SMART database (http://smart.embl-heidelberg.de). Protein sequences were aligned by using the CLUSTAL software at BCM and analyzed for phosphorylation sites by using PROSITE (http://us.expasy.org/prosite/). Open reading frames (ORFs) in the cloned DNA sequence were identified by using the ORF Finder at NCBI (http://www.ncbi.nlm.nih.gov/gorf/gorf.html).
Nucleotide sequence accession number. The nucleotide sequence obtained from the pMAR541 clone, which contains the whole carDSa gene, has been deposited at the DDBJ/EMBL/GenBank databases (accession number AJ536154).
| RESULTS |
|---|
|
|
|---|
|
The gene homologous to carD in S. aurantiaca, designated carDSa, was cloned by using phage genomic DNA libraries (see Materials and Methods). Both strands of a 1.2-kb stretch were sequenced, and potential ORFs were identified in computer-aided searches, taking into account the bias for G or C at the third codon position (owing to the high GC content of myxobacterial DNA [6, 12]). One ORF of 915 bp that would yield a gene product of 305 amino acids with high sequence similarity to the 316-residue-long CarD was identified. This would be in accord with the results from the immunoblot and Southern hybridization analyses described above (Fig. 1). The sequence alignment between CarD and its S. aurantiaca counterpart (Fig. 2) indicates that the characteristic HMGA-type domain is also present in CarDSa. It is, however, precisely in this domain where most of the differences in the primary structures of CarDSa and CarD appear. One significant difference is the presence of only three AT hooks (the RGRP DNA-binding motif) in CarDSa as opposed to the four in CarD; another is that only one consensus CKII phosphorylation site, S198, can be identified in the acidic region of CarDSa by using PROSITE, in contrast to the five predicted in CarD. In fact, other than Ser198, no other Ser or Thr appears in the entire acidic region of CarDSa. By contrast, the approximately 180-residue N-terminal segment that precedes the HMGA-like domain in CarDSa shows high sequence identity (86%) with that in CarD. This segment includes a region that may be involved in CarD dimerization (31).
|
CarD is essential for the expression of the light-inducible carQRS operon, a key regulatory gene cluster in M. xanthus carotenogenesis (29). A site containing two appropriately spaced AT-rich tracts upstream of the promoter of carQRS is required in vivo for its expression, and in vitro, CarD binds to this site with the characteristic HMGA minor-groove binding specificity (28, 31). Figure 3A shows that CarDSa binds to a 169-bp double-stranded DNA probe that includes the M. xanthus carQRS promoter region and the CarD-binding site. Moreover, the retarded band is observed in the presence of poly(dG-dC) as a nonspecific competitor (lane 4) but not with poly(dA-dT) or poly(dI-dC) (lanes 2 and 3). This behavior has been used to infer that HMGA1a binds to the minor groove of AT-rich tracts (43), and this is also observed with CarD (31). Therefore, our results show that CarDSa has the same AT hook DNA-binding specificity as CarD (or HMGA1a).
|
Complementation of a carD deletion by carDSa. The considerable similarity between CarD and CarDSa led us to examine if the latter would substitute for CarD function in M. xanthus. For this, we first constructed an M. xanthus strain (MR1900) with a precise in-frame deletion of carD (see Materials and Methods). As expected, this mutant strain was defective both in light-induced carotenogenesis and in fruiting body formation. To check for complementation by the carDSa gene, plasmid pMR2698 (Kmr) was introduced into strain MR1900 by electroporation. Electroporants should arise by integration of the plasmid into the M. xanthus chromosome by homologous recombination. The resulting Kmr colonies would therefore be merodiploids bearing both the carD deletion allele (carD3) and the carDSa gene. On plates, all of the Kmr electroporants showed light-induced carotenogenesis, as judged by the intense red color developed on illumination with blue light. Several of these electroporants were also picked on TPM and CF agar plates to monitor multicellular development, and the fruiting bodies were examined for the presence of spores as previously described (29). In all cases, development proceeded at the same rate as for the wild-type M. xanthus strain and yielded the same number of mature fruiting bodies. The behavior was equivalent to that observed in control experiments where MR1900 was electroporated with plasmid pMR2696 (Kmr), which has the carD' gene instead of carDSa. All of these results are summarized in Table 2. Restoration of the wild-type phenotype for both carotenogenesis and fruiting body development on introduction of the carDSa gene into the M. xanthus strain MR1900 clearly indicates that carD can be replaced by carDSa in vivo.
|
|
N) coding for the C-terminal segment spanning residues 179 to 316. We have previously shown that this CarD fragment is stably overexpressed in E. coli and that it retains the in vitro DNA-binding and CKII phosphorylation properties of CarD, but it appears to be largely monomeric, unlike the whole protein (31). Plasmid pMR2768 was electroporated into strain MR1900 to obtain merodiploids bearing the carD3 allele as well as the carD
N truncated form. These Kmr electroporants were deficient in both light-induced carotenogenesis and fruiting body formation (Table 2), even though the corresponding CarD C-terminal fragment was stably expressed, as verified by Western blots of whole cell extracts (data not shown). This demonstrates that the 180-residue N-terminal region that is so highly conserved in CarD and CarDSa (and absent in HMGA proteins) is essential for the protein to carry out its in vivo roles in carotenogenesis and fruiting body formation; the presence of the HMGA-like domain alone is not sufficient. Of the five predicted CKII phosphorylation sites in the acidic region of CarD, only S198 is present in CarDSa (Fig. 2), and the mutation S198A eliminates in vitro phosphorylation of CarDSa by CKII, as described earlier. We exploited this and the ability of CarDSa to substitute for CarD function to examine the role, if any, of phosphorylation by CKII-type kinases in vivo. For this, plasmid pMR2745, containing the carDSa(S198A) mutant gene, was introduced into strain MR1900. This restored the wild-type phenotype for carotenogenesis and multicellular development, just as do the plasmids pMR2698 (containing carDSa) and pMR2696 (containing carD') (Table 2). Moreover, light induction of PQRS by the CarDSa(S198A) mutant was essentially identical to that brought about by CarDSa (Fig. 4). These results suggest that phosphorylation of the HMGA domain by a CKII-type kinase is not a necessary regulatory event for normal carotenogenesis and fruiting body formation in M. xanthus.
carD-like genes in other bacteria. M. xanthus CarD had been the only protein with an HMGA-type domain identified in a prokaryote until the description in this study of the homologous CarDSa protein in S. aurantiaca. Based on 16S RNA analysis, the myxobacteria form a monophyletic grouping consisting of three distinct subgroups (Myxococcus, Chondromyces, and Nannocystis), and both M. xanthus and S. aurantiaca fall into the Myxococcus subgroup (42). To examine if the occurrence of HMGA-type proteins in myxobacteria is a general phenomenon, a DNA fragment corresponding to nucleotides 354 to 948 of carD (coding for its entire HMGA part) was used to probe under low-stringency conditions the genomic DNAs from M. coralloides and C. fuscus (Myxococcus subgroup), P. cellulosum (Chondromyces subgroup), and N. exedens (Nannocystis subgroup). As with M. xanthus and S. aurantiaca, a strong hybridization band was observed for M. coralloides and C. fuscus but not for P. cellulosum and N. exedens (Fig. 5A). Given that the hybridization was done under low-stringency conditions, the faint and diffuse signals observed for the last two strains could stem from the high GC content of the genomic DNA of myxobacteria. Consequently, these results suggest that the existence of HMGA domain-containing proteins in myxobacteria may be largely confined to the members of the Myxococcus subgroup.
|
subdivision of the proteobacteria (42). Therefore, we also performed a Southern hybridization analysis with the same probe used above and ClaI- or EcoRI-digested genomic DNA from Bdellovibrio sp. strain CP1. As shown in Fig. 5B, no hybridization signal was detected in this member of the same taxonomic subdivision as the myxobacteria. Whether bacteria other than myxobacteria contain proteins with HMGA-type domains was further addressed by similarity searches of the nonredundant protein sequence database and of the microbial genome protein database at NCBI (77 complete and 33 partial bacterial genome sequences at the time of this analysis, covering a wide range of taxonomic groups but no myxobacteria), using the gapped BLASTP program. Use of the AT hook segment of CarD or CarDSa alone as the query in a search for short, nearly exact matches in the nonredundant sequence database yielded as best hits eukaryotic HMGA or HMGA-like proteins. We also encountered one hit to a hypothetical protein in the bacterium Ralstonia metallidurans (accession no. ZP_00021411) which has one PGRP sequence, a less frequent core AT hook motif (3). That AT hooks are very rare among bacteria was further confirmed by searching the microbial genome database. Use of the AT hook region of CarD or CarDSa as the search query provided just one hypothetical protein in the bacterium Rhodopseudomonas palustris (accession no. ZP_00008414), which has two closely spaced RGRP sequences and one PGRP motif. However, neither the R. metallidurans hypothetical protein nor that in R. palustris contains a highly acidic region, which invariably lies adjacent to the AT hook segment in HMGA proteins. The SMART database with the AT hook (SMART accession number SM0384) as the query in a domain search of bacterial proteins lists 17 additional proteins. Interestingly, the majority of these are putative DNA-binding proteins, some of which could be implicated in transcriptional control, for example, transcriptional regulators of the TetR type in Caulobacter crescentus (accession no. AAK22321) and Streptomyces coelicolor (accession no. T36295) and of the WhiB type in Mycobacterium tuberculosis (accession no. NP_337818). Among the 17 proteins, a single RGRP AT hook is found in 13, one PGRP motif is found in another, and 2 have the very uncommon AGRP and LGRP motifs (3). The remaining protein has an RGRP AT hook and the very infrequent VGRP motif (Saccharopolyspora erythraea, accession no. AAL78056). Again, in none of these cases is the AT hook portion associated with an adjacent highly acidic region. Thus, based on presently available data, it appears that prokaryotic HMGA-type proteins are largely restricted to myxobacteria.
By contrast, 25 proteins found exclusively in bacteria showed significant similarity to the N-terminal domain of CarD or CarDSa in a search of the NCBI microbial genome protein database. All of these are hypothetical proteins grouped in the conserved-domain protein family pfam02559 (5), whose defining element is the CarD segment between residues 9 and 158. Figure 6 displays the multiple sequence alignment of the N-terminal portions of CarD and CarDSa and the eight proteins showing the highest similarity (sequence identities of 26 to 33%; similarities of 50 to 59%). It should be noted that the similarity to the 180-residue N-terminal segment of CarD and CarDSa extends over the entire length of the primary sequence of the homologs, whose sizes range from 153 to 165 residues, except for the 198-residue Corynebacterium glutamicum protein. Interestingly, a search of M. xanthus sequences in the Cereon microbial genome database allowed us to identify a 164-residue hypothetical protein with 35% identity and 58% similarity to the CarD/CarDSa N terminus, which is thus a new member of the pfam02559 protein family (Fig. 6). These results therefore uncover a distinct prokaryotic domain that exists as an independent module in various bacteria but which is linked to an HMGA-type DNA-binding domain in CarD and CarDSa.
|
| DISCUSSION |
|---|
|
|
|---|
The studies reported here demonstrate that the S. aurantiaca protein CarDSa shares several of the molecular properties of CarD in vitro and can replace the latter in vivo in its dual role in M. xanthus carotenogenesis and fruiting-body formation. This strongly suggests that carDSa is the S. aurantiaca ortholog of carD. The ability of CarDSa to mimic CarD in M. xanthus has also provided two additional insights into the molecular basis for CarD activity. The first is that CarD function is not impeded when there is one fewer AT hook, as in CarDSa. The second is on the significance of in vivo phosphorylation of CarD by a CKII-type kinase. Phosphorylation of eukaryotic HMGA proteins has been linked to their multifunctional roles in vivo, especially in the cell cycle, differentiation, and development (8, 27, 36, 41, 49). For instance, CKII phosphorylation of the acidic region appears to fine-tune the inherent structural plasticity of the randomly structured HMGA proteins, and this, as a consequence, modulates their intracellular stabilities as well as DNA-binding affinities. A similar scenario could have been envisaged for CarD given its role in M. xanthus cellular differentiation and development. However, as shown in this study, phosphorylation of CKII sites in the acidic region of the HMGA domain of CarD does not appear to play such an essential regulatory role, as can be inferred from the observed functional equivalence of CarDSa and its mutant form lacking the only target for CKII phosphorylation. Although CKII-type kinases have not yet been found in M. xanthus, a number of eukaryotic-type serine/threonine protein kinases have been reported, whose substrates remain to be identified (21). Whereas our results do not exclude the possibility that a CKII-type kinase activity exists in M. xanthus, they do suggest that CarD is not an obligatory substrate.
The ability of CarDSa to function in carotenogenesis and fruiting-body formation in M. xanthus suggests that it may be involved in similar roles in its natural context. S. aurantiaca is also capable of developing fruiting bodies, but these are considerably more elaborate than those in M. xanthus (12). Interestingly, efficient fruiting-body formation in S. aurantiaca is driven by light (34). Useful leads in defining the roles that CarDSa plays in these and other processes in S. aurantiaca may come from our comparison with the properties exhibited by M. xanthus CarD.
CarD-like proteins in other bacteria.
The results obtained in this study indicate that proteins containing HMGA-like domains in prokaryotes are found primarily in myxobacteria (order Myxococcales). In particular, such proteins appear to be confined to the suborder Cystobacterineae, which includes M. xanthus, S. aurantiaca, M. coralloides, and C. fuscus. We also have no evidence for the presence of HMGA proteins in the bdellovibrios, which belong to the taxonomic group closest to myxobacteria in the
subdivision of the proteobacteria. Our analysis of the many partial or complete available genome sequences for bacteria from other taxonomic groups yielded 19 proteins which have AT hook sequences. One of these has three AT hooks, one has two AT hooks, and the rest have no more than a single AT hook. However, a flanking highly acidic region, which is a hallmark of HMGA domains, is not apparent in any of these proteins. That proteins containing HMGA domains appear to be exclusive for one specific suborder of myxobacteria suggests two possible alternatives for the evolutionary origins of carD in myxobacteria: a de novo occurrence in the specific phylogenetic subgroup of myxobacteria or a consequence of horizontal acquisition of a eukaryotic HMGA-type gene. Nevertheless, when a protein is widely prevalent in eukaryotes but occurs in only one or very few bacterial species, horizontal transfer is the most parsimonious explanation (24). Horizontal gene transfer into the myxobacteria may be linked to their particular lifestyles, i.e., choice of habitat, feeding, and/or social behavior. In fact, the hypothesis that feeding habits may account for lateral transfer (13) has been invoked to account for the particular organization of the gene encoding the ß-1-4-endoglucanase in M. xanthus and S. aurantiaca, both of which are scavengers that prey on other soil microorganisms (35). Moreover, myxobacteria have generally been considered to be phylogenetically distinct and an evolutionarily advanced group of bacteria (42). The formation of fruiting bodies and spores in myxobacteria represent rather complex developmental and cellular differentiation cycles that involve elaborate signal transduction and gene regulatory cascades. These processes, which are generally not observed with most other bacteria, are more akin to those occurring in eukaryotes. It is therefore not surprising that a number of protein factors typically found in eukaryotes have analogs in myxobacteria. In M. xanthus these include, besides CarD, the serine/threonine protein kinases referred to above and at least one associated phosphatase (45) and CarF, a recently identified regulatory protein in carotenogenesis (15).
Our comparative sequence analysis of CarD and CarDSa revealed homologs to the N-terminal non-HMGA domain, but these were all entirely prokaryotic in origin. CarD and CarDSa could therefore exemplify interkingdom domain fusion between a preexisting bacterial domain and an acquired eukaryotic domain, a phenomenon that is particularly common in Actinomycetes (24). Furthermore, the additional presence of a distinct stand-alone version of the prokaryotic domain in M. xanthus is consistent with a probable two-step evolutionary process for the origin of CarD and CarDSa, that is, capture of the eukaryotic HMGA portion with subsequent fusion to the bacterial N-terminal part by recombination (24). Consistent with this two-domain makeup is our analysis of CarD, which revealed that the N-terminal and HMGA-type C-terminal parts are structurally distinct modules (31). Any selective advantage conferred by this particular domain organization would necessarily be speculative. The GC-rich nature of most myxobacteria would clearly narrow down the search for the specific AT-rich DNA-binding sites by the HMGA domain. Any speculation with regard to the N-terminal domain, on the other hand, must await future work aimed at identifying the exact nature of the part it plays in CarD/CarDSa functions. At least for CarD, as we have shown above, the association of the N-terminal region with the HMGA domain is an essential functional requirement in carotenogenesis and multicellular development.
| ACKNOWLEDGMENTS |
|---|
DASH and
gt11 S. aurantiaca DNA library, Victoriano Garre for technical advice with the library screening experiments, and J. A. Madrid for technical assistance. Christophe Hoor contributed to the construction of the carD deletion strain. We are grateful to F. Torrella and A. Sánchez-Amat (University of Murcia) for some of the bacterial species used in this work and to David Hodgson (University of Warwick, Warwick, United Kingdom) for plasmids pDAH217 and pDAH231. We are also grateful to the Monsanto Company for access to the M. xanthus sequence database (now available at the TIGR microbial database). We thank the anonymous referees for useful comments. This work was supported by the Spanish Ministerio de Ciencia y Tecnología (grant BMC2000-1006 to F.J.M. and Programa Ramón y Cajal to S.P.), Ministerio de Educación y Cultura (fellowship to M.L.C.), and Fundación Séneca (fellowship to M.P.-M.).
| FOOTNOTES |
|---|
Present address: Department of Immunology and Oncology, Centro Nacional de Biotecnología-CSIC, E-28049 Madrid, Spain. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles: