INTRODUCTION
As an exclusively human pathogen,
Streptococcuspyogenes, or group A streptococcus (GAS), has specifically adapted to evade host innate immunity and survive in multiple tissue niches. GAS primarily colonizes the epithelia of the nasopharynx and skin, but a number of clinically important serotypes are capable of invading deep tissue and the bloodstream, resulting in pathologies with a high mortality rate (
1,
2). The invasive phenotype in some of these strains has been linked to mutation of global transcriptional regulators, such as the CovRS/CsrRS two-component system (TCS), which regulates important virulence factors such as SpeB protease, streptokinase, streptodornase, and the hyaluronic acid capsule (
3–5). Transcriptional regulators of GAS are crucial to immune evasion, as observed for both stand-alone regulators like Mga (
6,
7) and TCSs like Ihk/Irr (
8,
9) and CovRS/CsrRS (
10). In the case of the latter TCS pathways, adaptation of these regulators following exposure of GAS to macrophages enhanced survival of GAS during reinfection of immune cells (
9). Thus, control of virulence factor expression through transcriptional regulator modulation has significant consequences for GAS-associated disease development.
Pathogenic streptococci exposed to host blood employ a variety of strategies to survive and propagate in this environment. The GAS cell surface specifically interacts with a high number of blood plasma proteins, and differences have been observed between the interaction profiles of invasive and noninvasive GAS strains (
11). The host coagulation cascade functions during bacterial infection of the bloodstream to entrap bacteria within a fibrin clot for killing. However, GAS can circumvent this mechanism by sequestering host plasminogen on the bacterial cell surface, converting it into active plasmin via streptokinase, and using it to degrade the fibrin network, thus liberating GAS from clots (reviewed in reference
12). Investigation of isogenic invasive and noninvasive GAS isolates found that expression of immunomodulatory factors (e.g., C5a peptidase, Gls24, and M protein) increased in response to growth in whole blood, particularly in invasive strains, and were necessary for survival in whole blood and neutrophil resistance (
13).
GAS is also able to overcome the metabolic constraints imposed on it by the blood environment. It possesses a metal transport operon (
mtsABC) (
14) as well as several dedicated hemoprotein and iron acquisition systems (
siaABC,
shr,
htsABC) (
15–17) that are controlled by multiple regulators (
mtsR and
perR; reviewed in reference
18) and are necessary for survival in blood (
19). Regulation of amino acid metabolism is required for both growth and immune evasion. For example, arginine catabolism through CodY-mediated control of the
arcABC operon (
20) contributes to resistance to inducible nitric oxide synthase (iNOS)-mediated production of nitric oxide (
21). Access to carbon sources directly impacts immune evasion, as recently demonstrated by loss of the
fruB-encoded 1-phosphopfrucokinase adversely affecting GAS survival in human blood and resistance to phagocytes (
22). Thus, it is important to understand how GAS adapts to specific host niches, metabolically and through immune evasion, to effectively address the rising incidence of invasive GAS strains (reviewed in reference
23).
In a transposon site hybridization (TraSH) forward genetic screen of M1T1 GAS 5448,
cpsY was found to be required for survival in whole human blood (
19). GAS CpsY displays 83% identity with MetR, a LysR family transcriptional regulator of
Streptococcus mutans shown to activate transcription of genes involved in methionine biosynthesis (
metE) and uptake (
atmB) in a homocysteine-dependent manner (
24). Homologs of CpsY found in
Streptococcus iniae and
Streptococcus agalactiae (MtaR) have been characterized as virulence determinants that regulate methionine transport (
25), amino acid metabolism (
26), and cell wall modifications necessary for resistance to neutrophil-mediated killing and survival
in vivo (
27–29). Thus, CpsY appears to play a crucial role as a transcriptional regulator in the development of streptococcal invasive disease through the promotion of immune cell resistance and host niche adaptation.
Here, we present evidence that the CpsY transcriptional regulator of GAS acts as a regulator of resistance to innate immune cells and survival in human blood in a manner distinct from those of other streptococcal pathogens. Loss of CpsY in GAS does not produce the methionine-dependent in vitro growth defects or the in vivo survival phenotypes observed with other streptococcal species in mouse models of infection. However, disruption of CpsY expression significantly impacts GAS resistance to human phagocytic cells, indicating that CpsY function in immune evasion is specific to the human host. Transcriptome analysis revealed that CpsY influences the expression of a small number of genes, both known and hypothetical, under our experimental conditions. The regulon includes genes encoding cell wall-modifying proteins and surface proteins, suggesting that CpsY-mediated resistance to host innate immunity may require alterations to the GAS cell envelope.
DISCUSSION
The
cpsY gene was initially identified as being required for GAS survival in a transposon site hybridization (TraSH) screen in whole human blood (
19). CpsY shares homology (83% identity) with a LysR family transcriptional regulator of
S. mutans (MetR), as well as to proteins in
S. iniae (CpsY) and
S. agalactiae (MtaR) that regulate methionine transport (
25), amino acid metabolism (
26), and resistance to neutrophil-mediated killing and survival
in vivo (
27–29). This led us to hypothesize that CpsY functions as a transcriptional regulator of GAS metabolism required for survival in human host blood. Unlike with the CpsY homologs in
S. iniae (
28) and
S. agalactiae (
25), loss of CpsY function by insertional inactivation of the encoding gene (
cpsY [
Spy0701]) did not affect growth of GAS in human plasma or rich medium (THY). CpsY was also not required for GAS growth in C medium or CDM, as seen for MtaR in
S. agalactiae, nor did supplementing media with methionine have any effect on the growth of either the (
ι)
cpsY mutant or the parental 5448 GAS strain. Finally, we found no significant differences in utilization of various carbohydrate sources by the 5448 strain and (
ι)
cpsY mutant. Together, these data suggest that the reduced fitness of a GAS mutant lacking CpsY in whole blood is not related to a broader metabolic defect.
We hypothesized that reduced survival of the (
ι)
cpsY mutant in whole blood is due to reduced resistance to immune cell activity, and/or decreased tolerance to cellular stresses encountered within these cells. Circulating neutrophils found abundantly in the bloodstream employ multiple mechanisms to kill invading bacterial pathogens, including the production of antimicrobial peptides, neutrophil extracellular traps (NETs), and phagocytosis of the pathogen. Within the phagocytic cell, engulfed pathogens must deal with production of reactive oxygen species, release of antibacterial granule components, and acidification (
31,
32). The production of virulence factors by GAS to resist these mechanisms of killing is well documented (
33), yet their coordination and regulation in different strain backgrounds is not fully understood (
18). When we tested the resistance to challenge by phagocytes in an opsonophagocytic killing assay (OKA), survival of the (
ι)
cpsY mutant strain was significantly reduced across PMNs, monocytes, and neutrophil-like differentiated HL60 cells. Inhibiting phagocytosis and blocking opsonization by treating phagocytes with cytochalasin D and heat inactivating the plasma used in the assay rescued survival of the (
ι)
cpsY mutant strain. This indicates that CpsY function is required for GAS to effectively resist opsonophagocytic killing by innate immune cells. The fact that inhibition of phagocytosis rescues survival of the (
ι)
cpsY mutant suggests that extracellular killing mechanisms of neutrophils such as NETs are likely not involved in the CpsY-associated phenotype, though we did not specifically test for this. It appears that disruption of CpsY in GAS makes the pathogen more susceptible to phagocytosis and/or less capable of intracellular growth/survival within neutrophils. This is likely due to altered expression of a CpsY regulon that may include antiphagocytic factors and/or proteins necessary for resisting cellular stresses to sustain bacterial growth within host immune cells.
The observed susceptibility to opsonophagocytic killing led us to postulate that CpsY is required for virulence in an
in vivo model of infection. We utilized both a subcutaneous (localized and disseminated) and an intraperitoneal (systemic) mouse model of infection to test for changes in virulence resulting from the loss of CpsY. The (
ι)
cpsY mutant strain did not show any decrease in virulence in either mouse model compared to both the 5448 wild-type parent and
cpsYR rescue strain. We examined whether the (
ι)
cpsY mutant strain also displayed reduced growth in whole mouse blood, and no differences were observed between the mutant and wild-type strains in the mouse blood Lancefield assay. These results lead us to conclude that the role of CpsY in resistance to opsonophagocytic killing involves a mechanism of the innate immune cells of its human host that is absent from murine immune cells. Host restriction of GAS has been associated with differences in host innate immunity. Examination of Toll-like receptor (TLR) pathways of the innate immune response to GAS found that the TLR13 response in mice is specifically triggered by GAS rRNA and that the TLR13 molecule is absent in human macrophages, as is an equivalent route for GAS RNA recognition (
34). Furthermore, there are other instances of GAS regulators that influence resistance to host immunity in a host-specific manner. FruR, the transcriptional repressor of the
fruRBA operon encoding the fructose-specific EII of the GAS phosphoenolpyruvate phosphotransferase system (PTS), was shown to mediate innate immune evasion in human host blood but not in murine models of infection (
22). Also, the human-specific antimicrobial cathelicidin LL-37 stimulates the CovRS/CsrRS two-component regulatory system-dependent expression of the capsule synthesis operon (
hasABC), the interleukin 8 (IL-8) protease PrtS/ScpC, and the integrin-like/IgG protease Mac/IdeS (
35). Induction of these virulence factors by the LL-37 peptide through CovRS/CsrRS resulted in a marked increase in GAS resistance to opsonophagocytic killing by human leukocytes.
Given its homology to transcriptional regulators of S. mutans (MetR), S. agalactiae (MtaR), and S. iniae, we also proposed that CpsY regulates expression of factors that, in the absence of CpsY, are not adequately expressed and thus fail to sustain survival of the (ι)cpsY strain in the presence of phagocytic immune cell challenge. We investigated whether there was any correspondence in the regulatory targets of GAS CpsY and its homologs in other streptococci, despite the clear differences in phenotypes that we observed. With the exception of a modest increase in artP expression during late logarithmic growth in enriched medium (THY), both RNA-seq and qRT-PCR analysis revealed that the transcript levels of known MetR/MtaR/CpsY regulon genes were unchanged in the (ι)cpsY GAS strain in all growth phases and under all medium conditions tested. These results support our finding that GAS CpsY activity does not reproduce the phenotypes observed for its homologs in other streptococci.
Moreover, our results indicate that CpsY is a pleiotropic regulator of GAS that influences expression of genes markedly different from those regulated by CpsY homologs of related streptococcal species. RNA-seq analysis of cultures grown in THY and RPMI 1640 supplemented with human plasma (OKA medium) revealed two small and distinct medium-associated sets of differentially expressed genes in the (
ι)
cpsY mutant. Genes involved in modifying the cell surface of GAS dominate the OKA medium-associated CpsY regulon. The reduced expression of genes encoding the cell wall-anchored M protein (
emm1) and cell wall-modifying enzymes (
mur1.2 and
sibA [
cdhA]) in the (
ι)
cpsY mutant suggests that CpsY induces expression of these factors during phagocyte challenge. M protein is a known antiphagocytic factor of GAS (
36–38) that is also involved in resistance to antimicrobial peptides (
39). Likewise, modification of the cell wall is important for immune cell evasion in multiple streptococcal species (
8,
40–42). Cell wall-modifying enzymes have been linked to the expression of antiphagocytic factors in GAS. Eliminating
d-alanine esterification of surface lipoteichoic acid resulted in reduced expression of
emm transcripts in 5448 and a related
covS mutant strain (8004) (
41). The underlying mechanism was not described, but Cox et al. speculated that unidentified regulatory pathways may be responsible. Whether CpsY-mediated direct or indirect regulation of
emm1,
mur1.2, and
sibA (
cdhA) is part of a mechanism linking cell wall modification and antiphagocytic factor expression in GAS is being investigated.
The manganese-specific efflux pump encoded by
mntE (
Spy0980) is another known immune evasion factor of GAS (
43) identified as differentially expressed in the absence of CpsY. Increased transcript levels of
mntE in the (
ι)
cpsY mutant suggest that CpsY represses expression of the efflux pump. Turner et al. showed that deletion of
mntE reduces GAS resistance to neutrophils and to oxidative stress (
43). Their results indicate that these phenotypes are due to loss of manganese homeostasis, which disrupts regulation of the peroxide-responsive repressor PerR and its respective regulon. However, expression of PerR or its regulon was unaffected according to our RNA-seq analysis (data not shown). Mn
2+-dependent regulation of PerR is sensitive to ion flux across the bacterial membrane, and multiple targets of PerR regulation exhibit PerR-independent expression, so it is possible that our RNA-seq analysis failed to capture CpsY-related effects on
mntE-dependent regulation of PerR.
RNA-seq analysis also identified increased transcript levels of
Spy0500 in the absence of CpsY. The
M5005_Spy0500 gene (here called
murA2) encodes a putative
N-acetylmuramoyl-
l-alanine amidase that is 78% identical to the transcript identified by Allen and Neely as regulated by CpsY and annotated as
murA in
S. iniae (
29). The gene annotated in M1T1 5448 as
murA was not differentially expressed in the (
ι)
cpsY mutant. In
S. iniae, deletion of
cpsY produces a strain (Δ
cpsY) with increased expression of
murA, a result that mirrors our observations on expression of the
Spy0500 transcript in the (
ι)
cpsY mutant. Loss of CpsY-dependent repression of
murA in
S. iniae reduced the proportion of cross-linked muropeptide and altered the cell wall structure of the bacterium (
29). However, overexpression of
murA did not replicate the muropeptide profile or the reduced neutrophil resistance phenotype observed in the Δ
cpsY strain. This result in
S. iniae and the variety of transcripts we have identified as differentially expressed in the (
ι)
cpsY strain suggest that CpsY-associated resistance to neutrophils could be due to pleiotropic effects of CpsY.
In contrast to case with the OKA medium-associated CpsY regulon, an entirely distinct set of transcripts was differentially expressed in the (
ι)
cpsY mutant strain when it was grown in rich medium (THY). Most relevant to immune cell resistance was the induced set of known virulence factors of GAS (i.e.,
speB,
spd,
nga [
spn],
prtS [
SpyCEP], and
sse). Expression of the SpeB cysteine protease promotes colonization and aids in immune cell evasion (
44–46), whereas downregulation of SpeB expression is linked to neutrophil-dependent selection of invasive GAS strains (
47–49). Streptococcal nucleases like the one encoded by
spd enhance GAS evasion of the innate immune response (
50) by degrading neutrophil extracellular traps (NETs) (
51). The NAD-glycohydrolase encoded by
nga (
spn) promotes survival and escape from the phagolysosome of phagocytic cells (
52,
53) and induces immune cell death (
54). Coordinate regulation of these factors is crucial to resistance to host immunity, and our analysis suggests that CpsY does participate in this regulation during
in vitro growth.
Another intriguing finding was the CpsY-associated differential expression of the fatty acid biosynthetic pathway. Recent studies have linked FabT-mediated regulation of fatty acid synthesis in GAS to virulence and intrahost genetic variation in invasive disease (
55). In the (
ι)
cpsY mutant strain, the majority of the
fabKDGFZ-accABCD operon was repressed (
Table 1), as was the enoyl coenzyme A (enoyl-CoA) hydratase-encoding
phaB. Conversely, the acetaldehyde-CoA/alcohol dehydrogenases encoded by
adhA (
adhP) and
adh2 (
adhE) showed increased expression. Our qRT-PCR analysis confirmed the differential expression of
adh2 (
adhE) in the (
ι)
cpsY mutant, but not that of
fabK and
phaB. The greater differential expression of
adh2 (
adhE) than of
fabK and
phaB, indicated by the log
2 fold change in expression obtained from RNA-seq analysis, suggests that this discrepancy could be due to the difference in sensitivity between RNA-seq and qRT-PCR analysis. The
fabKDGFZ-accABCD operon and
phaB constitute the majority of the unsaturated fatty acid biosynthesis pathway of GAS, whereas
adhA (
adhP) and
adh2 (
adhE) participate in fatty acid degradation. The potential influence of CpsY on fatty acid biosynthesis merits further study.
Together, the results of our investigation indicate that CpsY has a role in GAS resistance to the innate immune cells of its human host. Furthermore, GAS CpsY appears to regulate a small number of genes of varied functions whose expression CpsY influences fairly modestly. OKA medium represents the environment influencing the GAS transcriptome during immune cell challenge in our studies and thus we believe most accurately identifies the CpsY regulon relevant to resistance to opsonophagocytic killing. There was no overlap between the THY and OKA medium-associated regulons, and whether the CpsY regulon identified in THY reflects any CpsY-dependent gene regulation relevant to GAS pathogenesis is an interesting question to explore. Our investigation also indicates that the GAS CpsY regulon is markedly different from that of related streptococcal species, as the only overlap with the known MetR/MtaR/CpsY regulons of other streptococci consisted of just two genes (artP and murA2 [Spy0500]) under specific medium and growth conditions. Additionally, the role of GAS CpsY in resistance to immune cells is specific to human innate immune cells. The source of this specificity is unknown at this time and a research question we are actively pursuing.
Given that our data indicate that the CpsY regulon is medium dependent, CpsY potentially influences expression of its regulon in response to signaling by either a metabolic or immunogenic molecule. The CovRS/CsrRS-mediated induction of GAS virulence factors by the LL-37 antimicrobial peptide is an example of the latter (
35). LysR family transcriptional regulators like CpsY have a C-terminal cofactor-binding domain (reviewed in reference
56), which in various family members has been shown to bind metabolites like homocysteine in other streptococci (
57) or signaling molecules like 4-hydroxy-2-heptylquinone in
Pseudomonas aeruginosa (
58) or to sense redox changes in
Escherichia coli (
59). Another possible source of specificity of the CpsY regulon for human immune cells is that the CpsY regulon is required for optimal resistance to immune factors exclusively present in human immune cells or absent from immune cells of murine models of infection. For example, murine neutrophils do not express defensins, and the murine immune response lacks several signaling factors like CD4 on macrophages, Fc receptors like FcαR1 (CD89) on neutrophils, and FcγRIIA and FcγRIIC (
60). Investigation of the source of specificity for immune cell opsonophagocytic killing resistance of the CpsY regulon and the mechanism of CpsY-mediated gene regulation is an important avenue of research for understanding GAS pathogenesis.
MATERIALS AND METHODS
Bacterial strains and media.
Streptococcus pyogenes (GAS) strain 5448 (
61) is an M1T1 strain isolated from an invasive infection. The genome of strain 5448 was used as a reference genome, and the gene nomenclature used for this study corresponds to the MGAS5005 genome (
62). GAS bacteria were cultured in either Todd-Hewitt medium supplemented with 0.2% yeast extract (THY), C medium, or chemically defined medium supplemented with 0.5% glucose (CDM).
Escherichia coli strain DH5α (
hsdR17 recA1 gyrA endA1 relA1) was used as the host for plasmid constructions. All
E. coli strains were grown in Luria-Bertani (LB) broth. Antibiotics were used at the following concentrations: spectinomycin (Sp) at 100 μg ml
−1 for both
E. coli and GAS and kanamycin (Km) at 50 μg ml
−1 for
E. coli and 300 μg ml
−1 for GAS. Growth of GAS was assayed by measuring absorbance using a Klett-Summerson colorimeter (A filter) and expressed in Klett units. Alternatively, overnight cultures of GAS (10 ml) were grown in THY and adjusted to an optical density at 600 nm (OD
600) of 0.2 in saline, and 50-μl aliquots were added to individual wells of a 24-well plate (Corning/Costar) containing CDM with a given carbohydrate.
Construction of (ι)cpsY and rescue strains.
The (ι)
cpsY mutant strain and its rescue strain were constructed in the 5448 background using the pSinS/pHlpK mutagenesis system (
63) for stable plasmid integration into the GAS chromosome. A 345-nucleotide (nt) internal fragment of
cpsY (
M5005_Spy0701) was PCR amplified from GAS 5448 gDNA using the primers oIISpy0701F and oIISpy0701R (Table S1), digested with BamHI, and cloned into the pSinS suicide plasmid to obtain pSinS-CpsY. The pSinS-CpsY construct was introduced into pHlpK-containing GAS cells, and transformants were selected in the presence of Km and Sp at a permissive temperature (30°C) to allow replication of the two plasmids. Chromosomal integration of the mutagenic plasmid was achieved by culturing the clone at a nonpermissive temperature (37°C) in the presence of Sp only. GAS clones that were Sp resistant and Km sensitive were selected, and insertional inactivation of
cpsY by chromosomal integration of the pSinS-CpsY construct was confirmed by PCR analysis and sequencing using the primers oSISpy0701V1 and oSISpy0700V2 (Table S1). The rescue strain was obtained by reintroduction of the pHlpK plasmid into the (ι)
cpsY mutant and two overnight passages at 30°C without antibiotic selection to allow excision of the integrated plasmid. Following passage, Sp- and Km-sensitive clones were selected for at 37°C and reconstitution of the
cpsY wild-type allele was confirmed by PCR sequencing using primers oSISpy0701V1 and oSISpy0700V2 (Table S1).
Bacterial growth and antimicrobial sensitivity assays.
Growth and antimicrobial sensitivity assays employed mid-exponential-phase bacterial cultures (OD600 = 0.4) of either the wild type (5448) or mutant (ι)cpsY or rescue (cpsYR) strains of GAS. Bacteria were washed twice in phosphate-buffered saline (PBS), normalized to 1 × 103 CFU ml−1, and incubated at 37°C in 5% CO2 unless otherwise indicated. For the assessment of GAS survival in human plasma, normalized bacterial cultures were used to inoculate plasma isolated from heparinized blood of volunteer donors. Inoculated plasma was suspended 1:1 in RPMI 1640 cell culture medium (HyClone; SH30027) supplemented with 2.05 mM l-glutamine and incubated in 48-well flat-bottom plates; at the desired time points, samples were taken for plating on THY agar to obtain viable CFU counts. For the assessment of the (ι)cpsY strain response to methionine, CDM was supplemented with 50 μg ml−1 of methionine.
Carbohydrate metabolic profiles were determined using API 50 CH strips (bioMérieux). Strains were cultured overnight on Trypticase soy agar with 5% sheep blood plates (TSA II; Becton Dickinson Diagnostic), resuspended in 1 ml of saline, and vortexed for 3 min. Strains were then diluted to a final OD600 of 0.14 in 10 ml of API 50 CHL medium (bioMérieux) and cell suspension added to each of the 50 cupules, representing one carbon source each, in addition to a negative control. Strains were incubated covered at 37°C in 5% CO2 for 48 h. Utilization scores were determined at 24 h and 48 h. A “+” was given if the cupule changed from purple to yellow, indicating complete utilization. A “+/−” was given if there was a partial color change, indicating partial utilization. A “−” was given if there was no color change at all, indicating no utilization.
GAS sensitivity to antimicrobial compounds was determined as previously described (
29), with some modifications. Briefly, working concentrations of the following antimicrobials (Sigma) were prepared in THY broth: carbenicillin (3 μg ml
−1), bacitracin (4 μg ml
−1), and nisin (40 μg ml
−1). Aliquots of antimicrobial stock solutions (500 μl) were added to the last column of a 48-well flat-bottom plate and serially diluted 1:2 into 500 μl of THY broth containing normalized culture. For lysozyme treatments, wells were prepared with increasing concentrations of chicken egg white lysozyme (46,900 U mg
−1; Sigma), from 6.25 mg ml
−1 to 50 mg ml
−1 in a 500-μl total volume of THY broth. Plates were incubated overnight, and samples were taken for plating on THY agar at 5 h and 16 h to determine viable CFU counts. The minimum bactericidal concentration (MBC) was scored as the lowest antimicrobial concentration with no detectable colonies on THY agar after incubation for 5 h.
Bactericidal assays in whole blood.
The ability of GAS strains to survive in heparinized whole human blood was tested using the Lancefield blood bactericidal assay as previously described (
19,
30). Briefly, strains were grown to early mid-exponential phase (OD
600 ∼ 0.15) and serially diluted in saline. A 50-μl volume of a 10
−4 dilution (ca. 50 to 200 CFU) was added to 500 μl of fresh whole human blood and rotated for 3 h at 37°C. The multiplication factor (MF) was calculated by dividing CFU obtained from blood challenge by initial CFU inoculated. Data are presented as percent growth in blood corresponding to the MF of the mutant divided by the MF of the wild type × 100.
The bactericidal assay was adapted to use whole mouse blood as follows: 8- to 10-week-old female CD-1 mice (Charles River Laboratories) were anesthetized with ketamine, and a blood volume of ≥500 μl was withdrawn by terminal cardiac puncture and the blood immediately heparinized. A 30-μl volume of GAS diluted mid-exponential-phase culture (ca. 30 to 150 CFU) was added to 300 μl of fresh whole mouse blood and incubated as indicated above, and the MF was calculated. For the MF of each strain in murine blood obtained from multiple mice (n = 4 to 6), the MF was normalized to the MF of the wild-type strain (5448). All animal experiments were approved by the Institutional Animal Care and Use Committee (IACUC) at the University of Maryland, College Park (protocol R-16-05). Results of bactericidal assays were tested for significance by unpaired t test comparing the (ι)cpsY mutant or cpsYR rescue strain to the wild-type 5448 strain.
Cell culture.
Neutrophils (polymorphonuclear leukocytes [PMNs]) and monocytes were isolated from heparinized blood of nonimmune volunteer donors using Polymorphprep (Axis-Shield) and Ficoll Paque PLUS (GE Healthcare), respectively, as per the manufacturer's instructions. Contaminating red blood cells were removed by treatment with red blood cell lysis solution (Epicentre) and washing with Dulbecco's modified PBS (HyClone). Isolated PMNs and monocytes were both maintained in RPMI 1640 cell culture medium (HyClone) supplemented with 2.05 mM l-glutamine and 20% plasma from donor blood for the duration of the assay.
HL60 cells (Sigma) were maintained as indicated in the UAB-GBS-OPA protocol of Nahm and Burton (version E.02 [
64]). Briefly, HL60 cells line were cultured in RPMI 1640 medium supplemented with 2.05 mM
l-glutamine and 10% FBS (HyClone). Low-passage HL60 cells were differentiated into neutrophil-like cells for opsonophagocytic killing assays by supplementing culture medium with 1 μM all-trans retinoic acid (ATRA) in dimethyl sulfoxide (DMSO; 3-mg ml
−1 stock solution). PMNs, monocytes, and HL60 cells were maintained at 37°C in 5% CO
2.
OKAs.
Isolated PMNs or monocytes or differentiated HL60 cells were seeded at a density of 106 ml−1 in 24-well plates with RPMI 1640 medium supplemented with 2.05 mM l-glutamine. Wild-type GAS 5448 and mutant (ι)cpsY strains from overnight cultures were diluted into fresh THY and grown to mid-log phase (OD600 ∼ 0.4). For opsonophagocytic killing assays (OKAs), bacteria were opsonized prior to phagocyte challenge by suspension in donor plasma for 30 min at 37°C. GAS organisms were then added by centrifugation (500 × g for 5 min in a Sorvall/Heraeus 6445 rotor) to seeded phagocytes to the desired multiplicity of infection (MOI; an MOI of 0.1 was used unless otherwise indicated) at a final volume of 1 ml of RPMI 1640 medium (HyClone)–2.05 mM l-glutamine plus 20% donor plasma. Phagocyte-challenged bacteria were incubated at 37°C in 5% CO2 for 30 min or 2 h. For nonphagocytic killing assays, opsonophagocytosis of GAS was prevented by the addition of cytochalasin D (10 μg ml−1) to PMNs, monocytes, or HL60 cells 10 min prior to bacterial challenge, as well as resuspension of GAS in heat-inactivated donor plasma as required prior to addition to seeded phagocytes. GAS organisms were also incubated in RPMI 1640–2.05 mM l-glutamine plus 20% donor plasma in the absence of PMNs, monocytes, or differentiated HL60 cells for the purpose of survival comparison.
Following incubation, viable GAS were harvested by collecting well supernatants as well as neutrophils. Phagocytes were immediately lysed by resuspension in sterile H2O and the intracellular contents pelleted and recombined with corresponding well supernatants for plating on THY agar to obtain total viable bacterial counts after overnight incubation at 37°C in 5% CO2. Resistance of GAS to opsonophagocytic killing was assessed by comparing CFU obtained from plating of viable bacteria isolated from killing assays to CFU obtained from GAS incubation in cell culture media in the absence of phagocytes [(CFU obtained in PMNs/CFU obtained in media) × 100]. The percentages of viable bacteria recovered presented are normalized to that of wild-type 5448 bacteria, set to 100%, under the condition tested. Data presented are the results of at least 3 biological replicates, each performed in triplicate, and results were tested for significance by one-way analysis of variance (ANOVA) with Dunnett's multiple-comparison test for comparison to the wild-type strain (5448).
In vivo mouse models of GAS infection.
All animal experiments were approved by the IACUC at the University of Maryland, College Park (protocol R-16-05). An overnight GAS culture was used to inoculate 80 ml of THY and incubated statically at 37°C to late-logarithmic phase (Klett value, 90 to 100). A bacterial suspension in saline of approximately 4 × 109 CFU ml−1, as determined by microscope counts and verified by plating for viable colonies, was used to infect 5- to 6-week-old female CD-1 mice (Charles River Laboratories). For the systemic infection model, the mice were injected intraperitoneally (i.p.) with 100 μl of the bacterial suspension in saline (4 × 108 CFU/mouse). Mice were monitored three times daily over a period not exceeding 72 h for morbidity. For the subcutaneous infection model, the mice were anesthetized with ketamine, fur was removed from an ∼3-cm2 area of the haunch with Nair (Carter Products), and 100 μl of a cell suspension in saline (3 × 108 CFU/mouse) was injected under the skin. Mice were monitored twice daily for 7 days and were euthanized by CO2 asphyxiation upon signs of morbidity. Lesions were recorded photographically at 24, 48, and 72 h postinfection and measured using ImageJ software. Lesion size data were analyzed using GraphPad Prism (GraphPad Software) and tested for significance using an unpaired two-tailed t test with 99% confidence. Survival data were assessed by Kaplan-Meier survival analysis and tested for significance by log rank test. Data shown represent values for 10 mice for each strain.
RNA isolation and qRT-PCR analysis.
Quantitative real-time PCR (qRT-PCR) was performed as follows. Total RNA was isolated from GAS grown in enriched medium (THY) to early and mid- and late-logarithmic phase or to mid-logarithmic phase in THY (Klett value, 55 to 65), resuspended in RPMI 1640–2.05 mM
l-glutamine plus 20% plasma isolated from heparinized whole blood from human donors and incubated for 2 h as described above. RNA was isolated using the Direct-zol RNA MiniPrep kit (Zymo Research) with a modified procedure to improve GAS cell disruption, as described previously (
22). Briefly, cells were resuspended in 700 μl of TRIzol supplemented with ca. 300 mg of acid-washed glass beads (Sigma) and disrupted by vortexing for 5 min. Beads were collected by brief centrifugation, and the cell lysate was used for RNA purification as recommended by the manufacturer. Total RNA (5 μg) was subjected to DNase I treatment with the Turbo DNase-free kit (Life Technologies) to avoid gDNA contamination. Subsequently, 25 ng of DNase-treated total RNA was added to SYBR green master mix (Applied Biosystems) with 6.5 μl of each gene-specific real-time primer from a 20 nM stock (Table S1) using the one-step protocol on a Light Cycler 480 (Roche). Real-time primers were designed using the interactive tool Primer3 (
http://biotools.umassmed.edu/bioapps/primer3_www.cgi ). The levels shown represent ratios of the experimental/wild-type levels relative to
gyrA transcripts as the internal control. Standard error was calculated from three biological replicates, and differences over 2-fold in expression were considered significant.
RNA-seq and data analysis.
For RNA sequencing (RNA-seq), bacteria were grown in enriched medium (THY) either to late logarithmic phase (Klett value, 85 to 95) or to mid-logarithmic phase (Klett value, 55 to 65), resuspended in RPMI 1640–2.05 mM l-glutamine plus 20% plasma isolated from human blood donors and incubated for 2 h as described above. Following growth, cells were resuspended 1:2 in RNAprotect bacterial reagent (Qiagen) and incubated at room temperature for 5 min to stabilize RNA. Total RNA was extracted using the Direct-zol RNA miniprep kit (Zymo Research) as described above. RNA samples were treated with the Turbo DNase-free kit (Life Technologies) to avoid genomic DNA (gDNA) contamination. A total of 5 μg of DNase-treated RNA was subjected to rRNA removal using the Ribo-Zero Magnetic kit (Epicentre) for Gram-positive bacteria, and rRNA-depleted RNA was then purified with an RNAClean XP kit (Agencourt). Sample quality was assessed using a Bioanalyzer 2100 (Agilent), and sample quantity was determined using a NanoDrop 8000 spectrophotometer (Thermo Scientific).
RNA-seq directional libraries were generated using the ScriptSeq v2 RNA-seq library preparation kit (Illumina) according to the manufacturer's recommendations. Briefly, 45 ng of rRNA-depleted RNA was fragmented and used for reverse transcription with random primers containing a 5′-tagging sequence. The 5′-tagged cDNA was then modified at its 3′ end by a terminal-tagging reaction to generate di-tagged, single-stranded cDNA that was then purified using the AMPure system (Agencourt). The purified di-tagged cDNA was used as a template to generate second-strand cDNA containing Illumina adaptor sequences, to incorporate index barcodes, and to amplify the library by limited-cycle PCR. The resulting RNA-seq libraries were purified using the AMPure system (Agencourt), and RNA-seq library quality was verified as described above. Rapid-run 100-bp single-read DNA sequencing was then performed at the Institute for Bioscience and Biotechnology Research (IBBR) Sequencing Facility at the University of Maryland, College Park, using the Illumina HiSeq 1500 platform. Data were generated in the standard Sanger FastQ format.
Read quality was measured using FastQC (
65), filtered and trimmed using trimmomatic (
66), and mapped against the GAS 5448 genome at the National Center for Biotechnology Information (accession number
CP008776 ) using bowtie (
67), bowtie2 (
68), and tophat (
69), with options to allow one mismatch and randomly map multihit reads. The resulting alignments were converted to sorted BAM alignments (
70) and counted (
71) by coding and intergenic region. Initial visualizations of the sequencing mapping were performed using the Integrative Genomics Viewer (IGV) (
72). Differential expression analyses were performed following size factor and quantile normalization of read counts, and batch effect estimation was taken into account by including date in the Limma (
73) statistical model. The resulting metrics of expression were visualized using circos (
74) as well as tested for ontology enrichment using KEGG (
75), goseq (
76), clusterProfiler (
77), GOstats (
78), and topGO (
79). Correlation coefficients for RNA-seq were determined by plotting the log
2 value of the array on the
x axis to the log
2 value of the quantitative real-time PCR on the
y axis. Linear regression was used to determine the line of best fit, and the resulting
R2 value was calculated, which represented the fitness of the data.
Accession number(s).
Raw reads were deposited with the Sequence Read Archive (SRA) at the National Center for Biotechnology Information (Bioproject PJRNA351857).