A frameshift in Yersinia pestis rcsD alters canonical Rcs signalling to preserve flea-mammal plague transmission cycles

Abstract
Editor's evaluation
eLife digest
Introduction
Results
Discussion
Materials and methods
Data availability
References
Article and author information
Metrics

Abstract

Multiple genetic changes in the enteric pathogen Yersinia pseudotuberculosis have driven the emergence of Yesinia pestis, the arthropod-borne, etiological agent of plague. These include developing the capacity for biofilm-dependent blockage of the flea foregut to enable transmission by flea bite. Previously, we showed that pseudogenization of rcsA, encoding a component of the Rcs signalling pathway, is an important evolutionary step facilitating Y. pestis flea-borne transmission. Additionally, rcsD, another important gene in the Rcs system, harbours a frameshift mutation. Here, we demonstrated that this rcsD mutation resulted in production of a small protein composing the C-terminal RcsD histidine-phosphotransferase domain (designated RcsD-Hpt) and full-length RcsD. Genetic analysis revealed that the rcsD frameshift mutation followed the emergence of rcsA pseudogenization. It further altered the canonical Rcs phosphorylation signal cascade, fine-tuning biofilm production to be conducive with retention of the pgm locus in modern lineages of Y. pestis. Taken together, our findings suggest that a frameshift mutation in rcsD is an important evolutionary step that fine-tuned biofilm production to ensure perpetuation of flea-mammal plague transmission cycles.

Editor's evaluation

This is a valuable study that provides convincing evidence for fine-tuning of a signal transduction pathway in the emergence of the bacterial pathogen Yersinia pestis. The work advances our understanding of bacterial signal transduction, and will be of interest to those studying the evolution of Y. pestis and its adaptation to flea-borne transmission.

https://doi.org/10.7554/eLife.83946.sa0

eLife digest

Yersinia pestis, the agent responsible for the plague, emerged 6,000 to 7,000 years ago from Yersinia pseudotuberculosis, another type of bacteria which still exists today. Although they are highly similar genetically, these two species are strikingly different. While Y. pseudotuberculosis spreads via food and water and causes mild stomach distress, Y. pestis uses fleas to infect new hosts and has killed millions.

A small set of genetic changes has contributed to the emergence of Y. pestis by allowing it to thrive inside a flea and maximise its transmission. In particular, some of these mutations have led to the bacteria being able to come together to form a sticky layer that adheres to the gut of the insect, with this ‘biofilm’ stopping the flea from feeding on blood. The starving flea keeps trying to feed, and with each bite comes another opportunity for Y. pestis to jump host. However, it remains unclear exactly how the mutations have influenced biofilm formation to allow for this new transmission mechanism to take place.

To examine this phenomenon, Guo et al. focused on rcsD, a gene that codes for a component of the signalling system that controls biofilm creation. In Y. pestis this sequence has been mutated to become a ‘pseudogene’, a type of sequence which is often thought to be non-functional. However, the experiments showed that, in Y. pestis, rcsD could produce small amounts of a full-length RcsD protein similar to the one found in Y. pseudotuberculosis. However, the gene mostly produces a short ‘RcsD-Hpt’ protein that can, in turn, alter the expression of many genes, including those that decrease biofilm formation. This may prove to be beneficial for Y. pestis, for example when the bacteria switches from living in fleas to living in humans, where it does not require a biofilm.

Guo et al. further investigated the impact of rcsD becoming a pseudogene inY. pestis, showing that if normal amounts of the full-length RcsD protein are produced, the bacteria quickly lose the gene that allows them to form biofilm in fleas, and cause disease in humans. In fact, additional analyses revealed that all sequenced strains of ancient and modern Y. pestis bacteria can produce RcsD-Hpt, even if they do not carry the same exact rcsD mutation. Overall, these results indicate that rcsD turning into a pseudogene marked an important step in the emergence of Y. pestis strains that can cause lasting plague outbreaks. They also point towards pseudogenes having more important roles in evolution than previously thought.

Introduction

Approximately 6000–7000 years ago, Yersinia pestis evolved to be an arthropod-borne pathogen from its ancestor Yersinia pseudotuberculosis (Cui et al., 2013; Rasmussen et al., 2015; Spyrou et al., 2018). Despite their recent divergence, these species have markedly different life cycles. Y. pseudotuberculosis is transmitted by the faecal-oral route and usually causes a mild, self-limiting enteric disease in mammals (Putzker et al., 2001). Y. pestis, uniquely amongst enteric Gram-negative bacteria, is highly virulent and relies on flea-borne transmission (Hinnebusch, 2005). The co-occurrence of both ancestor and descendant Yersinia species provides an exemplary model to study the evolution of bacterial pathogens (Hinchliffe et al., 2003; Hinnebusch, 1997; Wren, 2003).

Y. pestis is transmitted by flea bites via a crude regurgitation mechanism (Hinnebusch and Erickson, 2008). After entering the flea gut, the planktonic Y. pestis form large aggregates and colonize the proventriculus. When infected fleas feed again as early as 1–3 days post infection, the bacterial mass in the proventriculus is regurgitated into the bite site, leading to a phenomenon referred to as early phase transmission (Bosio et al., 2020; Dewitte et al., 2020; Eisen et al., 2007; Hinnebusch and Erickson, 2008). Concurrently, Y. pestis continues to multiply in the flea digestive tract and forms HmsHFRS-dependent biofilms in the proventriculus, which blocks flea feeding (Abu Khweek et al., 2010). Continuous attempts to feed by the starved flea promotes reflux of bacteria-contaminated blood to the bite site, a phenomenon termed biofilm-dependent transmission (Bland et al., 2018). In contrast to the low efficiency of early phase transmission, which usually requires several infected fleas simultaneously feeding on a naïve host, biofilm-dependent late-stage transmission is highly efficient, and a single blocked flea has high potential for transmission (Bosio et al., 2020).

Y. pestis diverged from Y. pseudotuberculosis through a series of gene gains, gene losses, and genome rearrangements (Cao et al., 2022; Hinchliffe et al., 2003; McNally et al., 2016; Rascovan et al., 2019; Sun et al., 2014). Acquisition of the ymt gene enabled Y. pestis to survive and colonize the flea midgut to sustain flea-borne plague through expansion of its mammalian host range (Bland et al., 2021; Sun et al., 2014). Loss of the three genes, rcsA, pde2, and pde3, altered the c-di-GMP signalling pathway, which increased biofilm-forming capability (Sun et al., 2008; Sun et al., 2014). These four genetic changes enabled the Y. pseudotuberculosis progenitor strain to form biofilms in the proventriculus of fleas, promoting a flea-borne transmission modality (Sun et al., 2014).

The Rcs phosphorelay system, a non-orthodox two-component signal transduction system, consists of a hybrid sensor kinase RcsC, the phosphotransfer protein RcsD, and a response regulator RcsB (Guo and Sun, 2017; Wall et al., 2018). In Enterobacteriaceae, the outer membrane (OM) lipoprotein, RcsF, senses OM- and peptidoglycan-related stress (Smith et al., 2021; Tata et al., 2021). The inner membrane protein, IgaA, subsequently relays these signals to RcsD, which activates the Rcs phosphorelay system (Cho et al., 2014; Wall et al., 2020). Autophosphorylated RcsC transfers a phosphate group to a conserved histidine residue in the C-terminal histidine-phosphotransferase (HPt) domain of RcsD, which is finally transferred to RcsB (Takeda et al., 2001). Phosphorylated RcsB acts either alone or in combination with auxiliary proteins to regulate expression of target genes (Clarke, 2010; Huesa et al., 2021; Pannen et al., 2016). Functional RcsA works in concert with RcsB as a heterodimer to inhibit Y. pseudotuberculosis biofilm formation, in part by repressing expression of the c-di-GMP synthesis genes hmsT and hmsD (Bobrov et al., 2011; Fang et al., 2015; Guo et al., 2015; Sun et al., 2012; Sun et al., 2011). In Y. pestis, however, RcsA (RcsA_pe) was disrupted by acquiring a 30 bp repeat insertion sequence, leading to enhanced capacity for biofilm formation (Sun et al., 2008). rcsD in Y. pestis (rcsD_pe) is a putative pseudogene due to a 1 bp deletion, but it retains a limited ability to modulate biofilm formation in vitro (Sun et al., 2008).

Here, we investigated the functional consequences of the frameshift mutation in rcsD_pe and the mechanism by which it modulates activity of the Rcs signalling system. We found that the rcsD_pe variant produces a 103-amino acid protein containing the C-terminal HPt domain of RcsD, designated as RcsD-Hpt, and that this protein plays a dominant role in Rcs signalling in Y. pestis. Frameshifted rcsD alters Rcs signal transduction, subsequently modulating the expression of dozens of genes and the capacity for biofilm formation in Y. pestis. These evolutionary events may represent an important step in the emergence of ubiquitous branches of Y. pestis that can capably maintain plague outbreaks.

Results

rcsD_pe negatively regulates, while rcsD_pstb positively regulates, biofilm formation in Y. pestis

RcsD, an inner membrane protein, has an HPt domain at its C-terminus (Takeda et al., 2001; Wall et al., 2018). The rcsD gene in Y. pestis has undergone a frameshift after codon 642 due to a single nucleotide deletion (Figure 1A). However, rcsD_pe is still functional as deletion of the N-terminal region (1–1846 bp, ΔrcsD_N-term) of rcsD_pe reduced biofilm formation and Congo red (CR) pigmentation (Figure 1B and C; Sun et al., 2008), while replacement of rcsD_pe with rcsD_pstb, and rcsB deletion significantly increased biofilm formation (Figure 1B; Sun et al., 2008). To further characterize the differences between rcsD_pe and rcsD_pstb, biofilm formation was determined in Y. pestis wild type and rcsD N-terminal deletion strains carrying plasmid-encoded copies of each gene, respectively (Figure 1B and C). Consistent with our previous report (Sun et al., 2008), overexpression of rcsD_pstb significantly increased CR pigmentation and biofilm formation in vitro (Figure 1B and C). By contrast, overexpression of rcsD_pe decreased the CR phenotype and formation of biofilm (Figure 1B and C), indicating that rcsD_pe might play a different role to rcsD_pstb in the modulation of Rcs signalling.

Figure 1

Download asset Open asset

*rcsD*_pe negatively regulates biofilm formation, while *rcsD*_pstb positively regulates biofilm formation in Y. *pestis*.

(A) Schematic representation of the *rcsD* frameshift mutation that occurred during speciation of *Y. pestis* from its ancestor *Y. pseudotuberculosis. rcsD*_pstb: *rcsD* from *Y. pseudotuberculosis*; *rcsD*_pe: *rcsD* from *Y. pestis*; △*rcsB*, *rcsB* deletion; RcsD_N-term: N-terminal fragment of RcsD; RcsD-Hpt: C-terminal HPt domain of RcsD. Crystal violet (CV) binding assay (B) and Congo red (CR) pigmentation assay (C) using *Y. pestis* KIM6+ (WT), an N-terminal deletion mutant (△*rcsD*_N-term), and an *rcsD_pstb* substitution strain (*rcsD*_pe::*rcsD*_pstb) and their derivatives that carry plasmids harbouring *rcsD*_pe (prcsD_pe) or *rcsD*_pstb (prcsD_pstb). (D) Schematic representation of the *rcsD* mutations constructed in this study. CV binding assay (E) and CR pigmentation assay (F) using *Y. pestis* KIM6+ (WT) and its derivative strains harbouring plasmids expressing different truncations of RcsD. CV assays in panels B and E were performed together. Error bars represent ± SD from three independent experiments with three replicates. Statistical analysis was performed using one-way analysis of variance (ANOVA) with Dunnett’s multiple comparisons post-test. ns, not significant; *p<0.05, **p<0.01, ***p<0.001.

