Skip to main content
Advertisement
  • Loading metrics

Regulation of pneumococcal epigenetic and colony phases by multiple two-component regulatory systems

  • Juanjuan Wang,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing, China

  • Jing-Wen Li,

    Roles Formal analysis, Investigation, Methodology, Project administration, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

    Affiliation Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing, China

  • Jing Li,

    Roles Conceptualization, Funding acquisition, Investigation, Project administration, Supervision, Writing – review & editing

    Affiliation Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing, China

  • Yijia Huang,

    Roles Investigation, Methodology, Project administration, Writing – review & editing

    Affiliation Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing, China

  • Shaomeng Wang,

    Roles Conceptualization, Investigation, Project administration, Writing – review & editing

    Affiliation Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing, China

  • Jing-Ren Zhang

    Roles Conceptualization, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – review & editing

    zhanglab@tsinghua.edu.cn

    Affiliation Center for Infectious Disease Research, School of Medicine, Tsinghua University, Beijing, China

Abstract

Streptococcus pneumoniae is well known for phase variation between opaque (O) and transparent (T) colonies within clonal populations. While the O variant is specialized in invasive infection (with a thicker capsule and higher resistance to host clearance), the T counterpart possesses a relatively thinner capsule and thereby higher airway adherence and colonization. Our previous study found that phase variation is caused by reversible switches of the “opaque ON-or-OFF” methylomes or methylation patterns of pneumococcal genome, which is dominantly driven by the PsrA-catalyzed inversions of the DNA methyltransferase hsdS genes. This study revealed that switch frequency between the O and T variants is regulated by five transcriptional response regulators (rr) of the two-component systems (TCSs). The mutants of rr06, rr08, rr09, rr11 and rr14 produced significantly fewer O and more T colonies. Further mutagenesis revealed that RR06, RR08, RR09 and RR11 enrich the O variant by modulating the directions of the PsrA-catalyzed inversion reactions. In contrast, the impact of RR14 (RitR) on phase variation is independent of PsrA. Consistently, SMRT sequencing uncovered significantly diminished “opaque ON” methylome in the mutants of rr06, rr08, rr09 and rr11 but not that of rr14. Lastly, the phosphorylated form of RR11 was shown to activate the transcription of comW and two sugar utilization systems that are necessary for maintenance of the “opaque ON” genotype and phenotype. This work has thus uncovered multiple novel mechanisms that balance pneumococcal epigenetic status and physiology.

Author summary

Human pathogen Streptococcus pneumoniae undergoes reversible switches or phase variation between opaque (O) and transparent (T) colonies within clonal populations. Phase variation is caused by reversible DNA methylation patterns or epigenetic status of pneumococcal genome, which is driven by the recombinase PsrA-catalyzed site-specific recombinations. This work has discovered that switch frequency between the O and T variants is regulated by five two-component transcriptional regulators. While RR06, RR08, RR09 and RR11 act on the PsrA-catalyzed inversion reactions, RR14 impacts phase variation in an inversion-independent manner. Further analysis indicated that the phosphorylated RR11 regulator enforces the O colony phase by activating the transcription of multiple pneumococcal gene loci that are necessary for the O colony genotype and phenotype. Because two-component systems represent a fundamental mechanism for bacterial sensing and responding to environmental conditions, this work has thus uncovered a new level of molecular sophistication and complexity in bacterial fine-tuning of its biological activities and behaviors by linking the environmental sensing systems to the epigenetic switching machinery.

Introduction

Streptococcus pneumoniae (pneumococcus) is a commensal at human nasopharynx and also an important opportunistic pathogen that causes pneumonia, meningitis, septicemia, and otitis media [1]. This bacterium is well known for its strain-to-strain genetic variations, which are predominantly driven by natural transformation-mediated horizontal gene transfer. These variations are exemplified by acquisition of the exogenous genes conferring antimicrobial resistance, change of capsule types, and permutation of surface-exposed antigens. In the absence of foreign DNA, the pneumococcus is also capable of phase variation or reversible switch between the opaque (O) and transparent (T) colony variants in clonal populations on translucent agar plates [2, 3]. While the O colony variant is specialized in invasive infection (with a thicker capsule and higher resistance to host clearance), the transparent counterpart is characterized by a relatively thinner capsule, more cell wall teichoic acids, and thereby higher airway adherence [24]. In animal models, the O variant is more virulent in systemic infection while the T variants exhibit higher levels of nasopharyngeal colonization with relatively lower virulence [3, 4]. Recent studies have revealed that pneumococcal phase variation is epigenetically defined by reversible switches of the methylomes or genome methylation patterns [5, 6], which is controlled by invertase PsrA-catalyzed inversions of three homologous methyltransferase hsdS genes in the colony opacity determinant (cod) locus, also referred to as Spn556II locus [5, 7, 8].

The cod locus, a type-I restriction-modification (RM) system, consists of the hsdR (restriction endonuclease), hsdM (DNA methyltransferase), psrA (invertase), and three homologous hsdS (hsdSA, hsdSB and hsdSC) genes [5]. PsrA catalyzes extensive inversions between hsdSA and two downstream homologs (hsdSB and hsdSC) [9, 10]. As a typical sequence recognition or S subunit of the type-I RM DNA methyltransferase (MTase) [11], hsdSA encodes two target recognition domains (TRDs), each of which recognizes a half of the type-I RM methylation bipartite sequence [5]. Our recent study has shown that only hsdSA is actively transcribed and produces a functional S subunit for the MTase, whereas hsdSB and hsdSC merely serve as templates for replacement of one or two TRD domains in the HsdSA MTase by the PsrA-catalyzed inversions [5, 9]. The inversions generate six different hsdSA alleles, each of which drives methylation of a unique DNA motif in the genome of S. pneumoniae [5, 7]. This site-specific recombination mechanism enables a single pneumococcal cell to generate a mixture of progeny cells each carrying one of the six hsdSA alleles (A1-A6) and thereby producing six distinct methylomes [5]. Only the pneumococci carrying hsdSA1 form O colonies and those with the other five hsdSA alleles produce T colonies [5, 6]. Therefore, the reversible switch between O and T colony phases is epigenetically driven by reversible “ON-or-OFF” alteration of the hsdSA1 genotype in the cod locus [5], although the biochemical basis of the colony opacity remains to be characterized.

The two-component signal transduction systems (TCSs) are widely found in bacteria to regulate a variety of cellular processes in response to changes in environmental conditions, such as chemotaxis, capsule production, balance of osmolarity, photosynthesis, sporulation and transformation [12]. A typical TCS is comprised of a histidine kinase (HK) as a sensor and a cognate response regulator (RR) as an effector. The histidine kinase auto-phosphorylates itself at the conserved histidine residue when it senses extracellular changes and the cognate response regulator is subsequently activated by accepting the phosphoryl group from histidine to its conserved aspartic acid [13]. The DNA-binding output domain of activated regulator typically interacts with the promoter(s) of the target genes and activates/represses their transcription. S. pneumoniae has 13 complete TCSs and an orphan regulator (RR14, RitR) [1416]. Previous studies have shown that most of the TCSs are necessary for pneumococcal pathogenicity in animal models [15, 17], but the precise functions for most of the TCSs are largely unknown.

TCS05 (CiaR/H) and TCS12 (ComE/D) represent the best-characterized pneumococcal TCSs. TCS05 inhibits the competence state [1821], and promotes pneumococcal resistance to certain antibiotics and other stressful conditions [20, 2224]. TCS12 activates the genes associated with the competence state and natural transformation [2527]. TCS02 (WalR/K), the only TCS essential for pneumococcal viability, is associated with cell wall peptidoglycan metabolism [28, 29]; TCS03 (LiaR/S) senses cell wall damage and inhibits autolysis [30], TCS04 (PnpR/S) regulates uptake of inorganic phosphate [31, 32]; TCS06 (CbpR/S) activates expression of choline-binding protein A or PspC [33, 34]; TCS08 is involved in the cellobiose uptake [35], production of pilus-1 [16] and surface protein PavB [16]; TCS13 (BlpR/H) regulates the production of bacteriocins (pneumocins) [3638]. The orphan response regulator RR14 represses the transcription of the pneumococcal iron uptake (Piu) and other nutrient uptake loci in response to oxidative stress [39, 40]. The cellular functions of the remaining five TCSs (1, 7, 9, 10 and 11) are poorly characterized. TCS09 appears to regulate transcription of genes involved in sugar transport, competence and pilus-1 production in a strain-dependent manner [41, 42], whereas the protein structure of RR11 has been described without a known function [43]. In this study, we systematically determined the impact of 13 pneumococcal response regulators (except for the essential RR02) on phase variation in colony opacity. Five of those response regulators were found to enhance the O colony phase, mostly through modulation of PsrA-catalyzed inversions in the cod locus.

Results

Five response regulators modulate pneumococcal colony phases

