The FUR-like regulators PerRA and PerRB integrate a complex regulatory network that promotes mammalian host-adaptation and virulence of Leptospira interrogans

Leptospira interrogans, the causative agent of most cases of human leptospirosis, must respond to myriad environmental signals during its free-living and pathogenic lifestyles. Previously, we compared L. interrogans cultivated in vitro and in vivo using a dialysis membrane chamber (DMC) peritoneal implant model. From these studies emerged the importance of genes encoding the Peroxide responsive regulators PerRA and PerRB. First described in in Bacillus subtilis, PerRs are widespread in Gram-negative and -positive bacteria, where regulate the expression of gene products involved in detoxification of reactive oxygen species and virulence. Using perRA and perRB single and double mutants, we establish that L. interrogans requires at least one functional PerR for infectivity and renal colonization in a reservoir host. Our finding that the perRA/B double mutant survives at wild-type levels in DMCs is noteworthy as it demonstrates that the loss of virulence is not due to a metabolic lesion (i.e., metal starvation) but instead reflects dysregulation of virulence-related gene products. Comparative RNA-Seq analyses of perRA, perRB and perRA/B mutants cultivated within DMCs identified 106 genes that are dysregulated in the double mutant, including ligA, ligB and lvrA/B sensory histidine kinases. Decreased expression of LigA and LigB in the perRA/B mutant was not due to loss of LvrAB signaling. The majority of genes in the perRA and perRB single and double mutant DMC regulons were differentially expressed only in vivo, highlighting the importance of host signals for regulating gene expression in L. interrogans. Importantly, the PerRA, PerRB and PerRA/B DMC regulons each contain multiple genes related to environmental sensing and/or transcriptional regulation. Collectively, our data suggest that PerRA and PerRB are part of a complex regulatory network that promotes host adaptation by L. interrogans within mammals.


