Identification of a vegetative promoter in Myxococcus xanthus. A protein that has homology to histones

The evolution, evolvability and engineering of gene regulatory DNA

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness^1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces^4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using such models to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatory sequences and addressing fundamental questions in regulatory evolution.

nsd, a locus that affects the Myxococcus xanthus cellular response to nutrient concentration

Expression of the previously reported Tn5lac Omega4469 insertion in Myxococcus xanthus cells is regulated by the starvation response. Interested in learning more about the starvation response, we cloned and sequenced the region containing the insertion. Our analysis shows that the gene fusion is located in an open reading frame that we have designated nsd (nutrient sensing/utilizing defective) and that its expression is driven by a sigma70-like promoter. Sequence analysis of the nsd gene product provides no information on the potential structure or function of the encoded protein. In a further effort to learn about the role of nsd in the starvation response, we closely examined the phenotype of cells carrying the nsd::Tn5lac Omega4469 mutation. Our analysis showed that these cells initiate development on medium that contains nutrients sufficient to sustain vegetative growth of wild-type cells. Furthermore, in liquid media these same nutrient concentrations elicit a severe impairment of growth of nsd cells. The data suggest that the nsd cells launch a starvation response when there are enough nutrients to prevent one. In support of this hypothesis, we found that, when grown in these nutrient concentrations, nsd cells accumulate guanosine tetraphosphate, the cellular starvation signal. Therefore, we propose that nsd is used by cells to respond to available nutrient levels.

Identification of an activator protein required for the induction of fruA, a gene essential for fruiting body development in Myxococcus xanthus

Myxococcus xanthus exhibits social behavior and multicellular development. FruA is an essential transcription factor for fruiting body development in M. xanthus. In the present study, the upstream promoter region was found to be necessary for the induction of fruA expression during development. A cis-acting element required for the induction was identified and was located between nucleotides -154 and -107 with respect to the transcription initiation site. In addition, it was found that two binding sites exist within this element of the fruA promoter. By using DNA affinity column chromatography containing the cis-acting element, a fruA promoter-binding protein was purified. The purified protein was shown by N-terminal sequence analysis to be identical to MrpC, a protein identified previously by transposon insertion mutagenesis as an essential locus for fruiting body development [Sun, H. & Shi, W. (2001) J. Bacteriol. 183, 4786-4795]. Furthermore, fruA mRNA was not detectable in the mrpC::km strain, demonstrating that MrpC is essential for fruA expression. Moreover, mutational analysis of the binding sites for MrpC in the fruA promoter indicates that binding of MrpC activates transcription of fruA in vivo. This report provides evidence for a direct molecular interaction involved in temporally regulated gene expression in M. xanthus.

Transcriptional activation of a heat-shock gene, lonD, of Myxococcus xanthus by a two component histidine-aspartate phosphorelay system

In vitro transcription of lonD, a heat-shock gene from Myxococcus xanthus, was stimulated in the presence of extract from heat-shocked cells. For this stimulation the upstream promoter region of lonD was found to be essential. Activation of lonD transcription was also observed when extract from non-heat-shocked cells was heat treated in vitro at 42 degrees C for 10 min. A DNA binding assay and footprinting analysis revealed that a factor(s) binds to the upstream region from -122 to -107 with respect to the transcription initiation site. This region was required for heat-shock induction of lonD expression both in vitro and in vivo. The lonD promoter-binding protein named HsfA was purified, and its gene was cloned. Analysis of the DNA sequence reveals that HsfA is a response regulator of the two-component system and shows high sequence similarity to the NtrC family or the enhancer-binding proteins. Upstream of hsfA, a gene encoding a histidine kinase was identified and named hsfB. HsfB was found to be autophosphorylated and able to phosphorylate HsfA. HsfA with HsfB activated in vitro transcription of lonD in a manner dependent on RNA polymerase containing SigA, the housekeeping sigma factor of M. xanthus.

Biology and global distribution of myxobacteria in soils

