ORIGINAL RESEARCH article

Front. Microbiol., 02 August 2019

Sec. Physiology and Metabolism of Microorganisms

Volume 10 - 2019 | https://doi.org/10.3389/fmicb.2019.01772

Sequence Determinants Spanning −10 Motif and Spacer Region Implicated in Unique Ehrlichia chaffeensis Sigma 32-Dependent Promoter Activity of dnaK Gene

  • Center of Excellence for Vector-Borne Diseases, Department of Diagnostic Medicine/Pathobiology, College of Veterinary Medicine, Kansas State University, Manhattan, KS, United States

Abstract

Ehrlichia chaffeensis is an obligate intracellular tick-borne bacterium that causes human monocytic ehrlichiosis. Studying Ehrlichia gene regulation is challenge, as this and related rickettsiales lack natural plasmids and mutagenesis experiments are of a limited scope. E. chaffeensis contains only two sigma factors, σ32 and σ70. We previously developed Escherichia coli surrogate system to study transcriptional regulation from RNA polymerase (RNAP) containing Ehrlichia σ32 or σ70. We reported that RNAP binding motifs of E. chaffeensis genes recognized by σ32 or σ70 share extensive homology and that transcription may be initiated by either one of the sigma factors, although transcriptional efficiencies differ. In the current study, we investigated mapping the E. chaffeensis dnaK gene promoter using the pathogen σ32 expressed in E. coli lacking its native σ32. The E. coli surrogate system and our previously described in vitro transcription system aided in defining the unique −10 motif and spacer sequence of the dnaK promoter. We also mapped σ32 amino acids/domains engaged in its promoter regulation in E. chaffeensis. The data reported in this study demonstrate that the −10 and −35 motifs and spacer sequence located between the two motifs of dnaK promoter are critical for the RNAP function. Further, we mapped the importance of all six nucleotide positions of the −10 motif and identified critical determinants within it. In addition, we reported that the lack of C-rich sequence upstream to the −10 motif is unique in driving the pathogen-specific transcription by its σ32 from dnaK gene promoter. This is the first study in defining an E. chaffeensis σ32-dependent promoter and it offers insights about how this and other related rickettsial pathogens regulate stress response genes.

Introduction

Ehrlichia chaffeensis, a Gram-negative, obligate intracellular tick-borne rickettsial pathogen, causes human monocytic ehrlichiosis (HME) (). HME as an emerging infectious disease, first reported in the United States in 1987, has become one of the most prevalent tick borne diseases in the United States and is also described from several other parts of the world (; Yabsley, 2010). Further, E. chaffeensis infects several other vertebrates, including dogs, goats, coyotes and white-tailed deer (; ; ; ; ; ). The pathogen infection in people may result in an acute flu-like illness with symptoms ranging from persistent fever, headache, myalgia, anorexia and chills (). HME patients may exhibit leukopenia, thrombocytopenia, anemia, and upgraded levels of serum hepatic aminotransferases (). In addition, other related Anaplasmataceae family pathogens included in the genera Ehrlichia and Anaplasma have been established as causative agents of emerging diseases in people and various vertebrate animals in recent years (; ; ). Though some progress is made in establishing genetics in E. chaffeensis (; ). and similarly in other related Ehrlichia and Anaplasma (; ; ; Wood et al., 2014; ), the genetic tool kit and its application is still limited. For example, it is not possible to investigate regulation of gene expression by transforming this group of important pathogens, possibly also because the pathogens lack naturally existing extrachromosomal plasmids. This major impediment limits the understanding of molecular mechanisms used by the pathogens in regulating gene expression in support of their continued survival in vertebrate and tick hosts and in causing pathogenesis (; ; ). Several prior studies reported differences in gene expression of E. chaffeensis impacted by different host environments (; ). However, it is unclear how the organism regulates its gene expression in support of its adaptation to the hosts.