To determine potential impact of pneumococcal TCSs on phase variation in colony opacity, we generated unmarked deletion mutants of 13 response regulator (rr) genes in strain ST606, a streptomycin-resistant rpsL1 derivative of type-19F strain ST556, except for the essential gene rr02 [14, 44]. The rationale was that disruption of the rr gene in bacterial TCS is more likely to yield a functional phenotype [45]. The colony phenotype or ratio between O and T colonies in a single clonally derived population was characterized for each mutant on the catalase-TSA plates. Similar to the parental strain ST606, seven of the 13 rr mutants generated populations that were dominated by O colonies (rr01, rr03, rr04, rr07, rr10, rr12 and rr13) (Table 1). As an example, the rr10 mutant produced 75.2% O and 24.8% T colonies, which is comparable with the parental strain (79.8% O and 20.2% T). In contrast, deleting the six other regulators (rr05, rr06, rr08, rr09, rr11 and rr14) led to significantly reduced proportion of O colonies as compared with the parental strain (Fig 1A). The % O values of the rr05, rr06, rr08, rr09, rr11 and rr14 mutants are 60.9, 9.0, 14.5, 17.0, 22.2 and 27.8, respectively (Table 1). It should be noted that the transparent colonies of the rr05 mutant were morphologically different from the parental counterpart. ST606 produced transparent colonies with a large translucent center and a thin halo around it, but the rr05 counterparts showed a relatively smaller transparent center and a thicker ring (Fig 1A).

thumbnail
Table 1. The ratios of T and O colonies formed by the rr mutants.

https://doi.org/10.1371/journal.ppat.1008417.t001

thumbnail
Fig 1. Colony phenotypes of the rr mutants.

A. Representative colonies formed by selected rr mutants. ST606 (WT) and isogenic rr mutants were separately spread on catalase-TSA plates, incubated for 17 hrs before the colonies were photographed under a dissection microscope. Genotype (top) and identification (bottom) of each strain are marked. Representative colonies are indicated by red (opaque) and light blue (transparent) arrowheads, respectively. B. Ratio between opaque (O) and transparent (T) colonies of selected rr mutants and their revertants. O and T colonies on each plate were prepared as in Fig 1A to quantify O and T colonies formed by ST606 (WT), isogenic rr mutants and their revertants. The mean ± SEM of 3 values (from 3 plates) for the O (filled) and T (open) colonies of each strain is presented in a single bar. C. Ratio between O and T colonies derived from single O seeding colonies of the six selected rr mutants. Three well-separated O colonies were spread on three catalase-TSA plates to assess the relative ratio between O and T colonies for each strain as described in Fig 1B. D. Same as Fig 1C except for using T colonies as the seeding bacteria.

https://doi.org/10.1371/journal.ppat.1008417.g001

To verify the impact of these six regulators on pneumococcal colony opacity, we constructed genetic revertants of these mutants by replacing each deletion with the corresponding wild type gene in its native locus of the genome. The revertants of the rr06, rr08, rr09, rr11 and rr14 mutants displayed a comparable fraction of O colonies as the parental strain (Fig 1B; S1 Fig; S1 Table). However, the rr05 revertant strain produced 100.0% O colonies. Whole genome sequencing of the rr05 mutant and its revertant revealed a single nucleotide mutation from adenine to cytosine at the 539th position in the coding region of MYY890, which encodes the substrate-binding protein in the phosphate uptake system [31]. This nonsynonymous mutation led to a substitution of the 180th glutamine (CAG) with proline (CCG) of this protein. It is worth of mentioning that other pneumococcal strains lacking rr05 also tended to carry spontaneous mutations in our collection. These findings strongly suggested that deleting rr05 somehow renders the genome more prone to spontaneous mutations, which makes it difficult to properly interpret the colony morphology data. The TCS05 was thus excluded from further investigation. These results demonstrated that RR06, RR08, RR09, RR11 and RR14 regulate pneumococcal colony opacity.

We next determined the impact of the five regulators on reversibility between O and T colonies by re-streaking individual O and T colonies of each mutant on catalase-TSA plates. For each of the five mutants, the progeny colonies derived from both the O and T colony seeds showed a comparable ratio between the O and T phenotypes. As an example, the populations derived from the O (Fig 1C) and T (Fig 1D) seeding colonies of the rr06 mutant showed the most dramatic but similar reduction in the fraction of O colonies (1.9% for O seed, 1.6% for T seed). This finding is consistent with the result from the initial screening of the rr06 mutant (Fig 1B). As a negative control, progeny colonies of the rr10 mutant did not show significant alteration in colony phenotype, regardless the original phenotypes of the seeding colonies (Fig 1C and 1D). This result indicated that the rr06, rr08, rr09, rr11 and rr14 mutants still retained the capability of phase switching.

Colony phases are regulated by the PsrA-dependent and -independent mechanisms

Our previous study has found that the reversible switch between O and T colony phenotypes is controlled by the PsrA-catalyzed inversions of the three invertible regions in the cod locus [5, 9] (Fig 2A). Therefore, enzymatic inactive psrAY247A mutant (with a point mutation in the catalytic residue tyrosine 247) is locked in either O or T colony phase [9]. To determine whether the impact of RR06, RR08, RR09, RR11 and RR14 on phase variation requires PsrA-mediated inversions, we constructed deletion mutant in each of these genes in the O-locked psrAY247A strain, which was previously shown to produce only the O colonies [5] (Fig 2B; S2 Fig). Additional deletion of rr14 in this inversion-deficient mutant still resulted in significant changes in colony phenotype. To a less extent, the Δrr14-psrAY247A mutant formed 66.9% O colonies. This result suggested that RR14 (RitR) regulates pneumococcal colony opacity in a PsrA (hsdS inversion)-independent manner. However, deleting rr06, rr08, rr09 or rr11 in the psrAY247A background did not yield any obvious impact on the phenotype of the strain. Virtually all of the colonies formed by these four double mutants showed the O phenotype (Fig 2B; S2 Fig). This result indicated that regulatory impact of RR06, RR08, RR09 and RR11 on phase variation depends on the PsrA-mediated hsdS inversions.

thumbnail
Fig 2. Characteristics of the colonies produced by the psrAY247A-rr double mutants.

A. Schematic illustration of the gene organization in the cod locus. The genes encoding the restriction enzyme (hsdR), DNA methyltransferase (hsdM), sequence recognition proteins (hsdSA, hsdSB and hsdSC) and invertase (psrA) are depicted as thin arrows at the top. Three pairs of inverted repeats (IR1, IR2 and IR3) flanking the invertible regions that mediate DNA inversions are indicated as small arrowheads in the first row of the lower panel which depicts the six major hsdSA allelic configurations generated by PsrA-catalyzed inversions. B. Ratio between O and T colonies produced by the five selected rr mutants that were generated in the psrAY247A background. The colonies were prepared and processed as in Fig 1B.

https://doi.org/10.1371/journal.ppat.1008417.g002

The PsrA-dependent regulators modulate pneumococcal methylome

Because only the methylome specified by the HsdSA1 MTase is required for the formation of O colonies [5], the abovementioned changes of the rr mutants in colony phenotypes strongly suggested that these response regulators regulate pneumococcal methylome through PsrA-driven inversions. We thus compared the methylomes of strain ST606 and 13 isogenic rr mutants by single molecule real-time (SMRT) sequencing. This trial detected N6-methyladenine (6-mA) in virtually all 664 loci of the Spn556I recognition motif (5’-TCTAGm6A-3’, type II RM) in ST606 and its 13 derivatives (S2 Table), and thereby demonstrated that our sequencing setup was adequate for genome-wide detection of all 6-mA methylated sequences.

Consistent with frequent hsdS inversions in the cod locus [5, 7], SMRT sequencing revealed 6-mA methylation in the motifs recognized by four of the six hsdSA allelic variants (HsdSA1, HsdSA2, HsdSA3 and HsdSA4) of the cod locus in ST606 (Fig 3A; Table 2). Virtually all 2,058 loci of the HsdSA1 motif in the genome were methylated (99.4%) in ST606, but there were much lower methylation percentages for the other motifs (HsdSA2−38.7%, HsdSA3−69.5% and HsdSA4−7.0%); no 6-mA methylation was detected for any loci of the HsdSA5 and HsdSA6 motifs in ST606 or its 13 derivatives. This observation is consistent with the previous finding that single clonal population of ST606 is predominantly made up by the cells possessing the hsdSA1, hsdSA2 or hsdSA3 allele [5]. As exemplified with the rr10 mutant in Fig 3B, nine of the 13 rr mutants (rr01, rr03, rr04, rr05, rr07, rr10, rr12, rr13 and rr14) exhibited a similar methylome as the parental strain (Table 2, S2 and S3 Tables).

thumbnail
Fig 3. Relative methylation rate of the DNA motifs recognized by four HsdSA allelic variants in the rr mutants.

Relative methylation rate of each DNA motif recognized by HsdSA1 (5’-CRAm6AN8CTT-3’/3’-GYTTN8Gm6AA-5’, 2,058 loci), HsdSA2 (5’-CRA m6AN9TTC-3’/3’-GYTTN9m6AAG-5’, 2,054 loci), HsdSA3 (5’-CRAm6AN8CTG-3’/3’-GYTTN8G m6AC-5’, 1,468 loci) and HsdSA4 (5’-C m6ACN7CTG-3’/3’-GTGN7G m6AC-5’, 888 loci) was calculated in each strain by dividing the number of methylated chromosomal loci for each motif with the total loci of the motif in the ST556 genome. Only the values for the parental strain ST606 (A) and isogenic mutants of rr10 (B), rr06 (C), rr08 (D), rr09 (E), rr11 (F) and rr14 (G) are presented. The results for the other rr mutants are described in Table 2.

https://doi.org/10.1371/journal.ppat.1008417.g003

thumbnail
Table 2. Methylation sequences specified by the cod MTases1.

https://doi.org/10.1371/journal.ppat.1008417.t002