Next, we constructed two mutants, one lacking the entire rcsD_pe gene, and the other lacking a C-terminal fragment (ΔrcsD and rcsD_C-term) (Figure 1D). In contrast to the rcsD N-terminal deletion mutant (ΔrcsD_N-term), these two strains showed comparable phenotypes with increased CR adsorption and in vitro biofilm formation (Figure 1E and F). The opposing effects on biofilm formation in these mutants suggest that the C-terminal region may contribute to the inhibition of biofilm formation. To test this, we expressed rcsD_C-term, rcsD_N-term, and rcsD_pe in a rcsD deletion mutant. Overexpression of rcsD_C-term and rcsD_pe but not rcsD_N-term greatly decreased CR pigmentation and biofilm formation (Figure 1E and F). Collectively, these results demonstrate that RcsD_pstb positively regulates, while the frameshifted rcsD_pe negatively regulates, biofilm formation, and that this activity depends on the C-terminus.

rcsD_pe expresses intact RcsD and a small HPt-containing domain protein RcsD-Hpt

The C-terminal HPt domain of rcsD is phosphorylated by RcsC and subsequently transfers a phosphate group to the response regulator RcsB in vitro (Takeda et al., 2001). Given that frameshifted rcsD_pe is functional, and expression of the C-terminus of rcsD is sufficient to repress biofilm formation, we hypothesized that a small protein containing the HPt domain is expressed by rcsD_pe. To identify the putative open reading frame (ORF) of this latter protein, we analysed candidate translational initiation sites in the C-terminus of rcsD_pe by searching for AUG, GUG, and UUG codons, which account for approximately 80%, 12%, and 8% of start codons in bacterial genomes (Srivastava et al., 2016). We found two putative ORFs (encoding genes composing 462 and 573 nucleotides, respectively, Figure 2—figure supplement 1A) with UUG start codons, denoted as UUG^–462 and UUG^–573. Ectopic expression of these two putative ORFs in Y. pestis KIM6+ showed comparable biofilm formation to expression of full-length rcsD_pe (Figure 2—figure supplement 1B). Mutation of these two start codons alone (UUG^–462→UUA, UUG^–573→CUU) or together did not alter the function of rcsD_pe (Figure 2—figure supplement 1C), indicating they are not the start sites promoting expression of the functional protein.

AUU, CUG, and AUC are occasionally used as start codons in bacteria (Cao and Slavoff, 2020; D’Lima et al., 2017). We identified a putative ORF, encoding a 103-residue protein, initiated by an AUU codon, which also had a predicted RBS binding site (Figure 2A). Mutation of this codon or the predicted RBS region in rcsD_pe abolished its function, while replacement of AUU with a strong AUG start codon significantly enhanced its function (Figure 2B and C). Furthermore, ectopic expression of the 103-residue HPt domain with an AUU start codon, but not with a GGU codon, using a modified pBAD vector in the Y. pestis KIM6+ wild type strain, or rcsD deletion mutant, strongly repressed biofilm formation and CR pigmentation (Figure 2—figure supplement 1D), indicating that AUU is a functional start codon initiating translation of RcsD-Hpt.

Figure 2 with 2 supplements see all

Download asset Open asset

Expression of RcsD-Hpt is initiated by an uncommon AUU start codon and negatively regulates biofilm formation in Y. *pestis*.

(A) Schematic representation of mutations introduced to the putative RBS and start codon of RcsD-Hpt. ATT^–312→GGT or ATT^–312→ATG indicate that the predicted ATT start codon was mutated to GGT or ATG, respectively. RBS* indicates that the putative RBS of RcsD-Hpt was mutated. Congo red (CR) pigmentation assay (B) and crystal violet (CV) binding assay (C) using *Y. pestis* KIM6+ (WT) and its derivative strains. Error bars represent ± SD from three independent experiments with three replicates. Statistical analysis was performed using one-way analysis of variance (ANOVA) with Dunnett’s multiple comparisons post-test, comparing each construct with WT. *p<0.05, **p<0.01, ***p<0.001. (D) Expression of RcsD in *Y. pestis* and its derivatives were detected by western blot analysis using an anti-Flag antibody (see Figure 2—source data 1 details). The 3xFlag and His6 epitope tags (FH) encoding sequences were fused to the C-terminus of *rcsD*_pe and *rcsD*_pstb. (E) Structure of RcsD-Hpt (103 residues) predicted by AlphaFold2. Conserved His50 is located at the a3 helix of RcsD-Hpt. (F) Structure comparison of RcsD-Hpt (103 residues) by AlphaFold2 (yellow) and HptB of *Pseudomonas aeruginosa* PAO1 (PDB 7C1l, cyan).

Figure 2—source data 1 Raw source data for Figure 2D. The original file of the full raw unedited blots of Figure 2D and the uncropped blot with the relevant bands clearly labelled as Figure 2D.: https://cdn.elifesciences.org/articles/83946/elife-83946-fig2-data1-v2.zip
Download elife-83946-fig2-data1-v2.zip

Next, we introduced 3xFlag and 6xHis tags to rcsD_pe and rcsD_pstb and analysed their expression by western blotting. A band in the tagged RcsD_pstb strain was observed at the expected size (Figure 2D, Lane3), but no visible band was detected in the tagged rcsD_pe strain (Figure 2D, Lane2). To enhance detection sensitivity, the rcsD_pe lysates were enriched before western blotting. An approximately 15 kDa band, corresponding to the predicted protein containing the HPt domain, was detected in both rcsD_pe, rcsD_C-term, and rcsD_pstb expressing strains (Figure 2D and Figure 2—figure supplement 1E, see Figure 2—source data 1, Figure 2—figure supplement 1—source data 1 for details), but not detected in the putative start codon mutated rcsD_pe expressing strain (Figure 2—figure supplement 1E, see Figure 2—figure supplement 1—source data 1 for details). In addition, full-length RcsD was detected in rcsD_pe but not in rcsD_C-term and rcsD_pe-stop (a stop codon was introduced upstream of the frameshift site) expressing strains (Figure 2D and Figure 2—figure supplement 1F, see Figure 2—source data 1, Figure 2—figure supplement 1—source data 2 for details), indicating that intact RcsD is weakly expressed in wild type Y. pestis through translational recoding (Rodnina et al., 2020). This is in accordance with the translational lacZ reporter system, which detected 1% readthrough once frameshifted (Figure 2—figure supplement 1G). Consistent with these results, Y. pestis with rcsD_pe-stop showed similar biofilm and CR phenotype as the rcsD_N-term deletion strain (Figure 2C and data not shown). Taken together, these data suggest that a 103-amino acid protein, designated as RcsD-Hpt, is expressed by rcsD_pe and functions as a negative regulator of biofilm formation in Y. pestis.

RcsD is a phosphorelay protein that can transfer phosphate from the conserved His residue in its HPt domain to a conserved Asp in the receiver domain of RcsB and can also dephosphorylate this site (Ancona et al., 2015; Takeda et al., 2001). An H844A mutation in the Y. pestis wild type or rcsD_N-term deletion strain showed a comparable phenotype to a full rcsD deletion strain (Figure 2B and C and Figure 1E), while plasmid-borne expression of the mutated version of rcsD_pe or rcsD-hpt in Y. pestis wild type and its rcsD mutants displayed a similar phenotype to an empty vector control (Figure 2—figure supplement 2A). In addition, plasmid-borne expression of the mutated version of rcsD_pstb in the Y. pestis rcsD_N-term deletion strain showed significantly less biofilm than that of rcsD_pstb (Figure 2—figure supplement 2B), indicating H844 might be important for the dephosphorylation ability of RcsD. Next, we modelled the structure of RcsD-Hpt with AlphaFold2 (Fowler and Williamson, 2022; Figure 2E), which revealed a similar structure to HptB, an HPt orphan protein in P. aeruginosa (Figure 2F; Chen et al., 2020). Like HptB, the predicted structure of RcsD-Hpt shows an elongated bundle of four helices α2, α3, α4, and α5, covered by the short N-terminal α1 helix. The imidazole side chain of the conserved active-site histidine residue His50 (His844 in RcsD_pstb) is located near the middle of helix α3 and protrudes from the bundle where it is exposed, as His57 is in HptB. Taken together, RcsD-Hpt may function as a classical HPt orphan protein, and the conserved His residue is crucial for its function.

A frameshift in rcsD alters Rcs signalling in Y. pestis

RcsF and IgaA, which regulate environmental stress sensing (Cho et al., 2014; Guo and Sun, 2017; Wall et al., 2018), transfer signals to RcsD through interaction of IgaA with its periplasmic domain (Wall et al., 2020). Given that RcsD-Hpt does not encode a functional periplasmic domain, we hypothesized that the roles of RcsF and IgaA are dispensable in the Y. pestis Rcs signalling system. We constructed rcsF and igaA deletion mutants in the Y. pestis wild type and rcsD_pstb substitution strains, respectively. In the Y. pestis rcsD_pstb substitution strain, deletion of rcsF increased CR adsorption and biofilm formation (Figure 3A and B), while deletion of igaA completely abolished biofilm formation and CR binding (Figure 3C and D). Furthermore, overexpression of RcsF but not its C125S mutant (as Escherichia coli RcsF cysteine mutants are inactive Rogov et al., 2011), decreased biofilm formation and CR adsorption (Figure 3A and B), while expression of IgaA but not its C413S mutant (the conserved cysteine is essential for the function of IgaA in Salmonella enterica Pucciarelli et al., 2017) complemented the phenotype of the igaA deletion strain (Figure 3C and D). These results indicate that when receiving signals transduced from RcsF and IgaA, RcsD_pstb mediates decreased phosphorylation of RcsB, and thus Rcs signalling is switched off. Alternately, deletion or ectopic expression of rcsF or igaA does not regulate biofilm formation and CR pigmentation in wild type Y. pestis KIM6+ (Figure 3A–D), suggesting a decoupling of the requirement for RcsF and IgaA and the Rcs system after the transition from rcsD_pstb to rcsD_pe.

Figure 3

Download asset Open asset

The frameshift mutation in *rcsD* alters the Rcs signalling pathway in Y. *pestis*.

Congo red (CR) pigmentation assay (**A and C**) and crystal violet (CV) biofilm assay (**B and D**) the derivatives of *Y. pestis* KIM6+ (wild type [WT]) and the *rcsD*_pstb substitution strain. prcsF, plasmid expressing *rcsF*; prcsF (C125S), plasmid expressing *rcsF* with cysteine (C) to serine (S) substitution at position 125; pigaA, plasmid expressing *igaA*; pigaA (C413S), plasmid expressing *igaA* with cysteine (C) to serine (S) substitution at position 413. (E) CV binding assay using WT and a *rcsC* deletion mutant expressing RcsD-Hpt (*rcsD-hpt*). CR pigmentation assay (**F and G**) and CV biofilm assay (H) using two *Y. pestis rcsC* deletion mutants expressing different *rcsC* variants. CV assays in panels B, D, and H were performed together. Error bars represent ± SD from three independent experiments with three replicates. Statistical analysis was performed using one-way analysis of variance (ANOVA) with Dunnett’s multiple comparisons post-test. ns, not significant; *p<0.05, **p<0.01, ***p<0.001.