Introduction
The ability of pathogenic Leptospira spp. to sense and respond to environmental signals encountered within mammals is generally believed to be critical to sustain the bacterium within its zoonotic lifestyle. The majority of studies investigating gene regulation by L. interrogans have done so by manipulating in vitro growth conditions [19][20][21][22][23][24][25][26][27]. However, numerous studies using another enzootic spirochetal pathogen, Borrelia burgdorferi, have shown that cultivation in vitro under "mammalian host-like" conditions (i.e., increased temperature, increased pH, high osmolality) does not replicate the full range of environmental signals and physiological cues that spirochetes respond to in vivo [28][29][30][31][32][33][34][35][36][37][38][39][40]. Thus, to gain better insight into the transcriptomic and antigenic changes that L. interrogans undergoes within mammals, we developed an in vivo model in which leptospires are cultivated within dialysis membrane chambers (DMCs) implanted into the peritoneal cavities of rats, a natural reservoir host [41,42]. Leptospires within DMCs (6-8 kDa MWCO) are exposed to host-derived nutrients and environmental signals but are protected from the host's cellular and humoral immune responses. Importantly, the DMC model provides sufficient numbers of host-adapted organisms (~10 8 per ml) for genome-wide transcriptomics [41] and proteomics [43]. Using this model, we identified 166 genes (110 upregulated and 56 downregulated) differentiallyexpressed by L. interrogans serovar (sv.) Copenhageni strain Fiocruz L1-130 in response to host-specific signals [41]. Almost all of the genes upregulated by the Fiocruz L1-130 strain within DMCs were unique to pathogenic leptospires (i.e., not found in the genomes of saprophytic Leptospira species).
Not surprisingly, many of the genes upregulated by L. interrogans in DMCs encode functions related to environmental signaling and gene regulation [41], including LIMLP10155 (LIC12034), which encodes a member of the Ferric Uptake Regulator (FUR) superfamily [44]. The namesake of this highly diverse superfamily, Fur, functions as a global regulator of iron homeostasis in Gram-negative and -positive bacteria, controlling both the induction of iron uptake systems under iron limitation and the expression of iron storage proteins and iron-utilizing enzymes under iron sufficiency [45]. The FUR superfamily is diverse and includes regulatory sensors for zinc (Zur), manganese (Mur) and nickel (Nur) [46][47][48][49][50][51][52]. Operating under a divergent regulatory scheme, Iron response regulators (Irrs) sense Fe-heme and repress heme biosynthetic genes under iron-limiting conditions [53]. Unlike most FURs, PerRs, are not involved in metal homeostasis per se but instead sense intracellular peroxide and regulate genes involved in detoxification of ROS in a metal-dependent manner [45,[54][55][56][57][58]. While Fur and Zur regulators are widely distributed across both Gram-positive and Gram-negative bacteria, other FUR family regulators have more limited distribution. PerRs are found mainly in Gram-positive bacteria, and Irrs are limited to α-proteobacteria [53]. So far, Mur and Nur have been characterized in α-proteobacteria and actinomycetes, respectively, but their distribution within other taxonomic groups is still unclear. Typically, FURs act as repressors; inactivation of the regulator leads to 'constitutive de-repression' of target genes in the mutant. Although the mechanisms are not well understood, examples of FURs, including PerRs, acting as activators are well documented [59][60][61][62][63][64]. Beyond balancing metal homeostasis and toxicity, FURs also modulate intermediary metabolism, host colonization and virulence [52, [65][66][67][68][69][70].
In diverse bacteria, iron serves as an essential co-factor for many cellular processes, including energy generation via electron transport, intermediary metabolism and DNA biogenesis [68][69][70][71][72]. For many pathogens, the shift from a high-to low-iron environment is a key environmental signal for induction of expression of virulence genes [69]. Unlike other spirochetes, such as B. burgdorferi and Treponema pallidum, which require iron in trace amounts, if at all [73][74][75][76][77], Leptospira spp. require this metal for growth in vitro and, presumably, in the host [78]. Consequently, leptospires have evolved elaborate mechanisms for iron sensing, scavenging and utilization [79][80][81]. At the same time, leptospires must balance their physiological Phylogenetic analysis of L. interrogans FUR-like regulators LIMLP10155 (PDB:5NL9, PerRA), LIMLP05620 (PerRB), LIMLP04825 (Ffr1) and LIMLP18590 (Ffr2) with well-characterized FUR superfamily members from diverse bacteria. Phylogenetic analyses were performed as described in Materials and Methods. A midpoint rooted tree was generated using iTOL [200]. Fur family regulators represented in the tree: Bacillus subtilis PerR Residues confirmed to be involved in regulatory metal coordination (•) are highlighted in yellow, green or gray; position 103 is used to discriminate between PerR/Irrs (Asp, green) and Fur/ Zur/Mur/Nur regulators (Glu, gray). CxxC-motif residues (�) confirmed to be involved in structural metal coordination are highlighted in cyan. Residues in red are predicted but not confirmed by X-ray crystallography to be involved in regulatory or structural metal coordination. Asparagine (N) or arginine (R) residues (E) in blue, located in DNA binding helix H4, can be used to distinguish between PerR and Fur, respectively [97]. � , the PDB structure for E. coli Fur includes only the DNA binding domain. Numbers on the top correspond to residues positions in L. interrogans PerRA.
PerRA and PerRB clustered most closely with PerRs and iron-response regulators (Irrs), while LIMLP18590 and LIMLP04825 clustered with Zur/Nur and Fur/Mur regulators, respectively.
We next surveyed the amino acid sequences of the leptospiral FUR-like proteins for conserved regulatory and structural metal binding sites (MBS), which promote DNA binding and folding/dimerization, respectively, in other Fur family regulators [49,50,95]. As noted recently by Zavala-Alvarado et al. [93], PerRA and PerRB contain two PerR canonical amino acid residues (Asn60 and Asn68 in PerRA and PerRB, respectively) involved in peroxide sensitivity and DNA recognition (Asp103 and Asp112 in PerRA and PerRB, respectively) [96,97]. Based on these features and increased expression of perRB upon exposure of L. interrogans to peroxide, LIMLP05620 was named perRB [93]. Interestingly, as shown in Fig 1B, the aspartate of the PerR regulatory MBSs and the asparagine in the PerR DNA-binding helices (DBH) also are conserved in Irr proteins. As noted previously by Kebouchi et al. [92] and Zavala-Alvarado et al. [93,94], both PerRA and PerRB lack the C-terminal conserved CxxC motif(s) used for structural metal-dependent dimerization by many, but not all, FUR family regulators; this cysteinate motif also is absent in Irrs. Overall, the PerRA and PerRB DBHs are not highly conserved, raising the possibility that they recognize different upstream sequences. LIMLP04825, on the other hand, contains features conserved across Fur, Mur, Zur and Nur regulators, including a glutamic acid at position 103, one or possibly two CxxC motifs (residues 95-98 and [133][134][135][136], and an arginine (Arg60) within its DBH ( Fig 1B). Interestingly, LIMLP18590 contains features of both PerR (Asp at position 103) and Fur/Mur/Zur/Nur (Arg residue within its putative DBH). The regulatory metal binding site(s) for LIMLP18590 most closely resembles that of a Zur (Fig 1B), which includes two putative tetra-coordinated zinc binding sites rather than the single penta-coordinated site used by PerR. However, without additional data regarding the peroxide responsiveness and/or regulatory metal-binding properties of LIMLP04825 or LIMLP18590, it is not possible to discern their function(s). For this reason, we propose designating them as Fur family regulators 1 (Ffr1) and 2 (Ffr2), respectively.
We next assessed the conservation of FUR family regulators across pathogenic (P1 and P2) and saprophytic (S1 and S2) leptospiral subclades [98]. Orthologs for PerRA and Ffr1 were identified in all highly pathogenic (P1), some intermediate (P2) and all saprophytic strains (S1 and S2), whereas orthologs for PerRB and Ffr2 were found exclusively in pathogenic strains (Figs 2A, 2B and S1). Our analyses also identified two additional FUR family regulators, both of which were found only in saprophytic leptospires (Figs 2A, 2B and S1). The first, designated PerRC, contains features of a canonical PerR (two CxxC motifs, an aspartic acid residue within its regulatory MBS and an asparagine within its putative DBH). The second saprophyte-specific FUR family regulator, designated Ffr3, resembles a Fur/Mur/Nur-like regulator (two CxxC motifs, a glutamic acid residue within its regulatory MBS and an arginine within its putative DBH).

L. interrogans FUR family regulators are expressed at higher or comparable levels in DMCs compared to in vitro
Previously, we reported that expression of perRA in L. interrogans sv. Copenhageni strain Fiocruz L1-130 was induced 3.83-fold in response to mammalian host signals compared to in vitro [41]. Using qRT-PCR, we compared transcript levels for all four FUR-like regulators in L. interrogans sv. Manilae strain L495 [99,100] grown in vitro (EMJH at 30˚C) and following cultivation within DMCs. As shown in Fig 3, perRA (7.59-fold), ffr1 (3.20-fold), and ffr2 (5.70-fold) were upregulated significantly (p<0.05) in vivo. perRB was upregulated 1.63-fold in DMCs compared to in vitro, but the difference was not statistically significant (Fig 3).

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans

Inactivation of both perRA and perRB in L. interrogans results in loss of virulence in mice
Previously, Murray et al. [99] and Zavala-Alvarado et al. [94] independently reported that a L. interrogans Manilae perRA Tn mutant is virulent in hamsters. More recently, Zavala-Alvarado et al. [93] established that L. interrogans lacking PerRB also retain virulence in hamsters. Zavala-Alvarado and colleagues also generated a double mutant by insertional inactivation of perRA in the perRB Tn mutant; the resulting double mutant (perRA/B) was avirulent in hamsters [93]. Golden Syrian hamsters are exquisitely sensitive to L. interrogans and develop acute, fulminant, disseminated disease at doses as low as 10 1 [17,99,101,102]. Mice, on the other hand, are a natural reservoir for L. interrogans and relatively resistant to infection; at sublethal doses, susceptible mouse strains develop a self-resolving hematogenous dissemination phase (~1 week) followed by chronic, asymptomatic renal colonization marked by shedding large numbers of leptospires in urine [17,[103][104][105][106].
Given the differences in leptospiral disease progression and severity between hamsters and mice, we asked whether PerRA, PerRB, or both are required to establish infection and persistence within a reservoir host model. At the outset, we first established that our wild-type serovar Manilae parent (WT) is virulent in C3H/HeJ mice. Female 10-week old mice (n = 5 per group) were infected intraperitoneally with 5 × 10 6 , 1 × 10 6 , 1 × 10 5 and 1 × 10 4 leptospires. Mice were monitored daily for signs of disease (i.e., weight loss). Within 6 days, all mice in the 5 x 10 6 group and 3 of 5 mice in the 10 6 group succumbed to infection, while all others survived the entire 42-day experimental time course ( Fig 4A). Based on these virulence studies, harvested from all surviving mice infected with the WT parent, including the single urine-negative mouse from the 10 4 group, were culture-positive.
Prior to using the perRA and perRB single and double mutants for murine virulence studies, we first confirmed their genotypes by amplicon sequencing using primers listed in S6 Table  and immunoblot and established that loss of one regulator had no obvious effect on expression of the other in the corresponding single mutants (S2 Fig). We next compared infectivity of the WT, perRA, perRB or perRA/B strains in C3H/HeJ mice using a sublethal intraperitoneal inoculum (1 x 10 5 ). All of the mice inoculated with the WT parent and single mutants were infected, shedding comparable numbers of leptospires in their urine at 14-and 21-days p.i. (Fig 4C). In contrast, no leptospires were detected in urine from mice inoculated with the perRA/B double mutant. Consistent with data from urine, at day 28 p.i., all perRA/B-infected mice were negative for leptospires by both culture and qPCR ( Fig 4D). Lastly, in contrast to mice infected with the WT or single mutant strains, all of which generated robust serological responses against L. interrogans, none of the mice infected with the double mutant seroconverted ( Fig 4E).
perRA and perRB single and double mutants grow normally in rat peritoneal cavities PerR regulators have been linked to a wide range of physiological functions outside of oxidative stress, including metal homeostasis, metabolism and virulence [55, 107,108]. In Bacillus subtilis, inactivation of perR leads to increased expression of fur and iron starvation [82]. To examine whether the avirulent phenotype of the perRA/B double mutant could be due to an inability to grow in mammals, we took advantage of our DMC model, whereby leptospires are cultivated for 9-10 days within dialysis membrane chambers implanted in the peritoneal cavity of a rat [41,109]. Originally developed for B. burgdorferi, this model is able to separate genes related to physiological adaption (i.e., nutrient acquisition and metabolism) from those encoding virulence determinants, such as adhesins, motility and immune evasion. However, we saw no significant difference (p>0.05) in the mean number of leptospires for the wild-type (1.75 × 10 8 /ml) strain versus each mutant (perRA, 8.5 × 10 7 /ml; perRB, 3.57 × 10 8 /ml; and perRA/B, 2.70 × 10 8 /ml) recovered from DMCs 9 days post-implantation (3 biological replicates per strain). These data demonstrate that the virulence-defect observed with the double mutant is not due to a metabolic lesion (i.e., metal starvation).

Defining the PerRA and PerRB regulons in vivo by comparative RNA-Seq
Prototypical FUR family regulators, including PerR, modulate transcription by binding to DNA via one or more~19-bp inverted repeats ('boxes') located upstream of their target genes [110]. Kebouchi et al. [92] previously identified three potential PerR binding sites upstream of perRA in L. interrogans. However, searches of the Manilae genome using these sequences, as well as canonical Fur and PerR boxes [110, 111], did not identify additional hits [41,44,91]. Therefore, to identify genes controlled by PerRA, PerRB, or both, in response to host signals, we performed comparative RNA-Seq using WT, perRA, perRB and perRA/B strains cultivated in DMCs (3 biological replicates per strain); a summary of the raw Illumina read data is presented in S1 Table. Reads were mapped using EDGE-pro [112] and analyzed for differentiallyexpressed genes using DESeq2 [113]. Genes expressed at �3-fold higher/lower levels in the WT versus mutant with a False-discovery rate (FDR)-adjusted-p value (q) �0.05 were considered differentially expressed. Complete RNA-Seq datasets for all comparisons are presented in S2-S4 Tables. Raw read files have been deposited in the NCBI Sequence Read Archive (SRA) database (BioProject accession PRJNA659512).
Overview of the PerRA DMC regulon. The PerRA DMC regulon contained a total of 81 differentially expressed genes; 43 were expressed at higher levels (i.e., upregulated directly or indirectly by PerRA) in the WT parent compared to the perRA mutant, while 38 were expressed at lower levels (i.e., downregulated/repressed directly or indirectly by PerRA) (S2 Table). Notably, the PerRA DMC regulon is substantially larger than its in vitro counterpart (17 genes total but only 14 dysregulated >3-fold), recently reported by Zavala-Alvarado et al. [94]. Overlap between the PerRA DMC and in vitro regulons consists primarily of seven genes located in a single chromosomal locus (S3 Fig) containing LipL48 (LIMLP04280), a putative outer-membrane embedded TonB-dependent receptor (TBDR, LIMLP04270) and one of the 2-3 putative TonB/ExbD/ExbB transporters systems (LIMLP04245-04230) encoded by L. interrogans [44,91]. TonB-dependent transporters (TBDT) for iron typically are repressed by Fur [114]; thus, it was surprising that this system was upregulated by PerRA both in vitro and in DMCs. Interestingly, none of the prototypical oxidative stress-related genes identified by Zavala-Alvarado et al. [94] as being under PerRA control in vitro were dysregulated in DMCs.
More than half (55%) of genes in the PerRA DMC regulon encode proteins of unknown function (S5A Fig). The remaining genes are distributed over a wide range of functional categories (COGs) related to cellular homeostasis and metabolism. Most notably, the PerRA regulon includes five genes (two upregulated, three downregulated) involved in signaling and/or gene regulation (Figs 5A and S4A and S2 Table). The three upregulated genes (LIMLP02515, LIMLP05780 and LIMLP01845) encode putative DNA binding proteins, including a CsoR-like metal sensitive repressor, while the three downregulated genes encode a two-component system (TCS) histidine kinase with four Per-Arnt-Sim (PAS)-type sensor domains (LIMLP10140), a putative DNA binding protein (LIMLP00900) and a putative serine/threonine kinase with GAF domain (LIMLP11575). PAS domains are ubiquitous in bacteria and sense a wide range of ligands, including heme, FAD, fatty acids and divalent metals [115,116]. GAF domains share a similar fold to PAS domains and often regulate the catalytic activity of cyclic nucleotide phosphodiesterases [117]. Of note, none of these putative regulatory factors were dysregulated �3 fold by loss of PerRA in vitro [94].
Overview of the PerRB DMC regulon. Inactivation of perRB resulted in dysregulation of 200 genes (131 upregulated and 69 downregulated) within DMCs (S3 Table). In contrast, only 30 genes were dysregulated in the perRB mutant in vitro, with only one affected >3-fold [93]. Remarkably, we saw no overlap between the in vitro and DMC PerRB regulons. Overlap between the PerRA and PerRB DMC regulons was limited to genes within the TonB-dependent transporter locus described above (S3 Fig). Notably, none of the TonB-related genes were dysregulated �3 fold in the perRB mutant in vitro (S3 Fig). The implications of these data are two-fold; differences between the in vitro and DMC regulons for the perRB mutant imply that PerRB is not fully activated under normal growth conditions in vitro, while the minimal overlap between the PerRA and PerRB DMC regulons suggests that they recognize different upstream binding sites.
The majority (66%) of genes in the PerRB DMC regulon were upregulated; most of these encode proteins with unknown or poorly characterized functions (S4B Fig and S3 Table). Notably, however, the PerRB DMC regulon includes 17 genes (11 upregulated, 6 downregulated) related to signaling and/or gene regulation ( Fig 5B). The 11 upregulated genes include six related to signal transduction, three putative DNA binding proteins (LIMLP16420, LIMLP07225 and LIMLP15900), an ECF-type sigma factor (LIMLP14515) and a putative serine/threonine phosphatase with a PAS-type sensor domain (LIMLP06960) (Fig 5B). The six downregulated signaling genes include two additional putative DNA binding proteins (LIMLP07895 and LIML10055), a second ECF-type sigma factor (LIMLP14940), a putative anti-sigma factor antagonist (LIMLP04815), and an EAL-type phosphodiesterase (LIMLP04775) (Fig 5B). None of these putative regulators were affected in vitro by loss of PerRB [93].

Inactivation of both PerRA and PerRB results in a DMC regulon that differs dramatically from its single mutant counterparts
The PerRA/B DMC regulon contains 106 differentially expressed genes, 74 upregulated and 32 repressed (Tables 1, 2 and S4). Surprisingly, we saw limited overlap between the DMC regulons for the double and single mutants ( Fig 6A); all of the overlapping genes were located in

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans 1 Gene identifications and descriptions are based on L. interrogans sv. Manilae strain UP-MMC-NIID LP (accession numbers NZ_CP011931.1, NZ_CP011932.1 and NZ_CP011933.1). In some cases, annotations have been manually curated to conform with prior studies or bioinformatics.

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans 3 Localization is based on cumulative data from BUSCA [193], LipoP [192] and SpLip [191]. Abbreviations: CY, cytoplasmic; EC, extracellular; IM, integral inner membrane; OM, outer membrane-embedded; and LIPO, lipoprotein, which may be localized to the surface (OM insertion) or periplasm (OM or IM insertion).  Table 1). LIMLP14200 and LIMLP12135 contain domains found in lipases/esterases (IPR0002489) [118] and alkaline phosphatases and sulfatases (IPR000917), respectively. LIMLP18070 contains an ATP-NAD kinase domain (IPR022504), suggesting a role in maintaining NADP homeostasis and, by extension, NADPH-dependent reductive biosynthetic pathways. LIMLP02545, encoding one of the few gene products in the PerRA/B regulon related to oxidative stress, is a

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans putative methionine sulfoxide reductase, which catalyzes the reversible thioredoxin-dependent oxidation-reduction (repair) of Met-SO to Met [119,120]. Lastly, LIMLP15435 contains a domain found in glyoxalase/bleomycin resistance proteins; in bacteria, glyoxalases are used to detoxify methylglyoxal, a reduced derivative of pyruvate, as part of the glutathione-dependent glyoxalase system [121].
Eleven genes downregulated by PerRA/B in DMCs are involved in cellular homeostasis and metabolism (Fig 6B, Tables 2 and S4). Only three (LIMLP02795, LIMLP14175 and LIMLP08980), encoding a cytochrome c peroxidase, a thiol oxidoreductase and a Grx4 family monothiol glutaredoxin, respectively, are involved in oxidative stress adaptation.
Sensing and responding to the mammalian host environment. The PerRA/B regulon includes at least 17 genes related to environmental sensing, signaling and, potentially, host adaptation. Two (LIMLP02835 and LIMLP02840) encode methyl-accepting chemotaxis proteins co-regulated with LIMLP02845, encoding a small (62 aa) hypothetical protein of unknown function (Table 1). Nine, including lvrAB (discussed below), encode sensory histidine kinases, most of which contain PAS-type sensor domains (Fig 5C). One of the nine (LIMLP05830) encodes a regulator that contains both PAS and GAF domains ( Fig 5C). The PerRA/B DMC regulon includes six putative DNA binding proteins, four upregulated and two downregulated ( Fig 5C). Three upregulated genes belong to the TetR (LIMLP02080), Cro/C1λ (LIMLP04315) and CRP-like (LIMLP07225) repressor families, while the fourth (LIMLP00755) encodes a hypothetical protein containing a homeobox winged helix-like domain of unknown function (DUF433). A second TetR-like repressor (LIMLP19135) and a BolA-like regulator (LIMLP08975) were repressed by PerRA/PerRB in DMCs (Fig 5C). In E. coli, BolA has been linked to a range of adaptive responses, including biofilm formation and entry into stationary phase [122]. All but two of the regulatory proteins in the PerRA/B DMC regulon were dysregulated only in the double mutant; LIMLP06340, encoding a histidine kinase, and LIMLP07225, encoding a CRP-like DNA binding protein, also were upregulated by PerRB alone in DMCs (Fig 5B).
Although the vast majority (70%) of genes upregulated by PerRA/B encode proteins of unknown function (Fig 6B and Table 1), seven contain conserved domains potentially related to mammalian host adaptation and/or virulence. LIMLP08585 contains a PPM-type phosphatase domain (IPR001932); PPM domains are found in diverse regulatory proteins, including SpoIIE in B. subtilis [123]. LIMLP15425 contains a putative Lambda_Bor-like domain (PF06291), which in E. coli has been associated with increased serum survival [124,125]. LIMLP14585, annotated as a host attachment protein, contains a domain of unknown function (IPR019291) found in virulence-associated proteins from the plant pathogens Agrobacterium tumefaciens and Xanthomonas spp. [126,127]. LIMLP02040 contains a SRPBCC-like domain (cd07812), which forms a deep, hydrophobic ligand binding pocket capable of binding diverse ligands [128,129]. Three hypothetical proteins (LIMLP04635, LIMLP10965 and LIMLP16555) upregulated by PerRA/B are predicted to form β-propeller structures, which are associated with a wide range of functions, including ligand-binding, enzymatic activity, cell signaling, and protein-protein interactions [130]. Interestingly, Thibeaux et al. [131] previously noted that proteins with β-propeller repeats are enriched in highly virulent Leptospira spp. Six upregulated genes encode uncharacterized lipoproteins of unknown function (Table 1).
Eighteen (56%) genes downregulated by PerRA/B in DMCs encode proteins of unknown function ( Table 2). LIMP04970 and LIMLP11660, both predicted to encode lipoproteins, contain domains (pectin lyase-fold/IPR011050 and Ricin B lectin/IPR000772, respectively) potentially involved in binding to and/or cleavage of host-derived carbohydrates. LIMLP04765 contains an alpha/beta hydrolase domain shared by a wide range of hydrolytic enzymes. Lastly, LIMLP01455, encoding an inner membrane protein, contains a DoxX-like domain; in Mycobacterium tuberculosis, DoxX complexes with a thiosulfate sulfurtransferase (SseA) to promote resistance to agents that disrupt thiol homeostasis [132].
Known or putative virulence determinants. The upregulated portion of the PerRA/B regulon contains at least four virulence-associated genes (Table 1). Two, LIMLP15405/ligA and LIMLP15415/ligB, encode the pathogen-specific, multifunctional, Leptospiral Immunoglobulin-like repeat proteins LigA and LigB, respectively [133,134], while LIMLP08490 and LIMLP08485 encode the hybrid histidine kinases LvrA and LvrB, respectively [19]. Although tandemly located on the chromosome, ligA and ligB are not co-transcribed ( Fig 7A). They do, however, have identical upstream regions and respond similarly in vitro to conditions used to mimic the mammalian host milieu (e.g., high osmolality and increased temperature) [21][22][23][24]135]. Three genes located downstream of ligB, all encoding hypothetical proteins, also were upregulated (Fig 7A). Using antisera against the shared N-terminal repeats (Fig 7A), we compared expression of LigA and LigB in WT and mutant strains. As shown in Fig 7B, both LigA and LigB were completely absent in whole cell lysates prepared from the perRA/B double mutant cultivated within DMCs (Figs 7B, S5A and S5B); expression of both Ligs was restored to near wild-type levels by trans-complementation with perRB alone (Figs 7B, S5A and S5B). Interestingly, we saw a modest to substantial reduction in LigA/LigB in the perRA and perRB single mutants. (Figs 7B, S5A and S5B). Given that the upstream regions for ligA and ligB are identical, the molecular basis (e.g., transcriptional, post-transcriptional or both) for the difference between Lig levels in the perA and perRB mutants is unclear.
The downregulated portion of the PerRA/B DMC regulon contains at least one gene potentially related to virulence. LIMLP03665/colA, encoding a collagenase precursor [136], was expressed at~500-fold lower levels in the WT parent compared to the perRA/B mutant

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans (Table 2). While collagenase-mediated degradation of host tissues likely enhances dissemination of leptospires during early infection [137], once in the kidneys, repression of colA could help reduce pathogen-mediated damage to renal epithelial cells. Further transcriptional analysis of this gene is needed to establish its expression profile in different tissues over the course of infection.

Loss of LvrAB alone is not responsible for avirulence of the perRA/B double mutant in mice
As noted above, expression of lvrAB is disrupted only in the perRA/B double mutant (Table 1); similar results were obtained using leptospires grown in vitro [93]. Using LvrA-and LvrB-specific antisera, we confirmed our transcriptomic data at the protein level by immunoblot using whole cell lysates from WT, perRA, perRB and perRA/B mutant strains cultivated in DMCs (Figs 8A, S5A, S5C and S5D). Previously, Adhikarla et al. [19] reported that inactivation of lvrAB by transposon mutagenesis results in dysregulation of a large number of genes in vitro, including ligB. However, in our hands, we saw no decrease in LigA or LigB in the lvrAB mutant strain following cultivation in DMCs (Figs 7B, S5A and S5B). Adhikarla et al. [19] also reported that loss of either lvrAB or lvrB alone resulted in a significant loss of virulence in hamsters. To explore whether the avirulence of the perRA/B double mutant in mice (Fig 4C-4E) is due solely to loss of LvrAB, we assessed the ability of lvrAB and lvrB transposon mutants to colonize the kidneys of C3H/HeJ mice (5 mice per strain, per experiment). While mice infected with either the lvrAB or lvrB mutant shed~2-log 10 less leptospires in their urine compared to the WT controls, all of the urine samples collected from mice infected with either mutant were

PLOS PATHOGENS
PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans darkfield positive by day 21 (Fig 8B). At day 28 p.i., kidneys harvested from all mice infected with either the lvrAB or lvrB mutant were positive for leptospires by both culturing in EMJH and qPCR (S6A Fig). Mice infected with the lvrAB mutant also seroconverted (S6B Fig). Thus, while LvrAB signal transduction contributes to virulence, loss of lvrAB expression alone is not responsible for the complete loss of virulence observed with the perRA/B double mutant in mice.

Discussion
L. interrogans must sense and respond to diverse signals and threats during the free-living and reservoir host phases of its zoonotic cycle. Not surprisingly, L. interrogans encodes substantially more sensory and regulatory proteins than B. burgdorferi and T. pallidum [138], two pathogenic spirochetes with far more restrictive growth niches. However, the regulatory networks and gene products that sustain L. interrogans in nature remain poorly understood. To gain insight into the transcriptomic changes that leptospires undergo within the host, we previously compared L. interrogans sv. Copenhageni strain Fiocruz L1-130 cultivated in vitro and in mammals using our DMC peritoneal implant model [41,109]. From these studies emerged >100 genes that were differentially expressed in response to host-specific signals, including LIC12034, encoding the peroxide stress response regulator PerRA, which was upregulated 3.83-fold in DMCs. Herein, we confirmed these data using L. interrogans sv. Manilae strain L495 and also established that the three remaining FUR family regulators are transcribed at comparable (perRB) or higher (ffr1 and ffr2) levels in DMCs compared to in vitro. The importance of FUR family regulators for host adaptation was confirmed recently by Zavala-Alvarado et al. [93], who demonstrated that leptospires lacking both PerRA and PerRB are unable to infect hamsters. In our current study, we establish that these regulators also are required for renal colonization of C3H/HeJ mice. In both animal models, loss of virulence was observed only when both PerRA and PerRB were inactivated, suggesting that these regulators may serve redundant or overlapping functions in vivo. Our finding that the perRA/B double mutant survives at wild-type levels in DMCs is particularly noteworthy as it demonstrates that the avirulent phenotype observed for this mutant is not due to a metabolic lesion (i.e., metal starvation) but instead reflects dysregulation of one or more virulence-related genes. Transcriptomic analyses of perRA and perRB single and double mutants cultivated in DMCs brought to light a number of novel aspects of FUR-mediated regulation in L. interrogans. Most notably, the majority of genes in the PerRA, PerRB and PerRA/B regulons were differentially expressed only in DMCs, highlighting the importance of mammalian host-specific signals for PerRmediated regulation in L. interrogans. Remarkably, inactivation of both PerRA and PerRB resulted in a DMC regulon that differs substantially from those of either single mutant and includes a large cohort of genes involved in environmental sensing, signal transduction and transcriptional regulation.
Despite several attempts, Zavala-Alvarado et al.
[93] was unable to restore virulence to the perRA/B double mutant by trans-complementation with perRA or perRB alone. Consequently. we cannot rule out the possibility that the loss of virulence observed with the perRA/B double mutant is due to a spontaneous genetic defect outside of perRA or perRB. However, by comparative genomic sequencing of WT, perRA/B and perRB (the parental background for the perRA/B double mutant) strains, Zavala-Alvarado et al. [93] identified only two differences in the perRA/B strain. The first was a single nucleotide insertion in LIMLP11570, encoding a putative 3-oxoacyl ACP synthase related fatty acid synthesis. It is important to note that this same insertion also is observed in several L. interrogans isolates from human and animals and, as noted earlier, we saw no difference in the growth of the single and double mutants either in vitro or in DMCs. The second nonsynonymous difference is in LIMLP01895, encoding a putative hybrid histidine kinase; the corresponding polymorphism results in an alanine to valine substitution at amino acid 146; the position of this mutation is within an inter-domain region and, therefore, not likely to affect the protein's putative signal transduction function(s). Nonetheless, further investigation is necessary to establish the extent to which LIMLP01895 contributes to gene regulation and/or loss of virulence in the perRA/B double mutant.
The presence of multiple FUR family regulators in Leptospira spp. was noted previously by Louvel et al. [44], who identified five distinct orthologs between L. interrogans and L. biflexa. Phylogenetic analyses presented herein identified a sixth FUR family regulator and established that two (PerRA and Ffr1) are conserved within both pathogenic and saprophytic species and two each are unique to either pathogenic (PerRB and Ffr2) or saprophytic subclades (PerRC and Ffr3). As their designations suggest, based on sequence alignments, three are predicted to function as PerRs. Thus far, only B. licheniformis has been shown to encode multiple PerRs (PerR BL , PerR2 and PerR3), each of which displays a different level of sensitivity to H 2 O 2 (PerR2 > PerR BL > PerR3) [139]; the extent of regulatory overlap between these three PerRs has yet to be determined. The peroxide responsiveness and/or metal sensing properties of the remaining three leptospiral FUR family regulators cannot be predicted based on sequence alone. Our finding that almost all saprophytic and pathogenic Leptospira spp. encode closelyrelated PerRA and Ffr1 orthologs, however, implies that these two regulators could function outside of a host (e.g., within soil and/or water). The presence of a single PerR (PerRB) in all but one of the four P2 subclade species examined (L. wolfii) may contribute to the 'intermediate' virulence of these Leptospira spp. compared to highly virulent P1 subclade [6,140].
L. interrogans cultivated in DMCs express increased levels of catalase, AhpC-type peroxiredoxin and cytochrome c peroxidase [41], three enzymes typically associated with detoxification of reactive oxygen species (ROS) in bacteria [85,87,141]. These data also provide strong evidence that L. interrogans is exposed to ROS in vivo [41]. Consistent with this notion, catalase-deficient leptospires are more susceptible to H 2 O 2 in vitro and show reduced virulence in hamsters [142]. Host phagocytic cells, which generate oxygen radicals via a dedicated NADPH oxidase [143,144], are one likely source of exogenously-derived ROS in vivo. Leptospires within renal tubules, a highly oxygenated niche, also would be exposed to elevated levels of oxygen. Incomplete reduction of oxygen by iron-containing cytochromes is another potential source of endogenous ROS [85,141]. In bacteria, oxidative stress responses often are coordinated by two evolutionarily distinct master regulators-OxyR and PerR. OxyR, the more common of the two, belongs to the LysR family and functions primarily as an activator [145]. In its oxidized state, OxyR activates transcription of genes involved in the detoxification of H 2 O 2 (catalase and AhpC), the prevention or repair of DNA damage (Dps) and/or redox homeostasis (glutathione reductase, thioredoxin) [85]. PerR, first described in B. subtilis [146], typically represses rather than activates many of the same genes as OxyR and is released from DNA by peroxidation [49, 54, 90]. Although OxyR and PerR regulate transcription by different mechanisms, they react with H 2 O 2 at essentially the same rate constant (10 5 M -1 s -1 ) [89] and orchestrate highly similar responses. L. interrogans does not encode an OxyR homolog but, as noted above, encodes at least two PerR orthologs, PerRA and PerRB. Consistent with PerR functions in other bacteria, as shown here and elsewhere [91, 93, 94], L. interrogans perRA mutants show enhanced survival following exposure to lethal levels of H 2 O 2 in vitro and increased expression levels of catalase, AhpC and cytochrome c peroxidase. While inactivation of perRB had no effect on the ability of leptospires to withstand killing by H 2 O 2 , the perRB mutant showed increased tolerance to the superoxide-generating compound paraquat [93]. Moreover, no genes associated with ROS defenses were dysregulated in the perRB mutant in vitro [93]. In DMCs, only cytochrome c peroxidase, AhpC and a glutaredoxin were dysregulated in the double mutant. Interestingly, all three genes were expressed at higher levels in the WT compared to the mutant, suggesting that they are activated rather than repressed by PerRB. Moreover, expression of catalase was not significantly different in the WT vs. perRB or perRA/B DMC comparison.
The above data argue that while PerRA and PerRB may be 'activated' by ROS, the adaptive responses they control likely extend beyond oxidative stress. The prototypical PerR in B. subtilis (PerR Bs ) can coordinate either Mn 2+ or Fe 2+ .When co-factored with Fe 2+ , DNA binding by PerR:Fe is highly sensitive to H 2 O 2 due to irreversible iron-dependent oxidation of metalcoordinating histidine residues [54,96]. When cofactored with Mn 2+ , however, PerR Bs is able to bind DNA but is no longer peroxide sensitive [96]. In this way, PerR functions both as a peroxide responsive regulator and a ratiometric sensor for iron and manganese, altering its transcriptomic output based on intracellular metal availability and/or oxidative stress. Our finding that the PerRA/B DMC regulon contained only three genes related to oxidative stress and no genes related to iron homeostasis raises the possibility that PerRA and PerRB function in vivo may be regulated by metal availability rather than oxidative stress. Moreover, it is possible that both PerR:Fe and PerR:Mn regulate different cohorts of virulence genes, depending on the host milieu. Consistent with differences in the putative DNA-binding helices, we saw very little overlap between the PerRA and PerRB DMC regulons. Seven of the eight genes common to both regulons are located in a single locus encoding a TonB-dependent transporter (TBDT) system. In Gram-negative bacteria, TBDT systems promote the uptake of substrates, such as iron siderophores, heme, vitamin B12, and carbohydrates, that are either poorly transported by non-specific outer membrane porins or are present in the extracellular milieu at low concentration [114,147]. Substrate binding and uptake is mediated by high affinity, substrate-specific TonBdependent receptor (TBDR) proteins, which form outer membrane-embedded 22-stranded βbarrels [114]. The energy required for substrate transport is provided in the form of proton motive force, which is transduced from the inner to outer membrane by the TonB-ExbB-ExbD complex [114]. L. interrogans encodes 11 putative TonB-dependent receptors and at least two complete TonB-ExbB-ExbD transporters. None of the leptospiral TBDRs possess a N-terminal extension capable of interacting with anti-sigma factors, similar to that of the iron and heme TBDRs FecA and HasR in E. coli and Serratia marcescens, respectively [148,149]. Only one TBDT system, LIMLP04240-04270, is differentially regulated by PerRA and PerRB in DMCs. While the substrate(s) recognized by the TBDR (LIMLP04270) cannot be predicted based on sequence, mutagenesis studies on its ortholog in L. biflexa suggest that it is not essential for uptake of iron or heme in vitro [44]. Moreover, a L. interrogans transposon mutant containing an insertion in LIMLP04270 is virulent in hamsters [93]. However, given the large number of TBDRs in L interrogans, one of these may compensate for loss of LIMLP04270 in vitro and/or in vivo. Interestingly, in vitro, expression of LIMLP04240-04270 was dysregulated in the perRA and perRA/B mutants but not the perRB single mutant; only LIMLP04255 was upregulated 1.93-fold in the WT compared to the perRB mutant [93]. A second TBDR, encoded by LIMLP08410, was repressed by PerRA/B only within DMCs. Together, these data suggest that mammalian host signals play a key role in modulating TonB-dependent nutrient uptake in L. interrogans and, moreover, that the activity of PerRB is enhanced in vivo.
By comparative RNA-Seq, we identified four distinct PerR regulatory categories in L. interrogans (Fig 8). The first two include genes whose expression is controlled exclusively by a single PerR (PerRA only and PerRB only ). The most straightforward explanation for this category is that PerRA and PerRB recognize different upstream boxes. Given that the perRA and perRB single mutants are fully virulent in hamsters [93] and mice, genes in these two categories either are not required in mammals or encode redundant functions. The third category, PerRA and B, includes the TBDT locus, described above, which requires PerRA and PerRB for expression. Presumably the upstream regions for PerRA and B loci contain separate PerRA-and PerRB-specific boxes, both of which must be engaged for transcription. The fourth category, PerRA or B, contains genes that are regulated by both PerRs but require only one for expression. The most likely explanation for this category is that the upstream regions for these genes contain separate PerRA and PerRB boxes, only one of which needs to be engaged for expression. Although less likely, PerRA and PerRB also could recognize a single 'degenerate' PerR box. None of these scenarios, however, explains why the majority of PerRA only and PerRB only genes continue to be expressed in the perRA/B mutant in DMCs. We hypothesize that this unexpected regulatory scheme reflects a natural requirement for these gene products at points during the zoonotic cycle when both PerRA and PerRB are inactive. We envision two non-mutually exclusive explanations for this intriguing finding: (i) The PerRA/B DMC regulon contains at least two putative DNA binding proteins (DBPs) that are downregulated by PerRA or B in wild-type leptospires. Continued expression of these DBPs in the perRA/B double mutant could help sustain expression of PerRA only or PerRB only genes. Alternatively, loss of PerRA and PerRB could lead to a physiological state that enables L. interrogans' other FUR family regulators, Ffr1 and Ffr2, to "take over" expression of PerRA only and PerRB only genes. Examples of regulatory overlap between FURs in other bacteria are well documented [82, [150][151][152][153][154]. Further studies are needed to determine which, if any, of these scenarios are operative in L. interrogans. Surprisingly, the majority of genes controlled by PerRA and/or PerRB in DMCs were upregulated (i.e., expressed at lower levels in the single or double mutants compared to WT) rather than repressed in vivo. While only one prior study has demonstrated PerR-mediated activation [155], there are multiple examples of FUR family regulators acting as transcriptional activators [51,52]. In Vibrio vulnificus, apo-Fur positively regulates its own expression by binding upstream of the fur promoter [156]. In Helicobacter pylori and Salmonella enterica sv. Typhimurium, Fur activates expression by binding upstream of target gene and helping to recruit RNA polymerase [157,158]. In α-proteobacteria, Irrs (see below) act as positive and negative transcriptional regulators of genes related to heme homeostasis [59][60][61][62]. BosR, a FUR family regulator in the Lyme disease spirochete B. burgdorferi, activates transcription of the alternative sigma factor rpoS as part of a complex that includes the alternative sigma factor RpoN and the response regulator Rrp2 [159][160][161][162][163][164]. PerRA and/or PerRB also could activate transcription of target genes indirectly via repression of a regulatory small RNA (e.g., RyhB in E. coli) [165] or by preventing the binding of another repressor (i.e., anti-repression) [63,64].
Designation of PerRA and PerRB as peroxide stress regulators in L. interrogans is based largely on in vitro studies showing increased survival of perRA and perRB mutants following exposure to H 2 O 2 and paraquat, respectively [91,93,94]. Several lines of evidence, however, raise the possibility that these gene products function as iron response regulators (Irrs) rather than PerRs. Based on amino acid sequence alignments, PerRA and PerRB appear to be more closely related to Irrs than PerRs. In α-proteobacteria, Irrs and their regulatory partner, RirA, coordinate the expression of genes involved in heme biosynthesis with iron availability. Similar to FURs, RirA functions as metal-dependent transcriptional repressor but senses iron within Fe-S clusters rather than Fe 2+ . Interestingly, LIMLP06290 (LIC11283), annotated as a hypothetical protein, contains domains consistent with it being a RirA; the contribution of this putative RirA to iron homeostasis in Leptospira spp. has not been examined. At the sequence level, Irrs share a number of features with PerRs, including the presence of Asp and Arg residues in their regulatory metal sites and DNA binding helices, respectively. Irrs and PerRs also are responsive to similar levels of ROS, albeit by a different mechanism, and regulate many of the same effector genes (i.e., catalases and peroxidases) [62,166,167]. In some, but not all cases, irr mutants also show increased survival in vitro under high H 2 O 2 levels [62,168]. As noted above, only B. lichenformis is known to encode multiple PerRs. Numerous bacteria, on the other hand, encode two or more Irrs [166]. Variable affinity of Irrs for their target promoters enables them to modulate gene expression over a wider range of conditions than PerRs [169]. The autoregulatory sequences identified upstream of perRA [92] diverge significantly from canonical PerR and Fur boxes but show strong similarity to "Irr-boxes" [170,171]. Moreover, Irrs are known to act as activators as well as repressors [166,170]. Although typically associated with peroxide-sensitive regulation of iron/heme acquisition and utilization, Irrs have been shown to control diverse cellular processes, including virulence. Moreover, the vast majority of histidine kinases upregulated by PerRA/B in DMCs contain one or more PAS-type sensor domains, which have been shown to function as heme sensors [172]. Given the established importance of heme for survival of L. interrogans in mammals [173][174][175], our findings raise the possibility that heme sensing by PerRA and/or PerRB in mammals could serve as an important initiating event for host adaptation.
Studies presented here and elsewhere [93] demonstrate for the first time that both PerRA and PerRB are required for full transcription of the virulence-related genes ligA and ligB. These pathogen-specific surface lipoproteins have been studied extensively for their contributions to host-pathogen interactions [176,177], virulence [178] and potential use as vaccinogens [179][180][181]. Using a TALE-based transcriptional knockdown approach, Pappas and Picardeau [178] reported that both Ligs are required for virulence in hamsters. As noted earlier, ligA and ligB are not co-transcribed but instead share virtually identical upstream regions and, consequently, are co-regulated by the same environmental signals. Matsunaga, Haake and others previously reported that ligA and ligB are upregulated in response to physiological osmolarity (EMJH supplemented with 120 mM sodium chloride) [24] and increased temperature [182]. However, temperature-dependent regulation is mediated by a cis-acting RNA secondary structure that prevents translation at lower temperature and disruption of this cis element had minimal effect on osmoregulation [182]. Eshghi et al. [27] reported that inactivation of lb139 (LIMLP18410), encoding a putative anti-ECF sigma factor, resulted in~2.5-fold decreased expression of ligB in vitro. However, LIMLP18410 is not in the PerRA/B regulon either in vitro [93] or in DMCs (this study) and, therefore, is not responsible for dysregulation of ligA and ligB in the perRA/B double mutant. Additional studies are needed to establish whether PerRA and/or PerRB regulate expression of ligA and ligB directly by binding to the lig promoter region or indirectly via another effector protein. The presence of multiple sensory and regulatory effector proteins in the PerRA and PerRB DMC regulons argues that activation of PerRA and PerRB, presumably by oxidative stress, initiates a complex regulatory network capable of sensing and responding to a wide range of mammalian host-specific signals. Our finding that LvrAB-deficient leptospires express normal levels of LigA and LigB argues that at least two PerRA/B-dependent regulatory pathways (LvrAB-dependent and -independent) are operative in L. interrogans in vivo.

Ethics statement
All experiments involving animals conducted at UConn Health were performed in accordance with The Guide for the Care and Use of Laboratory Animals (8th Edition) (Guide for the Care and Use of Laboratory Animals, 1996) using protocols reviewed and approved by the UConn Health Institutional Animal Care and Use Committee [Animal Welfare Assurance (AWA) number A347-01].

Bacterial cultivation in vitro
L. interrogans strains are described in S5 Table. Leptospires were cultivated routinely in vitro in Ellinghausen, McCullough, Johnson and Harris medium (EMJH) [183,184] supplemented with 1% rabbit serum at 30˚C under static conditions. Mutants were maintained in EMJH under appropriate antibiotic selection (spectinomycin, 40 μg/ml and/or kanamycin, 40 μg/ml). Cultures were harvested at late logarithmic phase (1-5 × 10 8 per ml). Culture viability (i.e., motility and cell morphology) was evaluated by darkfield microscopy. Leptospires were enumerated using a Petroff-Hausser counting chamber (Hausser Scientific Co., Horsham, PA). Escherichia coli strains were maintained in Lysogeny broth (LB) or LB agar supplemented with the appropriate antibiotics (ampicillin, 100 μg/ml; spectinomycin, 100 μg/ml; and/or kanamycin, 100 μg/ml). The genotypes of L. interrogans mutants used in these studies were confirmed by PCR and amplicon sequencing using primers listed in S6 Table. Routine DNA manipulation and cloning  Table. PLOS PATHOGENS PerRA and PerRB control regulatory network required for virulence of Leptospira interrogans

SDS-PAGE and immunoblot analyses
To analyze L. interrogans whole cell lysates, equivalent amounts of cells (*10 8 leptospires per lane) were re-suspended and boiled in reducing Laemmli sample buffer (BioRad, Hercules, CA), separated through 10-12.5% separating polyacrylamide mini-gels and then visualized by GelCode Blue Stain Reagent (ThermoFisher Scientific, Grand Island, NY). Recombinant proteins expressed in E. coli were separated by SDS-PAGE and stained with GelCode Blue Stain Reagent (ThermoFisher). For immunoblotting, proteins were transferred to nitrocellulose membrane (GE Healthcare Life Sciences, Pittsburgh, PA) using Trans-Blot SD semi-dry transfer cell (BioRad, Hercules, CA). Membranes were blocked using milk block solution (MBS; 5% dry milk, 0.1% Tween 20, 5% fetal calf serum in PBS) for 1 h at room temperature. His-tagged recombinant proteins were detected using an HRP-conjugated anti-His monoclonal antibody (Sigma-Aldrich, St. Louis, MO) according to the manufacturer's instructions. Antisera against recombinant His-tagged leptospiral proteins were diluted 1:500 (PerRA and PerRB), 1:1000 (LvrA, LvrB and FlaB1), 1:10,000 (LigA/B repeat region) in MBS and incubated overnight at 4˚C. After washing with PBS containing 0.05% Tween 20 (PBST), bound antibody was detected with horseradish peroxidase-conjugated secondary antibody (Southern Biotechnology Associates, Birmingham, AL) diluted 1:30,000. After 1 hr at room temperature, membranes were washed at least five times with PBST and developed using the SuperSignal West Pico chemiluminescence substrate (Pierce, Rockford, IL).

Generation of host-adapted leptospires
To obtain mammalian host-adapted organisms, L. interrogans sv. Manilae strain L495 wild-type and mutant strains were cultivated in DMCs as previously described [41,42,109]. Briefly, DMCs were prepared with 9-10 mls of EMJH medium (supplemented with an additional 10% bovine serum albumin to maintain osmotic pressure) at a starting inoculum of 10 4 organisms per ml. Using strict aseptic technique, DMCs were implanted into the peritoneal cavity of an anesthetized female Sprague-Dawley rat. After nine days, animals were euthanized by CO2 narcosis and DMCs harvested. The viability and density of leptospires were evaluated by dark field microscopy using a Petroff-Hausser counting chamber (Hausser Scientific Co., Horsham, PA).

Murine infection experiments
To determine the lethal dose to 50% of mice (LD 50 ) for L. interrogans sv. Manilae strain L495, ten-week-old female C3H/HeJ mice (Jackson Laboratories, Bar Harbour, ME) were inoculated intraperitoneally (IP) with 200 μl of EMJH containing 5 × 10 6 , 10 6 , 10 5 or 10 4 leptospires (5 mice per group). Animals were monitored twice a day for signs of leptospirosis and, when moribund, were euthanized by anesthetic overdose. LD 50 was calculated using the Reed-Muench method [186]. For virulence studies, 10 5 of wild-type parent, mutant or complemented strains were used to infect C3H/HeJ mice (5 animals per group, per experiment). Beginning 14 days post-infection (p.i.), animals were monitored for the presence of leptospires in urine, collected in a metabolic chamber for~45 min following subcutaneous administration of furosemide (2-10 mg/kg, IP). Burdens in urine were assessed by darkfield microcopy using a Petroff-Hausser counting chamber. Twenty-eight or 42 days p.i. (virulence and LD 50 experiments, respectively), animals were euthanized by CO 2 narcosis and blood and kidneys were collected for serology, culturing in EMJH, and qPCR. Sera from individual mice were used to probe whole cell lysates (~10 8 leptospires per lane) prepared from the wild-type parent grown in vitro in EMJH at 30˚C.

qRT-PCR
Total RNA was isolated from leptospires (four biological replicates per condition) cultivated in vitro at 30˚C or following cultivation in DMCs as previously described [41]. cDNAs (+ and-RT) were assayed in quadruplicate in 25 μl reactions performed with SsoAdvanced Universal SYBR or Probe (lipL32) Super Mixes (Bio-Rad). Oligonucleotide primers used for qRT-PCR are provided in S6 Table. Copy numbers were calculated using internal standard curves (10 7 -10 1 copies) generated using purified amplicons for perRA, perRB, LIMLP18590 and LIMLP04825 and then normalized against lipL32 [179]. The standard curve for lipL32 was generated using a copy of the lipL32 amplicon cloned into pCR2.1-TOPO plasmid (Invitrogen). Normalized copy numbers were compared using an unpaired t test with two-tailed p values and 95% confidence interval (Prism v. 6, GraphPad Software).

Quantitation of burdens by qPCR
DNA was extracted from infected kidneys using the Qiagen DNeasy Blood & Tissue kit according to the manufacturer's recommendations. DNAs were analyzed by quantitative PCR (qPCR) using a TaqMan-based assay for lipL32 [179] in 25 μl reactions performed with SsoAdvanced Universal Probes Super Mix (Bio-Rad). Copy numbers for lipL32 were determined using an internal standard curve for the lipL32 amplicon cloned into pCR2.1 TOPO (Invitrogen). Average values for each strain were compared using an unpaired t test with two-tailed p values and 95% confidence interval (Prism v. 6, GraphPad Software).

RNA sequencing and comparative transcriptomics
Total RNA was prepared from leptospires cultivated in DMCs using TRIzol Reagent (Thermo-Fisher) (3 biological replicates per strain) and then treated twice with TURBO DNase (Ther-moFisher) followed by purification using RNeasy columns (Qiagen) as previously described [41]. Samples were eluted in RNAse-free water and purified RNA was analyzed using Qubit RNA HS Assay Kit (Thermo) and Agilent TapeStation 4200 (Agilent Technologies, Santa Clara, CA, USA) using the RNA High Sensitivity assay. Only samples with Ribosomal Integrity number (RINe) values >7.5 were used for library preparation. Stranded libraries were prepared from ribo-depleted RNA using Zymo-Seq RiboFree Total RNA Library Kit according to manufacturer's instructions. Libraries were validated for length and adapter dimer removal using the Agilent TapeStation 4200 D1000 high-sensitivity assay and then quantified and normalized using the double-stranded DNA (dsDNA) high-sensitivity assay for Qubit 3.0 (Life Technologies, Carlsbad, CA). Libraries were run on an Illumina High Output 75-cycle v2.5

Bioinformatics
Routine and comparative sequence analyses were performed using MacVector (version 17.5.4; MacVector, Inc., Apex, NC). Clusters of Orthologous Group (COG) classifications are based on MicroScope, an integrated platform for the annotation of bacterial gene function through genomic, pangenomic and metabolic comparative analysis [187]. Conserved domain searches were performed using Conserved Domain Database (CDD) Search [188], UniProt [189] and InterPro [190]. Candidate lipoproteins were identified based on Setubal et al. [191] and LipoP server [192]. Subcellular localization predictions were performed by BUSCA (Bologna Unified Subcellular Component Annotator) [193]. Multiple sequence alignments were generated by Clustal Omega [194] and MAFFT 7 [195]. Phylogenetic trees were generated using PhyML 3.0 [196] with LG substitution model chosen after an Akaike Information Criterion (AIC) model selection [197]. Tree improvement was done by subtree pruning and regrafting (SPR) method [198] with ten random starting trees. Robustness of branches was assessed by Approximate Likelihood-Ratio Test (aLRT-SH) [199]. The resulting trees were visualized and annotated using Interactive Tree of Life (iTOL, v 4.3) [200]. FUR domain-containing proteins in 26 Leptospira spp. genomes (10, 5, 6 and 5 species from subclade P1, P2, S1 and S2, respectively [98]) were identified using the Leptospira species name as a query in the Ferric-uptake regulator domain entry (IPR002481) in EMBL-EBI InterProScan [201]. Orthologs shared between L. interrogans sv. Manilae strain L495 and sv. Copenhageni strain Fiocruz L1-130 strains were identified using OrthoVenn 2.0 [202]. Burdens of leptospires in kidneys harvested from mice in Fig 8B. DNA samples from kidneys harvested 28 days post-inoculation were assessed (in quadruplicate) by qPCR using a Taqmanbased assay for lipL32. Bars represent the average and standard error of the mean. p-values were determined by comparing burdens in mice infected with wild-type (WT) and mutant strains at the same timepoint using a two-tailed t-test; we saw no significant difference (p>0.05) between burdens between the WT, lvrAB and lvrB strains. B. Immunoblot analysis of sera collected from C3H/HeJ mice 28-days following intraperitoneal inoculation with 10 5 wild-type or lvrAB mutant strains and then used to probe whole cell lysates of L. interrogans sv. Manilae strain L495 grown in EMJH at 30˚C. (TIF) S1  Table. Comparative RNA-Seq data for L. interrogans sv. Manilae L495 wild-type and perRA strains cultivated in dialysis membrane chambers (DMCs). The genome sequence of L. interrogans sv. Manilae strain UP-MMC-NIID LP (accession numbers NZ_CP011931.1, NZ_CP011932.1 and NZ_CP011933.1) was used for mapping and differential gene expression analysis. All non-coding RNAs and pseudogenes were removed before DESeq2 analysis. Column A: RefSeq locus tag. Column B: Locus tag. Column C: L. interrogans serovar Copenhageni strain Fiocruz L1-130 (accession numbers NC_005823.1 and NC_005824.1) orthologs identified using OrthoVenn 2.0 [202]. Dashes (-) indicate genes for which no clear ortholog was identified. Column D: Description of gene product, following genome annotation. Columns E, F: Clusters of Orthologous Group (COG) classifications based on MicroScope [187].  Table. Comparative RNA-Seq data for L. interrogans sv. Manilae L495 wild-type and perRA/B strains cultivated in dialysis membrane chambers (DMCs). The genome sequence of L. interrogans sv. Manilae strain UP-MMC-NIID LP (accession numbers NZ_CP011931.1, NZ_CP011932.1 and NZ_CP011933.1) was used for mapping and differential gene expression analysis. All non-coding RNAs and pseudogenes were removed before DESeq2 analysis. Column A: RefSeq locus tag. Column B: Locus tag. Column C: L. interrogans sv. Copenhageni strain Fiocruz L1-130 (accession numbers NC_005823.1 and NC_005824.1) orthologs identified using OrthoVenn 2.0. [202]. Dashes (-) indicate genes for which no clear ortholog was identified. Column D: Description of gene product, following genome annotation. Column E. Identification of conserved domain(s) within the corresponding gene product based on search of the Interpro database [190,201]. The domain identification for each gene is followed by