In agreement with significant decrease in the fraction of O colonies in the rr06, rr08, rr09 and rr11 mutants, these strains showed dramatic alteration in the methylome specified by the “opaque ON” HsdSA1 MTase. The rr06 mutant showed the most striking overhaul in the methylome (Fig 3C). Compared with 99.4% methylation rate of the HsdSA1 motif in ST606 (Fig 3A), none of these sequences were detected as a methylated form in this mutant. The complete loss of methylated HsdSA1 motif loci in this strain was accompanied with increase in methylation rates of the HsdSA2 (2,054 copies) and HsdSA3 (1,468 copies) motifs. Identification of 6-mA methylation in the HsdSA2 and HsdSA3 motifs demonstrated that the HsdM subunit of the cod MTase was functional in the rr06 mutant. The complete loss of HsdSA1 motif methylation in the rr06 mutant was consistent with dramatic reduction in the fraction of the O colonies in this strain (Fig 1; Table 1). Similar to their common characteristics in the colony phenotype (Fig 1; Table 1), the mutants of rr08 and rr09 showed a similar pattern in methylome. Methylation of the HsdSA1 motif loci fell to 63.1% and 44.5% in the rr08 (Fig 3D) and rr09 (Fig 3E) mutants, respectively. Both mutants also displayed a moderate increase in methylation of the HsdSA3 motif loci. Interestingly, none of the 2,054 HsdSA2 motif copies were methylated in either rr08 or rr09 (Table 2). In sharp contrast, nearly all of the HsdSA2 motif loci were methylated in the rr11 mutant (94.5%) (Fig 3F). The methylation rates of the other HsdSA motifs were extremely low (5.5% for HsdSA1 and 6.3% for HsdSA3) or undetectable (HsdSA4-6) in this strain. The methylation level of various HsdSA motifs in the rr11 mutant is consistent with its colony phenotype (Fig 1). In summary, the SMRT sequencing confirmed the observation with the psrAY247A-rr double mutants that RR06, RR08, RR09 and RR11 regulate pneumococcal methylome through the PsrA-driven inversions.

In contrary to its T colony-dominant phenotype, the methylome of the rr14 mutant was virtually identical to that of the parental strain (Table 2; S2 and S3 Tables). In particular, virtually all 2,058 copies of the HsdSA1 motif were methylated in the rr14 mutant (Fig 3G). This result is consistent with their lack of association with PsrA-catalyzed inversions (Fig 2), and further validated our conclusion that RR14 impacts pneumococcal colony morphology through a non-epigenetic mechanism. DNA methylation detection also revealed certain interesting features associated with the Spn556III locus, the other functional type-I RM system in ST556 and many other pneumococcal strains [7, 46, 47]. A recent study reports that the two hsdS genes in this locus undergo excision and reintegration recombinations in other pneumococcal strains [47]. In addition to a previously identified motif (5’-GATm6AN7TCA-3’) for this system [46], the parental (5’-GGm6AN7TGA-3’) and Δrr13 (5’-GGm6AN7TCA-3’) strains each showed a new type-I RM MTase recognition sequence (S3 Table). Based on the published specificities of the HsdS variants for this system [47], it is apparent that these sequences were methylated by two new HsdS alleles generated by DNA excisions in the Spn556III locus. We designated these HsdS alleles as the HsdS1 (published), HsdS2 (new in ST606) and HsdS3 (new in Δrr13) (S3 Table). Complete shift of the methylation activity from the HsdS1 MTase to HsdS3 MTase in Δrr13 was caused by the formation of a hybrid hsdS3 allele in the Spn556III locus as reflected in the SMRT sequencing data. This observation suggested that RR13 influences the hsdS recombinations in the Spn556III locus. However, our preliminary mutagenesis trial did not reveal obvious connection between the Spn556III locus and colony phases. The mutants with unmarked deletion of the entire Spn556III locus displayed a similar methylation rate of the cod HsdSA1 motif to the parental strain.

The PsrA-dependent regulators modulate the colony phases in multiple pneumococcal strains

Based on the significant impact of RR06, RR08, RR09 and RR11 on the directions of phase variation in the ST556 background, we tested the role of these four regulators in P384 (serotype 6A) and ST877 (serotype 35B), two strains that exhibited typical hsdSA1-dependent phase variation in colony opacity [5]. In agreement with our previous study [5], strain P384 produced 76.2% O colonies on catalase-TSA plates. However, deleting rr06, rr08, rr09 or rr11 led to dramatic reduction in the fraction of O colonies within clonal populations (Fig 4A). The mutants of rr06, rr08, rr09 and rr11 produced 3.9%, 14.6%, 7.1% and 13.9% O colonies, respectively. As a negative control, the rr10 mutant retained a similar percentage of O colonies (80.9%) as the parental strain. In a similar manner, genetic deletions of rr06, rr08, rr09 or rr11 in strain ST877 also resulted in remarkable loss of O colonies (Fig 4B). As compared with the 75.3% O colonies produced by the parental strain, the mutants of rr06, rr08, rr09 and rr11 only formed 17.7%, 24.6%, 12.9% and 2.3% O colonies, respectively. In contrast, the isogenic rr10 mutant of ST877 produced a similar level of O colonies (71.7%) as the parental strain. Together, these results demonstrated that RR06, RR08, RR09 and RR11 are required for the O phenotype in multiple pneumococcal strains.

thumbnail
Fig 4. The colony characteristics of the rr06, rr08, rr09 and rr11 mutants in the P384 and ST877 strain backgrounds.

A. Representative colonies formed by P384 (6A) and its isogenic rr mutants. The colonies were generated, photographed and marked as in Fig 1A. The ratio between the O and T colonies (bottom) is displayed as in Fig 1B except for different strains. B. Representative colonies produced by ST877 (35B) and its isogenic rr mutants. Same as in (A) with the exception of different strains.

https://doi.org/10.1371/journal.ppat.1008417.g004

Phosphorylated form of RR11 drives pneumococci toward the O phase

We further characterized the molecular mechanism behind the action of RR11 in regulation of pneumococcal phase variation. We first tested whether RR11 depends on phosphorylation state of the 53rd aspartic acid residue (D53) for its role in phase variation by making a D-to-A point mutation at this site. D53 has been identified as an amino acid residue to receive phosphoryl group for activation of RR11 [43]. Similar to the rr11 deletion mutant (TH9063), the rr11D53A strain (TH11425) displayed a “hypo-O” phenotype, in which the proportion of O colonies decreased to 14.5% (Fig 5A). The effect of this point mutation on colony phenotype was successfully rescued by in situ replacing the rr11D53A allele with the wildtype rr11, suggesting that phosphorylation of RR11 is required for its modulation of phase variation. Previous studies have found that the substitution of the conserved aspartic acid residue (D) for phosphorylation in bacterial response regulators with glutamic acid (E) can trigger proteins in an activated or “constitutively phosphorylated” state in the absence of the cognate sensing kinase [48, 49]. We thus constructed a ST606 derivative carrying a rr11D53E allele by replacing the endogenous rr11 in the original chromosomal locus (TH13757). Consistent with the predicted “constitutive phosphorylation” status of the RR11 protein, the rr11D53E strain showed a “hyper-O” phenotype with 93.4% of O colonies in a clonally derived population (Fig 5A). Furthermore, reverting the rr11 deletion with the wildtype rr11 allele fully restored the colony phenotype of the mutant to that of the parental strain. The revertant produced 79.7% O colonies. These genetic and phenotypic characterizations indicated that the phosphorylated form of RR11 is required for the generation and/or maintenance of an O colony-dominated population.

thumbnail
Fig 5. Colony opacity characteristics of the tcs11 mutants.

A. Colonies characteristics of the rr11 mutants. The colony phenotypes of ST606 derivatives lacking rr11 or carrying the D53A, D53E or wildtype rr11 allele were characterized and presented essentially as in Fig 1. B. Colonies characteristics of the hk11 mutants. Same as in A except for using ST606 derivatives with various hk11 alleles. C. Colonies characteristics of the tcs11 mutants. Same as in A except for using ST606 derivative lacking both the rr11 and hk11 of the 11th two-component system (Δtcs11) and tcs11 revertant.

https://doi.org/10.1371/journal.ppat.1008417.g005

We next tested the role of HK11, the sensor histidine kinase partner of RR11, in modulation of pneumococcal colony phenotype by unmarked deletion of hk11 in ST606. The Δhk11 strain formed 41.4% O colonies, a significant reduction from that of the parental strain (80.3%), although it still produced more O colonies than the Δrr11 strain (23.2%) (Fig 5B). This result suggested that the HK11 is involved in phase variation. To verify this finding, we constructed a H-to-A mutation at the 184th histidine residue of the HK11 because this represents the highly conserved site for auto-phosphorylation of sensing kinases in bacterial TCSs [13]. The hk11H184A mutant showed only a modest reduction in proportion of O colonies (with 67.8% O colonies) although this reduction was fully restored in the hk11 revertant (Fig 5B). The milder phenotype of the hk11H184A strain than that of rr11D53A suggested that RR11 is phosphorylated by an alternative phosphoryl group donor(s) besides HK11. Hentrich et al. reported that intracellular acetyl phosphate generated by pyruvate kinase SpxB is capable of phosphorylating the RR05 (CiaR) response regulator [50]. To elaborate on the function of HK11 in phase variation, we constructed a “hyper-active” mutant of HK11 by making a T-to-P mutation in the 190th threonine residue. Previous study has demonstrated that the same mutation in the corresponding threonine residue at the 230th position of HK05 (CiaH) renders the protein constitutively active as the sensor kinase [20]. The “constitutively active” hk11T190P strain displayed a “hyper-O” phenotype (with 92.2% O colonies), which resembles the phenotype of the “constitutively phosphorylated” RR11 (rr11D53E) strain. These lines of evidence showed that the kinase activity of the sensing kinase HK11 is necessary for phosphorylation of RR11 in generation and/or maintenance of an O colony-dominated population.