RcsC, a bifunctional histidine kinase and phosphatase, phosphorylates and dephosphorylates RcsD, which subsequently, via RcsB, activates or represses expression of its target genes (Guo and Sun, 2017; Wall et al., 2018). Deletion of RcsC resulted in increased biofilm formation and CR pigmentation (Sun et al., 2008; Figure 3E), indicating that RcsC is involved in Rcs signalling in Y. pestis. In addition, high expression of RcsD-Hpt in the rcsC deletion mutant did not affect biofilm formation and pigmentation (Figure 3E), indicating that RcsC is required for RcsD-Hpt-mediated biofilm repression in Y. pestis. His489 and Asp885 in RcsC are the predicted autophosphorylation and subsequent transfer receipt sites, respectively (Clarke et al., 2002; Latasa et al., 2012). A D885A mutation in RcsC abolished its function in an rcsD_pstb substitute or wild type strain (Figure 3F–H), indicating that Asp885 was crucial for phosphate transfer in both conditions. An H489A mutation in RcsC did not alter its function in the wild type strain (Figure 3F and H), indicating that His489 is not important for the phosphorylation of RcsD-Hpt. Surprisingly, mutation of His489 resulted in decreased biofilm formation and CR adsorption in the rcsD_pstb substitution strain (Figure 3G and H), indicating that His489 might function as a phosphate reservoir involved in dephosphorylation of RcsD. RcsC T903A constitutively activates Rcs in S. enterica (García-Calderón et al., 2005; Latasa et al., 2012). The corresponding T913A mutation in RcsC strongly decreased biofilm formation and CR pigmentation in the RcsD_pstb substitution strain (Figure 3F–H), indicating that an RcsC T913A mutant stimulates Rcs in a similar manner to S. enterica. The same mutation in RcsC resulted in enhanced biofilm formation in Y. pestis, indicating that an RcsC T913A mutant has an impaired ability to phosphorylate RcsD-Hpt. Taken together, these data indicate that RcsC remains a crucial component of the Rcs phosphorelay system in Y. pestis.

RcsB is a phosphoacceptor in the Rcs system and contains a conserved Asp site in its receiver domain which can be phosphorylated by RcsD (Fredericks et al., 2006; Takeda et al., 2001). Phosphorylated RcsB regulates expression of biofilm-related genes (hmsT, hmsD, hmsP, and hmsHFRS) in Y. pestis (Fang et al., 2015; Guo et al., 2015; Sun et al., 2012). Expression of rcsD_pe or rcsD_pstb only conferred altered biofilm formation in the presence of RcsB (Figure 4A), indicating that RcsD_pe and RcsD_pstb might differentially modulate the phosphorylation of RcsB. To test this hypothesis, we detected the phosphorylation of RcsB using Phos-tag SDS-PAGE gels and western blotting (Madec et al., 2014). As expected, replacement of rcsD_pe by rcsD_pstb in Y. pestis resulted in decreased phosphorylation of RcsB, while mutation of the conserved Asp residue (D56Q) in RcsB abolished phosphorylation (Figure 4B, see Figure 4—source data 1 for details). Consistent with the phosphorylation status of RcsB, the transcription and expression of HmsT were differentially regulated by RcsD_pe and RcsD_pstb (Figure 4C and D, see Figure 4—source data 2 for details).

Figure 4

Download asset Open asset

The frameshift mutation in *rcsD* increases the phosphorylation of RcsB and represses the expression of HmsT.

(A) Crystal violet (CV) biofilm assays using a wild type (WT) (white), Δ*rcsB* (grey), or *rcsB* (D56Q) (black) mutant strain, each harbouring pUC19 vectors expressing *rcsD*_pe or *rcsD*_pstb. (B) Phosphorylation analysis of RcsB in the *Y. pestis* KIM6+ strain (WT) harbouring an RcsB expression plasmid (prcsB), a plasmid expressing RcsB fused with a 3xflag tag (prcsB-Flag), or a modified prcsB-Flag expression plasmid in which the conserved phosphorylation site Asp50 was mutated to Gln, and the *Y. pestis* rcsD_pstb substitution strain harbouring the prcsB-Flag expression plasmid (see Figure 4—source data 1 for details). (C) Quantification of HmsT expression using a β-galactosidase assay. The *lacZ* reporter gene was fused with the *hmsT* promoter in plasmid pGD926. (D) Expression of HmsT was analysed by western blotting using an anti-Flag antibody (see Figure 4—source data 2 for details). Error bars represent ± SD from three independent experiments with three replicates. Statistical analysis was performed using one-way analysis of variance (ANOVA) with Dunnett’s multiple comparisons post-test. ns, not significant; *p<0.05, **p<0.01, ***p<0.001.

Figure 4—source data 1 Raw source data for Figure 4B. The original file of the full raw unedited blots and the uncropped blot with the relevant bands clearly labelled as Figure 4B.: https://cdn.elifesciences.org/articles/83946/elife-83946-fig4-data1-v2.zip
Download elife-83946-fig4-data1-v2.zip
Figure 4—source data 2 Raw source data for Figure 4D. The original file of the full raw unedited blots and the uncropped blot with the relevant bands clearly labelled as Figure 4D.: https://cdn.elifesciences.org/articles/83946/elife-83946-fig4-data2-v2.zip
Download elife-83946-fig4-data2-v2.zip

The frameshift mutation in rcsD promotes retention of the pgm locus during Y. pestis flea infection

Loss of function in rcsA is a crucial step for Y. pestis to establish flea-borne transmission (Sun et al., 2008; Sun et al., 2014). We therefore speculated that mutation of rcsD might also play a role in the adaptation of Y. pestis to the flea. We therefore infected the Oriental rat flea, Xenopsylla cheopis, with Y. pestis wild type (KIM6+), the rcsD_pe::rcsD_pstb and the rcsD_N-term deletion strain. Bacterial burdens in infected fleas at 0, 7, and 28 day (s) post infection were not significantly different between any strain combinations (Figure 5A). In addition, we did not observe significant differences in flea blockage, despite different in vitro capacities to form biofilms in these strains (Figure 5B). Taken together, these results suggest that the frameshift mutation in rcsD does not alter the infection, persistence, and blockage-forming capacity of Y. pestis in fleas.

Figure 5

Download asset Open asset

The frameshift mutation in *rcsD* stabilizes the *pgm* locus in Y. *pestis*.

(A) Bacterial burdens in fleas infected with *Y. pestis* wild type (WT), rcsD_pe::*rcsD*_pstb and *rcsD*_N-term strains after 0, 7, and 14 days of infection. (B) Cumulative blockage of fleas after 4 weeks of infection with *Y. pestis* WT, rcsD_pe::*rcsD*_pstb and ΔrcsD_N-term strains. Two independent infection experiments are shown. (C) Percent of *pgm* locus loss in fleas infected with *Y. pestis* WT and rcsD_pe::*rcsD*_pstb after 4 weeks of infection. Ten infected fleas were used for this assay. Statistical analysis was performed using a Fisher’s exact test. (D) Percent of *pgm* locus loss *in vitro* with *Y. pestis* KIM6+ and mutants. Two-way analysis of variance (ANOVA) with Dunnett’s multiple comparisons were performed for statistical analysis of mutants with WT strain KIM6+. Error bars represent ± SD from three independent experiments with six replicates. ns, not significant; *p<0.05, **p<0.01, ***p<0.001.

We fortuitously observed that a Y. pestis rcsD_pstb strain, but not the wild type, displayed a pgm- phenotype after 4 weeks of flea infection. The pigmentation phenotype (pgm+) is defined by Y. pestis absorption of exogenous hemin or CR to form pigmented colonies. A spontaneous deletion of the 102 kb pgm locus imposed by the instability of two flanking IS elements results in a pgm- phenotype (Fetherston et al., 1992; Tong et al., 2005). Interestingly, an accumulation in pgm locus mutations was previously reported for strains exhibiting enhanced biofilm formation (Fetherston et al., 1992; Kirillina et al., 2004; Silva-Rohwer et al., 2021). Prior to the infectious blood meal, no colonies exhibited a pgm- phenotype when the wild type and rcsD_pstb strains were plated on CR agar plates. Twenty-eight days post infection, we observed the pgm- phenotype in bacteria isolated from all fleas infected with the rcsD_pe::rcsD_pstb strain. The mean percentage of isolates displaying the pgm- phenotype was 16.3%, ranging from 0.1% to 47.9%, whereas only one out of ten fleas infected with the wild type strain displayed this phenotype with a frequency of ~1.0% (Figure 5C). We further analysed the pgm mutation rate of Y. pestis in liquid medium. To mimic the environment in the flea, Y. pestis were grown in liquid medium for several days, where they remained in stationary phase before reinoculation. Consistent with the flea infection data, Y. pestis encoding rcsD_pstb displayed a significantly higher pgm- mutation rate relative to the wild type strain (Figure 5D). Strains with deletion of rcsD-Hpt or mutated rcsD_pe also exhibited increased pgm mutation rates (Figure 5D), indicating that RcsD-Hpt is important for stability of the pgm locus. This data agrees with previous reports of strains showing increased pgm locus loss concomitant with enhanced biofilm production levels (Fetherston et al., 1992; Kirillina et al., 2004; Silva-Rohwer et al., 2021). PCR analysis confirmed that the pgm- phenotype was caused by the spontaneous deletion of the pgm locus (data not shown). Taken together, these results suggest that the frameshift mutation in rcsD promoted stable maintenance of the pgm locus of Y. pestis KIM6+ in infected fleas.

Genome-wide identification of genes regulated by the Rcs phosphorelay system in Y. pestis

Rcs has been reported to modulate virulence in other bacteria (Li et al., 2015; Wall et al., 2018; Wang et al., 2012). KIM6+ is an avirulent derivative of the fully virulent KIM strain, which was cured of the pCD1 plasmid. To investigate the role of Rcs on pathogen virulence in a mammalian host, we took advantage of the Y. pestis biovar Microtus strain 201 (Zhang et al., 2014), a human-avirulent but rodent-virulent strain, isolated from a natural reservoir, the Brandt’s vole (Microtus brandti) (Song et al., 2004). Mutation of rcsB or rcsD in the Y. pestis biovar Microtus strain 201 showed a similar CR absorption and in vitro biofilm phenotype as the Y. pestis KIM6+ strain (data not shown). To our surprise, deletion of rcsB or replacement of rcsD_pe by rcsD_pstb in this strain did not significantly affect their virulence when mice were subcutaneously infected with different doses of bacteria (Figure 6A–C).

Figure 6 with 1 supplement see all

Download asset Open asset

Genome-wide identification of genes regulated by the Rcs phosphorelay system in *Y. pestis*.

(**A–C**) Survival of C57BL/6 mice infected with *Y. pestis* Microtus strain 201 and its derivatives using an infectious dose of 40, 160, and 800 colony-forming units (CFU). (D) Venn diagram of upregulated and downregulated genes in *Y. pestis* strains growing at temperatures (26°C or 37°C). (E) qPCR analysis of *hmsT* and differentially expressed genes (DEGs) identified by RNA-seq. The screening threshold for DEGs was defined as |logFC|≥1, and p≤0.05. Error bars indicate SD from at least three samples. (F) Heatmaps showing the differential expression of genes identified through clusters of orthologous group (COG) analysis.