This review presents an overview of the present status of the biology of the myxobacteria, including the molecular biology of the systems that control and regulate myxobacterial gliding movement and morphogenesis. The present status of myxobacterial taxonomy and phylogeny is described. The evolutionary biology of the myxobacteria is emphasized with respect to their social behavior and the molecular basis of their signal chains. Most important within the metabolic physiology are the biologically active secondary metabolites of myxobacteria and their molecular mechanisms of action. The global distribution of myxobacteria in soils is described on the basis of data given in the literature as well as of comprehensive analyses of 1398 soil samples from 64 countries of all continents. The results are analyzed with respect to the spectrum and number of species depending on ecological and habitat-specific factors. The myxobacterial floras of different climate zones are compared. Included are myxobacterial species adapted to extreme biotopes. The efficiency of different methods used presently for isolation of myxobacteria is compared.

Identification of the Omega4499 regulatory region controlling developmental expression of a Myxococcus xanthus cytochrome P-450 system

Omega4499 is the site of a Tn5 lac insertion in the Myxococcus xanthus chromosome that fuses lacZ expression to a developmentally regulated promoter. Cell-cell interactions that occur during development, including C signaling, are required for normal expression of Tn5 lac Omega4499. The DNA upstream of the Omega4499 insertion has been cloned, and the promoter has been localized. Analysis of the DNA sequence downstream of the promoter revealed one complete open reading frame and a second partial open reading frame that is interrupted by Tn5 lac Omega4499. The predicted products of these open reading frames are highly similar to reductase and oxidase components of bacterial cytochrome P-450 systems, which allow catabolism or anabolism of unusual compounds. However, the function of the gene products of the Omega4499 locus remains unclear because M. xanthus containing Tn5 lac Omega4499 exhibits no apparent defect in growth, developmental aggregation, fruiting body formation, or sporulation. Deletion analysis of the Omega4499 regulatory region showed that multiple DNA elements spanning more than 500 bp upstream of the transcriptional start site contribute to developmental promoter activity. At least two DNA elements, one downstream of -49 bp and one between -49 and -218 bp, boosted activity of the promoter in response to intercellular C signaling. Three sequences in the Omega4499 promoter region, centered at -55, -33, and -1 bp, nearly match a 7-bp sequence found in other C signal-dependent promoters. We propose that these sequences, matching the consensus sequence 5'-CAYYCCY-3', be called C box sequences, and we speculate that these sequences are cis-acting regulatory elements important for the expression of M. xanthus genes that depend upon intercellular C signaling during development.

A new sigma factor, SigD, essential for stationary phase is also required for multicellular differentiation in Myxococcus xanthus

BACKGROUND

Myxococcus xanthus is a gram-negative bacterium that undergoes spectacular development to form multicellular fruiting bodies under nutrient deprivation. Inside a fruiting body, vegetative cells differentiate into spores. A number of sigma factors have been shown to play roles in the regulation of gene expression in the M. xanthus life cycle. Additional sigma factors were searched to further explore the M. xanthus life cycle.

RESULTS

A new sigma factor was identified, SigD, which consists of 297 amino acid residues. Two transcription initiation sites for the sigD gene were detected by primer extension analysis using total RNA from the vegetative and developmental cells, one of which was specific for development. The characterization of sigD-lacZ fusion strains demonstrated that sigD expression increased during entry into stationary phase of vegetative growth and during early development. A deletion mutant of sigD exhibited growth defects during the late-log phase and stationary phase, with dramatically reduced cell viability. The patterns of protein synthesis at late log phase of vegetative growth and at early development on CF agar plates were significantly different between the deletion mutant and the wild-type strain. The deletion mutant was more sensitive to various stresses when compared with the wild-type strain and did not accumulate trehalose in response to osmotic stress. The deletion mutant also showed a significant delay in fruiting body formation and sporulation and yielded fewer spores than the wild-type strain.

CONCLUSIONS

SigD shows characteristic features of the stationary phase sigma factors and also plays important roles in multicellular differentiation of M. xanthus.

Translation initiation factor IF2 of the myxobacterium Stigmatella aurantiaca: presence of a single species with an unusual N-terminal sequence