Finally, we tested the phenotypic impact of the entire tcs11 locus by simultaneous deletion of both the rr11 and hk11 genes in ST606. In a dramatic fashion, the tcs11 mutant lost the ability to produce O colonies (0.4%) and produced only T colonies (Fig 5C). This deficiency was fully complemented to the parental level (81.6%) in the revertant. Taken together, the combinations of gene deletion and complementation unequivocally have demonstrated that the TCS11 system regulates the direction of phase variation in favor of the O phase by phosphorylation of the RR11 response regulator.

Phosphorylated RR11 stabilizes the “opaque ON” genotype in a PsrA-dependent manner

Because only the hsdSA1 allele, one of the six functional hsdSA alleles derived from the PsrA-driven inversions in the cod locus, enables pneumococci to produce O colonies [5, 6], the requirement of TCS11 for the O colony-dominant phenotype implied that this two-component system enriches the fraction of the hsdSA1-carrying variant in the population. We thus determined relative abundance of the hsdSA1-carrying variant in the populations of ST606 and various TCS11 mutants by measuring the relative mRNA abundance of the hsdSA1 allele in the context of the six hsdSA allelic variants. The quantitative reverse transcriptase PCR (qRT-PCR) analysis showed 41.9% hsdSA1 allele out of the six potential hsdSA alleles in the parental strain ST606 (Fig 6A). However, the hsdSA1 transcript was reduced to 17.1% and 16.5% in the rr11-null and rr11D53A mutants, respectively, but the significantly decreased representation of the hsdSA1 allele in the RR11-deficient populations was restored to the parental level in the rr11 revertant (rr11 rev). Consistent with its “opaque ON” phenotype, the “hyper-phosphorylated” rr11D53E mutant showed significantly increased hsdSA1 transcript (65.4%). This result indicated that phosphorylated RR11 enriches the hsdSA1 variant out of the six hsdSA allelic variants.

thumbnail
Fig 6. Impact of the rr11 mutations on the hsdSA1 allelic configuration.

A. Relative abundance of the hsdSA1 mRNA in the rr11 mutants. The mRNA levels of the hsdSA1 allele in the clonal populations of ST606 (WT), isogenic Δrr11 and psrAY247A mutants were detected by qRT-PCR. The transcripts of the 5’ non-invertible segment of hsdSA were similarly detected in all strains as a reference to calculate the relative CT values. The relative CT values of the ST606 and Δrr11 strains were then normalized to those of the psrAY247A strain that has the locked hsdSA1 allele (100%). The data are shown as mean ± SEM of a representative experiment. Each experiment was replicated at least twice. B. Colony opacity of the psrAY247A-rr11 double mutants. ST606 (WT) and its derivatives with either the psrAY247A allele alone or additional rr11 allelic modifications were processed for colony enumeration, photographing and data presentation as in Fig 1B. C. Relative abundance of the hsdSA1 mRNA in the psrAY247A-rr11 double mutants. Same as in Fig 6A except for different strains.

https://doi.org/10.1371/journal.ppat.1008417.g006

We next determined if the function of RR11 in regulation of pneumococcal colony phases involves PsrA-catalyzed inversions in the cod locus by constructing multiple rr11 mutations in the invertase-negative psrAY247A strain. In sharp contrast to remarkable impact of RR11 on the colony opacity phenotype in the wild type pneumococci (Fig 5), all genetic manipulations of rr11 in the psrAY247A strain did not show obvious impact on the opaque-locked phenotype. All of the rr11 mutants (Δrr11, rr11D53A and rr11D53E) in the psrAY247A background showed the “O” phenotype (Fig 6B; S3 Fig), indicating that RR11 is no longer able to impact the colony phases in the absence of the hsdS inversions. To verify the importance of PsrA in the regulation of RR11 on pneumococcal colony phases, we determined the fraction of the hsdSA1-carrying cells in various genetic backgrounds of psrA and rr11 by measuring the relative abundance of the hsdSA1 mRNA molecules in these strains. As shown in Fig 6C, all of the three rr11 mutants (Δrr11, rr11D53A and rr11D53E) generated in the psrAY247A background had a uniformly increased level of the hsdSA1 mRNAs as compared with the psrA+ ST606 strain. Taken together, these results have demonstrated that RR11 collaborates with the PsrA invertase in regulating the hsdSA allelic configurations in the cod locus.

Multiple RR11-regulated genes are required for stabilization of the O phase

To determine how RR11 regulates pneumococcal colony phases, we compared the transcriptomes between ST606 and Δrr11 mutant by RNA sequencing (RNA-seq). This trial identified a total of 23 genes with significant reduction (>2-fold change) in transcription between ST606 and Δrr11 (Table 3). This result thus revealed that RR11 acts as a transcriptional activator under these conditions. Except for comW, MYY691 and MYY1747, all the RR11-activated genes are concentrated in six apparent operons (Table 3; S4A Fig). While comW encodes the DNA-binding regulator ComW, a transcriptional activator of the competence genes [51, 52], virtually all of the other down-regulated genes in the rr11 mutant are associated with substrate uptake and metabolism. However, the rr11 mutant showed a similar level of the psrA mRNA as the parental strain (S4 Table). This result indicated that RR11 does not regulate pneumococcal colony phases through modulating the transcription of psrA.

thumbnail
Table 3. Differentially expressed genes in the rr11 mutant.

https://doi.org/10.1371/journal.ppat.1008417.t003

We next determined potential contribution of the RR11-regulated genes to the colony phases. The qRT-PCR result confirmed significant down-regulation of 24 genes within the RR11-regulated loci identified by the RNA-seq (S4B Fig). Further deletion analysis revealed significant shift in the colony phenotype in the mutants of three RR11-regulated gene loci (Fig 7A). The fraction of the O colonies in the comW mutant was dramatically diminished to 9.3% from 80.7% in the parental strain (Fig 7B). To a less extent, unmarked deletion of the genes in the two sugar phosphotransferase systems (PTS) also led to significant reduction in the fraction of the O colonies (Fig 7A). The mutants of the bgaC (MYY134-139) and ugl (MYY403-408) loci showed 42.4% and 31.9% O colonies, respectively (Fig 7B). While the bgaC locus is implicated in pneumococcal uptake and utilization of galactose and N-acetylgalactosamine [53, 54], the ugl gene cluster is necessary for uptake and metabolism of hyaluronic acid [55, 56]. The mutagenesis did not detect obvious changes in the mutants of the other RR11-regulated loci as exemplified by the ΔMYY1793-1796 (putative ABC transporter, gtfA locus), ΔMYY1924-1925 (putative ABC transporter, rr11 locus) and ΔMYY2067-2068 (arginine metabolism, arcA locus) strains (Fig 7B). This result revealed that the competence regulator ComW and two sugar utilization systems are associated with pneumococcal phase variation.

thumbnail
Fig 7. Impact of the RR11-regulated genes on the colony phases.

A. Colony opacity characteristics of the RR11-regulated gene mutants. The O and T colonies were prepared and photographed as in Fig 1A. B. Relative composition of the O and T colony types in each strain were enumerated and presented as in Fig 1B. C. Relative abundance of the hsdSA1 mRNA in the populations of the RR11-regulated gene mutants. The mRNA levels of the hsdSA1 allele in ST606 (WT) and its mutant lacking RR11-regulated genes were detected, analyzed and presented in the same manner as in Fig 6A.

https://doi.org/10.1371/journal.ppat.1008417.g007

To define how the RR11-regulated genes impact the colony phenotype, we tested if comW and the two sugar utilization loci are involved in the dominance of the hsdSA1-carrying cells in the ST606 population. Relative abundance of the hsdSA1 allele (out of the six potential hsdSA allelic variants in the cod locus) was compared between ST606 and isogenic mutants of the RR11-regulated gene loci by testing the levels of the hsdSA1 mRNA using qRT-PCR. As compared with the invertase-negative psrAY247A strain (an hsdSA1-locked strain producing 100% hsdSA1 mRNA), the parental strain ST606 produced 44.5% hsdSA1 mRNA (Fig 7C). However, the level was reduced to 3.5% and 1.9% in the mutants of the comW and ugl loci, respectively, indicating that comW and the hyaluronate utilization locus are essential for stabilization of hsdSA1 allelic configuration in the cod locus. In sharp contrast to its significant change in colony phenotype (Fig 7B), the mutant of the bgaC locus had a comparable level of the hsdSA1 allele as the parental strain (Fig 7C). This result suggested that the galactose utilization system impacts pneumococcal colony phases in an inversion-independent manner. As represented by the ΔMYY1793-1796 (gtfA locus), ΔMYY1924-1925 (rr11 locus) and ΔMYY2067-2068 (arcA locus) mutants, deleting other RR11-regulated genes did not lead to significant change in the allelic dominance of hsdSA1 in ST606 (Fig 7C). Taken together, these data showed that three of the RR11-regulated gene loci are necessary for stabilization of the O phase. While the comW and hyaluronate utilization loci achieve this function by enforcing the hsdSA1 allelic configuration in the cod locus, the galactose utilization system impacts pneumococcal colony phases in an inversion-independent manner.

Discussion