Rcs has been reported to modulate the expression of many genes in response to environmental stress. To characterize the genes regulated by the Rcs system in Y. pestis, we performed RNA-seq on total RNA isolated from Y. pestis biovar Microtus strain 201, rcsB deletion and rcsD substitution strains cultured at 26°C and 37°C. A total of 139 genes (43 upregulated and 96 downregulated) and 49 genes (18 upregulated and 31 downregulated) were significantly differentially expressed (|logFC|≥1, p≤0.05) in the rcsB deletion and rcsD substitution strain, respectively, when compared with the wild type strain at 26°C (Figure 6D, Figure 6—figure supplement 1A–1B , and Supplementary file 4). At mammalian temperature (37°C), 37 genes (26 upregulated and 11 downregulated) and 177 genes (97 upregulated and 80 downregulated) were significantly differentially expressed (Figure 6D, Figure 6—figure supplement 1C–1D and Supplementary file 4). Several differentially expressed genes (DEGs) identified by RNA-seq were verified by quantitative real-time PCR (qRT-PCR), indicating comparable patterns of expression (Figure 6E and Supplementary file 4). Substitution of rcsD_pstb had an opposing effect on gene expression to deletion of rcsB (Figure 6F and Supplementary file 4), indicating that the frameshift mutation in rcsD lessens the regulatory function of Rcs. Furthermore, Rcs positively regulated genes such as those encoding the type 6 secretion system, and those related to biosynthesis of yersiniabactin and ABC transporter genes (Figure 6F). These observations suggest that Rcs might play an important role in the environmental fitness and virulence of Y. pestis, and thus could be required for the flea-mammalian host transmission cycle in the wild.

A frameshift mutation in rcsD is an evolutionary step present in modern Y. pestis lineages

To investigate the role of the rcsD mutation in the evolution of Y. pestis, we analysed the evolutionary changes that occurred during the divergence of Y. pestis from Y. pseudotuberculosis (Figure 7 and Supplementary file 5). Mutation of rcsD is present in all Y. pestis that harbour five genetic changes (pde3’, ymt, rcsA, pde2, and ureD) required for flea colonization (Figure 7), except for the 0.PE7 branch (Figure 7). An IS element in the ompC gene is one of two IS elements driving instability of the pgm locus (Figure 7 and Supplementary file 5; Fetherston et al., 1992; Tong et al., 2005). The IS element in ompC was present in most modern Y. pestis branches, but not in RT5, an ancient Y. pestis strain isolated from a bubonic plague patient during the Bronze Age (Spyrou et al., 2018). This indicates that emergence of the rcsD mutation is likely not due to loss of the pgm locus enabled by IS elements. Although multiple rcsD mutations are present across the phylogeny (Figure 7 and Supplementary file 5), the HPt encoding region is present in all sequenced Y. pestis isolates, indicating an important role of RcsD-Hpt in refining stable blockage-mediated flea-borne transmission of Y. pestis.

Figure 7

Download asset Open asset

Genetic changes in Rcs genes during speciation of *Y. pestis*.

Acquisition events are indicated by green lightning symbols, loss of genetic material by a red cross, and mutation by blue circles. All comparisons shown are relative to *Y. pseudotuberculosis*. Increasingly darker shades of grey represent *Y. pseudotuberculosis*, and ancient strains from the Iron Age and the LNBA lineage, respectively. Increasingly darker shades of green represent branches one to four as annotated in the figure. Rcs-related genes and IS_ompC were considered as intact (dark blue), mutated (light blue), or absent (white). *: not analyzed. 1, present in all *Y. pestis* except branch 0 strain 0.PE4b; 2, *rcsD*_N-term and *rcsD*_C-term were located in a different genome site in the Nairobi (1.ANT), Angola (0.PE3), and Algeria3 (ORI) strains due to chromosome rearrangement; 3, indels are present in *rcsD*_N-term in some strains I-3086 (0.PE4m) and A-1825 (2.MED1). MRCA, the most recent common ancestor. Figure 7 and nomenclature are adapted from Figure 1 of Demeure et al., 2019.

Discussion

Y. pestis and its ancestor Y. pseudotuberculosis have historically been studied as models for pathogen evolution, and have helped to shape our understanding of the evolutionary processes driving niche adaptation, transmission, and pathogenesis (Wren, 2003). A recent paleogenomic study has clarified major steps driving evolution of Y. pestis (Rasmussen et al., 2015; Spyrou et al., 2018). An ancestral Y. pseudotuberculosis, which has a mutation in the promoter region of pde3 (Sun et al., 2014), acquired two plasmids, pPla and pMT1, in addition to other genetic changes, to evolve into the ancient virulent Y. pestis (Bearden et al., 2009; Cui et al., 2013; Sodeinde et al., 1992; Zimbler et al., 2015). At this point, Y. pestis may still have been prevalent in the environment, where it could be transmitted to humans and animals by the faecal-oral route and occasionally by flea bites through early phase transmission (Figure 8—figure supplement 1). Later, other genetic changes, including acquisition of ymt, and mutations in rcsA, pde2, and ureD, occurred in the ancient Y. pestis, converting the pathogen into full competency for a flea-borne transmission modality (Chain et al., 2004; Cui et al., 2013; Hinnebusch, 2005; Figure 8—figure supplement 1).

Compared to its ancestor, which faced multiple changing environments, the establishment of a flea-mammalian host transmission cycle limited environmental exposure of Y. pestis. The progenitor lineage of flea-borne Y. pestis still required a series of genetic changes to repurpose its environmental signal sensing and transduction systems to adapt to its new lifestyle and niche (Hinnebusch et al., 2017; Yang et al., 2023). Loss of redundant genes and response pathways may have contributed to the fitness of Y. pestis. For example, mutation of flagella-related genes in Y. pestis occurred in parallel to its colonization of the flea (Minnich and Rohde, 2007). Analysis of a broad set of Yersinia genomes demonstrated that deletion of one thymine in rcsD occurred after the ancient Y. pestis acquired the major genetic changes required for flea-borne transmission (Figures 7 and 8). The frameshift mutation in rcsD might compensate for the fitness cost imposed by loss of rcsA function by dampening the subsequent drastic changes in gene expression.

Figure 8 with 1 supplement see all

Download asset Open asset

Predicted model of Rcs signal transduction and phosphoryl transfer in *Y. pseudotuberculosis* and *Y. pestis*.

In *Y. pseudotuberculosis*, full-length RcsD plays a dominant role while the lowly expressed RcsD-Hpt might play a moonlighting function. When Rcs is in ‘off’ state, RcsF is not activated leading to release from repression of IgaA, which interacts with RcsD. In this situation, RcsB, which might receive the phosphoryl group from another source such acetyl phosphate, is dephosphorylated by intact RcsD, which is then dephosphorylated by RcsC. When Rcs is in an ‘on’ state, activated RcsF interacts with IgaA, releasing the interaction of RcsD and IgaA. In this situation, RcsC transfers the phosphoryl group to RcsD, and then to RcsB, leading to an activated Rcs system. In *Y. pestis*, RcsD-Hpt might play a dominant role. Rcs system might be only slightly regulated by IgaA and RcsF pathway. In this situation, RcsC transfers the phosphoryl group to RcsD-Hpt, and then to RcsB, leading to an activated Rcs system. The arrows indicate direction of flow of phosphate and weight of the arrows correlates with the flow of amount of phosphate. Weight of the red inhibitor lines correlates with magnitude of inhibition. Phosphate is denoted as small circles with the letter P.

The frameshift present in rcsD_pe leads to the expression of two functional proteins: RcsD-Hpt and intact RcsD. Low levels of intact RcsD may be expressed by translational readthrough, while RcsD-Hpt is expressed from a rare AUU start codon. Intact RcsD and RcsD-Hpt might have different functions in Y. pestis (Figure 8 and Figure 8—figure supplement 1). Under normal conditions, RcsD dephosphorylates RcsB, while RcsD-Hpt phosphorylates RcsB (Figure 8). This subsequently promotes different capacities for biofilm formation and likely multiple other phenotypes. The periplasmic domain of RcsD receives environmental signals sensed by RcsF and IgaA, which in turn regulates the phosphorylation of RcsD by RcsC (Wall et al., 2018; Wall et al., 2020). RcsD-Hpt lost the ability to respond to environmental signal transduction by RcsF and RcsD but could still receive a phosphate group from RcsC (Figure 8). Although RcsD-Hpt and intact RcsD were expressed in Y. pestis, RcsD-Hpt appears to play a dominant role in regulation of the Rcs pathway (Figure 8). This hypothesis is supported by two observations: (1) expression of rcsD_pe conferred a similar phenotype as expression of rcsD-hpt, and (2) RcsF and IgaA modulate Rcs signalling in the rcsD_pstb substitution strain but not in wild type Y. pestis.

Although rcsD_pstb has the same RBS and start codon, only a very small amount of RcsD-Hpt relative to full-length RcsD was detected by western analysis, indicating intact RcsD plays a major role in the background of rcsD_pstb. Sequence analysis indicated that a putative start codon and RBS are present in rcsD in many organisms (Supplementary file 6). This indicates that RcsD-Hpt may play a moonlighting role in Rcs signalling in Y. pseudotuberculosis (Figure 8). Indeed, wild type RcsD in E. coli produces low levels of a short phosphotransfer protein (Rogov et al., 2004; Wall et al., 2020). HPt orphan proteins function as phosphate transfer components in multiple phosphorelay systems in numerous prokaryotes and eukaryotes (Hérivaux et al., 2018; Kennedy et al., 2016; Mohanan et al., 2017; Valentini et al., 2016), and have evolved from larger phosphotransferase proteins containing multiple domains. The frameshift present in rcsD of Y. pestis may represent an ongoing evolutionary process generating an orphan HPt protein and consequently a new regulatory pathway.

The rcsD frameshift alters the Rcs signalling pathway, which in turn decreases Y. pestis biofilm formation. The rcsD frameshift in Y. pestis does not significantly affect mammalian virulence and flea colonization in our study, but it may promote bacterial fitness during successive flea-mammal host transmission cycles. In agreement with previous findings that spontaneous deletion of the pgm locus is increased with enhanced biofilm formation (Fetherston et al., 1992; Podladchikova et al., 2002; Silva-Rohwer et al., 2021), our work shows that increased biofilm formation in an rcsD mutant promotes biofilm formation with a subsequent increase in pgm loss. The pgm locus harbours the ybt operon, which is involved in iron acquisition and is required for virulence of Y. pestis (Fetherston et al., 2010; Sebbane et al., 2010). Indeed, Y. pestis strains lacking the pgm locus are avirulent and have been used as live plague vaccines in some countries (Bearden and Perry, 1999; Podladchikova et al., 2002). It is notable that prolonged exposure to biofilm-stimulating conditions in the flea gut and during in vitro growth resulted in exacerbated loss of the pgm locus (Figure 5C and D). This may explain why no effect on virulence was noted for the rcsD_pe::rcsD_pstb strain in our mouse infection studies which utilized strains cultured at 37°C overnight. Enhanced biofilm production through the rcsA mutation in the absence of accompanying rcsD mutation in ancient Y. pestis strains likely caused a high rate of pgm locus loss during flea infection thus conferring vaccine-like protection when transferred to the mammalian host. Eventually this would break the transmission cycle between flea and mammalian host. In modern lineages the rcsD mutation therefore serves to dampen biofilm production without obvious compromise to flea blockage and infection rates. Thus, stable maintenance of the pgm locus required to amplify infection in the mammalian host and perpetuate the flea-mammal transmission cycle, and intensity of plague outbreaks, is promoted (Figure 7). Finally, the frameshift mutation of rcsD might represent an important step in the emergence of extant ubiquitous lineages of Y. pestis.

Materials and methods

Bacterial strains and plasmids

Request a detailed protocol