The structural gene for translation initiation factor IF2 (infB) was isolated from the myxobacterium Stigmatella aurantiaca on a 5.18-kb BamHI genomic restriction fragment. The infB gene (ca. 3.16 kb) encodes a 1,054-residue polypeptide with extensive homology within its G domain and C terminus with the equivalent regions of IF2s from Escherichia coli, Bacillus subtilis, Bacillus stearothermophilus, and Streptococcus faecium. The N-terminal region does not display any significant homology to other known proteins. The S. aurantiaca infB gene encodes a single protein which cross-reacted with antiserum to E. coli IF2 and was able to complement an E. coli infB mutant. The S. aurantiaca IF2 is distinguished from all other IF2s by a sequence of 160 residues near the N terminus that has an unusual composition, made up essentially of alanine, proline, valine, and glutamic acid. Within this sequence, the pattern PXXXAP is repeated nine times. Complete deletion of this sequence did not affect the factor's function in initiation of translation and even increased its capacity to complement the E. coli infB mutant.

Photolyase of Myxococcus xanthus, a Gram-negative eubacterium, is more similar to photolyases found in Archaea and "higher" eukaryotes than to photolyases of other eubacteria

We report the identification of the gene encoding a DNA photolyase (phrA) from the Gram-negative eubacterium Myxococcus xanthus. The deduced amino acid sequence of M. xanthus photolyase indicates that the protein contains 401 amino acids (Mr 45,071). By comparison of the amino acid and DNA sequences with those of other known photolyases, it has been found that it is more similar to the deduced amino acid sequences of the photolyases of "higher" eukaryotes than to the photolyases of other eubacteria. Recombinant plasmids carrying M. xanthus phrA rescue the photoreactivation activity of an irradiated strain of Escherichia coli with a deletion in phrA. This rescue is light-dependent.

Suppression of a signaling defect during Myxococcus xanthus development

The csgA gene encodes an extracellular protein that is essential for cell-cell communication (C-signaling) during fruiting body development of Myxococcus xanthus. Two transposon insertions in the socABC operon, soc-560 and socC559, restore development to csgA null mutants. Mixing soc-560 csgA cells or socC559 csgA cells with csgA cells at a ratio of 1:1 stimulated the development of csgA cells, suggesting that soc mutations allow cells to produce the C-signal or a similar molecule via a csgA-independent mechanism. The socABC operon contains the following three genes: socA, a member of the short-chain alcohol dehydrogenase gene family; socB, a gene encoding a putative membrane anchoring protein; and socC, a negative autoregulator of socABC operon expression. Both suppressor mutations inactivate socC, leading to a 30- to 100-fold increase in socA transcription; socA expression in suppressor strains is at least 100-fold higher than csgA expression during all stages of development. The amino acid sequence of SocA has 28% identity and 51% similarity with that of CsgA. We suggest that CsgA suppression is due to overproduction of SocA, which can substitute for CsgA. These results raise the possibility that a cell surface dehydrogenase plays a role in C-signaling.

A cluster of structural and regulatory genes for light-induced carotenogenesis in Myxococcus xanthus

In the bacterium Myxococcus xanthus, several genes for carotenoid synthesis lie together at the carA-carB chromosomal locus and are co-ordinately activated by blue light. A 12-kb DNA stretch from wild-type M. xanthus has been sequenced that includes the entire carA-carB gene cluster. According to sequence analysis, the cluster contains 11 different genes. Intergenic distances are very short or nil (implying translational coupling), giving further support to previous evidence indicating that most (or all) of the genes in the cluster form a single operon. At the promoter region, a potential -35 site for the binding of sigma factors is found. However, the -10 region shows little similarity with analogous sites in other bacterial promoters. Five (possibly six) genes in the carA-carB operon code for enzymes acting on early or late steps of the pathway for carotenoid synthesis. Other genes in the operon show no overall similarity with previously known genes. However, peptide stretches in the predicted products of two genes exhibit strong similarity with the DNA binding domain of the MerR family of transcriptional regulators. At least one of the predicted DNA-binding domains is altered in a mutant strain affected in light-regulation of the car genes.

An early A-signal-dependent gene in Myxococcus xanthus has a sigma 54-like promoter