Phase variation in colony opacity has been well characterized as an important strategy for pneumococcal adaption to various host environmental conditions [3, 4]. However, it is completely unknown whether the directions of the reversible switches between O and T colonies are re-balanced in response to environmental conditions. Our systematic screening of the 13 two-component response regulators, for the first time, has demonstrated that the balance between the O and T phases of S. pneumoniae is subjective to modulation by five TCSs in the hsdS inversion-dependent (RR06, RR08, RR09 and RR11) and -independent (RR14) manners. Although the precise environmental signal(s) sensed by each of the regulators remains to be identified, this study uncovers a new level of molecular sophistication and complexity in bacterial fine tuning of its biological activities and behaviors by linking the two-component environmental sensing systems to the epigenetic switching machinery (illustrated in Fig 8). Furthermore, the colony phenotypes of the TCS mutants identified in this work are highly valuable for systematic characterization of pneumococcal two-component systems in the future because their biological contributions remain largely unknown due to the lack of phenotypic hints.

thumbnail
Fig 8. Working model depicting modulation of pneumococcal colony phases by the two-component systems.

In response to unspecified environmental cue(s) sensed by HK11, the 53rd conserved aspartic acid residue of RR11 is phosphorylated by HK11 or an alternative donor(s) of phosphoryl group, which in turn activates the comW and hyaluronate utilization locus. By the undefined mechanism(s), the RR11-activated genes act to drive the directions of the hsdS inversions toward the hsdSA1 allelic configuration and eventually the O colony phase. Likewise, RR06, RR08 and RR09 promote the hsdSA1 allelic configuration and O colony phase by the uncharacterized mechanisms. RR14 modulates colony phases in a non-epigenetic manner.

https://doi.org/10.1371/journal.ppat.1008417.g008

This study has shown that, once being phosphorylated by its cognate sensing kinase HK11 and/or other phosphoryl group donor(s), RR11 indirectly enforces the O colony phase of S. pneumoniae by transcriptional activation of its target genes. The importance of RR11 phosphorylation is supported by the similarly diminished O phase of the rr11-null and phosphorylation deficient rr11D53A strains. Consistently, the RR11 hyper-phosphorylation rr11D53E strain displayed a hyper-O colony phenotype. Mutations in the key amino acid residues of the HK11 have also corroborated the significance of RR11 phosphorylation in favoring the O phase. Lastly, complete loss of the O colonies in the TCS11-deficient mutant has provided the striking phenotypic evidence for the essential role of this two-component system in modulating the directions of PsrA-mediated inversions and phase variation. More severe phenotype of the tcs11-null mutant (0.4% O colonies) than the rr11 counterpart (reduced O colonies) suggests that HK11 may cross-phosphorylate other response regulators (e.g. RR06, RR08 and RR09) that are involved in modulating the directions of the inversion reactions in the cod locus. Cross-phosphorylation of bacterial response regulators has been well documented [57].

RR11 activates the transcription of comW and two sugar utilization gene (bgaC and ugl) loci that are necessary for the maintenance of the O colony phase. Among the 23 RR11-regulated genes identified by RNA-seq, comW and two sugar utilization gene (bgaC and ugl) loci are shown to promote the O phase. The comW and hyaluronate utilization (ugl) loci accomplish this function by modulation of the PsrA-catalyzed inversions in the cod locus because deleting comW or the entire ugl locus not only led to significant loss of the O colonies but also diminished the dominance of the bacterial cells carrying the “opaque ON” hsdSA1 allele in the clonal populations. In contrast, loss of the galactose utilization (bgaC) system only shrank the representation of the O colonies without significant change in the status of hsdSA1 in the mutant population. This information points to the conclusion that the bgaC locus impacts pneumococcal colony phases in an inversion-independent manner. While it remains to be determined if all of the genes in the ugl and bgaC sugar utilization loci are involved in modulation of the colony phases, the existing information has provided more mechanistic hint for the action of ComW.

ComW may regulate hsdS inversions through its physical interactions with unknown enhancer sequence(s) and/or PsrA. ComW is known as a competence-specific protein that stabilizes and activates the alternative sigma factor ComX for transcriptional activation of the natural transformation genes [51, 52]. A recent study shows that ComW binds DNA independent of DNA sequence, and its activity in activating pneumococcal transformation depends on the DNA binding activity [58], but this protein has not been associated with any non-competence function. In the context, ComW could act like the HU and/or Fis DNA-binding proteins of Salmonella enterica serovar Typhimurium. HU and Fis (factor for inversion stimulation) are the nucleoid-associated proteins that promote the Hin recombinase-catalyzed inversions, the best-characterized DNA inversion system in prokaryotes [59, 60]. The Hin-mediated inversions of a ~1 kb sequence flanked by the two hix sites (inverted repeat sequences recognized by the Hin recombinase) lead to alternative expression of the FljB and FliC flagellins and thereby drives phase (antigenic) variation in two different flagellar serotypes [59, 61]. While HU loops the invertible sequence [60, 62, 63], Fis enhances assembly of the supercoiling-dependent invertasome by binding to a 65-bp enhancer sequence (within the invertible sequence) [64]. The Fis/enhancer complex also determines the direction and rotation of the Hin synaptic complex [65]. It is possible that ComW enhances the formation and activity of the PsrA invertasome through its DNA-binding activity. Alternatively, binding of ComW to its target sequence(s) may regulate the hsdS inversions by altering the supercoiling state of the local sequences. It has been shown that negative DNA supercoiling is essential for the Hin-mediated inversions [60, 66].

Multiple lines of experimental evidence have also revealed that RR06, RR08 and RR09 promote the O colony phase by modulating the direction of the PsrA-catalyzed hsdS inversions. The first set of the proof came from the impact of the psrA gene on the colony phenotypes of the rr06, rr08 and rr09 mutants. Although these mutants generated in the psrA+ genetic background showed significantly reduced proportions of O colonies, the phase-locked O counterpart derived from the PsrA-deficient strain uniformly produced O colonies. Consistent with their reduced fractions of O colonies, the rr06, rr08 and rr09 mutants also exhibited diminished methylation for the DNA motif recognized by the “opaque ON” HsdSA1 MTase. Although the current data cannot explain how the pneumococci utilize multiple response regulators for regulation of hsdS inversion reactions, these factors must act in a highly coordinated manner in the context of environmental/cellular conditions. As exemplified with TCS11, further identification of the genes regulated by RR06, RR08 and RR09 under these conditions will be necessary for elucidation of the molecular mechanisms governing the actions of these TCSs in pneumococcal phase variation.

The existing information strongly suggests that pneumococcal phase variation in colony opacity is a gross reflection of pneumococcal physiological conditions, such as redox and metabolism. The linkage between pneumococcal colony phases and redox conditions is supported by significant impact of the orphan response regulator RR14 (RitR) on the colony phases. RR14 represses the transcription of the pneumococcal iron uptake (Piu) and nutrient uptake loci [39]. A recent study has shown that RR14 serves as a principal regulator in response to oxidative stress via the cysteine-mediated dimerization, which enhances pneumococcal tolerance to hydrogen peroxide and other oxidants [40]. S. pneumoniae is famous for its production of millimolar range of hydrogen peroxide as a metabolic byproduct of pyruvate oxidase SpxB [67, 68]. Consistently, SpxB is under-expressed in the opaque colony variant [69]. Potential association of pneumococcal metabolism with colony phases agrees with the known regulation of the nutrient uptake systems by RR14 [39]. Moreover, significant impact of the bgaC and ugl sugar utilization systems on the O phase could be explained by their roles in cellular metabolism. Lastly, broad impact of the five TCSs on the colony phases implies that the appearance of pneumococcal colonies is defined by vastly different mechanisms/pathways.

Materials and methods

Bacterial strains and reagents

The bacterial strains used in this work are listed in S5 Table. S. pneumoniae was grown in Todd-Hewitt broth with 0.5% yeast extract (THY) or on tryptic soy agar (TSA) plate as described [70]. Streptomycin (150 μg/ml) and kanamycin (400 μg/ml) were added in the medium when necessary. All chemical reagents used in this study were purchased from Sigma (Shanghai, China) unless otherwise indicated. DNA processing enzymes were purchased from New England Biolabs (Beijing, China). All primers were synthesized by Synbio Tech (Beijing, China) and are listed in S6 Table. All Sanger sequencing data were obtained from Ruibiotech company (Beijing, China).

Construction of pneumococcal mutants

Pneumococcal mutants were constructed in streptomycin-resistant strains ST606 (derivative of ST556, serotype 19F), TH6671 (derivative of P384, serotype 6A) and TH6675 (derivative of ST877, serotype 35B) essentially as described [5]. DNA templates, primers, restriction enzymes, resulting strains, genotypes and other details associated with mutant construction are listed in S7 Table. Briefly, JC1 (a modified Janus cassette used in our previous studies) replacement of the target sequences was generated as followed: the up- and down-stream arms of each target sequence were PCR amplified, digested with XbaI and XhoI, and ligated to a XbaI/XhoI-digested JC1 before being transformed into target pneumococcal strains. JC1 was amplified by primers Pr9840 and Pr1098 from genomic DNA of strain TH7919 (TH5445ΔbgaA::JC1) [5]. Unmarked deletions were subsequently constructed in the JC1 replacement strains by transformation with either ligation or fusion PCR products of the up- and down-stream target sequences.