This study utilized Y. pestis strain KIM6+, which is derived from the sequenced strain KIM strain (Deng et al., 2002), but is cured of the pCD1/pYV plasmid required for mammalian virulence, and is competent for flea blockage (Hinnebusch et al., 1996) and biofilm formation (Darby et al., 2002; Sun et al., 2008). Studies in mice utilized the biovar Microtus strain 201 which was avirulent to human but fully virulent to mice (Zhang et al., 2014).

Deletion of rcsC, rcsF, and igaA was achieved using a one-step method to integrate PCR products into the chromosome with pKD46, as previously described (Datsenko and Wanner, 2000). Double mutants were made by the sequential application of pKD46-mediated deletion. CRISPR-Cas12a-assisted recombineering were used to introduce point mutations, deletions, insertions, and gene replacements in this study (Yan et al., 2017). All strains were verified by PCR, DNA sequencing, and plasmid complementation.

For construction of the hmsT::lacZ reporter, 350 bp of hmsT upstream sequence, together with the first seven codons of the ORF, were amplified by PCR using KIM6+ chromosomal DNA as the template. The DNA fragments were digested with HindIII and BamHI restriction enzymes and cloned into pGD926 (Ditta et al., 1985; Vieille and Elmerich, 1990), generating plasmids pYC593 and pYC287. Plasmids expressing rcsC (D885A), rcsC (T913A), rcsD (H844A), rcsD::rcsD-3xflag, and rcsD-3xflag::rcsD-3xflag-his6 were generated by overlapping PCR as described previously (Li et al., 2015).

For inducible rcsD expression, the gene was cloned downstream of the arabinose-inducible promoter of plasmid pBAD/Myc-His (Invitrogen). The plasmid used for determining readthrough was generated by cloning a partial sequence of rcsD (159 bp for rcsD_pstb and 158 bp for rcsD_pe) containing the 8T (frameshifted region) into pMal-lacZ, generating plasmids pMal-rcsD_pstb-lacZ and pMal-rcsD_pe-lacZ. An additional stop codon was introduced into the 158 bp rcsD_pe sequence of pMal-rcsD_pe-lacZ, generating plasmid pMal-rcsD_pe-stop-lacZ.

All strains and plasmids used in this study are shown in Supplementary file 1 and oligonucleotides used in this study are shown in Supplementary file 3.

In vitro biofilms

Request a detailed protocol

Microtiter plate biofilm assays were performed as previously described (Sun et al., 2012). Briefly, bacteria were cultured overnight in LB broth supplemented with 4 mmol CaCl₂ and 4 mmol MgCl₂. Cultures were subsequently diluted into 96-well plates and incubated with shaking for 24 hr at 26°C. The wells were washed, and the adherent biofilm was stained with crystal violet, solubilized with 80% ethanol and 20% acetone, and measured by A₆₀₀. Results are from three independent experiments with three technical replicates per experiment.

β-Galactosidase assays

Request a detailed protocol

β-Galactosidase activities were measured as previously described (Guo et al., 2015; Sun et al., 2012). Briefly, overnight cultures of Y. pestis harbouring lacZ reporters were diluted to an OD₆₀₀ of 0.05 and grown in LB broth at room temperature to an OD₆₀₀ of 1.5. The active β-galactosidase, encoded by the lacZ reporter gene in Y. pestis strains, can cleave o-nitrophenyl-β-D-galactopyranoside (ONPG) substrate to a bright yellow product. The cells were lysed and ONPG solution were added. After incubation at 37°C, the reaction was stopped by adding 1 M Na₂CO₃, then absorbance was measured at 420 nm. Results were normalized against cell density and incubation time, and shown in Miller units (Miller, 1972). At least two independent experiments with technical triplicates were performed.

Quantitative real-time PCR

Request a detailed protocol

qRT-PCR was carried out as previously described (Sun et al., 2011). Briefly, cells were first grown in LB broth overnight before diluting to an OD₆₀₀ of 0.05 in LB and incubating at room temperature to an OD₆₀₀ of 0.8. Total RNA was isolated using the Rneasy Mini Kit (Qiagen). Residual DNA was removed by treatment with rDNase I (Ambion) and confirmed by PCR. cDNA was synthesized from the RNA and used for quantitative PCR on an Applied Biosystems unit (Quant Studio 5). The quantity of mRNA was normalized relative to the reference gene 16sRNA (YP_r1). The relative mRNA expression levels in each strain were normalized to the wild type samples. Primers and probe sets used in this study are listed in Supplementary file 3. Results from three independent experiments performed in technical triplicate were analyzed by one-way analysis of variance (ANOVA) with Bonferroni’s test.

Western blotting

Request a detailed protocol

Western blotting was performed as previously described (Ren et al., 2017). For detection of enriched RcsD_pe, Y. pestis strains expressing rcsD fused with 3xFlag-His6 were grown at 26°C to stationary phase. Cells were harvested by centrifugation and disrupted by sonication. The protein was enriched by Ni-nitrilotriacetic acid His resin for western blot analysis. For detection of enriched RcsD and RcsD-Hpt in Figure 2—figure supplement 1E and F, Y. pestis strains expressing rcsD fused with 3xFlag were grown at 26°C to stationary phase. Cells were harvested by centrifugation and disrupted by sonication. The protein was enriched by flag resin for western blot analysis. For detection of HmsT, Y. pestis strains harbouring a plasmid expressing HmsT with 3xFlag were grown at 26°C to stationary phase. Cells were harvested by centrifugation and disrupted by sonication. Approximately 20 ng proteins were loaded for detection of HmsT. The proteins were separated on 10% SDS-PAGE gels transferred to PVDF membranes (Millipore), analyzed by immunoblotting with an anti-Flag antibody produced by Invitrogen (Catalog number: MA1-91878-HRP; RRID: AB_2537626), and detected with ECL Western Detection Reagents (Bio-Rad). Resulting bands were quantitated by densitometry using NIH ImageJ (Gallo-Oller et al., 2018).

Phos-tag SDS-PAGE

Request a detailed protocol

For detection of protein phosphorylation, acrylamide gel was mixed with 25 μM Phos-tag acrylamide (AAL-107, Wako) and 25 μM MnCl₂ (Madec et al., 2014). Cells grown to an OD₆₀₀ of 0.8 were centrifuged and pellets were resuspended in PBS containing Protease Inhibitor Cocktail (Roche) and then lysed with a sonicator. Samples were quickly loaded onto gels containing 25 µM Phos-tag acrylamide and 25 µM MnCl₂ and run. After a 10 min wash with WB transfer buffer supplied with 1 mM EDTA, followed by a 10 min wash with transfer buffer without EDTA. Proteins were transferred to a PVDF membrane and blotted using an anti-Flag antibody produced by Invitrogen (Catalog number: MA1-91878-HRP; RRID: AB_2537626).

RNA-seq

Request a detailed protocol

Total RNA was extracted using the RNeasy Kit (Qiagen). RNA-seq and expression quantification were performed by Genewiz. Gene expression levels were further normalized using the fragments per kilobase of transcript per million mapped reads method to eliminate the influence of different gene lengths and sequencing depth (Wang et al., 2009). The edgeR package was used to identify DEGs across samples with fold changes ≥2 and a false discovery rate-adjusted p (p-value)≤0.05 (Anders and Huber, 2010; Robinson et al., 2010). DEGs were then subjected to an enrichment analysis of GO function and KEGG pathways (Harris et al., 2004; Kanehisa and Goto, 2000).

Murine infection

Request a detailed protocol

Animals were handled in strict accordance with the Guidelines for the Welfare and Ethics of Laboratory Animals of China and all the animal experiments were approved by the Institutional Animal Care Committee of Military Medical Sciences. Bacterial cultures at 37°C were washed twice with PBS (pH 7.2) and then subjected to serial 10-fold dilutions with PBS. Dilutions were plated onto brain heart infusion (BHI) agar plates to calculate the numbers of colony-forming units (CFU). For each strain, different doses of bacterial suspension were inoculated subcutaneously at the inguinal region of 10 female BALB/c mice (aged 6–8 weeks), which were obtained from Charles River Laboratories (Beijing, China). Survival was monitored at regular intervals, and a survival curve was generated with GraphPad Prism 5.0. p-Values were determined using the log-rank (Mantel–Cox) test and the Gehan–Breslow–Wilcoxon test; p<0.01 was considered statistically significant.

Flea blockage

Request a detailed protocol

Flea infections and blockage analysis were carried out as previously described (Silva-Rohwer et al., 2021; Sun et al., 2014). Y. pestis strains were grown overnight in 3 mL HIB at 26°C with shaking, then diluted into 100 mL HIB to cultivate at 37°C without shaking. The following day harvested cells were suspended in sterile PBS, and optical density at absorbance of 600 nm was determined. A commercial preparation of heparinized mouse blood (BioIVT, New York) was inoculated with Y. pestis to a final concentration of CFU/mL of ~5 × 10⁸ to 1 × 10⁹. X. cheopis fleas were allowed to feed on the infected blood through a mouse skin membrane. Studies with mice were performed in strict accordance with the U.S. National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals (National Research Council Committee for the Update of the Guide for the Care and Use of Laboratory Animals et al., 2011) and as approved by the Washington State University Institutional Animal Care and Use Committee.