A-signaling plays an essential role in the early stages of Myxococcus xanthus fruiting body development. Expression of the 452I gene, which is regulated at the level of RNA accumulation, depends on starvation and on A-signaling. To identify the cis-acting regulatory elements which allow gene 4521 to respond to the nutritional and A-factor signals, the 4521 transcription start site was mapped. The region just upstream of the start site showed sequence similarity to the sigma 54 family of promoters and to the developmentally regulated mbhA promoter of M. xanthus. A mutational analysis of this region established that the bases which were conserved between the sigma 54 consensus, mbhA, and 4521 promoters were also important for 4521 promoter activity. Changes which altered the spacing between two conserved regions centered around positions -14 and -24 abolished promoter activity. In contrast, mutations in a putative -10 region for a sigma 70-like promoter had little effect on expression of 4521. Despite their similar promoter regions, the regulation of the 4521 and mbhA genes was shown to differ with respect to timing of expression and requirement for a solid surface and extracellular signals. This suggests a model in which different activator proteins may be responsible for regulating expression of these two genes.

Genetic and molecular analysis of the tRNA-tufB operon of the myxobacterium Stigmatella aurantiaca

The tufB gene, encoding elongation factor Tu (EF-Tu), from the myxobacterium Stigmatella aurantiaca was cloned and sequenced. It is preceded by four tRNA genes, the first ever described in myxobacteria. The tRNA synthesized from these genes and the general organization of the locus seem identical to that of Escherichia coli, but differences of potential importance were found in the tRNA sequences and in the intergenic regions. The primary structure of EF-Tu was deduced from the tufB DNA sequence. The factor is composed of 396 amino acids, with a predicted molecular mass of 43.4 kDa, which was confirmed by expression of tufB in maxicells. Sequence comparisons between S.aurantiaca EF-Tu and other bacterial homologues from E.coli, Salmonella typhimurium and Thermus thermophilus displayed extensive homologies (75.9%). Among the variable positions, two Cys residues probably involved in the temperature sensitivity of E.coli and S.typhimurium EF-Tu are replaced in T.thermophilus and S.aurantiaca EF-Tu. Since two or even three tuf genes have been described in other bacterial species, the presence of multiple tuf genes was sought for. Southern and Northern analysis are consistent with two tuf genes in the genome of S.aurantiaca. Primer extension experiments indicate that the four tRNA genes and tufB are organized in a single operon.

Light-induced carotenogenesis in Myxococcus xanthus: DNA sequence analysis of the carR region

The carR region encodes a light-inducible promoter, a negative regulator of the promoter and a trans-acting activator that controls the light-inducible Myxococcus xanthus carotenoid biosynthesis regulon. DNA sequence analysis revealed, downstream of the promoter, three translationally coupled genes, carQ, carR and carS. Sequencing of mutations demonstrated that carR encoded the negative regulator and was an integral membrane protein. Mutant construction and sequencing revealed that carS was the trans-acting activator and that carQ was a positive regulator of the promoter. Neither gene encodes proteins with known sequence-specific DNA-binding motifs. The sequence of the light-inducible promoter region, identified by primer extension analysis, showed similarity to the consensus sequence of the Escherichia coli stress response ('heat-shock') promoters.

Oar, a 115-kilodalton membrane protein required for development of Myxococcus xanthus

Myxococcus xanthus is a developmental gram-negative bacterium which forms multicellular fruiting bodies upon nutrient starvation. This bacterium was found to contain a 115-kDa membrane protein which separated with the inner membrane fraction by sucrose density gradient centrifugation. The gene for this protein was cloned, and its DNA sequence was determined. The deduced amino acid sequence consists of 1,061 residues. This protein contains a putative signal sequence and many short segments, found scattered throughout the entire protein, that have sequence similarities with OmpA, a major outer membrane protein of Escherichia coli. Thus, the gene was designated oar (OmpA-related protein). A second open reading frame was found 36 bases downstream of the oar termination codon. This open reading frame encodes a protein of 236 residues and contains a putative lipoprotein signal sequence. An aor disruption mutation (delta oar) showed no effect on vegetative growth but caused abnormal morphogenesis during development and reduced myxospore formation. When examined with a light microscope, delta oar cells were unable to aggregate on developmental agar, indicating that Oar is required for cellular adhesiveness during development.