The genetic revertants were constructed by transformation of JC1 replacement strains with amplicon of the corresponding wildtype target genes from genomic DNA of ST606. Gene deletions in the psrAY247A background were similarly constructed in strain TH6552 [5]. The point mutations in rr11 and hk11 were generated by introducing the mutations in the primers used for fusion PCR. Specifically, the codon of D53 (GAT) in rr11 was converted to GCG (A) and GAG (E), respectively; the codons of H184 (CAT) and T190 (ACC) in hk11 to GCC (A) and CCT (P), respectively. The bgaC, ugl, gtfA, rr11 and arcA locus deletion mutants were constructed by unmarked deletion of MYY134-139, MYY403-408, MYY1793-1796, MYY1924-1925 and MYY2067-2068 as indicated in S7 Table, respectively. Mutated sequence in each strain was confirmed by PCR amplification and DNA sequencing.

Microscopic quantification of O and T colonies

Pneumococcal colony opacity was microscopically assessed with colonies grown on catalase-TSA plates (6,000 units catalase/9-cm diameter plate) under 37 °C, 5% CO2 for 17 hrs as described [71]. Each inoculum was adjusted according to its OD620nm value to yield ~200 colonies per plate for quantification of the O and T colonies for each strain. The ratio between the two types of colonies for each strain was obtained with triplicate plates each time, and subsequently repeated at least twice. The ratio of the progeny colonies generated from individual O and T seeding colonies was determined as previously described [5].

SMRT sequencing

Methylation of pneumococcal genome was assessed by SMRT sequencing essentially as described [5]. Genomic DNA was isolated from colonies grown on catalase-TSA plates using the HiPure Bacterial DNA Kit (Magen, Beijing, China) according to the manufacturer’s protocol. Briefly, pneumococci were grown on the plates under 37 °C, 5% CO2 for 17 hrs before the colonies were scraped off with a glass spreader and washed once with pre-chilled Ringer’s solution. Bacterial lysis was achieved in 250 μl Buffer STE supplemented with 10 μl of 10% sodium deoxycholate (DOC) solution before the lysates were processed to remove RNA and proteins. The final DNA extracts were eluted with 35 μl of deionized distilled water and quantified with the NanoDrop 2000 spectrophotometer (Thermo Scientific, USA). 10Kb SMRT Bell library was constructed and sequenced using the PacBio RSII sequencing platform at the Novogene Bioinformatics Technology (Beijing, China). The low-quality reads were filtered by the SMRT Link v5.0.1 and the filtered reads was subsequently assembled by SMRT portal (Version 2.3.0) to generate one contig without gaps. The raw data for the results presented in this work are available at the NCBI database under the following accessions: SRR10083119 (ST606, wild type), SRR10083118 (TH9048 Δrr01), SRR10083113 (TH7009 Δrr03), SRR10083112 (TH9054 Δrr04), SRR10083111 (TH10784 Δrr05), SRR10083110 (TH9164 Δrr06), SRR10083109 (TH9057 Δrr07), SRR10083108 (TH9181 Δrr08), SRR10083107 (TH8468 Δrr09), SRR10083106 (TH9060 Δrr10), SRR10083117 (TH9063 Δrr11), SRR10083116 (TH9259 Δrr12), SRR10083115 (TH9066 Δrr13) and SRR10083114 (TH9167 Δrr14).

RNA-seq

RNA sequencing (RNA-seq) was carried out as described with minor modifications [19]. Specifically, pneumococci were cultured on catalase-TSA plates for 17 hrs (a time for the routine colony photographing) before colonies on multiple plates were collected in pre-chilled Ringer’s solution and pooled as described for genomic DNA purification. Bacteria in suspensions were pelleted by centrifugation at 4 °C, frozen in liquid nitrogen and stored at -80 °C. Total RNA was extracted from the frozen samples with the Purelink RNA Mini Kit (Invitrogen, USA) and further purified with Qiagen RNeasy MinElute spin columns according to the manufacturer’s protocol (Qiagen, Germany). RNA-seq was performed at the Novogene Bioinformatics Technology (Beijing, China). Trimmed reads were mapped to the genome of S. pneumoniae ST556 (CP003357.2) using Bowtie 2.3.1 and Tophat 2.1.1. Significant difference was defined by an at least 2-fold change and a P value of < 0.001. All of the raw RNA-seq data presented in this work are available in NCBI’s Gene Expression Omnibus (GEO) database (accession GSE137447). The result of each sample represents the means of two independent experiments.

qRT-PCR

Pneumococcal mRNAs were quantified by quantitative real-time reverse transcriptase PCR (qRT-PCR) as described [72]. Briefly, total RNA extracts were prepared as described for RNA-seq and used to prepare cDNA pools with iScript cDNA Synthesis Kit (Bio-Rad, USA). The 367-bp hsdSA1 allele-specific sequence was amplified with primers Pr16174/Pr16175. As an internal reference for PCR, the 5’ non-invertible region shared by the six hsdSA alleles was also amplified using primers Pr16178/Pr16179. The relative abundance of hsdSA1 mRNA was obtained by a two-tier approach. The ΔCT value for each strain was firstly calculated by subtracting the average CT value of the non-invertible hsdSA reactions from the CT values of the hsdSA1-specific reactions. Because the hsdSA1-specific sequence (367 bp) was substantially longer than the common hsdSA region (267 bp), the ΔCT value of each strain was further normalized to the counterpart of psrAY247A by subtracting the average ΔCT of psrAY247A from that of each strain. Our previous study showed that the loss-of-function mutation in psrA made the psrAY247A strain genetically locked in the hsdSA1 allelic state [5], and should produce only the hsdSA1 mRNA. The relative abundance of hsdSA1 mRNA of each strain is presented as (2-ΔΔCT) % given that the relative abundance of hsdSA1 mRNA in psrAY247A is 100%. The transcriptional levels of RR11-regulated genes were detected by qRT-PCR with the era gene as an internal control and primers listed in S8 Table. The era gene was amplified with primers Pr7932/Pr7933, which is commonly used as an internal control [72]. The relative gene expression was calculated according to the comparative 2-ΔΔCT method [73] and the ΔΔCT was calculated using the following equation: ΔΔCT = (CT gene of interest -

CT era) mutant-(CT gene of interest-CT era) ST606. The data from one representative experiment are presented as mean value of triplicate samples ± SEM (standard error of mean) for each strain. Each experiment was repeated at least twice.

Statistical analysis

The colony ratio data was statistically analyzed by two-sided Chi-square test (means); qRT-PCR, relative abundance of hsdSA1 mRNA data by two-tailed unpaired Student’s t test. The relevant data are presented as mean ± SEM. Significant differences are defined by P values of < 0.05 (*), < 0.01 (**), < 0.001(***) and < 0.0001 (****).

Supporting information

S1 Table. The opacity ratio of six regulator mutants and their isogenic revertants.

https://doi.org/10.1371/journal.ppat.1008417.s001

(DOCX)

S2 Table. Methylation sequences specified by the Spn556I MTase.

https://doi.org/10.1371/journal.ppat.1008417.s002

(DOCX)

S3 Table. Methylation sequences specified by the Spn556III MTases.

https://doi.org/10.1371/journal.ppat.1008417.s003

(DOCX)

S4 Table. The transcripts of Δrr11 mutant.

https://doi.org/10.1371/journal.ppat.1008417.s004

(XLSX)

S5 Table. Bacterial strains used in this study.

https://doi.org/10.1371/journal.ppat.1008417.s005

(DOCX)

S7 Table. PCR amplifications used for pneumococcal mutagenesis in this study.

https://doi.org/10.1371/journal.ppat.1008417.s007

(DOCX)

S8 Table. The qRT-PCR settings in this study.

https://doi.org/10.1371/journal.ppat.1008417.s008

(DOCX)

S1 Fig. Colony opacity of the isogenic revertants of rr mutants.

ST606 rr isogenic revertants were grown and processed for photographing of the colonies, and marked as in Fig 1A.

https://doi.org/10.1371/journal.ppat.1008417.s009

(TIF)

S2 Fig. Colony opacity of the psrAY247A-rr double mutants.

ST606 derivatives with either the inactive psrAY247A allele alone (TH6552) or both the psrAY247A allele and unmarked deletion of a single rr gene were grown and processed for photographing of the colonies, and marked as in Fig 1A.

https://doi.org/10.1371/journal.ppat.1008417.s010

(TIF)

S3 Fig. Colony opacity of the psrAY247A-rr11 double mutants.

ST606 derivatives with either the inactive psrAY247A allele alone (TH6552) or both the psrAY247A allele and rr11 mutants were grown and processed for photographing of the colonies, and marked as in Fig 1A.

https://doi.org/10.1371/journal.ppat.1008417.s011

(TIF)

S4 Fig. The genetic organization and transcriptional expression of RR11-regulated gene loci.

A. The genetic organization of the six RR11-regulated gene loci. The translational orientations of the genes in the six RR11-regulated gene loci are indicated with arrowheads; each gene identified with its functional or genomic names below; the number of nucleotides between two adjacent genes marked at the intergenic region. The genes deleted for mutagenesis were shown in bold. B. Transcription of the RR11-regulated genes in the Δrr11 mutant. Transcriptions of comW and the genes in the loci of MYY134-139, MYY403-408, MYY1791-1796, MYY1923-1925 and MYY2067-2068 in the ST606 (WT) and Δrr11 strains were detected by qRT-PCR. Relative transcriptional difference of each gene in the rr11 mutant is calculated by normalizing the CT value of each gene to that of the parental strain.

https://doi.org/10.1371/journal.ppat.1008417.s012

(TIF)

Acknowledgments

We are grateful for the help of Dr. Jiaxing Cui for the analysis of RNA sequencing data.

