ANR is relatively conserved among diarrheagenic pathogens. Over the last five years, massive sequencing of new bacterial genomes has identified hundreds of new ANR members in multiple pathogens. ANR is widely distributed in at least 26 Gram-negative bacterial species25,27. Phylogenetic analysis of the amino acid sequence of ANR members from clinically relevant diarrheagenic bacteria, such as pathogenic E. coli, Salmonella enterica, and Vibrio sp., revealed divergence in ANR cognates that fall in at least three clades (termed 1 to 3) (Fig. 1A). The archetype ANR, Aar (from EAEC) and Cnr (from ETEC) fall in Clade-1. ANR homologs from AE pathogens such as EPEC, C. rodentium (Rnr), and Salmonella enterica (ANRSe) fall in Clade-2. ANRs from Vibrio (ANRVibrio) are grouped in Clade-3 (Fig. 1A).
All predicted ANR members have a low molecular mass (4.36–9.54 kDa), and they exhibit 25–100% identity with Aar (Fig. 1B). In general terms, ANRs from Clade-1 are smaller (~ 7.233 kDa) than ANRs from Clade-2 (8.903 kDa) and display the most significant amino acidic discrepancy between the three Clades. ANRs from Clade-2 are highly conserved among AE pathogens (100% coverage, 93.5% identity) (Fig. 1B). However, Rnr and ANRVibrio are the most distantly related ANRs to the archetype Aar, sharing only 25% identity.
Characterization of the Rnr regulon in EPEC. Since enteric pathogens use distinct mechanisms to colonize and invade their hosts, and the fact that the amino acidic identity between ANR members differs significantly among pathogens (Fig. 1B), it is uncertain whether ANR accomplishes the same regulatory function in different enteric pathogens. We sought to determine this gap in knowledge by dissecting the regulon and biological role of Rnr in attaching and effacing pathogens, which shares only ~ 25% of amino acid identity with the archetype Aar. Accordingly, the rnr gene was deleted in the prototype EPEC O127:H6 strain E2348/69, and RNA-seq determined its transcriptome. For these experiments, the wild-type (WT) EPEC, its isogenic EPECrnr mutant, and EPECrnr complemented with rnr in-trans [EPECrnr(pRnr)] were grown in DMEM-high glucose for six hours to activate the expression of Rnr. Subsequently, total RNA from all strains was extracted and processed for cDNA synthesis, library construction, and DNA sequencing by CD genomics (NY, USA), as indicated in the material and methods. Bioinformatics analysis revealed approximately 500 genes that were differentially expressed (DEGs) in the Rnr regulon (+/- 1.5 fold, p < 0.05) (Fig. 2 panels C-F). The majority of Rnr-regulated genes were located in the chromosome of EPEC E2348/69 (Fig. 2A) and associated with six major functional categories: genes involved in metabolism (46%), protein transport (11%), regulation (5%), ribosomal activities (7%), virulence (including bacterial adherence and motility) (3%), and other functions (28%) (Fig. 2, panels B, G-L) and Supplementary data 1). 8% of Rnr-regulated genes were encoded in the LEE pathogenicity island located in the chromosome and 6% in the pMAR2 plasmid (GenBank FM180569.1) (Fig. 2A).
Ler and HNS were identified in the Rnr-regulon with 26 other regulatory proteins (Fig. 2G and Supplementary Fig. S1). Most of these regulators belong to the AraC/XylS family, including PerA, EutR, MelR, AdiY, YdiP, and GadX. Interestingly, we previously found that Aar also regulates GadX in response to the acid environment in EAEC 26,27. Other important regulators under Rnr control are BssG, FimG, and CsgD, associated with biofilms or bacterial adherence (Fig. 2G).
In addition, we observed that a large number of genes regulated by Rnr are involved in metabolism (~ 200 genes), including the Rut operon (Supplementary Fig. S1), intricate in the degradation of exogenous pyrimidines as the sole nitrogen source, and the arginine succinyltransferase pathway which uses arginine as a source of carbon and nitrogen. Numerous genes of these operons are also under the NtrC control28,29.
Rnr negatively regulates the locus of enterocyte effacement pathogenicity island (LEE-PAI). The LEE-PAI (~ 36 kb) is composed of 42 genes and seven operons encoding the T3SS (Fig. 3A). Our transcriptomic data shows increased expression of LEE genes in EPECrnr, whose complementation in trans with pRnr plasmid restored the expression of genes to comparable wild-type levels (Fig. 3B). Among those, twenty genes encode for proteins that form the core of the T3SS were ~ 2 fold upregulated, including proteins of the basal body (EscC, EscD, EscJ); inner membrane machinery (EscV, EscR, EscS, EscT, EscU); needle tip and translocon (EspD, EspB, EspA) and the EscN ATPse. The intimin gene (eaeA) and its receptor tir were also upregulated (~ 5 fold), which are associated with EPEC intimate adherence to host cells (Fig. 3B). Several LEE regulators (Ler, GrlA and GrlR) located in the LEE-PAI were upregulated (~ 2 to 50 fold) in EPECrnr as compared to WT and complemented EPECrnr strain (Fig. 3B).
To validate our transcriptome dataset, nine genes were selected based on their relevance in the virulence of EPEC for qRT-PCR analysis. We analyzed genes encoding structural proteins of the T3SS apparatus (espA, espB, and espD), regulatory proteins (hns, ler, and perA), and proteins involved in adherence (tir, eaeA, and bfpA). For qRT-PCR experiments, WT EPEC, EPECrnr, and EPECrnr(pRnr) strains were grown in DMEM-high glucose for 2, 3, and 4 h, and total RNA was isolated and prepared for gene expression analysis. Our data showed that 8 out of 9 analyzed genes exhibited higher levels of expression in the EPECrnr mutant (~ 3 fold) compared to WT after 4 h of growth (middle log phase) (Fig. 4).
We previously found that Aar increases the expression of AAF fimbrial genes in early EAEC growth stages by decreasing the expression of HNS, which acts as a repressor of AAF expression; however, when Aar is increased, it acts as a negative regulator of AAF by inactivating AggR, the positive AraC/XyS regulator of AAF30. Similar findings were observed in EPEC with bfpA and eaeA (Fig. 4G and 4I), which are regulated by PerA and HNS, respectively, and these, in turn, are regulated by Rnr (Fig. 4A and 4H).
We sought to determine if changes observed at the transcriptional level also correlated with changes at the protein level. Accordingly, EPEC derivatives were grown in DMEM until the late exponential growth phase (OD600 ηm of 1.0), and the abundance of T3SS structural proteins (EspA and EspB) (Fig. 5A) were evaluated in whole-cell (Fig. 5B) and supernatant proteins (Fig. 5C) by western blot using specific polyclonal antibodies for EspA, and EspB. GroEL was used as an internal loading control for whole bacterial preps and as an indicator of cytoplasmic protein contamination in the secreted proteins fraction. In agreement with our transcriptomic data, we found that deletion of rnr correlated with an increased amount of EspA and EspB in whole bacterial and supernatant preps (Fig. 5B, C, D). In contrast, complementation of EPECrnr with either pAar or pRnr plasmids drastically reduced the expression of EspA and EspB in whole-cell preps and supernatants (Fig. 5B, C, D). Taken together, our findings indicate that members of the ANR family (Rnr and Aar) can regulate a variety of AraC/XylS and HNS regulators, including those controlling the LEE PAI in AE pathogens.
Rnr protein directly interacts with HNS and Ler proteins. We previously showed that Aar interacts with HNS global repressor affecting its regulatory activity 26. Therefore, we sought to determine if Rnr can interact with HNS and Ler, both members of the HNS family. The Bacterial-two hybrid (BACTH) system is broadly used to scrutinize protein interactions between regulatory proteins31, and we have successfully used this approach to examine interactions between Aar-AggR and Aar-HNS26,27. Thus, we used the BACTH system to investigate interactions between Rnr, HNS, and Ler. Accordingly, rnr, hns, and ler genes were fused to T25 and T18 fragments of the catalytic domain of Bordetella pertussis adenylate cyclase, expressed in plasmids pKNT25 and pUT18, respectively (Fig. 6A)31. The resulting plasmids were co-transformed in different combinations of pUT18 and pKNT25 derivatives into the reporter strain E. coli BTH101. Remarkably, we observed protein-protein interactions of Rnr with members of HNS family; HNS, and Ler in the BACTH system manifested by the appearance of a moderate to intense blue color on agar plates (Fig. 6B) and by quantification of the β-galactosidase activity (Fig. 6C). Taken together our findings suggest that Rnr is regulating gene expression of the LEE-PAI by direct interaction with Ler and HNS global regulators.
Rnr and Aar are interchangeable for regulating the T3SS in EPEC and fimbriae in EAEC. Since the heterologous expression of Aar in EPEC is capable of downregulating the expression of T3SS in EPEC in our previous experiment (Fig. 5B and 5C), we sought to determine whether the heterologous expression of Rnr in EAEC042 downregulates the expression of the AggR-regulated AAF fimbriae, the main virulence factor of EAEC associated with host-interactions. Accordingly, we analyzed whole-cell proteins from EAEC derivatives expressing Aar and Rnr by SDS-PAGE and Western blot (Fig. 7). As expected, we found that Rnr was able to downregulate the expression of the major AAF fimbria subunit, AafA, in EAEC (Fig. 7), suggesting that despite the low homology between Aar and Rnr, they may possess structural features that allow function conservation between distantly related ANR members.
Aar and Rnr negatively impact intestinal colonization in their respective pathogens. Despite the extensive molecular characterization of ANR in EAEC, its role in bacterial pathogenesis is not entirely understood, partly due to the lack of adequate animal models for E. coli pathogens. Since Aar downregulates the expression of AggR-regulated AAF fimbria in EAEC (Fig. 7)25, and Rnr downregulates PerA-regulated BfpA and genes associated with intimate adherence mediated by T3SS in EPEC (Fig. 4), we sought to determine the impact of ANR regulation in bacterial adherence and intestinal colonization.
Human intestinal organoids have become the gold standard for studying host-pathogen interactions and have been successfully used to investigate essential features of EAEC and EPEC pathogenesis32,33. We, therefore, used this relevant intestinal model to examine the role of ANR in bacterial colonization. For these experiments, human intestinal colonoid monolayers were infected with parental 042, 042aar, and 042aar(pAar) at 37oC for 6 h, and bacterial adherence was analyzed by confocal microscopy (Fig. 8A-P). We observed that the deletion of aar significantly increases bacterial colonization in human colonoids compared to the parental strain (Fig. 8C, 8G, 8K, and 8O). In agreement with the negative role of Aar, complementation of 042aar with the pAar plasmid drastically reduced biofilm formation (Figs. 8D, 8H, 8L, and 8P) as judged by the enumeration of bacterial cells on colonoids (Fig. 8R). Moreover, microscopic examination of bacterial biofilms revealed increased bacterial aggregation in colonoids infected with 042aar than parental 042 strain (Fig. 8N, O).
We next determined whether Rnr impacts EPEC intestinal colonization and the formation of AE lesions in the human intestine. Accordingly, human cell monolayers were infected with 106 CFU of EPEC E2348/69, EPECrnr, and EPECrnr(pRnr) for 6 h at 37oC. Subsequently, infected cells and uninfected controls were analyzed for EPEC adherence and formation of AE lesions by confocal microscopy (Fig. 9). The confocal images were pixel-quantified as previously reported (Fig. 9M)32. We observed a more significant number of adhered EPECrnr strain on intestinal cell monolayers and which correlated with a greater number of AE lesions compared to the parental EPEC strain (Fig. 9A, 9B, 9D, 9E, 9G, 9H and 9M). Complementation of EPECrnr with the pRnr plasmid drastically reduced bacterial adherence (Fig. 9C, 9F and 9I) and the number of AE lesions on intestinal cell monolayers as judged by actin polymerization beneath the adherent bacteria (Fig. 9L). Although AE lesions were observed on cells infected with all EPEC strains, cells infected with EPECrnr(pRnr), which overexpresses Rnr, exhibit smaller actin pedestals than the WT or EPECrnr strains as judged by qualitative analysis of confocal images (Fig. 9J, 9K and 9L). Taken together, our data suggest a central role of ANR in modulating intestinal colonization by diarrheagenic pathogens (Figs. 8 and 9).