Mice: A CD-1 mouse breeding colony originally sourced from Envigo (https://www.envigo.com/model/hsd-icr-cd-1) is maintained at WSU. Males and females are used. Neonates between the ages of 2–6 days are used for feeding fleas for breeding and maintenance of infected fleas.

Fleas: X. cheopis fleas are maintained at WSU since 2010 in Dr. Vadyvaloo’s lab. These fleas were originally sourced from Dr Joseph Hinnebusch’s lab at the NIH. Males and females were used for experiments.

Analysis of pigmentation

Request a detailed protocol

For in vivo pigmentation phenotype detection, 10 fleas infected with Y. pestis wild type or the rcsD_pe::rcsD_pstb strain were collected at the end of the infection period (T=28 days). Fleas were individually triturated in sterile PBS, and the fractions were serially diluted and cultivated on a CR plate. Two days later, CFU enumeration relevant to the pigmentation phenotype was performed. For in vitro pigmentation assays, fresh Y. pestis KIM6+ strains were inoculated into BHI medium and grown at 26°C for 24 hr. Cultures were diluted 1:100 into fresh BHI medium every 3 days. The cultures at 0, 1, 4, and 7 days post inoculation were plated on CR plates for analysis. The loss of the pgm locus was confirmed by PCR using two sets of primers described in Supplementary file 3. One set of primers targets the hmsS gene, which can obtain PCR products from the wild type strain but not the pgm locus deletion mutant. Another set of primers target the upstream and downstream regions of the pgm locus, and can be used to obtain PCR products from the pgm locus deletion mutant but not the wild type strain.

Quantification and statistical analysis

Request a detailed protocol

Figure legends detail the quantification and statistical analyses methods. We conducted the statistical analyses by GraphPad Prism.

Data availability

All data is available within the paper, its Supporting Information files, and the NCBI GenBank. RNA-seq sequencing data can be accessed in NCBI GenBank using BioProject ID: PRJNA876755. Source data files have been provided for Figures 2D, 4B, 4D, Figure 2—figure supplement 1E and F.

The following data sets were generated

1. Guo X-P
(2022) NCBI BioProject
ID PRJNA876755. Differential expression of Y. pestis genes in Y. pestis mutant strains growing at 26 degrees C and 37 degrees C.

https://www.ncbi.nlm.nih.gov/bioproject/PRJNA876755

References

(2010) Analysis of HmsH and its role in plague biofilm formation
Microbiology 156:1424–1438.

https://doi.org/10.1099/mic.0.036640-0
- Google Scholar
(2015) Conserved aspartate and lysine residues of RcsB are required for amylovoran biosynthesis, virulence, and DNA binding in Erwinia amylovora
Molecular Genetics and Genomics 290:1265–1276.

https://doi.org/10.1007/s00438-015-0988-8
- PubMed
- Google Scholar
1. Anders S
2. Huber W
(2010) Differential expression analysis for sequence count data
Genome Biology 11:R106.

https://doi.org/10.1186/gb-2010-11-10-r106
- PubMed
- Google Scholar
1. Bearden SW
2. Perry RD
(1999) The Yfe system of Yersinia pestis transports iron and manganese and is required for full virulence of plague
Molecular Microbiology 32:403–414.

https://doi.org/10.1046/j.1365-2958.1999.01360.x
- PubMed
- Google Scholar
1. Bearden SW
2. Sexton C
3. Pare J
4. Fowler JM
5. Arvidson CG
6. Yerman L
7. Viola RE
8. Brubaker RR
(2009) Attenuated enzootic (pestoides) isolates of Yersinia pestis express active aspartase
Microbiology 155:198–209.

https://doi.org/10.1099/mic.0.021170-0
- PubMed
- Google Scholar
(2018) Infectious blood source alters early foregut infection and regurgitative transmission of Yersinia pestis by rodent fleas
PLOS Pathogens 14:e1006859.

https://doi.org/10.1371/journal.ppat.1006859
- Google Scholar
(2021) Acquisition of Yersinia murine toxin enabled Yersinia pestis to expand the range of mammalian hosts that sustain flea-borne plague
PLOS Pathogens 17:e1009995.

https://doi.org/10.1371/journal.ppat.1009995
- Google Scholar
1. Bobrov AG
2. Kirillina O
3. Ryjenkov DA
4. Waters CM
5. Price PA
6. Fetherston JD
7. Mack D
8. Goldman WE
9. Gomelsky M
10. Perry RD
(2011) Systematic analysis of cyclic di-GMP signalling enzymes and their role in biofilm formation and virulence in Yersinia pestis
Molecular Microbiology 79:533–551.

https://doi.org/10.1111/j.1365-2958.2010.07470.x
- PubMed
- Google Scholar
1. Bosio CF
2. Jarrett CO
3. Scott DP
4. Fintzi J
5. Hinnebusch BJ
6. Monack DM
(2020) Comparison of the transmission efficiency and plague progression dynamics associated with two mechanisms by which fleas transmit Yersinia pestis
PLOS Pathogens 16:e1009092.

https://doi.org/10.1371/journal.ppat.1009092
- Google Scholar
1. Cao X
2. Slavoff SA
(2020) Non-AUG start codons: expanding and regulating the small and alternative ORFeome
Experimental Cell Research 391:111973.

https://doi.org/10.1016/j.yexcr.2020.111973
- PubMed
- Google Scholar
1. Cao S
2. Jiao Y
3. Jiang W
4. Wu Y
5. Qin S
6. Ren Y
7. You Y
8. Tan Y
9. Guo X
10. Chen H
11. Zhang Y
12. Wu G
13. Wang T
14. Zhou Y
15. Song Y
16. Cui Y
17. Shao F
18. Yang R
19. Du Z
(2022) Subversion of GBP-mediated host defense by E3 ligases acquired during Yersinia pestis evolution
Nature Communications 13:4526.

https://doi.org/10.1038/s41467-022-32218-y
- Google Scholar
1. Chain PSG
2. Carniel E
3. Larimer FW
4. Lamerdin J
5. Stoutland PO
6. Regala WM
7. Georgescu AM
8. Vergez LM
9. Land ML
10. Motin VL
11. Brubaker RR
12. Fowler J
13. Hinnebusch J
14. Marceau M
15. Medigue C
16. Simonet M
17. Chenal-Francisque V
18. Souza B
19. Dacheux D
20. Elliott JM
21. Derbise A
22. Hauser LJ
23. Garcia E
(2004) Insights into the evolution of Yersinia pestis through whole-genome comparison with Yersinia pseudotuberculosis
PNAS 101:13826–13831.

https://doi.org/10.1073/pnas.0404012101
- PubMed
- Google Scholar
1. Chen SK
2. Guan HH
3. Wu PH
4. Lin LT
5. Wu MC
6. Chang HY
7. Chen NC
8. Lin CC
9. Chuankhayan P
10. Huang YC
11. Lin PJ
12. Chen CJ
(2020) Structural insights into the histidine-containing phosphotransfer protein and receiver domain of sensor histidine kinase suggest a complex model in the two-component regulatory system in Pseudomonas aeruginosa
IUCrJ 7:934–948.

https://doi.org/10.1107/S2052252520009665
- Google Scholar
1. Cho SH
2. Szewczyk J
3. Pesavento C
4. Zietek M
5. Banzhaf M
6. Roszczenko P
7. Asmar A
8. Laloux G
9. Hov AK
10. Leverrier P
11. Van der Henst C
12. Vertommen D
13. Typas A
14. Collet JF
(2014) Detecting envelope stress by monitoring β-barrel assembly
Cell 159:1652–1664.

https://doi.org/10.1016/j.cell.2014.11.045
- PubMed
- Google Scholar
1. Clarke DJ
2. Joyce SA
3. Toutain CM
4. Jacq A
5. Holland IB
(2002) Genetic analysis of the RcsC sensor kinase from Escherichia coli K-12
Journal of Bacteriology 184:1204–1208.

https://doi.org/10.1128/jb.184.4.1204-1208.2002
- Google Scholar
1. Clarke DJ
(2010) The Rcs phosphorelay: more than just a two-component pathway
Future Microbiology 5:1173–1184.

https://doi.org/10.2217/fmb.10.83
- PubMed
- Google Scholar
Book
(2011) Guide for the Care and Use of Laboratory Animals
Washington, D.C: National Academies Press.

https://doi.org/10.17226/25801
- Google Scholar
1. Cui Y
2. Yu C
3. Yan Y
4. Li D
5. Li Y
6. Jombart T
7. Weinert LA
8. Wang Z
9. Guo Z
10. Xu L
11. Zhang Y
12. Zheng H
13. Qin N
14. Xiao X
15. Wu M
16. Wang X
17. Zhou D
18. Qi Z
19. Du Z
20. Wu H
21. Yang X
22. Cao H
23. Wang H
24. Wang J
25. Yao S
26. Rakin A
27. Li Y
28. Falush D
29. Balloux F
30. Achtman M
31. Song Y
32. Wang J
33. Yang R
(2013) Historical variations in mutation rate in an epidemic pathogen, Yersinia pestis
PNAS 110:577–582.

https://doi.org/10.1073/pnas.1205750110
- PubMed
- Google Scholar
1. Darby C
2. Hsu JW
3. Ghori N
4. Falkow S
(2002) Caenorhabditis elegans: plague bacteria biofilm blocks food intake
Nature 417:243–244.

https://doi.org/10.1038/417243a
- PubMed
- Google Scholar
1. Datsenko KA
2. Wanner BL
(2000) One-step inactivation of chromosomal genes in Escherichia coli K-12 using PCR products
PNAS 97:6640–6645.

https://doi.org/10.1073/pnas.120163297
- PubMed
- Google Scholar
(2019) Yersinia pestis and plague: an updated view on evolution, virulence determinants, immune subversion, vaccination, and diagnostics
Genes and Immunity 20:357–370.

https://doi.org/10.1038/s41435-019-0065-0
- PubMed
- Google Scholar
1. Deng W
2. Burland V
3. Plunkett G
4. Boutin A
5. Mayhew GF
6. Liss P
7. Perna NT
8. Rose DJ
9. Mau B
10. Zhou S
11. Schwartz DC
12. Fetherston JD
13. Lindler LE
14. Brubaker RR
15. Plano GV
16. Straley SC
17. McDonough KA
18. Nilles ML
19. Matson JS
20. Blattner FR
21. Perry RD
(2002) Genome sequence of Yersinia pestis KIM
Journal of Bacteriology 184:4601–4611.

https://doi.org/10.1128/JB.184.16.4601-4611.2002
- PubMed
- Google Scholar
1. Dewitte A
2. Bouvenot T
3. Pierre F
4. Ricard I
5. Pradel E
6. Barois N
7. Hujeux A
8. Bontemps-Gallo S
9. Sebbane F
10. Mecsas J
(2020) A refined model of how Yersinia pestis produces a transmissible infection in its flea vector
PLOS Pathogens 16:e1008440.

https://doi.org/10.1371/journal.ppat.1008440
- Google Scholar
1. Ditta G
2. Schmidhauser T
3. Yakobson E
4. Lu P
5. Liang XW
6. Finlay DR
7. Guiney D
8. Helinski DR
(1985) Plasmids related to the broad host range vector, pRK290, useful for gene cloning and for monitoring gene expression
Plasmid 13:149–153.

https://doi.org/10.1016/0147-619x(85)90068-x
- PubMed
- Google Scholar
1. D’Lima NG
2. Khitun A
3. Rosenbloom AD
4. Yuan P
5. Gassaway BM
6. Barber KW
7. Rinehart J
8. Slavoff SA
(2017) Comparative proteomics enables identification of nonannotated cold shock proteins in E. coli
Journal of Proteome Research 16:3722–3731.

https://doi.org/10.1021/acs.jproteome.7b00419
- PubMed
- Google Scholar
(2007) Early-Phase transmission of Yersinia pestis by unblocked Xenopsylla cheopis (siphonaptera: pulicidae) is as efficient as transmission by blocked fleas
Journal of Medical Entomology 44:678–682.

https://doi.org/10.1093/jmedent/44.4.678
- Google Scholar
1. Fang N
2. Yang H
3. Fang H
4. Liu L
5. Zhang Y
6. Wang L
7. Han Y
8. Zhou D
9. Yang R
(2015) RcsAB is a major repressor of Yersinia biofilm development through directly acting on hmsCDE, hmsT and hmsHFRS
Scientific Reports 5:9566.

https://doi.org/10.1038/srep09566
- Google Scholar
(1992) Loss of the pigmentation phenotype in Yersinia pestis is due to the spontaneous deletion of 102 kb of chromosomal DNA which is flanked by a repetitive element
Molecular Microbiology 6:2693–2704.

https://doi.org/10.1111/j.1365-2958.1992.tb01446.x
- PubMed
- Google Scholar
(2010) The yersiniabactin transport system is critical for the pathogenesis of bubonic and pneumonic plague
Infection and Immunity 78:2045–2052.

https://doi.org/10.1128/IAI.01236-09
- PubMed
- Google Scholar
1. Fowler NJ
2. Williamson MP
(2022) The accuracy of protein structures in solution determined by alphaFold and NMR
Structure 30:925–933.

https://doi.org/10.1016/j.str.2022.04.005
- PubMed
- Google Scholar
(2006) Acetyl phosphate-sensitive regulation of flagellar biogenesis and capsular biosynthesis depends on the Rcs phosphorelay
Molecular Microbiology 61:734–747.

https://doi.org/10.1111/j.1365-2958.2006.05260.x
- PubMed
- Google Scholar
(2018) A new background subtraction method for Western blot densitometry band quantification through image analysis software
Journal of Immunological Methods 457:1–5.

https://doi.org/10.1016/j.jim.2018.03.004
- PubMed
- Google Scholar
(2005) Virulence attenuation in Salmonella enterica rcsC mutants with constitutive activation of the Rcs system
Microbiology 151:579–588.

https://doi.org/10.1099/mic.0.27520-0
- Google Scholar
1. Guo XP
2. Ren GX
3. Zhu H
4. Mao XJ
5. Sun YC
(2015) Differential regulation of the hmsCDE operon in Yersinia pestis and Yersinia pseudotuberculosis by the Rcs phosphorelay system
Scientific Reports 5:8412.

https://doi.org/10.1038/srep08412
- PubMed
- Google Scholar
1. Guo XP
2. Sun YC
(2017) New insights into the non-orthodox two component Rcs phosphorelay system
Frontiers in Microbiology 8:2014.

https://doi.org/10.3389/fmicb.2017.02014
- Google Scholar
1. Harris MA
2. Clark J
3. Ireland A
4. Lomax J
5. Ashburner M
6. Foulger R
7. Eilbeck K
8. Lewis S
9. Marshall B
10. Mungall C
11. Richter J
12. Rubin GM
13. Blake JA
14. Bult C
15. Dolan M
16. Drabkin H
17. Eppig JT
18. Hill DP
19. Ni L
20. Ringwald M
21. Balakrishnan R
22. Cherry JM
23. Christie KR
24. Costanzo MC
25. Dwight SS
26. Engel S
27. Fisk DG
28. Hirschman JE
29. Hong EL
30. Nash RS
31. Sethuraman A
32. Theesfeld CL
33. Botstein D
34. Dolinski K
35. Feierbach B
36. Berardini T
37. Mundodi S
38. Rhee SY
39. Apweiler R
40. Barrell D
41. Camon E
42. Dimmer E
43. Lee V
44. Chisholm R
45. Gaudet P
46. Kibbe W
47. Kishore R
48. Schwarz EM
49. Sternberg P
50. Gwinn M
51. Hannick L
52. Wortman J
53. Berriman M
54. Wood V
55. de la Cruz N
56. Tonellato P
57. Jaiswal P
58. Seigfried T
59. White R
60. Gene Ontology Consortium
(2004) The Gene Ontology (GO) database and informatics resource
Nucleic Acids Research 32:D258–D261.

https://doi.org/10.1093/nar/gkh036
- PubMed
- Google Scholar
(2018) Progressive loss of hybrid histidine kinase genes during the evolution of budding yeasts (Saccharomycotina)
Current Genetics 64:841–851.

https://doi.org/10.1007/s00294-017-0797-1
- PubMed
- Google Scholar
1. Hinchliffe SJ
2. Isherwood KE
3. Stabler RA
4. Prentice MB
5. Rakin A
6. Nichols RA
7. Oyston PCF
8. Hinds J
9. Titball RW
10. Wren BW
(2003) Application of DNA microarrays to study the evolutionary genomics of Yersinia pestis and Yersinia pseudotuberculosis
Genome Research 13:2018–2029.

https://doi.org/10.1101/gr.1507303
- PubMed
- Google Scholar
(1996) Role of the Yersinia pestis hemin storage (hms) locus in the transmission of plague by fleas
Science 273:367–370.

https://doi.org/10.1126/science.273.5273.367
- PubMed
- Google Scholar
1. Hinnebusch BJ
(1997) Bubonic plague: a molecular genetic case history of the emergence of an infectious disease
Journal of Molecular Medicine 75:645–652.

https://doi.org/10.1007/s001090050148
- PubMed
- Google Scholar
1. Hinnebusch BJ
(2005)
The evolution of flea-borne transmission in Yersinia pestis

Current Issues in Molecular Biology 7:197–212.
- PubMed
- Google Scholar
1. Hinnebusch BJ
2. Erickson DL
(2008) Yersinia pestis biofilm in the flea vector and its role in the transmission of plague
Current Topics in Microbiology and Immunology 322:229–248.

https://doi.org/10.1007/978-3-540-75418-3_11
- PubMed
- Google Scholar
(2017) “Fleaing” the plague: adaptations of Yersinia pestis to its insect vector that lead to transmission
Annual Review of Microbiology 71:215–232.

https://doi.org/10.1146/annurev-micro-090816-093521
- PubMed
- Google Scholar
(2021) Structure-based analyses of Salmonella RcsB variants unravel new features of the Rcs regulon
Nucleic Acids Research 49:2357–2374.

https://doi.org/10.1093/nar/gkab060
- PubMed
- Google Scholar
1. Kanehisa M
2. Goto S
(2000) KEGG: kyoto encyclopedia of genes and genomes
Nucleic Acids Research 28:27–30.

https://doi.org/10.1093/nar/28.1.27
- PubMed
- Google Scholar
(2016) Extended N-terminal region of the essential phosphorelay signaling protein Ypd1 from Cryptococcus neoformans contributes to structural stability, phosphostability and binding of calcium ions
FEMS Yeast Research 16:fow068.

https://doi.org/10.1093/femsyr/fow068
- PubMed
- Google Scholar
(2004) HmsP, a putative phosphodiesterase, and HmsT, a putative diguanylate cyclase, control Hms-dependent biofilm formation in Yersinia pestis
Molecular Microbiology 54:75–88.

https://doi.org/10.1111/j.1365-2958.2004.04253.x
- PubMed
- Google Scholar
(2012) Salmonella biofilm development depends on the phosphorylation status of RcsB
Journal of Bacteriology 194:3708–3722.

https://doi.org/10.1128/JB.00361-12
- Google Scholar
1. Li Y
2. Hu Y
3. Francis MS
4. Chen S
(2015) RcsB positively regulates the Yersinia Ysc-Yop type III secretion system by activating expression of the master transcriptional regulator LcrF
Environmental Microbiology 17:1219–1233.

https://doi.org/10.1111/1462-2920.12556
- PubMed
- Google Scholar
(2014) Increased phosphorylation of the RcsB regulator of the RcsCDB phosphorelay in strains of Dickeya dadantii devoid of osmoregulated periplasmic glucans revealed by Phos-Tag gel analysis
Microbiology 160:2763–2770.

https://doi.org/10.1099/mic.0.081273-0
- Google Scholar
1. McNally A
2. Thomson NR
3. Reuter S
4. Wren BW
(2016) “Add, stir and reduce”: Yersinia spp. as model bacteria for pathogen evolution
Nature Reviews. Microbiology 14:177–190.

https://doi.org/10.1038/nrmicro.2015.29
- PubMed
- Google Scholar
Book
1. Miller JH
(1972)
Experiments in Molecular Genetics

New York: Cold Spring Harbor Laboratory Press.
- Google Scholar
1. Minnich SA
2. Rohde HN
(2007) A rationale for repression and/or loss of motility by pathogenic Yersinia in the mammalian host
Advances in Experimental Medicine and Biology 603:298–310.

https://doi.org/10.1007/978-0-387-72124-8_27
- PubMed
- Google Scholar
(2017) Fungal Histidine Phosphotransferase plays a crucial role in Photomorphogenesis and pathogenesis in Magnaporthe oryzae
Frontiers in Chemistry 5:31.

https://doi.org/10.3389/fchem.2017.00031
- PubMed
- Google Scholar
(2016) Interaction of the RcsB response regulator with auxiliary transcription regulators in Escherichia coli
The Journal of Biological Chemistry 291:2357–2370.

https://doi.org/10.1074/jbc.M115.696815
- PubMed
- Google Scholar
(2002)
Study of PGM mutation mechanism in Yersinia pestis (plague pathogen) vaccine strain EV76

Molekuliarnaia Genetika, Mikrobiologiia i Virusologiia 1:14–19.
- PubMed
- Google Scholar
(2017) A disulfide bond in the membrane protein IgaA is essential for repression of the RcsCDB system
Frontiers in Microbiology 8:2605.

https://doi.org/10.3389/fmicb.2017.02605
- PubMed
- Google Scholar
1. Putzker M
2. Sauer H
3. Sobe D
(2001)
Plague and other human infections caused by Yersinia species

Clinical Laboratory 47:453–466.
- PubMed
- Google Scholar
(2019) Emergence and spread of basal lineages of Yersinia pestis during the Neolithic decline
Cell 176:295–305.

https://doi.org/10.1016/j.cell.2018.11.005
- PubMed
- Google Scholar
1. Rasmussen S
2. Allentoft ME
3. Nielsen K
4. Orlando L
5. Sikora M
6. Sjögren K-G
7. Pedersen AG
8. Schubert M
9. Van Dam A
10. Kapel CMO
11. Nielsen HB
12. Brunak S
13. Avetisyan P
14. Epimakhov A
15. Khalyapin MV
16. Gnuni A
17. Kriiska A
18. Lasak I
19. Metspalu M
20. Moiseyev V
21. Gromov A
22. Pokutta D
23. Saag L
24. Varul L
25. Yepiskoposyan L
26. Sicheritz-Pontén T
27. Foley RA
28. Lahr MM
29. Nielsen R
30. Kristiansen K
31. Willerslev E
(2015) Early divergent strains of Yersinia pestis in Eurasia 5,000 years ago
Cell 163:571–582.

https://doi.org/10.1016/j.cell.2015.10.009
- PubMed
- Google Scholar
1. Ren GX
2. Guo XP
3. Sun YC
(2017) HmsC controls Yersinia pestis biofilm formation in response to redox environment
Frontiers in Cellular and Infection Microbiology 7:355.

https://doi.org/10.3389/fcimb.2017.00355
- PubMed
- Google Scholar
(2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data
Bioinformatics 26:139–140.

https://doi.org/10.1093/bioinformatics/btp616
- Google Scholar
1. Rodnina MV
2. Korniy N
3. Klimova M
4. Karki P
5. Peng BZ
6. Senyushkina T
7. Belardinelli R
8. Maracci C
9. Wohlgemuth I
10. Samatova E
11. Peske F
(2020) Translational recoding: canonical translation mechanisms reinterpreted
Nucleic Acids Research 48:1056–1067.

https://doi.org/10.1093/nar/gkz783
- Google Scholar
1. Rogov VV
2. Bernhard F
3. Löhr F
4. Dötsch V
(2004) Solution structure of the Escherichia coli YojN histidine-phosphotransferase domain and its interaction with cognate phosphoryl receiver domains
Journal of Molecular Biology 343:1035–1048.

https://doi.org/10.1016/j.jmb.2004.08.096
- Google Scholar
1. Rogov VV
2. Rogova NYu
3. Bernhard F
4. Löhr F
5. Dötsch V
(2011) A disulfide bridge network within the soluble periplasmic domain determines structure and function of the outer membrane protein RcsF
Journal of Biological Chemistry 286:18775–18783.

https://doi.org/10.1074/jbc.M111.230185
- Google Scholar
(2010) Role of the Yersinia pestis yersiniabactin iron acquisition system in the incidence of flea-borne plague
PLOS ONE 5:e14379.

https://doi.org/10.1371/journal.pone.0014379
- PubMed
- Google Scholar
(2021) CsrA enhances cyclic-di-GMP biosynthesis and Yersinia pestis biofilm blockage of the flea foregut by alleviating Hfq-dependent repression of the hmsT mRNA
MBio 12:e0135821.

https://doi.org/10.1128/mBio.01358-21
- PubMed
- Google Scholar
1. Smith LM
2. Jackson SA
3. Malone LM
4. Ussher JE
5. Gardner PP
6. Fineran PC
(2021) The Rcs stress response inversely controls surface and CRISPR-Cas adaptive immunity to discriminate plasmids and phages
Nature Microbiology 6:162–172.

https://doi.org/10.1038/s41564-020-00822-7
- PubMed
- Google Scholar
1. Sodeinde OA
2. Subrahmanyam YVBK
3. Stark K
4. Quan T
5. Bao Y
6. Goguen JD
(1992) A surface protease and the invasive character of plague
Science 258:1004–1007.

https://doi.org/10.1126/science.1439793
- Google Scholar
1. Song Y
2. Tong Z
3. Wang J
4. Wang L
5. Guo Z
6. Han Y
7. Zhang J
8. Pei D
9. Zhou D
10. Qin H
11. Pang X
12. Han Y
13. Zhai J
14. Li M
15. Cui B
16. Qi Z
17. Jin L
18. Dai R
19. Chen F
20. Li S
21. Ye C
22. Du Z
23. Lin W
24. Wang J
25. Yu J
26. Yang H
27. Wang J
28. Huang P
29. Yang R
(2004) Complete genome sequence of Yersinia pestis strain 91001, an isolate avirulent to humans
DNA Research 11:179–197.

https://doi.org/10.1093/dnares/11.3.179
- PubMed
- Google Scholar
1. Spyrou MA
2. Tukhbatova RI
3. Wang C-C
4. Valtueña AA
5. Lankapalli AK
6. Kondrashin VV
7. Tsybin VA
8. Khokhlov A
9. Kühnert D
10. Herbig A
11. Bos KI
12. Krause J
(2018) Analysis of 3800-year-old Yersinia pestis genomes suggests Bronze Age origin for bubonic plague
Nature Communications 9:2234.

https://doi.org/10.1038/s41467-018-04550-9
- Google Scholar
1. Srivastava A
2. Gogoi P
3. Deka B
4. Goswami S
5. Kanaujia SP
(2016) In silico analysis of 5’-UTRs highlights the prevalence of Shine-Dalgarno and leaderless-dependent mechanisms of translation initiation in bacteria and archaea, respectively
Journal of Theoretical Biology 402:54–61.

https://doi.org/10.1016/j.jtbi.2016.05.005
- PubMed
- Google Scholar
(2008) Experimental evidence for negative selection in the evolution of a Yersinia pestis pseudogene
PNAS 105:8097–8101.

https://doi.org/10.1073/pnas.0803525105
- PubMed
- Google Scholar
(2011) Differential control of Yersinia pestis biofilm formation in vitro and in the flea vector by two c-di-GMP diguanylate cyclases
PLOS ONE 6:e19267.

https://doi.org/10.1371/journal.pone.0019267
- PubMed
- Google Scholar
1. Sun YC
2. Guo XP
3. Hinnebusch BJ
4. Darby C
(2012) The Yersinia pestis Rcs phosphorelay inhibits biofilm formation by repressing transcription of the diguanylate cyclase gene hmsT
Journal of Bacteriology 194:2020–2026.

https://doi.org/10.1128/JB.06243-11
- PubMed
- Google Scholar
(2014) Retracing the evolutionary path that led to flea-borne transmission of Yersinia pestis
Cell Host & Microbe 15:578–586.

https://doi.org/10.1016/j.chom.2014.04.003
- PubMed
- Google Scholar
1. Takeda S
2. Fujisawa Y
3. Matsubara M
4. Aiba H
5. Mizuno T
(2001) A novel feature of the multistep phosphorelay in Escherichia coli: a revised model of the RcsC -- > YojN -- > RcsB signalling pathway implicated in capsular synthesis and swarming behaviour
Molecular Microbiology 40:440–450.

https://doi.org/10.1046/j.1365-2958.2001.02393.x
- PubMed
- Google Scholar
1. Tata M
2. Kumar S
3. Lach SR
4. Saha S
5. Hart EM
6. Konovalova A
(2021) High-throughput suppressor screen demonstrates that RcsF monitors outer membrane integrity and not Bam complex function
PNAS 118:e2100369118.

https://doi.org/10.1073/pnas.2100369118
- PubMed
- Google Scholar
1. Tong ZZ
2. Zhou DS
3. Song YJ
4. Zhang L
5. Pei D
6. Han YP
7. Pang X
8. Li M
9. Cui BZ
10. Wang J
11. Guo ZB
12. Qi ZZ
13. Jin LX
14. Zhai JH
15. Du ZM
16. Wang XY
17. Wang J
18. Huang PT
19. Yang HM
20. Yang RF
(2005) Genetic variations in the pgm locus among natural isolates of Yersinia pestis
The Journal of General and Applied Microbiology 51:11–19.

https://doi.org/10.2323/jgam.51.11
- PubMed
- Google Scholar
(2016) The diguanylate cyclase HsbD intersects with the HptB regulatory cascade to control Pseudomonas aeruginosa biofilm and motility
PLOS Genetics 12:e1006354.

https://doi.org/10.1371/journal.pgen.1006354
- PubMed
- Google Scholar
1. Vieille C
2. Elmerich C
(1990) Characterization of two Azospirillum brasilense Sp7 plasmid genes homologous to Rhizobium meliloti nodPQ
Molecular Plant-Microbe Interactions 3:389–400.

https://doi.org/10.1094/mpmi-3-389
- PubMed
- Google Scholar
(2018) The complex Rcs regulatory cascade
Annual Review of Microbiology 72:111–139.

https://doi.org/10.1146/annurev-micro-090817-062640
- PubMed
- Google Scholar
(2020) IgaA negatively regulates the Rcs phosphorelay via contact with the RcsD phosphotransfer protein
PLOS Genetics 16:e1008610.

https://doi.org/10.1371/journal.pgen.1008610
- PubMed
- Google Scholar
(2009) RNA-Seq: a revolutionary tool for transcriptomics
Nature Reviews. Genetics 10:57–63.

https://doi.org/10.1038/nrg2484
- PubMed
- Google Scholar
1. Wang D
2. Qi M
3. Calla B
4. Korban SS
5. Clough SJ
6. Cock PJA
7. Sundin GW
8. Toth I
9. Zhao Y
(2012) Genome-Wide identification of genes regulated by the Rcs phosphorelay system in Erwinia amylovora
Molecular Plant-Microbe Interactions 25:6–cs17.

https://doi.org/10.1094/MPMI-08-11-0207
- PubMed
- Google Scholar
1. Wren BW
(2003) The yersiniae--a model genus to study the rapid evolution of bacterial pathogens
Nature Reviews. Microbiology 1:55–64.

https://doi.org/10.1038/nrmicro730
- PubMed
- Google Scholar
1. Yan MY
2. Yan HQ
3. Ren GX
4. Zhao JP
5. Guo XP
6. Sun YC
(2017) CRISPR-Cas12a-assisted recombineering in bacteria
Applied and Environmental Microbiology 83:e00947-17.

https://doi.org/10.1128/AEM.00947-17
- PubMed
- Google Scholar
1. Yang R
2. Atkinson S
3. Chen Z
4. Cui Y
5. Du Z
6. Han Y
7. Sebbane F
8. Slavin P
9. Song Y
10. Yan Y
11. Wu Y
12. Xu L
13. Zhang C
14. Zhang Y
15. Hinnebusch BJ
16. Stenseth NC
17. Motin VL
(2023) Yersinia pestis and plague: some knowns and unknowns
Zoonoses 3:2022-0040.

https://doi.org/10.15212/ZOONOSES-2022-0040
- Google Scholar
1. Zhang Q
2. Wang Q
3. Tian G
4. Qi Z
5. Zhang X
6. Wu X
7. Qiu Y
8. Bi Y
9. Yang X
10. Xin Y
11. He J
12. Zhou J
13. Zeng L
14. Yang R
15. Wang X
(2014) Yersinia pestis biovar microtus strain 201, an avirulent strain to humans, provides protection against bubonic plague in Rhesus macaques
Human Vaccines & Immunotherapeutics 10:368–377.

https://doi.org/10.4161/hv.27060
- PubMed
- Google Scholar
(2015) Early emergence of Yersinia pestis as a severe respiratory pathogen
Nature Communications 6:7487.

https://doi.org/10.1038/ncomms8487
- PubMed
- Google Scholar

Article and author information

Author details

Xiao-Peng Guo

NHC key laboratory of Systems Biology of Pathogens, Institute of Pathogen Biology, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China

Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing

Contributed equally with
Hai-Qin Yan and Wenhui Yang

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0001-5745-2866
Hai-Qin Yan
1. Department of Basic Medical Sciences, Anhui Key Laboratory of Infection and Immunity, Bengbu Medical College, Bengbu, China
2. Paul G. Allen School for Global Health, Washington State University, Pullman, United States
Contribution
Conceptualization, Data curation, Formal analysis, Funding acquisition, Validation, Investigation, Visualization, Writing – original draft, Writing – review and editing

Contributed equally with
Xiao-Peng Guo and Wenhui Yang

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-1063-3840
Wenhui Yang

State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China

Contribution
Conceptualization, Data curation, Formal analysis, Validation, Investigation, Visualization

Contributed equally with
Xiao-Peng Guo and Hai-Qin Yan

Competing interests
No competing interests declared
Zhe Yin

State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China

Contribution
Formal analysis

Competing interests
No competing interests declared
Viveka Vadyvaloo

Paul G. Allen School for Global Health, Washington State University, Pullman, United States

Contribution
Conceptualization, Resources, Supervision, Funding acquisition, Writing – review and editing

For correspondence
vvadyvaloo@wsu.edu

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0003-4842-0525
Dongsheng Zhou

State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Microbiology and Epidemiology, Beijing, China

Contribution
Conceptualization, Resources, Supervision, Project administration, Writing – review and editing

For correspondence
zhouds@bmi.ac.cn

Competing interests
No competing interests declared
Yi-Cheng Sun

NHC key laboratory of Systems Biology of Pathogens, Institute of Pathogen Biology, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China

Contribution
Conceptualization, Resources, Supervision, Funding acquisition, Writing – original draft, Project administration, Writing – review and editing

For correspondence
sunyc@ipbcams.ac.cn

Competing interests
No competing interests declared

"This ORCID iD identifies the author of this article:" 0000-0002-5790-7071

Funding

National Key R&D Program of China (2022YFC2303200)

Xiao-Peng Guo

National Natural Science Foundation of China (31700072)

Xiao-Peng Guo

National Natural Science Foundation of China (31670139)

Yi-Cheng Sun

National Natural Science Foundation of China (31800120)

Hai-Qin Yan

Chinese Academy of Medical Sciences (Non-profit Central Research Institute Fund 2019HY310001)

Yi-Cheng Sun

Chinese Academy of Medical Sciences (Innovation Fund for Medical Sciences (CIFMS) 2021-I2M-1-043)

Yi-Cheng Sun

Fundamental Research Funds for the Central Universities (3332021092)

Yi-Cheng Sun

National Institutes of Health (Research Project Grant R01AI117016-01A1)

Viveka Vadyvaloo

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

We thank Dr. Xiaoyun Pang at Institute of Biophysics, Chinese Academy of Sciences (Beijing, China) for help to predict the structure of RcsD-Hpt by AlphaFold2.

Ethics

The animal study of flea blockage, related to Figure 5A, 5B and 5C, was performed in strict accordance with the U.S. National Institutes of Health (NIH) Guide for the Care and Use of Laboratory Animals (National Research Council Committee for the Update of the Guide for the Care and Use of Laboratory, 2011) and as approved by the Washington State University Institutional Animal Care and Use Committee, under the Animal Subject Approval Form (ASAF) 6641 and 6396. The animal study of murine infection, related to Figure 6A, 6B and 6C, was performed in strict accordance to the Guidelines for the Welfare and Ethics of Laboratory Animals of China and all the animal experiments were approved by the Institutional Animal Care and Use Committee (IACUC) of Academy of Military Medical Sciences (AMMS), ethical approval number IACUC-DWZX-2021-057.

Copyright

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.