Refrences

  1. 1. Bogaert D, de Groot R, Hermans PWM. Streptococcus pneumoniae colonisation: the key to pneumococcal disease. Lancet Infect Dis. 2004;4(3):144–54. pmid:14998500
  2. 2. Weiser JN. Phase variation in colony opacity by Streptococcus pneumoniae. Microb Drug Resist. 1998;4(2):129–35. pmid:9651000.
  3. 3. Weiser JN, Austrian R, Sreenivasan PK, Masure HR. Phase variation in pneumococcal opacity: relationship between colonial morphology and nasopharyngeal colonization. Infect Immun. 1994;62(6):2582–9. pmid:8188381.
  4. 4. Kim JO, Weiser JN. Association of intrastrain phase variation in quantity of capsular polysaccharide and teichoic acid with the virulence of Streptococcus pneumoniae. J Infect Dis. 1998;177(2):368–77. pmid:9466523
  5. 5. Li J, Li JW, Feng Z, Wang J, An H, Liu Y, et al. Epigenetic switch driven by DNA inversions dictates phase variation in Streptococcus pneumoniae. PLoS Pathog. 2016;12(7):e1005762. pmid:27427949.
  6. 6. Li J, Zhang JR. Phase variation of Streptococcus pneumoniae. Microbiol Spectr. 2019;7(1):GPP3-0005-2018. pmid:30737916.
  7. 7. Manso AS, Chai MH, Atack JM, Furi L, De Ste Croix M, Haigh R, et al. A random six-phase switch regulates pneumococcal virulence via global epigenetic changes. Nat Commun. 2014;5:5055. pmid:25268848.
  8. 8. De Ste Croix M, Vacca I, Kwun MJ, Ralph JD, Bentley SD, Haigh R, et al. Phase-variable methylation and epigenetic regulation by type I restriction-modification systems. FEMS Microbiol Rev. 2017;41(Supp_1):S3–S15. pmid:28830092.
  9. 9. Li JW, Li J, Wang J, Li C, Zhang JR. Molecular mechanisms of hsdS inversions in the cod locus of Streptococcus pneumoniae. J Bacteriol. 2019;201(6):e00581–18. pmid:30617241.
  10. 10. De Ste Croix M, Chen KY, Vacca I, Manso AS, Johnston C, Polard P, et al. Recombination of the phase-variable spnIII locus is independent of all known pneumococcal site-specific recombinases. J Bacteriol. 2019;201(15)e00233–19. pmid:31085693.
  11. 11. Loenen WAM, Dryden DTF, Raleigh EA, Wilson GG. Type I restriction enzymes and their relatives. Nucleic Acids Res. 2014;42(1):20–44. pmid:24068554
  12. 12. Stock AM, Robinson VL, Goudreau PN. Two-component signal transduction. Annu Rev Biochem. 2000;69:183–215. pmid:10966457.
  13. 13. Stock J, Park P, Surette M, Levit M. Two-component signal transduction system: structure-function relationship and mechanisms of catalysis. 1995; p. 25–51. In: Hoch J, Silhavy T (ed), Two-Component Signal Transduction. ASM Press, Washington, DC. https://doi.org/10.1128/9781555818319.ch3
  14. 14. Lange R, Wagner C, de Saizieu A, Flint N, Molnos J, Stieger M, et al. Domain organization and molecular characterization of 13 two-component systems identified by genome sequencing of Streptococcus pneumoniae. Gene. 1999;237(1):223–34. pmid:10524254.
  15. 15. Throup JP, Koretke KK, Bryant AP, Ingraham KA, Chalker AF, Ge Y, et al. A genomic analysis of two-component signal transduction in Streptococcus pneumoniae. Mol Microbiol. 2000;35(3):566–76. pmid:10672179.
  16. 16. Gomez-Mejia A, Gamez G, Hammerschmidt S. Streptococcus pneumoniae two-component regulatory systems: The interplay of the pneumococcus with its environment. Int J Med Microbiol. 2018;308(6):722–37. pmid:29221986.
  17. 17. Trihn M, Ge X, Dobson A, Kitten T, Munro CL, Xu P. Two-component system response regulators involved in virulence of Streptococcus pneumoniae TIGR4 in infective endocarditis. PLoS One. 2013;8(1):e54320. pmid:23342132.
  18. 18. Schnorpfeil A, Kranz M, Kovacs M, Kirsch C, Gartmann J, Brunner I, et al. Target evaluation of the non-coding csRNAs reveals a link of the two-component regulatory system CiaRH to competence control in Streptococcus pneumoniae R6. Mol Microbiol. 2013;89(2):334–49. pmid:23710838.
  19. 19. Liu Y, Zeng Y, Huang Y, Gu L, Wang S, Li C, et al. HtrA-mediated selective degradation of DNA uptake apparatus accelerates termination of pneumococcal transformation. Mol Microbiol. 2019. pmid:31396996.
  20. 20. Guenzi E, Gasc AM, Sicard MA, Hakenbeck R. A two-component signal-transducing system is involved in competence and penicillin susceptibility in laboratory mutants of Streptococcus pneumoniae. Mol Microbiol. 1994;12(3):505–15. pmid:8065267.
  21. 21. Sebert ME, Patel KP, Plotnick M, Weiser JN. Pneumococcal HtrA protease mediates inhibition of competence by the CiaRH two-component signaling system. J Bacteriol. 2005;187(12):3969–79. pmid:15937159
  22. 22. Pinas GE, Cortes PR, Orio AG, Echenique J. Acidic stress induces autolysis by a CSP-independent ComE pathway in Streptococcus pneumoniae. Microbiology. 2008;154(Pt 5):1300–8. pmid:18451038.
  23. 23. Mascher T, Heintz M, Zahner D, Merai M, Hakenbeck R. The CiaRH system of Streptococcus pneumoniae prevents lysis during stress induced by treatment with cell wall inhibitors and by mutations in pbp2x involved in beta-lactam resistance. J Bacteriol. 2006;188(5):1959–68. pmid:16484208.
  24. 24. Dagkessamanskaia A, Moscoso M, Henard V, Guiral S, Overweg K, Reuter M, et al. Interconnection of competence, stress and CiaR regulons in Streptococcus pneumoniae: competence triggers stationary phase autolysis of ciaR mutant cells. Mol Microbiol. 2004;51(4):1071–86. pmid:14763981
  25. 25. Salvadori G, Junges R, Morrison DA, Petersen FC. Competence in Streptococcus pneumoniae and close commensal relatives: mechanisms and implications. Front Cell Infect Microbiol. 2019;9:94. pmid:31001492.
  26. 26. Pestova EV, Havarstein LS, Morrison DA. Regulation of competence for genetic transformation in Streptococcus pneumoniae by an auto-induced peptide pheromone and a two-component regulatory system. Mol Microbiol. 1996;21(4):853–62. pmid:8878046.
  27. 27. Havarstein LS, Gaustad P, Nes IF, Morrison DA. Identification of the streptococcal competence-pheromone receptor. Mol Microbiol. 1996;21(4):863–9. pmid:8878047
  28. 28. Gutu AD, Wayne KJ, Sham LT, Winkler ME. Kinetic characterization of the WalRKSpn (VicRK) two-component system of Streptococcus pneumoniae: dependence of WalKSpn (VicK) phosphatase activity on its PAS domain. J Bacteriol. 2010;192(9):2346–58. pmid:20190050.
  29. 29. Wayne KJ, Sham LT, Tsui HC, Gutu AD, Barendt SM, Keen SK, et al. Localization and cellular amounts of the WalRKJ (VicRKX) two-component regulatory system proteins in serotype 2 Streptococcus pneumoniae. J Bacteriol. 2010;192(17):4388–94. pmid:20622066.
  30. 30. Eldholm V, Gutt B, Johnsborg O, Bruckner R, Maurer P, Hakenbeck R, et al. The pneumococcal cell envelope stress-sensing system LiaFSR is activated by murein hydrolases and lipid II-interacting antibiotics. J Bacteriol. 2010;192(7):1761–73. pmid:20118250.
  31. 31. Zheng JJ, Sinha D, Wayne KJ, Winkler ME. Physiological roles of the dual phosphate transporter systems in low and high phosphate conditions and in capsule maintenance of Streptococcus pneumoniae D39. Front Cell Infect Microbiol. 2016;6:63. pmid:27379215.
  32. 32. Novak R, Cauwels A, Charpentier E, Tuomanen E. Identification of a Streptococcus pneumoniae gene locus encoding proteins of an ABC phosphate transporter and a two-component regulatory system. J Bacteriol. 1999;181(4):1126–33. pmid:9973337
  33. 33. Standish AJ, Stroeher UH, Paton JC. The two-component signal transduction system RR06/HK06 regulates expression of cbpA in Streptococcus pneumoniae. Proc Natl Acad Sci U S A. 2005;102(21):7701–6. pmid:15897461.
  34. 34. Ma Z, Zhang JR. RR06 activates transcription of spr1996 and cbpA in Streptococcus pneumoniae. J Bacteriol. 2007;189(6):2497–509. pmid:17220227.
  35. 35. McKessar SJ, Hakenbeck R. The two-component regulatory system TCS08 is involved in cellobiose metabolism of Streptococcus pneumoniae R6. J Bacteriol. 2007;189(4):1342–50. pmid:17028271
  36. 36. de Saizieu A, Gardes C, Flint N, Wagner C, Kamber M, Mitchell TJ, et al. Microarray-based identification of a novel Streptococcus pneumoniae regulon controlled by an autoinduced peptide. J Bacteriol. 2000;182(17):4696–703. pmid:10940007.
  37. 37. Reichmann P, Hakenbeck R. Allelic variation in a peptide-inducible two-component system of Streptococcus pneumoniae. Fems Microbiol Lett. 2000;190(2):231–6. pmid:11034284
  38. 38. Kjos M, Miller E, Slager J, Lake FB, Gericke O, Roberts IS, et al. Expression of Streptococcus pneumoniae bacteriocins is induced by antibiotics via regulatory interplay with the competence system. PLoS Pathog. 2016;12(2):e1005422. pmid:26840404.
  39. 39. Ulijasz AT, Andes DR, Glasner JD, Weisblum B. Regulation of iron transport in Streptococcus pneumoniae by RitR, an orphan response regulator. J Bacteriol. 2004;186(23):8123–36. pmid:15547286.
  40. 40. Glanville DG, Han L, Maule AF, Woodacre A, Thanki D, Abdullah IT, et al. RitR is an archetype for a novel family of redox sensors in the streptococci that has evolved from two-component response regulators and is required for pneumococcal colonization. PLoS Pathog. 2018;14(5):e1007052. pmid:29750817.
  41. 41. Blue CE, Mitchell TJ. Contribution of a response regulator to the virulence of Streptococcus pneumoniae is strain dependent. Infect Immun. 2003;71(8):4405–13. pmid:12874319
  42. 42. Hendriksen WT, Silva N, Bootsma HJ, Blue CE, Paterson GK, Kerr AR, et al. Regulation of gene expression in Streptococcus pneumoniae by response regulator 09 is strain dependent. J Bacteriol. 2007;189(4):1382–9. pmid:17085554
  43. 43. Park AK, Lee JH, Chi YM, Park H. Structural characterization of the full-length response regulator spr1814 in complex with a phosphate analogue reveals a novel conformational plasticity of the linker region. Biochem Biophys Res Commun. 2016;473(2):625–9. pmid:27038544.
  44. 44. Tsui HC, Zheng JJ, Magallon AN, Ryan JD, Yunck R, Rued BE, et al. Suppression of a deletion mutation in the gene encoding essential PBP2b reveals a new lytic transglycosylase involved in peripheral peptidoglycan synthesis in Streptococcus pneumoniae D39. Mol Microbiol. 2016;100(6):1039–65. pmid:26933838.
  45. 45. Lukat GS, McCleary WR, Stock AM, Stock JB. Phosphorylation of bacterial response regulator proteins by low molecular weight phospho-donors. Proc Natl Acad Sci U S A. 1992;89(2):718–22. pmid:1731345.
  46. 46. Feng Z, Li J, Zhang JR, Zhang X. qDNAmod: a statistical model-based tool to reveal intercellular heterogeneity of DNA modification from SMRT sequencing data. Nucleic Acids Res. 2014;42(22):13488–99. pmid:25404133.
  47. 47. Kwun MJ, Oggioni MR, De Ste Croix M, Bentley SD, Croucher NJ. Excision-reintegration at a pneumococcal phase-variable restriction-modification locus drives within- and between-strain epigenetic differentiation and inhibits gene acquisition. Nucleic Acids Res. 2018;46(21):11438–53. pmid:30321375
  48. 48. Smith JG, Latiolais JA, Guanga GP, Pennington JD, Silversmith RE, Bourret RB. A search for amino acid substitutions that universally activate response regulators. Mol Microbiol. 2004;51(3):887–901. pmid:14731287.
  49. 49. Martin B, Soulet AL, Mirouze N, Prudhomme M, Mortier-Barriere I, Granadel C, et al. ComE/ComE~P interplay dictates activation or extinction status of pneumococcal X-state (competence). Mol Microbiol. 2013;87(2):394–411. pmid:23216914.
  50. 50. Hentrich K, Lofling J, Pathak A, Nizet V, Varki A, Henriques-Normark B. Streptococcus pneumoniae senses a human-like sialic acid profile via the response regulator CiaR. Cell Host Microbe. 2016;20(3):307–17. pmid:27593514.
  51. 51. Sung CK, Morrison DA. Two distinct functions of ComW in stabilization and activation of the alternative sigma factor ComX in Streptococcus pneumoniae. J Bacteriol. 2005;187(9):3052–61. pmid:15838032.
  52. 52. Luo P, Li H, Morrison DA. Identification of ComW as a new component in the regulation of genetic transformation in Streptococcus pneumoniae. Mol Microbiol. 2004;54(1):172–83. pmid:15458414.
  53. 53. Jeong JK, Kwon O, Lee YM, Oh DB, Lee JM, Kim S, et al. Characterization of the Streptococcus pneumoniae BgaC protein as a novel surface beta-galactosidase with specific hydrolysis activity for the Galbeta1-3GlcNAc moiety of oligosaccharides. J Bacteriol. 2009;191(9):3011–23. pmid:19270088.
  54. 54. Afzal M, Shafeeq S, Ahmed H, Kuipers OP. N-acetylgalatosamine-mediated regulation of the aga operon by AgaR in Streptococcus pneumoniae. Front Cell Infect Microbiol. 2016;6:101. pmid:27672623.
  55. 55. Maruyama Y, Nakamichi Y, Itoh T, Mikami B, Hashimoto W, Murata K. Substrate specificity of streptococcal unsaturated glucuronyl hydrolases for sulfated glycosaminoglycan. J Biol Chem. 2009;284(27):18059–69. pmid:19416976.
  56. 56. Marion C, Stewart JM, Tazi MF, Burnaugh AM, Linke CM, Woodiga SA, et al. Streptococcus pneumoniae can utilize multiple sources of hyaluronic acid for growth. Infect Immun. 2012;80(4):1390–8. pmid:22311922
  57. 57. Agrawal R, Sahoo BK, Saini DK. Cross-talk and specificity in two-component signal transduction pathways. Future Microbiol. 2016;11(5):685–97. pmid:27159035
  58. 58. Inniss NL, Prehna G, Morrison DA. The pneumococcal sigma(X) activator, ComW, is a DNA-binding protein critical for natural transformation. J Biol Chem. 2019;294(29):11101–18. pmid:31160340.
  59. 59. Johnson RC. Site-specific DNA inversion by serine recombinases. Microbiol Spectr. 2015;3(1):MDNA3-0047-2014. pmid:26104558.
  60. 60. Johnson RC, Bruist MF, Simon MI. Host protein requirements for in vitro site-specific DNA inversion. Cell. 1986;46(4):531–9. pmid:3524854.
  61. 61. Zieg J, Silverman M, Hilmen M, Simon M. Recombinational switch for gene expression. Science. 1977;196(4286):170–2. pmid:322276.
  62. 62. Hillyard DR, Edlund M, Hughes KT, Marsh M, Higgins NP. Subunit-specific phenotypes of Salmonella typhimurium HU mutants. J Bacteriol. 1990;172(9):5402–7. pmid:2168381.
  63. 63. Haykinson MJ, Johnson RC. DNA looping and the helical repeat in vitro and in vivo: effect of HU protein and enhancer location on Hin invertasome assembly. EMBO J. 1993;12(6):2503–12. pmid:8508775.
  64. 64. McLean MM, Chang Y, Dhar G, Heiss JK, Johnson RC. Multiple interfaces between a serine recombinase and an enhancer control site-specific DNA inversion. Elife. 2013;2:e01211. pmid:24151546
  65. 65. Dhar G, Heiss JK, Johnson RC. Mechanical constraints on Hin subunit rotation imposed by the Fis/enhancer system and DNA supercoiling during site-specific recombination. Mol Cell. 2009;34(6):746–59. pmid:19560425.
  66. 66. Johnson RC, Bruist MB, Glaccum MB, Simon MI. In vitro analysis of Hin-mediated site-specific recombination. Cold Spring Harb Symp Quant Biol. 1984;49:751–60. pmid:6099257.
  67. 67. Pericone CD, Park S, Imlay JA, Weiser JN. Factors contributing to hydrogen peroxide resistance in Streptococcus pneumoniae include pyruvate oxidase (SpxB) and avoidance of the toxic effects of the fenton reaction. J Bacteriol. 2003;185(23):6815–25. pmid:14617646.
  68. 68. Lisher JP, Tsui HT, Ramos-Montanez S, Hentchel KL, Martin JE, Trinidad JC, et al. Biological and chemical adaptation to endogenous hydrogen peroxide production in Streptococcus pneumoniae D39. mSphere. 2017;2(1). pmid:28070562.
  69. 69. Overweg K, Pericone CD, Verhoef GG, Weiser JN, Meiring HD, De Jong AP, et al. Differential protein expression in phenotypic variants of Streptococcus pneumoniae. Infect Immun. 2000;68(8):4604–10. pmid:10899862.
  70. 70. Chen H, Ma Y, Yang J, O’Brien CJ, Lee SL, Mazurkiewicz JE, et al. Genetic requirement for pneumococcal ear infection. Plos One. 2008;3(8):e2950. pmid:18670623
  71. 71. Li J, Wang J, Jiao F, Zhang J-R. Observation of pneumococcal phase variation in colony morphology. Bio-protocol. 2017;7(15):e2434.
  72. 72. Liu X, Li JW, Feng Z, Luo Y, Veening JW, Zhang JR. Transcriptional repressor PtvR regulates phenotypic tolerance to vancomycin in Streptococcus pneumoniae. J Bacteriol. 2017;199(14). pmid:28484041.
  73. 73. Schmittgen T. D. Livak K. J. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8. pmid:18546601