Comparison of beta-galactosidase production by two inducible promoters in Myxococcus xanthus

The inducibility of two promoter systems, one heterologous and one homologous, has been assessed in the Gram-negative bacterium Myxococcus xanthus. The heterologous system involved the hybrid tac promoter and the presence of lacIq, the lac repressor from Escherichia coli. This system is inducible in its natural host with isopropyl-beta-D-thiogalactopyranoside (IPTG). The homologous promoter system involves the light-inducible carQRS promoter, which is normally involved in the expression of the regulators of the light-inducible light-protective carotenoid synthesis regulon in M. xanthus. In each case, promoter activity and strength was assayed using the E. coli gene lacZ. In our constructs, which were present in a single copy in the M. xanthus chromosome, the carQRS promoter yielded at least a 47-fold increase in beta-galactosidase production upon light induction, whilst IPTG increased by 8-fold the amount of enzyme produced under the control of the ptac-lacIq system. Regulation by the latter was significantly higher than that obtained with the unmodified lacZ promoter.

Cloning and DNA sequence of the gene coding for the major sigma factor from Myxococcus xanthus

The gene for a sigma factor (rpoD) was cloned from Myxococcus xanthus, a soil bacterium which differentiates to form fruiting bodies upon starvation for nutrients. The DNA sequence of the gene was determined, and an open reading frame encoding a polypeptide of 708 amino acid residues (Mr = 80,391) was identified. Except for the amino-terminal sequence consisting of 100 residues, the M. xanthus sigma factor (sigma-80) showed extensive similarity with Escherichia coli sigma-70 as well as Bacillus subtilis sigma-43. In particular, the carboxy-terminal sequence of 242 residues that is known to be required for promoter recognition and core recognition showed 78 and 72% amino acid sequence identity with the E. coli and B. subtilis sigma factors, respectively. The putative RpoD protein was detected at the position of an apparent molecular weight of 86,000 by Western blot (immunoblot) analysis by using antiserum against B. subtilis sigma-43, which agreed well with the position of a vegetative sigma factor of M. xanthus previously identified by Rudd and Zusman (K. Rudd and D. R. Zusman, J. Bacteriol. 151:89-105, 1982).

Nucleotide sequence and transcriptional products of the csg locus of Myxococcus xanthus

The csg locus of Myxococcus xanthus appears to control the production of an intercellular signal that is essential for development. The complete nucleotide sequence of a clone containing the csg locus was determined by the dideoxy-chain termination method. Pattern recognition analyses of the DNA sequence revealed the presence of two protein-coding regions that are convergently oriented and separated by only 8 nucleotides. Tn5 lac insertions into this clone detected two transcriptional units that are transcribed in a convergent fashion and whose expression increases during development. The two genes represented by these protein-coding regions and transcriptional units have been designated csgA and fprA. Northern (RNA) blot analyses detected an 800-nucleotide RNA specific to the csgA gene and a 900-nucleotide RNA specific to the fprA gene. Our results, along with mutational studies, identify csgA as the gene involved in cell communication. The function of the fprA gene is described in an accompanying paper (L. J. Shimkets, J. Bacteriol. 172:24-30, 1990).

Reverse transcriptase associated with the biosynthesis of the branched RNA-linked msDNA in Myxococcus xanthus

Myxobacteria have been shown to produce a peculiar RNA-DNA complex called msDNA, in which a single-stranded DNA is branched out from a RNA molecule (msdRNA) by a 2',5' phosphodiester linkage. It has been predicted that reverse transcriptase is required for msDNA biosynthesis. We identified a gene for reverse transcriptase in M. xanthus in the region that has been demonstrated to code for a cis- or transacting element for msDNA synthesis. This gene is located immediately downstream of the msdRNA coding region, and codes for a polypeptide of 485 amino acid residues. The polypeptide shows sequence similarity with retroviral reverse transcriptases. This fact, together with the mode of msDNA synthesis, suggests a possible relationship between retroviruses and the msDNA system. The analysis of the gene and the distribution of the msDNA system in independent isolates of M. xanthus indicate that the element is as old as other essential genes in M. xanthus and that it was not recently acquired into the genome.