Regulation of gene expression in bacteria is primarily controlled at the transcription. An RNA polymerase (RNAP) core enzyme with a sigma (σ) factor offers a simple and valid mechanism for bacteria to rapidly accommodate to diverse environmental changes by suitably modifying the transcriptional profiles (; ; ). Typically, an RNAP holoenzyme is a multi-subunit complex consists of a core enzyme containing two alpha (α), a beta (β), a beta′ (β′), and a omega (ω) subunits and then the inclusion of a σ factor (). A σ factor enables a core enzyme in specific binding to the promoter region of a gene for initiating transcription. Numbers of σ factors differ depending on the genome size variations and the environmental diversification of a bacterium (). For example, Escherichia coli possesses seven σ factors, while 109 σ factors are identified in Sorangium cellulosum (; ). Obligate intracellular bacteria generally tend to have reduced genomes and consequently also have fewer σ factors (). For example, E. chaffeensis genome of 1,176 kb has only two σ factor genes; rpoD (ECH_0760) (the primary housekeeping σ70 gene) and rpoH (ECH_0655) (the alternate σ32 gene) () (GenBank # NC_007799.1).

To study gene regulation in E. chaffeensis, we previously described in vitro transcription system and an E. coli surrogate system that is valuable in investigating gene regulation driven from its σ70 (; , ). Earlier, we also defined the promoters of several pathogen genes by utilizing in vitro transcription assays where E. coli RNAP core enzyme is reconstituted with the recombinant E. chaffeensis σ factors (; ). Our studies demonstrated that the RNAP binding motifs (−10 and −35 regions) of E. chaffeensis gene promoters share extensive homology and that they are recognizable by RNAP with either one of its only two sigma factors; σ32 or σ70, although affinities vary for different gene promoters (). In E. coli, gene expression of heat shock proteins; Dnak-DnaJ-GrpE and GroES-GroEL chaperone complexes, are controlled by its σ32 (). Similarly, we discovered that the E. chaffeensis chaperon protein gene (Ech_0471) encoding for DnaK protein is transcribed primarily by σ32 (). Genes regulated by σ32 are known to induce cellular responses under varieties of stresses confronted during bacterial growth and are likely important for E. chaffeensis survival in its hostile host environments and that they may contribute to pathogenicity similar to other Gram negative bacteria (; ; ; ; ). For example, DnaK in Vibrio cholerae, the gene expression regulated by its σ32, is involved in causing virulence in a host (; ; ).

Prior research on E. coli (a γ-proteobacteria) provides abundant knowledge regarding gene regulation from its σ32- and σ70-bound RNAP, while such knowledge for both σ32-and σ70 bound RNAP in other Gram-negative bacteria, particularly for α-proteobacteria, including for pathogenic organisms is very limited. Importantly, it is unclear how intracellular pathogens, such as E. chaffeensis, regulate gene expression to overcome the host stress and adapt to host environmental changes within its arthropod (tick) and vertebrate hosts. To extend knowledge on how E. chaffeensis regulates its gene expression, we continue investigations in defining the functions of its RNAP holoenzyme comprising σ32 or σ70. Transcription derived by an RNAP typically implicates in recognizing and binding to DNA sequence motifs of a promoter; −10 and −35 regions, and the spacer sequences located between the two motifs of a gene promoter (; ).

In the current study, we described the mapping of dnaK gene promoter recognized primarily by the E. chaffeensis RNAP containing σ32 using the previously developed E. coli surrogate system in the strain, CAG57101 (; ). In E. coli CAG57101, its endogenous rpoH gene (encoding for σ32) is inactivated () and in its place, we expressed the E. chaffeensis σ32 from a plasmid in defining dnaK promoter mapping, as direct gene mapping studies are not possible in this and other related intracellular rickettsials. We also investigated the functional domains of E. chaffeensis σ32 likely important for the RNAP function and in its interactions with the −10 motif and the spacer sequence of dnaK.

Materials and Methods

E. coli Strains and Plasmids

Escherichia coli strains used in this study were TOP10 (Invitrogen Technologies, Carlsbad, CA, United States), BL21(DE3) (Novagen, San Diego, CA, United States), and CAG57101 (). Several plasmid constructs used in this study were obtained from commercial sources or modified from one or more of the existing plasmids. They include the derivatives of pSAKT-Eco_rpoH (previously known as pSAKT32) (; ), pQF50K-Ech_dnaK () and pMT504 (). Genetic makeup of plasmids described in this study were included in Supplementary Table S1, except those obtained from a commercial source. The plasmid pSAKT-Eco_rpoH containing a p15A origin of replication and an ampicillin resistance gene has E. coli rpoH gene under the control of IPTG inducible Plac promoter (; ). The E. coli rpoH from this plasmid was removed by using Q5 Site-Directed Mutagenesis Kit (New England Biolab, Inc., Ipswich, MA, United States) and the modified plasmid was referred as the pSAKT. The pSAKT-Ech_rpoH (previously known as pSAKT32- Ech_rpoH) contained E. chaffeensis rpoH (). E. chaffeensis rpoH variants with substitutions within the 2.3, 2.4, and 3.0 regions of σ32 were constructed by mutagenesis using Q5 Site-Directed Mutagenesis Kit (New England Biolab, Inc., Ipswich, MA, United States). The names of the modified pSAKT-Ech_rpoH are provided in Supplementary Table S1.

The pQF50K-Ech_dnaK plasmid and the pQF50K-Ech_dnaK with deletion of −35 motif, which contains the β-galactosidase coding sequence (lacZ) with a pMB1 origin of replication and with a kanamycin resistance gene, was reported earlier (). The modified plasmid with deletion of −10 motif was generated from pQF50K-Ech_dnaK plasmid using Q5 Site-Directed Mutagenesis Kit (New England Biolab, Inc., Ipswich, MA, United States). Site directed mutagenesis constructs with mutations at every nucleotide of the −10 motif were also generated from the pQF50K-Ech_dnaK plasmid using Q5 Site-Directed Mutagenesis Kit. The mutants to modify the AT rich spacer sequence of the dnaK promoter were generated similarly by modifying the pQF50K-Ech_dnaK plasmid. The expression plasmids of E. chaffeensis wildtype σ32 was constructed and used for preparing purified recombinant proteins of σ32 as in early reports ().

For in vitro transcription analysis, pMT504-Ech_dnaK as transcription template was prepared and reported earlier (). Constructs with various mutations at −10 motif for the dnaK promoter for in vitro transcription assays were similarly prepared from this plasmid using Q5 Site-Directed Mutagenesis Kit. The lengths of transcripts for the various promoter segments of dnaK gene are 162 nucleotides. Integrity of all cloned segments in the plasmid constructs was confirmed by DNA sequence analysis using CEQ 8000 Genetic Analysis System (Beckman Coulter, Fullerton, CA, United States). The names of all engineered plasmids were listed in Supplementary Table S1. Mutagenic oligonucleotides were described in the Supplementary Table S2.

E. coli Growth Conditions and β-Galactosidase Assays

The E. coli strain CAG57101 transformed with the recombinant plasmids were grown as in early reports (; ). Briefly, cultures were grown at 30°C in Luria–Bertani (LB) medium with chloramphenicol (30 μg/ml) and spectinomycin (50 μg/ml) in support of the strain’s growth, and by ampicillin (100 μg/ml) for maintaining the pSAKT-derived plasmids. To assess the functions and impact of various mutations within the promoter regions of genes encoding dnaK, pQF50K-derived plasmid containing the promoter segments were also maintained by growing E. coli cultures with the addition of kanamycin (50 μg/ml). E. coli cultures of CAG57101 in LB medium were grown overnight with appropriate antibiotic supplements which were diluted 1:100 into a fresh medium containing appropriate antibiotics and the growth was continued for 2 h. Subsequently, cultures were then induced with 1 mM IPTG for 3 h before harvesting, when OD at 600 nm reached between 0.6 and 0.8. Lysates were prepared and used to measure β-galactosidase enzyme activity using a β-Gal Assay Kit (Invitrogen Technologies, Carlsbad, CA, United States). All experiments were performed three independent times with independently grown cultures; specific activity of β-galactosidase was calculated as outlined in the kit protocol.

In vitro Transcription Assays

In vitro transcription reactions were performed in 10 μl reaction mixture containing 0.1 picomoles each of the supercoiled plasmid DNA as the template and using RNAP holoenzyme containing recombinant E. chaffeensis σ32 (). The holoenzyme was prepared by mixing 0.5 μl of 1:10 diluted stock of E. coli core enzyme (Epicentre, Madison, WI, United States) mixed with 10-fold molar excess of purified recombinant E. chaffeensis σ32 and kept in ice for 30 min prior to using for the reactions. The transcription reactions were performed at 37°C for 20 min, and the reactions were terminated by adding 7 μl of stop solution (95% formamide, 20 mM EDTA, 0.05% bromophenol blue and 0.05% xylene cyanol). Six microliters each of the samples were resolved on a 6% polyacrylamide sequencing gel with 7 M urea, then gels were transferred to a Whatman paper, dried and 162 nucleotide transcripts were visualized by exposing an X-ray film to the gels. The transcripts were quantified using ImageJ software1.

Bioinformatics

Multiple DNA alignments were done using the programs Clustal X version 2.0 with default parameters ().

Statistical Analysis

Statistical analyses were performed using Student’s t-test, and a P-value < 0.05 was considered significant with a single asterisk.

Results

The −10 Motif Is Needed for the E. chaffeensis dnak Gene Transcription by RNAP Holoenzyme Constituting of Its σ32

Our prior studies demonstrated that −35 motif, but not −10 motif, is required for the σ70-bound RNAP holoenzyme transcription from σ70-dependent promoters in E. chaffeensis (). We also reported earlier that −35 motif is similarly required for the dnaK gene transcription by E. chaffeensis σ32–bound RNAP holoenzyme (). To test whether or not the −10 motif of dnaK promoter is required for the σ32-dependent gene regulation, plasmid constructs lacking −10 or −35 motifs of the promoter were cloned upstream to the β-galactosidase coding sequence in an E. coli mutant deficient for its σ32 expression that is functionally complemented with the E. chaffeensis σ32 (Figure 1). The −10 motif deletion and similarly the −35 motif deletion resulted in a significant reduction of β-galactosidase activity (to 11 and 21%, respectively; p ≤ 0.0005) compared to the dnaK wildtype (WT). The reduction of the β-galactosidase activity for the −10 motif deletion was similar to the negative control where the promoter segment was absent (NP).

FIGURE 1

Identifying the Critical Sequence Determinants of −10 Motif of E. chaffeensis dnaK

The consensus sequence of −35 motif for E. chaffeensis promoters recognized by its σ32 (TTGTAT) is similar to the consensus −35 motif of σ32-dependent promoters in E. coli (TTGAAA) and similarly it shares extensive homology to −35 motif for the genes recognized by σ70 (TTGNTT) (; ). The consensus −10 motif of E. chaffeensis promoters recognized by its σ32 (TATATN) is also similar to its consensus −10 motif recognized by its σ70 (TATTNT), however, it differs significantly from the consensus −10 motif of E. coli σ32-dependent promoters (CCCCATWT) (; ; ). While deletion of −10 motif from σ70-dependent promoters has no impact on promoter activities in E. chaffeensis (), such deletion from the σ32-dependent dnaK resulted in significant and 90% reduction in the promoter activity (Figure 1). These novel data suggest that, contrary to σ70-bound RNAP (), the −10 motif plays a critical role for σ32-bound RNAP in E. chaffeensis. We therefore performed detailed point mutation experiments to define the critical sequence determinants of the −10 motif for the dnaK promoter activity (Figure 2). We have made substitution mutations at each base of the six-nucleotide motif (TATATC) and evaluated the impact of each mutation by measuring changes in β-galactosidase expression in CAG57101 E. coli functionally complemented with E. chaffeensis σ32. A specific substitution mutation was indicated by combination of letters and numbers. For example, T1A indicates a change from T to A transversion at the first position in the −10 motif. One or more substitutions at all six nucleotide positions resulted in significant decline in the promoter activity of dnaK. Substitutions at the first five nucleotides to any other nucleotide resulted in significant promoter activity drop. Substitutions in the first position from T1G or T1C resulted in the near complete loss of promoter activity, while T1A resulted in the reduction of promoter activity to 39%. Similarly, in the second position, A2C and A2T mutations caused significant loss of promoter activity (reduced to 17% and 7% compared to the wildtype), whereas A2G mutation caused decline of two thirds of activity similar to T1A substitution. In the third and fourth positions, changes to any other nucleotide had the greatest impact in promoter activity decline (78–98% reduction). Mutations in the fifth position to T5G or T5C had a greater impact (a near 90% decline), while reduction was less apparent for the T5A that is similar to the mutations in the first position T. In the sixth position, only C6G mutation resulted in significant decline in the promoter activity to 35%.

FIGURE 2

In vitro Transcription Assays to Verify Sequence Determinants of −10 Motif Mapped in E. coli CAG57101

To validate the results in defining the −10 motif in the E. coli surrogate system, we performed in vitro transcription assays using several randomly selected promoter mutation constructs; the assays were performed using RNAP holoenzyme reconstituted with the recombinant E. chaffeensis σ32. We selected five dnaK −10 motif mutants for this experiment and compared the results with the wildtype promoter. Both wildtype and mutated versions of dnaK promoter segments were cloned into the G-less cassette and used as templates in the in vitro transcription assays (Figures 3A,B). Consistent with the results recorded with the E. coli CAG57101 system, mutants TIA, T1G and A2G produced lesser transcripts compared to the wildtype dnaK. There was no significant difference for the T5A mutation compared to the wildtype, which is also consistent with the results observed in the E. coli system. In the sixth position, C6T mutation caused an increase of in vitro transcript level, which is also similar to the enhanced promoter activity observed in the E. coli CAG57101experiments.

FIGURE 3

The Spacer Sequences Affect Promoter Activity

Previous studies in E. coli demonstrate that spacer sequences located between −10 and −35 motifs contribute to promoter activities (; ; , ; ). In particular, nucleotides present in the spacer sequence and its length play critical role for a promoter activity. Further, a short C-rich region upstream to −10 motif in E. coli and in other γ-proteobacteria is identified as important extended −10 motif required for efficient transcription by σ32-containing RNAP holoenzyme (; ; ; ), while such C-rich sequence does not exist in E. chaffeensis gene promoters, including in the dnaK promoter. Our previous studies for σ70-dependent promoter genes suggested that changes to the spacer sequence impact a promoter activity (). We, therefore, investigated the importance of dnaK spacer sequence, including in determining about how the lack of C-rich sequence impact the promoter activity. Nine spacer mutants were prepared where nucleotides within the spacer sequence were modified; they included replacing the spacer sequence with its complementary sequence (CP), or with a high GC content spacer while keeping the spacer length constant (GC), or by increasing the spacer sequence lengths from 17 bp to 18, 19, or 20 bps or by decreasing it to 16, 15, or 14 bp and finally by deleting the spacer sequence completely. These different spacer mutant constructs were depicted in Figure 4A. The β-galactosidase expression was then assessed for all these modified spacer promoter segments and compared to wildtype (WT) promoter construct in the E. coli surrogate system (CAG57101) (Figure 4B). The CP mutant caused a minor, non-significant increase in the promoter activity (30% increase). The spacer substitution with GC resulted in a significant 50% decline in the β-galactosidase expression. Increasing the spacer length to 18 bp or decreasing to 14 or 15 bp caused a major decline in the promoter activity, although the greatest decline was observed with the 14 bp spacer (96% drop), while reducing the length to 16 bp had no impact. Increasing the spacer lengths to 19 or 20 bp resulted in much higher enhancement of the promoter activity (586% and 249%, respectively). Deletion of the entire spacer sequence had no impact on the promoter activity compared to the WT promoter. We detected the presence of another 15 bp spacer like sequence and an alternative −35 motif sequence in the complete deletion spacer construct; thus, it is highly likely that these sequences served as alternate spacer and −35 motif for the RNAP (Supplementary Figure S1).

FIGURE 4

Substitution Mutations in E. chaffeensis σ32

Prior studies in E. coli revealed that the −10 motif of σ32-dependent promoters is recognized by few amino acids within the regions 2.3, 2.4, and 3.0 of σ32 protein (; ). To determine what amino acids spanning these regions of E. chaffeensis σ32 would contribute to the promoter activity, we have made substitution mutations at 6 different amino acid positions likely alter the functional domains of the pathogen σ32. The amino acids for substitution mutations were selected based on their homology with the E. coli σ32 (Supplementary Figure S2). Tryptophan (W) 108 of E. coli σ32 within the region 2.4 is identified as important for recognition of −13C in a promoter recognized by it (; ). In E. chaffeensis σ32, W is located at 106. Glutamic acid (E) at position 112 in E. coli σ32 is also implicated for its contribution to −13C recognition and interaction with its σ32-dependent promoters (). Phenylalanine (F) at position 110 of E. chaffeensis σ32 is in the homologous position to E112 in E. coli σ32. As E. chaffeensis dnaK lacks C-rich region upstream to its −10 motif, we expected that amino acid substitution mutations at these two positions in E. chaffeensis σ32 would not have any impact. Lysine (K) at position 130 in region 3.0 of E. coli σ32 is similarly implicated for the recognition of the C-rich region for the σ32-promoters; groE and grpE (). Glutamin (Q) at position 128 in region 3.0 of E. chaffeensis σ32 is the amino acid at the homologous position to K130 of E. coli σ32. Substitution of W106 to A caused significant reduction of the β-galactosidase expression (80% decline), while F110 substitution to A or E (alanine or glutamic acid) in E. chaffeensis σ32 did not significantly alter the promoter function (Figure 5). Similarly, substitution of Q128 to A did not significantly impact E. chaffeensis dnaK promoter activity. The mutational data with the exception of W106 are consistent with the lack of C-rich region in E. chaffeensis. However, W106 may be critical for the promoter activity independent of the C-rich sequence, at least in E. chaffeensis. Phenylalanine (F) at position 104 within the 2.3 region of E. coli σ32 is identified as critical for its structural integrity and activity of σ32 (). In E. chaffeensis σ32, a polar amino acid (at position 102) {tyrosine (Y)} is present at the position homologous to F104. Similarly, A111 of E. coli σ32 is implicated for its binding to core RNAP (). Y102 to A caused a significant decline of E. chaffeensis σ32 activity (88% reduction). Likewise, substitution of A109 to glutamine (Q) in E. chaffeensis σ32 that is homologous to A111of E. coli σ32 resulted in a significant decline of activity (about 31% of WT level). Previous studies in E. coli demonstrate that substitution of F136A within the region 3.0 reduces the interaction between core RNAP and σ32 thus leading to 80% decline in promoter activity for its groE gene (; ; ). A similar substitution in E. chaffeensis σ32; F134A also caused similar decline of its activity for the dnaK promoter (reduced to 28% activity compared to WT E. chaffeensis σ32). Together, E. chaffeensis σ32 substitution mutation experiments allowed the identification of critical functional domains engaged in σ32-bound RNAP interactions with the dnaK promoter, including in confirming that the C-rich region is not critical for its function.

FIGURE 5

Differences in the −10 Motifs Between E. coli and E. chaffeensis Are Sufficient in Having Differential σ32 Functions

Unlike E. coli and other γ- proteobacterial gene promoters, −10 motif of E. chaffeensis σ32-dependent dnaK promoter (TATATN) is distinct in lacking a C-rich sequence upstream to −10 motif (; ; ). Results presented in the previous section suggest that the C-rich region is indeed not required for E. chaffeensis σ32 dependent dnaK promoter function. To further map how the variations in −10 motif and spacer sequence in E. chaffeensis dnaK make it unique for this intracellular pathogen promoter function, we prepared two modified constructs where two or four nucleotides spanning between the spacer sequence and −10 motif (TT or TATT, respectively) were replaced with either two Cs or four Cs (Figure 6A). Importantly, these sequence modifications change the E. chaffeensis −10 motif to be more similar to the E. coli consensus −10 motif (CCCCATWT) (Figure 6A). The modified constructs having 2Cs and 4Cs were then assessed in the E. coli surrogate system expressing either E. chaffeensis σ32 (Ech- σ32) or E. coli σ32 (Eco-σ32). Compared to the wildtype dnaK promoter, both the 2C and 4C mutants caused drastic decline in the β-galactosidase expression by as much as 15-fold for E. chaffeensis σ32 (Figure 6B). On the contrary, while wildtype dnaK promoter had a minimal β-galactosidase expression with the E. coli σ32, the 2Cs and 4Cs substitutions caused a significant and a step-wise increase of β-galactosidase expression to 2.8-fold and 14.5-fold, respectively. Together, these results suggest that the inclusion of C-rich sequences is sufficient in altering the promoter specificities of E. chaffeensis to be similar to E. coli σ32-dependent RNAP.

FIGURE 6

Discussion

Because obligate intracellular bacteria reside within the cytosol or a phagosome of an infected host cell, they encounter minimal environmental changes, possibly compared to free-living bacteria, such as E. coli. Obligate intracellular bacteria typically have reduced genomes and consequently their genomes also have limited numbers of sigma factors. For example, E. chaffeensis genome has only two sigma factor genes encoding for a constitutive and an alternative sigma factor, σ70 and σ32, respectively (). On the contrary, a non-pathogenic bacterium, such as E. coli, contains seven sigma factors, including σ70 and σ32 (; ). Alternate sigma factor; σ32 is known to regulate genes involved in overcoming environmental stresses (Zhao et al., 2005; ). Prior studies in E. coli, V. cholerae, Neisseria gonorrhoeae, and Francisella tularensis suggest that σ32 regulon involves several heat-shock proteins, also known as chaperon proteins; Hsp40, HspG, Dnak, DnaJ, GroES, and GroEL. These proteins are involved in temperature-induced stress control of the organisms (Zhao et al., 2005; ; ; ). These proteins play critical roles in both preventing denaturation of proteins and also to promote renaturation of aggregated proteins (; ; ). Some of the proteins involved in regulating stress response, such as ClpB, HspG, and DnaK, are also considered important for the intracellular survival of a pathogenic bacterium, F. tularensis (; ; ; ). The sigma factor, σ32, is identified as critical for promoting transcription of genes engaged in overcoming stressful environments for bacteria to promote the reduction of accumulation of misfolded and aggregated proteins (; ; ).

Ehrlichia chaffeensis and related Anaplasmataceae family pathogens have the ability adapt to vertebrate and tick hosts. However, it is unclear how tick transmitted pathogens regulate gene expression in tick and vertebrate hosts. Studying gene regulation will aid in understanding how rickettsial pathogens adapt to dual hosts and sense nutrient, starving, temperature, and other stressful environments within an infected host cell. Previous studies revealed that E. chaffeensis has two morphological forms, dense core cell as the pathogen’s infectious form and reticulate cell, which replicates within a phagosome of an infected host cell (Zhang et al., 2007; ). We recently reported that a stress response protein, ClpB, transcripts are higher during the replicative stage of the pathogen, while gene expression of another heat shock protein, DnaK, and the stress response sigma factor, RpoH, remained as constitutively expressed throughout the replicating stage (; Zhang et al., 2013). Further, we also reported that E. chaffeensis sigma factors, σ32 and σ70, function cooperatively in transcribing pathogen genes, but with varying affinities (). The dnaK has higher affinity for RNAP containing E. chaffeensis σ32 compared to σ70 (). The current study is the first in defining E. chaffeensis σ32-dependent gene promoter region of dnaK.

Despite our recent advances in establishing tractable genetics for E. chaffeensis (), studying gene regulation remains a challenge due to lack of a well-established methods to maintain extracellular plasmids, as the pathogen and other related Anaplasmataceae pathogens lack natural plasmids. To overcome this challenge, in the current study, we utilized the E. coli surrogate system to map promoter-binding domains required for gene expression in E. chaffeensis for a σ32-dependent gene promoter. We validated the results from the E. coli surrogate system using in vitro transcription assays which we developed earlier (, ). In vitro transcription system is applied well in understanding intracellular bacterial gene regulation for pathogens, such as Chlamydia trachomatis, for which a tractable genetic system is equally not well developed, which can support gene regulation studies (; ; ; Yu and Tan, 2003; ; ; ; , ). Hence, in view of the technical challenges, methods described in the current study are innovative in mapping E. chaffeensis gene regulation, and that the study will lead the way for similar investigations in other related Anaplasmataceae family pathogens.

RNA polymerase holoenzyme containing σ70 of E. coli transcribes housekeeping genes by recognizing two highly conserved motifs of a gene promoter; referred as −10 and −35 motifs (). We reported earlier that the consensus −10 and −35 motifs for E. chaffeensis σ70-dependent promoters (TATTNT and TTGNTT, respectively) () are similar to E. coli consensus −10 and −35 (TATAAT and TTGACA, respectively) (; ). We reported that the −35 motif and the AT-rich spacer sequences are important for a gene-specific regulation by σ70-dependent promoters (). The predicted E. chaffeensis consensus −35 motif for the alternative sigma factor (σ32) is also similar to E. coli consensus −35 motif (; ,; ). The −10 motif of E. chaffeensis σ32-dependent promoters (TATATN), however, differs substantially from the consensus E. coli−10 motif (CCCCATNT) (; ,; ). The consensus −35 motif is also extensively conserved among σ32-dependent promoters in proteobacteria (Supplementary Table S3), while −10 motif of σ32-dependent promoters vary considerably among various classes of proteobacteria. For example, many γ-proteobacteria contain a C-rich sequence upstream to the −10 motif and is implicated in regulating σ32-dependent gene regulation (; ; ; ). Such sequence, however, is absent in E. chaffeensis σ32-dependent gene promoters (). Thus, we hypothesized that the pathogen is evolved to possess a distinctive −10 motif and that it plays a critical role in σ32-dependent promoter regulation. Previous studies suggest that the length of a spacer sequence between the −10 and −35 motifs and the composition of specific nucleotides present within a spacer affect promoter activity (; ; ; ; ). In particular, spacer sequences influence the three-dimensional structure of a promoter and any modifications to a spacer sequence, therefore, affect the efficiency of interactions between the RNAP holoenzyme with a promoter sequence either positively or negatively (; ; ). Consistent with prior studies, our current investigation similarly revealed the importance of spacer sequence in contributing to the promoter activity for the E. chaffeensis dnaK gene transcribed by its σ32. We previously reported that modifications to spacer sequences similarly impact E. chaffeensis genes transcribed by σ70 genes and that changes to spacer sequences alter the curvature of a promoter region ().

In summary, we mapped E. chaffeensis dnaK−10 and −35 motifs and the spacer sequence upstream to it by performing detailed mutational analysis. Furthermore, this study confirmed that the C-rich region–specific interactions between a gene promoter and σ32 protein domains, known to be critical for E. coli, are absent in E. chaffeensis, which makes the pathogen gene regulation distinct, as judged from characterizing the dnaK gene promoter.

Statements

Data availability statement

The datasets generated for this study are available on request to the corresponding author.

Author contributions

RG conceived and directed the research design. HL contributed to the research plan and performed the experiments. HL and RG evaluated the data and prepared the manuscript.

Funding

This work was supported by the PHS grant # AI070908 from the National Institute of Allergy and Infectious Diseases, National Institutes of Health, United States. This manuscript is a contribution from the Kansas Agricultural Experiment Station (Contribution Number 19-291-J). The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgments

We thank Ms. Mal Rooks Hoover for her help in preparing the figures.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.01772/full#supplementary-material

References

Summary

Keywords

gene regulation, intracellular bacteria, Anaplasmataceae, sigma factor, Ehrlichia chaffeensis, RNA polymerase

Citation

Liu H and Ganta RR (2019) Sequence Determinants Spanning −10 Motif and Spacer Region Implicated in Unique Ehrlichia chaffeensis Sigma 32-Dependent Promoter Activity of dnaK Gene. Front. Microbiol. 10:1772. doi: 10.3389/fmicb.2019.01772

Received

10 May 2019

Accepted

18 July 2019

Published

02 August 2019

Volume

10 - 2019

Edited by

Hari S. Misra, Bhabha Atomic Research Centre (BARC), India

Reviewed by

Jere W. McBride, The University of Texas Medical Branch at Galveston, United States; Bing Zhu, The University of Texas Medical Branch at Galveston, United States; Tian Luo, The University of Texas Medical Branch at Galveston, United States

Updates

Copyright

*Correspondence: Roman R. Ganta,

This article was submitted to Microbial Physiology and Metabolism, a section of the journal Frontiers in Microbiology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics