Analysis of the Genomics and Mouse Virulence of an Emergent Clone of Streptococcus dysgalactiae Subspecies equisimilis

Streptococcus dysgalactiae subsp. equisimilis is a bacterial pathogen that is increasingly recognized as a cause of severe human infections. Much less is known about the genomics and infection pathogenesis of S. dysgalactiae subsp. equisimilis strains compared to the closely related bacterium Streptococcus pyogenes. To address these knowledge deficits, we sequenced to closure the genomes of seven S. dysgalactiae subsp. equisimilis human isolates, including six that were emm type stG62647. Recently, for unknown reasons, strains of this emm type have emerged and caused an increasing number of severe human infections in several countries. The genomes of these seven strains vary between 2.15 and 2.21 Mbp. The core chromosomes of these six S. dysgalactiae subsp. equisimilis stG62647 strains are closely related, differing on average by only 495 single-nucleotide polymorphisms, consistent with a recent descent from a common progenitor. The largest source of genetic diversity among these seven isolates is differences in putative mobile genetic elements, both chromosomal and extrachromosomal. Consistent with the epidemiological observations of increased frequency and severity of infections, both stG62647 strains studied were significantly more virulent than a strain of emm type stC74a in a mouse model of necrotizing myositis, as assessed by bacterial CFU burden, lesion size, and survival curves. Taken together, our genomic and pathogenesis data show the strains of emm type stG62647 we studied are closely genetically related and have enhanced virulence in a mouse model of severe invasive disease. Our findings underscore the need for expanded study of the genomics and molecular pathogenesis of S. dysgalactiae subsp. equisimilis strains causing human infections. IMPORTANCE Our studies addressed a critical knowledge gap in understanding the genomics and virulence of the bacterial pathogen Streptococcus dysgalactiae subsp. equisimilis. S. dysgalactiae subsp. equisimilis strains are responsible for a recent increase in severe human infections in some countries. We determined that certain S. dysgalactiae subsp. equisimilis strains are genetically descended from a common ancestor and that these strains can cause severe infections in a mouse model of necrotizing myositis. Our findings highlight the need for expanded studies on the genomics and pathogenic mechanisms of this understudied subspecies of the Streptococcus family.

Genetic heterogeneity of the Spy1336/R28-Spy1337 virulence axis in Streptococcus pyogenes and effect on gene transcript levels and pathogenesis

Streptococcus pyogenes is a strict human pathogen responsible for more than 700 million infections annually worldwide. Strains of serotype M28 S. pyogenes are typically among the five more abundant types causing invasive infections and pharyngitis in adults and children. Type M28 strains also have an unusual propensity to cause puerperal sepsis and neonatal disease. We recently discovered that a one-nucleotide indel in an intergenic homopolymeric tract located between genes Spy1336/R28 and Spy1337 altered virulence in a mouse model of infection. In the present study, we analyzed size variation in this homopolymeric tract and determined the extent of heterogeneity in the number of tandemly-repeated 79-amino acid domains in the coding region of Spy1336/R28 in large samples of strains recovered from humans with invasive infections. Both repeat sequence elements are highly polymorphic in natural populations of M28 strains. Variation in the homopolymeric tract results in (i) changes in transcript levels of Spy1336/R28 and Spy1337 in vitro, (ii) differences in virulence in a mouse model of necrotizing myositis, and (iii) global transcriptome changes as shown by RNAseq analysis of isogenic mutant strains. Variation in the number of tandem repeats in the coding sequence of Spy1336/R28 is responsible for size variation of R28 protein in natural populations. Isogenic mutant strains in which genes encoding R28 or transcriptional regulator Spy1337 are inactivated are significantly less virulent in a nonhuman primate model of necrotizing myositis. Our findings provide impetus for additional studies addressing the role of R28 and Spy1337 variation in pathogen-host interactions.

New Pathogenesis Mechanisms and Translational Leads Identified by Multidimensional Analysis of Necrotizing Myositis in Primates

A fundamental goal of contemporary biomedical research is to understand the molecular basis of disease pathogenesis and exploit this information to develop targeted and more-effective therapies. Necrotizing myositis caused by the bacterial pathogen Streptococcus pyogenes is a devastating human infection with a high mortality rate and few successful therapeutic options. We used dual transcriptome sequencing (RNA-seq) to analyze the transcriptomes of S. pyogenes and host skeletal muscle recovered contemporaneously from infected nonhuman primates. The in vivo bacterial transcriptome was strikingly remodeled compared to organisms grown in vitro, with significant upregulation of genes contributing to virulence and altered regulation of metabolic genes. The transcriptome of muscle tissue from infected nonhuman primates (NHPs) differed significantly from that of mock-infected animals, due in part to substantial changes in genes contributing to inflammation and host defense processes. We discovered significant positive correlations between group A streptococcus (GAS) virulence factor transcripts and genes involved in the host immune response and inflammation. We also discovered significant correlations between the magnitude of bacterial virulence gene expression in vivo and pathogen fitness, as assessed by previously conducted genome-wide transposon-directed insertion site sequencing (TraDIS). By integrating the bacterial RNA-seq data with the fitness data generated by TraDIS, we discovered five new pathogen genes, namely, S. pyogenes 0281 (Spy0281 [dahA]), ihk-irr, slr, isp, and ciaH, that contribute to necrotizing myositis and confirmed these findings using isogenic deletion-mutant strains. Taken together, our study results provide rich new information about the molecular events occurring in severe invasive infection of primate skeletal muscle that has extensive translational research implications.IMPORTANCE Necrotizing myositis caused by Streptococcus pyogenes has high morbidity and mortality rates and relatively few successful therapeutic options. In addition, there is no licensed human S. pyogenes vaccine. To gain enhanced understanding of the molecular basis of this infection, we employed a multidimensional analysis strategy that included dual RNA-seq and other data derived from experimental infection of nonhuman primates. The data were used to target five streptococcal genes for pathogenesis research, resulting in the unambiguous demonstration that these genes contribute to pathogen-host molecular interactions in necrotizing infections. We exploited fitness data derived from a recently conducted genome-wide transposon mutagenesis study to discover significant correlation between the magnitude of bacterial virulence gene expression in vivo and pathogen fitness. Collectively, our findings have significant implications for translational research, potentially including vaccine efforts.

Polymorphisms in Regulator of Cov Contribute to the Molecular Pathogenesis of Serotype M28 Group A Streptococcus

Two-component systems (TCSs) are signal transduction proteins that enable bacteria to respond to external stimuli by altering the global transcriptome. Accessory proteins interact with TCSs to fine-tune their activity. In group A Streptococcus (GAS), regulator of Cov (RocA) is an accessory protein that functions with the control of virulence regulator/sensor TCS, which regulates approximately 15% of the GAS transcriptome. Whole-genome sequencing analysis of serotype M28 GAS strains collected from invasive infections in humans identified a higher number of missense (amino acid-altering) and nonsense (protein-truncating) polymorphisms in rocA than expected. We hypothesized that polymorphisms in RocA alter the global transcriptome and virulence of serotype M28 GAS. We used naturally occurring clinical isolates with rocA polymorphisms (n = 48), an isogenic rocA deletion mutant strain, and five isogenic rocA polymorphism mutant strains to perform genome-wide transcript analysis (RNA sequencing), in vitro virulence factor assays, and mouse and nonhuman primate pathogenesis studies to test this hypothesis. Results demonstrated that polymorphisms in rocA result in either a subtle transcriptome change, causing a wild-type-like virulence phenotype, or a substantial transcriptome change, leading to a significantly increased virulence phenotype. Each polymorphism had a unique effect on the global GAS transcriptome. Taken together, our data show that naturally occurring polymorphisms in one gene encoding an accessory protein can significantly alter the global transcriptome and virulence phenotype of GAS, an important human pathogen.

Gene fitness landscape of group A streptococcus during necrotizing myositis

Necrotizing fasciitis and myositis are devastating infections characterized by high mortality. Group A streptococcus (GAS) is a common cause of these infections, but the molecular pathogenesis is poorly understood. We report a genome-wide analysis using serotype M1 and M28 strains that identified GAS genes contributing to necrotizing myositis in nonhuman primates (NHP), a clinically relevant model. Using transposon-directed insertion-site sequencing (TraDIS), we identified 126 and 116 GAS genes required for infection by serotype M1 and M28 organisms, respectively. For both M1 and M28 strains, more than 25% of the GAS genes required for necrotizing myositis encode known or putative transporters. Thirteen GAS transporters contributed to both M1 and M28 strain fitness in NHP myositis, including putative importers for amino acids, carbohydrates, and vitamins and exporters for toxins, quorum-sensing peptides, and uncharacterized molecules. Targeted deletion of genes encoding 5 transporters confirmed that each isogenic mutant strain was significantly (P < 0.05) impaired in causing necrotizing myositis in NHPs. Quantitative reverse-transcriptase PCR (qRT-PCR) analysis showed that these 5 genes are expressed in infected NHP and human skeletal muscle. Certain substrate-binding lipoproteins of these transporters, such as Spy0271 and Spy1728, were previously documented to be surface exposed, suggesting that our findings have translational research implications.

Postgenomic adventures with Rhodobacter sphaeroides

This review describes some of the recent highlights taken from the studies of Rhodobacter sphaeroides 2.4.1. The review is not intended to be comprehensive, but to reflect the bias of the authors as to how the availability of a sequenced and annotated genome, a gene-chip, and proteomic profile as well as comparative genomic analyses can direct the progress of future research in this system.

Combining microarray and genomic data to predict DNA binding motifs

The ability to detect regulatory elements within genome sequences is important in understanding how gene expression is controlled in biological systems. In this work, microarray data analysis is combined with genome sequence analysis to predict DNA sequences in the photosynthetic bacterium Rhodobacter sphaeroides that bind the regulators PrrA, PpsR and FnrL. These predictions were made by using hierarchical clustering to detect genes that share similar expression patterns. The DNA sequences upstream of these genes were then searched for possible transcription factor recognition motifs that may be involved in their co-regulation. The approach used promises to be widely applicable for the prediction of cis-acting DNA binding elements. Using this method the authors were independently able to detect and extend the previously described consensus sequences that have been suggested to bind FnrL and PpsR. In addition, sequences that may be recognized by the global regulator PrrA were predicted. The results support the earlier suggestions that the DNA binding sequence of PrrA may have a variable-sized gap between its conserved block elements. Using the predicted DNA binding sequences, a whole-genome-scale analysis was performed to determine the relative importance of the interplay between the three regulators PpsR, FnrL and PrrA. Results of this analysis showed that, compared to the regulation by PpsR and FnrL, a much larger number of genes are candidates to be regulated by PrrA. The study demonstrates by example that integration of multiple data types can be a powerful approach for inferring transcriptional regulatory patterns in microbial systems, and it allowed the detection of photosynthesis-related regulatory patterns in R. sphaeroides.

Interacting regulatory circuits involved in orderly control of photosynthesis gene expression in Rhodobacter sphaeroides 2.4.1

FnrL, the homolog of the global anaerobic regulator Fnr, is required for the induction of the photosynthetic apparatus in Rhodobacter sphaeroides 2.4.1. Thus, the precise role of FnrL in photosynthesis (PS) gene expression and its interaction(s) with other regulators of PS gene expression are of considerable importance to our understanding of the regulatory circuitry governing spectral complex formation. Using a CcoP and FnrL double mutant strain, we obtained results which suggested that FnrL is not involved in the transduction of the inhibitory signal, by which PS gene expression is "silenced," emanating from the cbb(3) oxidase encoded by the ccoNOQP operon under aerobic conditions. The dominant effect of the ccoP mutation in the FnrL mutant strain with respect to spectral complex formation under aerobic conditions and restoration of a PS-positive phenotype suggested that inactivation of the cbb(3) oxidase to some extent bypasses the requirement for FnrL in the formation of spectral complexes. Additional analyses revealed that anaerobic induction of the bchE, hemN, and hemZ genes, which are involved in the tetrapyrrole biosynthetic pathways, requires FnrL. Thus, FnrL appears to be involved at multiple loci involved in the regulation of PS gene expression. Additionally, bchE was also shown to be regulated by the PrrBA two-component system, in conjunction with hemN and hemZ. These and other results to be discussed permit us to more accurately describe the role of FnrL as well as the interactions between the FnrL, PrrBA, and other regulatory circuits in the regulation of PS gene expression.

From redox flow to gene regulation: role of the PrrC protein of Rhodobacter sphaeroides 2.4.1

Activation of photosynthesis (PS) gene expression by the PrrBA two-component activation system in Rhodobacter sphaeroides 2.4.1 results from the interruption of an inhibitory signal originating from the cbb(3) cytochrome c oxidase via its interaction with oxygen, in conjunction with the Rdx redox proteins. The CcoQ protein, encoded by the ccoNOQP operon, which encodes the cbb(3) cytochrome c oxidase, was shown to act as a "transponder" that conveys the signal derived from reductant flow through cbb(3) to oxygen, to the Prr system. To further define the elements comprising this signal transduction pathway we considered the prrC gene product, which to date possessed no definable role in this signal transduction pathway despite its being part of the prrBCA gene cluster. Similar to mutations in cbb(3) and rdx, suitably constructed prrC deletion mutations lead to PS gene expression in the presence of high oxygen. Unlike mutations that remove cbb(3) terminal oxidase activity or Rdx function, the PrrC deletion mutant shows no effect upon cbb(3) activity, nor does it affect the ratio of the carotenoid (Crt) spheroidene (SE) to spheroidenone (SO). Thus, the PrrC deletion mutant behaves identically to the CcoQ deletion mutant. Taking these and previous results together, we suggest that PrrC is located upstream of the two-component PrrBA activation system in the signal transduction pathway but downstream of the cbb(3) cytochrome c oxidase and its "transponder" CcoQ. The PrrC deletion mutant was also shown to lead to an increase in the DorA protein under aerobic conditions as was shown earlier for the cbb(3) mutant. Finally, PrrC is a member of a highly conserved family of proteins found in both prokaryotes and eukaryotes, and this appears to be the first instance in which a direct regulatory role has been ascribed to a member of this protein family.

A redox-responsive pathway for aerobic regulation of photosynthesis gene expression in Rhodobacter sphaeroides 2.4.1

To further understand the proposed signal transduction pathway involving the presumed redox proteins RdxBH and cbb3 cytochrome oxidase in Rhodobacter sphaeroides 2.4.1, a series of mutants lacking components of both the Prr two-component activation system and the cbb3-type cytochrome oxidase or RdxBH were constructed. We report that under highly aerobic conditions, aberrant photosynthesis gene expression and spectral complex formation typical of cbb3- or RdxBH-deficient mutants were no longer observed when either prrA (encoding the response regulator of the Prr system) or prrB (encoding the presumed sensor kinase) was also deleted. These double-mutant strains are phenotypically identical to single-mutant PrrA and PrrB strains, suggesting that the signal(s) originating from the cbb3 terminal oxidase affects downstream puc and puf operon expression by acting exclusively through the Prr system. When the same double-mutant strains were examined under anaerobic dark dimethyl sulfoxide growth conditions, photosynthesis gene expression was obligatorily linked to the two-component activation system. However, photosynthesis gene expression under the same growth conditions was significantly higher in the cbb3 mutant strain when compared to that in the wild type. Similarly, under anaerobic photosynthetic conditions the high levels of the oxidized carotenoid, spheroidenone, which accumulate in cbb3-deficient mutants were nearly restored to normal in a PrrB- CcoP- double mutant. This observation, together with previously published results, suggests that the regulation of the CrtA-catalyzed reaction possesses both transcriptional and posttranscriptional regulatory effectors. We propose that the cbb3 cytochrome oxidase, which by definition can interact with external oxygen, serves to control the activity of the Prr two-component activation system under both aerobic and anaerobic conditions. Although independent from the cbb3 oxidase, the RdxBH proteins are also required for normal functioning of the Prr two-component activation system and are therefore believed to lie between the cbb3 oxidase in this oxygen-sensing, redox signaling pathway and the Prr activation system.

Complex regulatory activities associated with the histidine kinase PrrB in expression of photosynthesis genes in Rhodobacter sphaeroides 2.4.1

Rhodobacter sphaeroides 2.4.1 synthesizes a specialized photosynthetic membrane upon reduction of the O2 tension below threshold levels. The genes prrB and prrA encode a sensor kinase and a response regulator, respectively, of a two-component regulatory system that presumably is involved in transduction of the signal(s) that monitors alterations in oxygen levels. A third gene, prrC, is also involved in this cascade of events. Previously, we described a mutant form of PrrB, namely, PrrB78 (J. M. Eraso and S. Kaplan, J. Bacteriol. 177:2695-2706, 1995), which results in aerobic expression of the photosynthetic apparatus. Here we examine three mutated forms of the prrB gene that have the potential to encode truncated polypeptides containing the N-terminal 6, 63, or 163 amino acids, respectively. The resulting mutant strains showed residual levels of the light-harvesting spectral complexes and had diminished photosynthetic growth rates at high light intensities with no discernible growth under intermediate or low light conditions. When either lacZ transcriptional fusions or direct mRNA determinations were used to monitor specific photosynthesis gene expression, all the mutant strains showed unexpectedly high levels of gene expression when compared to mutant strains affected in prrA. Conversely, when translational fusions were used to monitor photosynthesis gene expression in these mutant strains, expression of both puc and puf operons was reduced, especially puf expression. In light of these studies and those of the PrrB78 mutant, the data suggest that PrrA can be activated in situ by something other than PrrB, and it also appears that PrrB can function as a negative regulator acting through PrrA. Finally, we consider the role of the Prr regulatory system in the posttranscriptional control of photosynthesis gene expression.

Differential carotenoid composition of the B875 and B800-850 photosynthetic antenna complexes in Rhodobacter sphaeroides 2.4.1: involvement of spheroidene and spheroidenone in adaptation to changes in light intensity and oxygen availability

Rhodobacter sphaeroides 2.4.1 is a member of the nonsulfur purple facultative photosynthetic proteobacteria, capable of growth under a variety of cultivation conditions. In addition to the structural polypeptides and bacteriochlorophyll, the two major antenna complexes, B875 and B800-850, contain a variety of carotenoids which are an important structural and functional component of the membrane-bound photosynthetic complexes of this bacterium. Two major carotenoids, spheroidene and its keto derivative, spheroidenone, are differentially synthesized by R. sphaeroides, depending on the growth conditions. Spheroidene prevails during growth under anaerobic conditions and low light intensities, whereas spheroidenone is predominant in semiaerobically grown cells or during anaerobic growth at high light intensities. In this study, we demonstrate that in wild-type cells, spheroidene is predominantly associated with the B800-850 photosynthetic antenna complex and spheroidenone is more abundant in the B875 complex. Exploiting mutants defective in the biosynthesis of either the B875 or B800-850 light-harvesting complex, we demonstrate an association between the formation of either the B875 or B800-850 complex, on the one hand, and the accumulation of spheroidenone or spheroidene, on the other. The possible involvement of the conversion of spheroidene to spheroidenone as a significant control mechanism involved in the adaptation of R. sphaeroides to changes in light intensity and oxygen tension is discussed.

Increased production of colicin E1 in stationary phase

The synthesis of colicin E1 is known to be regulated by the SOS response, anaerobiosis, and catabolite repression. The expression of cea-lacZ fusions was also found to be stimulated when cells reached stationary phase. This increase in expression was determined to be due to depletion of nutrients from the medium, since the addition of fresh medium reversed the effect. Expression of the fusion increased when cells were starved in 10 mM MgSO4 and when they were grown in conditioned medium in which cells had been grown previously. The stimulation of expression occurred when the cea-lacZ fusion was present in single-copy as well as in multicopy plasmids. Finally, the data were consistent with this increase being independent of the SOS response, anaerobiosis, catabolite repression, and integration host factor as well as the stationary-phase regulators encoded by rpoS and lrp.

Oxygen-insensitive synthesis of the photosynthetic membranes of Rhodobacter sphaeroides: a mutant histidine kinase

Two new loci, prrB and prrC, involved in the positive regulation of photosynthesis gene expression in response to anaerobiosis, have been identified in Rhodobacter sphaeroides. prrB encodes a sensor histidine kinase that is responsive to the removal of oxygen and functions through the response regulator PrrA. Inactivation of prrB results in a substantial reduction of photosynthetic spectral complexes as well as in the inability of cells to grow photosynthetically at low to medium light intensities. Together, prrB and prrA provide the major signal involved in synthesis of the specialized intracytoplasmic membrane (ICM), harboring components essential to the light reactions of photosynthesis. Previously, J. K. Lee and S. Kaplan (J. Bacteriol. 174:1158-1171, 1992) identified a mutant which resulted in high-level expression of the puc operon, encoding the apoproteins giving rise to the B800-850 spectral complex, in the presence of oxygen as well as in the synthesis of the ICM under conditions of high oxygenation. This mutation is shown to reside in prrB, resulting in a leucine-to-proline change at position 78 in mutant PrrB (PRRB78). Measurements of mRNA levels in cells containing the prrB78 mutation support the idea that prrB is a global regulator of photosynthesis gene expression. Two additional mutants, PRRB1 and PRRB2, which make two truncated forms of the PrrB protein, possess substantially reduced amounts of spectral complexes. Although the precise role of prrC remains to be determined, evidence suggests that it too is involved in the regulatory cascade involving prrB and prrA. The genetic organization of the photosynthesis response regulatory (PRR) region is discussed.

prrA, a putative response regulator involved in oxygen regulation of photosynthesis gene expression in Rhodobacter sphaeroides

A new locus, prrA, involved in the regulation of photosynthesis gene expression in response to oxygen, has been identified in Rhodobacter sphaeroides. Inactivation of prrA results in the absence of photosynthetic spectral complexes. The prrA gene product has strong homology to response regulators associated with signal transduction in other prokaryotes. When prrA is present in multiple copies, cells produce light-harvesting complexes under aerobic growth conditions, suggesting that prrA affects photosynthesis gene expression positively in response to oxygen deprivation. Analysis of the expression of puc::lacZ fusions in wild-type and PrrA- cells revealed a substantial decrease in LacZ expression in the absence of prrA under all conditions of growth, especially when cells were grown anaerobically in the dark in the presence of dimethyl sulfoxide. Northern (RNA) and slot blot hybridizations confirmed the beta-galactoside results for puc and revealed additional positive regulation of puf, puhA, and cycA by PrrA. The effect of truncated PrrA on photosynthesis gene expression in the presence of low oxygen levels can be explained by assuming that PrrA may be effective as a multimer. PrrA was found to act on the downstream regulatory sequences (J. K. Lee and S. Kaplan, J. Bacteriol. 174:1146-1157, 1992) of the puc operon regulatory region. Finally, two spontaneous prrA mutations that abolish prrA function by changing amino acids in the amino-terminal domain of the protein were isolated.

Transcriptional regulation of puc operon expression in Rhodobacter sphaeroides. Involvement of an integration host factor-binding sequence

The putative overlapping consensus sequences (-129 to -105) for binding of fumarate nitrate reductase regulator- and integration host factor (IHF)-like proteins to puc operon upstream DNA of Rhodobacter sphaeroides was protected from DNase I digestion by purified Escherichia coli IHF. The binding of E. coli IHF to the purported IHF-binding site in the puc upstream DNA is highly sequence-specific. The recorded binding affinity was significantly lower than that of E. coli IHF to the lambda attP site. Employing site-directed changes in the DNA sequence within the -129 to -105 region, a loss in IHF binding, as monitored through gel retardation analysis, was correlated with alterations in puc operon expression monitored through the use of puc::lacZ transcriptional fusions. These results suggest that the IHF-binding site is involved in repression of puc operon transcription by oxygen as well as modulation of puc operon transcription levels by incident light intensity. Mutations specific to the upstream half of the putative fumarate nitrate reductase regulator-binding site of the puc upstream DNA did not show any physiological effects under the experimental conditions employed. Taken together, these studies reveal that the DNA sequence between -129 to -105 may involve facilitation of the interaction between upstream and downstream cis-acting regulatory sequences involved in puc operon expression.

Anaerobic control of colicin E1 production

Expression of the cea gene, which is carried by the ColE1 plasmid and which encodes colicin E1, was found to be greatly increased when the cells were grown anaerobically. By using cea-lacZ fusions to quantitate expression, aerobic levels were found to be only a few percent of the anaerobic levels. The anaerobic increase in expression was observed both in protein and in operon fusions, indicating that its regulation occurred at the level of transcription. It was also found to require a functional fnr gene and to occur when the cea-lacZ fusion was present as a single copy in the bacterial chromosome instead of in the multicopy ColE1 plasmid. Anaerobic expression was regulated by the SOS response and catabolite repression as is aerobic expression. The start site of the mRNA produced under anaerobic conditions was mapped by primer extension and found to be the same as the start for mRNA produced under aerobic conditions. These observations show that the cea gene is anaerobically regulated and that the Fnr protein is a positive regulator of transcription of this gene.