Evaluation of Incompatibility Group I1 (IncI1) Plasmid-Containing Salmonella enterica and Assessment of the Plasmids in Bacteriocin Production and Biofilm Development

Mobile genetic elements, such as plasmids, can potentially increase the ability of bacteria to infect and persist in vertebrate host cells. IncI1 plasmids are widely distributed in Salmonella from food animal sources and associated with clinically important strains. These plasmids often encode antimicrobial resistance; however, little is known about their impact on the virulence of Salmonella strains. To assess the potential impact of the plasmids on virulence, 43 IncI1-positive Salmonella isolates from human and animal sources were subjected to whole genome sequence (WGS) analyses and evaluated for their abilities to invade and persist for 48 h in Caco-2 human intestinal epithelial cells, form biofilms and encode bacteriocins. Draft WGS data were submitted to predict the presence of virulence and antimicrobial resistance genes, plasmid replicon types present, conduct plasmid multilocus sequence typing (pMLST), and core genome MLST (cgMLST) in the isolates. Caco-2 cells were infected with Salmonella strains and incubated for both one and 48 h for the invasion and persistence assays, respectively. Additionally, Salmonella isolates and IncI1 plasmid carrying transconjugants (n = 12) generated in Escherichia coli were assessed for their ability to produce biofilms and bacteriocin inhibition of growth of other bacteria. All Salmonella isolates infected Caco-2 cells and persisted in the cells at 48 hrs. Persistent cell counts were observed to be significantly higher than invasion assay cell counts in 26% of the isolates. Among the IncI1 plasmids, there were 18 pMLST types. Nearly 35% (n = 15) of Salmonella isolates produced biofilms; however, none of the IncI1-positive transconjugants produced increased biofilms compared to the recipient. Approximately 65% (n = 28) of isolates and 67% (n = 8) of IncI1-positive transconjugants were able to inhibit growth of at least one E. coli strain; however, none inhibited the growth of strains from species other than E. coli. The study characterized IncI1 positive Salmonella isolates and provided evidence about the potential contributions of IncI1 plasmids virulence phenotypes and areas where they do not. These findings should allow for more focused efforts to assess the impact of plasmids on bacterial pathophysiology and human health.

Mobile genetic elements, such as plasmids, can potentially increase the ability of bacteria to infect and persist in vertebrate host cells. IncI1 plasmids are widely distributed in Salmonella from food animal sources and associated with clinically important strains. These plasmids often encode antimicrobial resistance; however, little is known about their impact on the virulence of Salmonella strains. To assess the potential impact of the plasmids on virulence, 43 IncI1-positive Salmonella isolates from human and animal sources were subjected to whole genome sequence (WGS) analyses and evaluated for their abilities to invade and persist for 48 h in Caco-2 human intestinal epithelial cells, form biofilms and encode bacteriocins. Draft WGS data were submitted to predict the presence of virulence and antimicrobial resistance genes, plasmid replicon types present, conduct plasmid multilocus sequence typing (pMLST), and core genome MLST (cgMLST) in the isolates. Caco-2 cells were infected with Salmonella strains and incubated for both one and 48 h for the invasion and persistence assays, respectively. Additionally, Salmonella isolates and IncI1 plasmid carrying transconjugants (n = 12) generated in Escherichia coli were assessed for their ability to produce biofilms and bacteriocin inhibition of growth of other bacteria. All Salmonella isolates infected Caco-2 cells and persisted in the cells at 48 hrs. Persistent cell counts were observed to be significantly higher than invasion assay cell counts in 26% of the isolates. Among the IncI1 plasmids, there were 18 pMLST types. Nearly 35% (n = 15) of Salmonella isolates produced biofilms; however, none of the IncI1-positive transconjugants produced increased biofilms compared to the recipient. Approximately 65% (n = 28) of isolates and 67% (n = 8) of IncI1-positive transconjugants were able to inhibit growth of at least one E. coli strain; however, none inhibited the growth of strains from species other than E. coli.

INTRODUCTION
Salmonella enterica is estimated to cause more than one million infections resulting in 400 deaths per year in the U.S. (1). Most Salmonella infections are self-limiting and resolve with symptomatic treatment (2). In some cases, such as co-infection with other bacteria, infection of an immunocompromised host and/or infection by highly virulent strains of Salmonella can lead to more severe, invasive and sometimes fatal infections (3). In such cases, patients typically need antimicrobial treatment to recover from disease. The annual economic impact of all Salmonella infections has been estimated to be up to 4.4 billion dollars due to the cost of treatment, loss of wages, and quality of life (4). Most Salmonella infections are enteric in nature and foodborne salmonellosis is an important economic and public health concern (5, 6). The origin of foodborne salmonellosis can be eggs, poultry products, meat, and fresh produce (5)(6)(7)(8)(9)(10). Salmonella is relatively widespread among poultry and other food animals (11,12).
Salmonella maintains its high presence in diverse hosts through genetic plasticity (13), which allows Salmonella to alter its genetic composition to adapt to changing environmental conditions. This plasticity can be achieved with the help of mobile genetic elements such as plasmids. Plasmids that have been characterized in Salmonella can carry genes associated with increased antimicrobial resistance and virulence for their hosts (5,14,15). Plasmid transfer between bacteria leads to the potential for rapid horizontal spread of genes among bacteria. If these genes influence phenotypic characteristics, such as antimicrobial resistance and virulence, then a better understanding of genetic determinants could be vital to the management of foodborne illnesses arising from Salmonella (5).
Plasmids can be grouped based on their incompatibility to co-exist in the same cell (16). Plasmid incompatibility assays for typing are based on phenomena that prevent coexistence of plasmids with the same replication and division mechanisms in a bacterium (17). Incompatibility group I1 (IncI1) plasmids are commonly found in enteric bacteria from food animal sources and are associated with clinically relevant strains. They are known for their potential to carry and disseminate antimicrobial resistance genes among enteric pathogens (18,19). For example, many of the plasmids associated with the dissemination of genes encoding resistance to ceftriaxone, an antimicrobial agent used in management of severe Salmonella infections, have been reported to be IncI1 plasmids (20).
Similarly, the spread of virulence-associated genes via plasmids could lead to Salmonella-related illnesses difficult to manage, which raises a potential public health concern. Bacterial virulence genes encode factors/molecules that aid a bacterium in colonization of host niche; entry into, survival within and exit from the host cell; evasion or suppression of the host immune response to the bacterium; and/or obtaining limited nutrients from the host environment (21,22). These factors can include biofilm formation that aid in niche colonization, bacteriocin production that can limit niche competition, nutrient acquisition, secretion systems that facilitate bacterial uptake into host cells and improve intracellular survival, and mechanisms that regulate the expression of virulence mechanisms (22). Genes potentially associated with virulence have been identified on IncI1 plasmids (23); however very few, if any studies have been conducted to directly evaluate IncI1 plasmids and their virulence potential in Salmonella. Thus, the objectives of this study were to characterize Salmonella strains that harbor IncI1 plasmids, both for the presence of antimicrobial resistance and virulence genes and the potential to invade and persist in human intestinal epithelial cells; and to specifically assess the impact of IncI1 plasmids on the ability to produce biofilms and encode bacteriocins that could provide a selective colonization advantage for the strains that carry them.

Biofilm Production Assay
The ability of Salmonella strains, recipient E. coli J53 and their respective transconjugants to produce biofilms was evaluated using the polystyrene microtiter plate method (36). Briefly, each strain was grown overnight in tryptic soy broth (TSB) at 37 • C. The following day, the cells were diluted 1:100 in TSB and 200 µL of each bacterial suspension were inoculated into three wells of sterile 96-well flat-bottomed polystyrene plates. Negative control wells containing only TSB were also inoculated in triplicate. The plates were covered and incubated statically for 24 h at 37 • C. After incubation, the growth medium was removed, and the wells were washed three times with 250 µL of phosphate buffered saline (PBS). The attached bacteria were then fixed with 200 µL of methanol (Sigma-Aldrich, St. Louis, MO) per well for 15 min, after which the plates were emptied and dried at room temperature. The plates were stained with 50 µL per well of crystal violet stain (0.41%; Fisher Diagnostics, Middletown, VA) for 5 min. The stain was removed, and the plate was washed under running tap water and subsequently airdried. Next, 200 µL of 33% glacial acetic acid (Sigma-Aldrich) was added to the wells and the optical density (OD) of each well was measured at 550 nm in a SpectrMax Absorbance Reader (Molecular Devices, San Jose, CA). For analysis, the absorbance of the blank was subtracted from the wells to compensate for background absorbance. The experiments were repeated and the arithmetic mean value and standard deviations of the OD values for the six replicates was calculated. The results were further analyzed using Excel (Microsoft, Redmond, WA) using a t-test for significant difference (p < 0.05) in biofilm formation relative to the negative control and/or E. coli J53 for the transconjugants.

Bacteriocin Inhibition Assay
The ability of Salmonella and transconjugant strains to produce colicins was evaluated by assessing growth inhibition of E. coli J53 (24) and the P. aeruginosa, E. faecalis, S. aureus, S. enterica, K. pneumoniae, E. cloacae, and the eight additional E. coli strains described above. Each of these isolates was grown overnight on tryptic soy agar with 5% sheep's blood and incubated overnight at 37 • C. The following day bacterial growth was removed from the plates with a sterile swab and suspended in sterile demineralized water and the bacterial cell concentration adjusted to approximately 0.5 McFarland standard. The suspension was swabbed for confluence on blood agar plates and 5 µl of a Salmonella or E. coli transconjugant inoculum were prepared by suspending the equivalent of 5-6 bacterial colonies in 200 µl of sterile water, was spotted onto swabbed plate. The plates were incubated at 37 • C for 16 to 18 hrs and then examined for growth inhibition of the bacterial lawn adjacent to the Salmonella or transconjugant growth. Inhibition was identified by a complete zone of growth inhibition around the Salmonella or transconjugant spots.

Bacteriocin-Associated Gene Detection
For the wildtype Salmonella strains, bacteriocin-associated genes were identified from whole genome sequencing data using PATRIC and BLAST searching against the colicin genes. For the transconjugants, PCR reactions were carried out to detect the presence of 10 different bacteriocin-related genes (cia, cib, colA, colD, colE1, colM, cvaA, cvaB, cvaC, and imm). The primer sequences are listed in Table 2. Three to four colonies were collected from the plates and suspended in 200 µl of sterile water. The cells were lysed using a boiling method to prepare DNA template (37). PCR reactions included 2X Master Mix (Promega, Madison, WI), primer pairs (10 pmol), 2.5 µl of boiled template and sterile water to the final 25 µl volume. The genes were amplified with the following steps: denaturation at 95 • C for 5 min, 30 cycles of denaturation at 95 • C for 30 secs, annealing at either 52 • C (cvaA, cvaB, and cvaC), 58 • C (cib and imm) or 60 • C (cia, colA, colD, colE1, and colM) for 30 secs, and extension at 72 • C for 90 s, and a final extension at 72 • C for 7 min. The resulting PCR products were separated on 2% E-gels (Invitrogen, Carlsbad, CA) and visualized under UV-light using a Gel-Doc XR system (Bio-Rad, Hercules, CA).
pMLST For the wildtype Salmonella strains, pMLST profiles were identified from whole genome sequencing data. For the transconjugants, DNA template were prepared as previously described above for PCR and the repI, ardA, trbA, sogS, and pilL gene loci were amplified following the methods outlined previously (38). The resulting PCR products were separated on 2% E-gels, visualized under UV-light and the positive PCR reactions purified using the QIAquick PCR Purification Kit (Qiagen). The purified PCR products were sequenced at the University of Arkansas for Medical Sciences Sequencing Core Facility. The fasta files from WGS and individual allelespecific sequencing reads were submitted for identification of Frontiers in Veterinary Science | www.frontiersin.org

Invasion Assay
Caco-2 cells were then infected with a 10 times greater number of Salmonella (i.e., multiplicity of infection 10:1). After infection, the cells were incubated at 37 • C in 5% CO 2 for 1 h. Gentamicin (200 µg/ml) was added to each well, and the cells were incubated for another hour. Cells were subsequently washed three times with PBS and then lysed with 0.1% 4 • C Triton-X and suspensions were serially diluted followed by plating on LB agar plates.
After overnight incubation at 37 • C, the cells on the plates were enumerated to quantify the number of bacteria that invaded the Caco-2 cells (39).

Persistence Assay
After an initial 1-h incubation, gentamicin (100 µg/ml) was added to each well and the plates were incubated for 48 h. After which, the cells were washed, lysed and plated as described for the invasion assay.

Statistical Analyses
For each of the invasion and persistence experiments, statistical analyses were conducted using Microsoft Excel (ver. 2016, Redmond, WA). To normalize the impact of inoculum differences across the experiments, each bacterial cell count was divided by respective inoculum dose to obtain the relative invasion or persistence ratios. The average of the ratios of the six replicates was calculated along with the standard deviation for each set. Differences between the cell survival of invasion and persistence were evaluated to facilitate an understanding the contribution of the different plasmids using two-tailed paired T-test and p < 0.05 was considered as a significant difference between two groups compared.

RESULTS
Forty-three IncI1 positive S. enterica isolates, representing eight different serovars were examined in the present study (Tables 1, 3). From these strains, 12 transconjugant containing IncI1 plasmids, but lacking other known plasmids were produced in recipient E. coli J53 and were further evaluated for the impact of IncI1 plasmids on bacterial physiology ( Table 4). Each of the Salmonella isolates underwent WGS analyses and the sequencing assemblies were deposited in GenBank under the accession numbers shown in Table 1. The assembled sequencing contigs were submitted to the ResFinder program to identify predicted antimicrobial resistance genes (Supplemental Table 1). All isolates were positive for aac(6 ′ )-Iaa, which is not associated with a resistance phenotype in Salmonella (40) and not included in the table. All but three isolates (856, 891, and N865), were positive for at least one additional antimicrobial resistance gene. The plasmid content of each strain was assessed using PlasmidFinder and the results presented in Table 1. All strains carried the IncI1 replicon and most carried additional identified replicon types, including IncA/C (42%, n = 18), IncHI2 (33%, n = 14), ColpVC (26%, n = 11), IncX1 (23%, n = 10), and IncFIB (16%, n = 7) ( Table 1). In addition, the WGS data was analyzed using cgMLST and PATRIC to characterize the bacterial genomes. PATRIC identified putative virulence factors by comparison to their PATRIC-VF database that matched

For the Bacteriocin Inhibition Assay: "+" indicates complete inhibition and "-" indicates no apparent inhibition. For the Bacteriocin Gene Detection: "+" indicates the presence and "-" indicates a absence of the gene in WGS data.
Frontiers in Veterinary Science | www.frontiersin.org  sequence annotations to the predicted Salmonella virulence factors. Supplemental Table 2 indicates the putative virulence factors predicted in the strains. The result of phylogenetic analysis based on the presence of the virulence factors is presented in Figure 1. Figure 1 also contains the data for the average ratio of bacterial persistence counts to the inoculum in Caco-2 cells for the respective isolates, which can serve as a phenotypic marker of potential virulence. The summary results of the cgMLST analyses are shown in Table 1, there were 34 different cgMLST types, with the largest groups being ST 18855 and 43318 (each 9%, n = 4). Supplemental Figure 1 shows a cladogram based on the cgMLST SNP analyses indicating a representation of the genetic relatedness of the strains. The results of the IncI1 pMLST experiments of the Salmonella isolates identified 18 different allele profiles ( Table 1). The most common STs identified were 26 (26%, n = 11), 12 (23%, n = 10), and an undefined ST with the allele profile (3, 2, 17, 3, 3; Table 1, 9%, n = 4). For seven of the isolates not all five alleles were detected; however, three (7%) and two (5%) of these isolates had shared overall profiles ( Table 1). Among the 12 IncI1 transconjugants evaluate, there were six different allele profiles, including three that did not have a previously identified STs. The most common STs were 26 (50%, n = 6) and 12 (17%, n = 2) ( Table 4).
The results of the bacteriocin inhibition assays and PCR results for the Salmonella isolates are shown in Table 1. None of the IncI1-positive Salmonella isolates inhibited the non-E. coli strains, including P. aeruginosa, E. faecalis, S. aureus, S. enterica, K. pneumoniae, and E. cloacae. As noted previously (24), the majority of isolates (65%, 28/43) were able to inhibit the growth of E. coli 53 (27). For the other seven E. coli strains tested, the results were varied, none of the isolates inhibited E. coli ATCC 25922 and veterinary clinical strains 164, 524, 550, and 586. E. coli O157:H7 isolate 734, and veterinary E. coli isolates 542 and 590 were inhibited by Salmonella isolates N134 and N136 ( Table 3). The detection of the bacteriocin-associated genes in the Salmonella isolates are shown in Table 3. All of the isolates were positive for the colicin receptor genes cirA and omr and the colicin v-associated gene cvpA, with the exception of isolate 67 for cvpA. Among the predicted bacteriocin-encoding gene families, colicin lysis protein (clp: n = 1, 2%), channel-forming colicins (cib: n = 32, 74% and cia: n = 2, 5%), colicin production (cvpA: n = 42, 98%, cvaC: n = 3, 7%, colE: n = 1, 2%) were detected. Two different colicin immunity proteins were predicted among the strains (NCBI protein ID: Q2VNY0: n = 26, 60%) and (protein ID: V1L9U1: n = 4, 9%). All strains that were able to inhibit E. coli J53 growth carried either cib (n = 27, 96%) and/or cia (n = 2, 7%), and with the exception of isolates 67 and 76 carried the colicin b immunity (imm) gene ( Table 3). Isolate 76 was negative for cib, but positive for cia. Five strains were stains were positive for cib but did not inhibit E. coli J53, each of these lacked the colicin immunity gene (Q2VNY0). Isolates N134 and N136, which inhibited the greatest number of E. coli were positive for cib, cvpA, cvaC, and imm. Isolate 991, was positive for cib, cvpA, and cvaC, but negative for imm and subsequently there was no observed inhibition of any of the bacterial strains.  Table 2. In addition, the calculation of the average ratio of persistence cell counts to the inoculum levels are shown in the right-hand column. The color shading of the values is based on the ratio, with green at the lower levels of persistence through red at the higher levels.
Of the 12 IncI1 transconjugants, eight (67%) were able to inhibit the growth of E. coli J53, but did not inhibit the growth of any of the other strains tested, including all other E. coli ( Table 4). All eight of the transconjugants that inhibited E. coli J53 carried the cib and imm genes, while those that did not inhibit growth were negative for each of the bacteriocin genes tested ( Table 4). Each of the strains that inhibited the growth of E. coli J53 were either ST26 (n = 5), ST12 (n = 2) or an unclassified ST (n = 1).
The biofilm experiments demonstrated that there was variability of the wild-type strains to form biofilms. Fifteen of the 43 (35%) isolates had significantly increased absorbance versus the negative control, indicating an ability to form biofilms. All 12 transconjugants appeared to have biofilm formation ability; however only two (17%) reached the level of significance. The recipient strain E. coli J53 also had non-significantly higher absorbance (p = 0.11). When the biofilm forming ability of the transconjugants was compared to E. coli J53, none had significantly greater biofilm ability than J53 alone, indicating that the potential contribution of the plasmid(s) to biofilm formation is minimal.
The results of the invasion and persistence assays for the Salmonella isolates are displayed in Figure 2. All isolates tested could invade Caco-2 cells. Twelve of the forty-three (28%) isolates demonstrated significantly higher bacterial cell counts after 48 hrs of infection compared to their initial inoculum levels, which is indicative of intracellular growth. An additional 12 isolates (26%) exhibited increased cell numbers after 48 h that did not reach statistical significance. Seven of the isolates (16%) examined, yielded persistent cell populations significantly lower than those observed in the invasion assays (p < 0.05). An additional 12 (28%) showed decreased persistence that did not reach statistical significance.
When the ratio of persistence to the initial inoculum were examined based on serotype (Figure 1), the isolates from S. Newport demonstrated a significantly lower survival in host cells than those of serovars Typhimurium (p = 0.00015) and Heidelberg (p = 0.0036), while S. Kentucky strains had significantly lower persistence than those of S. Typhimurium (p = 0.0098), but not S. Heidelberg (p = 0.089). Similarly, when the sources of the isolates were compared, isolates collected from swine showed a significantly lower persistence in Caco-2 cells than those isolated from chickens (p = 0.007), turkeys (p = 0.015), and human patients (p = 0.017). There were no significant differences between isolates containing the three most common pMLST types (noted above) and increased persistence within the Caco-2 cells.

DISCUSSION
Salmonella enterica can infect several different hosts, with their transmission being primarily via the fecal-oral route (5). Several Salmonella strains have the ability to form biofilms that allow groups of cells to be more resilient to killing and persist in the environment, potentially facilitating contamination of food products, animal production environments and within host organisms (36). After ingestion, Salmonella enters the intestinal tract and can cause gastroenteritis, the most common illness caused by Salmonella infection. To cause infection, Salmonella must compete with intestinal bacteria and colonize the intestinal epithelia. The production of bacteriocins can provide a competitive advantage for the Salmonella against members of the microbiota, such as E. coli (41,42). Intestinal epithelial cells provide Salmonella a niche for colonization in their host (43). A type 3 secretion system (T3SS) encoded by Salmonella pathogenicity island 1 (SPI-1) facilitates the uptake of the Salmonella into the host cells (44). After entering epithelial cells Salmonella become engulfed into Salmonella-containing vacuoles (SCV). Within the SCV, Salmonella can move toward the basal side of the epithelium where they can be taken up by macrophages or dendritic cells (45). Through these immune cells, Salmonella can enter the lymphatic and systemic circulation leading to systemic infection (11). Therefore, intestinal epithelial cells play an important role in the pathogenesis of Salmonella infection. The Caco-2 cell line is derived from human intestinal epithelial cells and has been commonly used to study Salmonella pathogenesis (39,43,46). Caco-2 cells in combination with the WGS analyses were used in this study to gain a better understanding of IncI1 plasmid-positive Salmonella virulence, using invasion, persistence and intracellular proliferation in intestinal epithelial cells as the model system.
In the virulence assay, we observed that persistent cell populations after 48 hrs of infection were higher than the invasion cell counts in many of the isolates examined (Figure 2). All of the isolates in the study contained the genes encoding the SPI-1 T3SS apparatus (Supplemental Table 2), which likely facilitated the uptake of the Salmonella into the host cells (44). Among the different serovars tested, isolates from serovar Typhimurium appeared to have increased persistence capacity by showing higher average persistent cell counts, in general. Salmonella Typhimurium is known to be able to persist in the SCV, as well as the host cell cytosol (47,48). Survival in the SCV is aided by a second T3SS encoded on SPI-2, whose effectors are important for intracellular survival. The presence of the SPI-2 encoded genes were identified in the strains and can contribute to the survival and proliferation within the host cells (49). Differences in the SPI-2 regulation and expression could give S. Typhimurium an advantage for persistence within intestinal epithelial cells, potentially allowing for more efficient intracellular multiplication than other serovars, however their exact contributions need to be further quantified. While statistical significance was observed, these numbers need to be evaluated in the context that there were uneven distributions of isolates from different sources and serotypes, with S. Newport and swine being on the lower end of representation, thus the impact of each isolate examined could disproportionately alter the overall results.
To attempt to get a better handle on the potential specific contribution that the IncI1 plasmids may have had on persistence and other characteristics examined, we used pMLST to subtype the IncI1 plasmids in the Salmonella isolates and the transconjugants. If there were particular IncI1 plasmid types that were associated with the more persistent phenotypes, pMLST might help discern these differences. Based on the results of the Caco-2 tissue culture experiments and corresponding results of pMLST, there did not appear to a connection between plasmid type and persistence in the epithelial cells. Of the strains that had significantly increased and decreased persistence or intracellular multiplication compared to the initial inoculum, there were overlaps in plasmid types, for example, isolates 706 and 93 had the same allele profiles (ST 26), yet opposite persistence profiles, with strain 706 having the highest ratio of cell numbers to inoculum at 48 hrs and 93 the lowest (Figure 1). Ideally, these invasion and persistence experiments would have been completed in the transconjugants or in plasmid-cured Salmonella; however, attempts to cure the plasmids or transfer the plasmids into a well-characterized Salmonella strain to conduct the invasion and persistence experiments were unsuccessful in the present study. Also, the E. coli transconjugants generated were not ideal for the comparative assessment of invasion and persistence because they lack the T3SS important for entry into intestinal epithelial cells. Our ongoing efforts in the lab are focused on optimization of methods for plasmid curing and transfer of plasmids into Salmonella recipients which facilitate further genetic characterization of these and other plasmids. Because isolates with similar plasmid types spanned the spectrum of persistence, from significantly reduced persistence to significantly increased cell numbers, the results indicate that the IncI1 plasmids may not be major contributors to persistence in Salmonella. Likewise, the IncI1 plasmids did not play a significant role in the ability to form biofilms. While some of the strains formed biofilms, there was no consistent trend across the IncI1-positive isolates (Figure 3). Over half of the strains tested did not have significantly more biofilm formation compared to the blank wells which had no cells. Furthermore, in the experiments with transconjugants, the addition of the IncI1 plasmids to the recipient strain did not significantly increase the biofilm producing ability of E. coli J53 (Figure 4).
In addition to the biofilm and persistence studies, this study examined the production of bacteriocins and their potential ability to inhibit multiple species including enteric organisms such as S. enterica, K. pneumoniae, E. cloacae, and E. coli. Our earlier work indicated that several of the Salmonella isolates and transconjugants generated from these strains were able to inhibit E. coli J53 (24), thus in the current study we examined multiple E. coli isolates, the other enteric organisms noted above, and potential pathogens such as P. aeruginosa, E. faecalis, and S. aureus to further explore the potential spectrum of activity. When bacteriocin-associated genes were evaluated, some isolates were positive for multiple bacteriocin genes. Even with the multiple genes, the spectrum of activity for most isolates was narrow, only inhibiting a limited number of E. coli. Two isolates N134 and N136 were able to inhibit four different E. coli, including an E. coli O157:H7 strain. These two isolates contained multiple bacteriocin-associated genes, such as cib, which is commonly associated with IncI1 plasmids (50) and cvaC that is typically associated with IncFIB plasmids (51). Both isolates were positive for both IncI1 and IncFIB plasmids, however neither plasmid appeared to be conjugative for either isolate so were not able to evaluate their impact in the transconjugants.
The development of the transconjugants in E. coli generated using the well-characterized E. coli J53 strain (27,28), helped to determine the impact of the IncI1 plasmids on bacterial inhibition. The majority (67%) of transconjugants that carried IncI1 plasmid were able to inhibit the growth of E. coli J53. All of these transconjugants received the colicin immunity gene, imm and cib bacteriocin gene. The IncI1 plasmids from each of these strains fell into ST 26, 12 or a unique profile (isolate 482). One Salmonella isolate (100) with ST 26 IncI1 plasmid lacked cib or imm and was not able to inhibit the grow of any strains tested, which was contrary to all other ST 26 plasmid-positive isolates. Taken together, these results demonstrate the importance of specific IncI1 plasmids for the dissemination of bacteriocins that may provide a selective advantage in an environment containing certain E. coli populations.
IncI1 plasmids may persist in Salmonella populations due to their contribution to antimicrobial resistance and bacteriocin production but at the same time appear to have minimal metabolic costs on the host strains. Johnson et al. (52) showed that the acquisition of IncI1 plasmids did not significantly add to the fitness cost of a host bacterium and in some cases led to a negative fitness cost (i.e., increased fitness) to the strains that acquired them. In strains carrying IncA/C plasmids, the fitness cost of also acquiring an IncI1 plasmid was no greater than carrying the IncA/C plasmid alone or in some cases a negative cost (52). This may help explain why many of the isolates with the greatest ability to multiply and persist in the epithelial cells (i.e., highest persistence/invasion ratios) carried multiple large plasmids, including IncA/C, IncHI2, IncFIB, and IncX (Figure 2 and Table 1).
Additionally, there is evidence that intestinal inflammation and iron limitation may impact bacteriocin function due the expression of the enterobactin siderophore (iron chelator) receptors on the bacterial cell surface, which can also serve as a receptor for colicin 1b (encoded by cib) (41,53). Results from the PATRIC-VF analyses show that each of the strains have the siderophore receptors. With the increased expression of the colicin receptors by commensal E. coli, it may allow the Salmonella to impair their competitors for nutrients and colonization that facilitates uptake into the intestinal epithelial cells (42). The presence of plasmid-associated genes within the host may also impact the susceptibility to bacteriocins, expression of plasmid-associated siderophores, salmochelin, and aerobactin, by E. coli have been shown to limit the impact of Salmonella in the gastrointestinal tract (54). Interestingly, the two isolates that inhibited the greatest number of E. coli (N134 and N136) carried an IncFIB plasmid in addition to IncI1 and the veterinary E. coli strains that they inhibited (542 and 590) lacked the aerobactin operon, while the four (164, 524, 550, and 586) not inhibited were positive for the aerobactin operon (29), which is commonly carried by IncFIB plasmids (39). The reason for this dichotomy could either be due to E. coli carrying an IncFIB-encoded immunity protein or having the aerobactin machinery to acquire iron without needing to highly express the enterobactin receptor. Thus, there is a continued need to be decipher the biology of plasmids to help understand factors that contribute to the public health challenges faced by consumers.

CONCLUSIONS
Salmonella are important foodborne pathogens in humans and are major concerns for food animal production. Increased virulence of Salmonella and the capability to disseminate virulence determinants rapidly on plasmids, make them a concern for public health. This study was undertaken to characterize IncI1 plasmid carrying Salmonella isolates and identify potential contributions of IncI1 plasmids to virulence. The WGS analyses indicated that the IncI1 plasmids are spread across a generically diverse range of Salmonella including multiple serotypes, cgMLST profiles and having diverse virulence gene profiles. Additionally, there was diversity among the IncI1 plasmids evaluated in the study as evidenced by the multiple pMLST subtypes among the Salmonella isolates evaluated in the study. While these plasmids were not apparent significant contributors to increased biofilm formation or persistence in intestinal epithelial cells, the study does indicate a role of specific IncI1 plasmids to inhibit the growth of other bacteria and potentially contribute to virulence. Bacteriocins can provide a colonization advantage to Salmonella and subsequently facilitate invasion of the intestinal epithelial cells. The study demonstrated that Salmonella containing IncI1 plasmids frequently carry multiple antimicrobial resistance genes and can persist and often multiply with human intestinal epithelial cells. Overall, the study shows that it is important to begin to decipher areas of virulence where IncI1 plasmids likely contribute, while at the same time describing other areas where they do not. These results will allow us to hone in on specific plasmidencoded factors that impact bacterial pathophysiology and human health.

DATA AVAILABILITY
The datasets generated for this study can be found in GenBank, Accession numbers are listed in Table 1.

AUTHOR CONTRIBUTIONS
SR and SF secured the funds for the project. SF, SR, PK, NA, and JH designed the studies. PK, AC, NA, JD, and BK conducted the laboratory studies. SF, PK, SR, NA, YS, and JH conducted the data analyses. PK, AC, SF, and SR wrote the final manuscript. All contributed to the editing and refinement of the finished manuscript.

ACKNOWLEDGMENTS
We would like to thank Dr. Rajesh Nayak, Regulatory Health Scientist, Regulatory Compliance and Risk Management, Food and Drug Administration, National Center for Toxicological Research, for providing some of the isolates used in this study. PK's graduate assistantship was provided by the Center for Advanced Surface Engineering, under the National Science Foundation Grant No. OIA-1457888 and the Arkansas EPSCoR Program, ASSET III. AC and NA are supported through fellowships administered through the Oak Ridge Institute for Science and Engineering. The information in this manuscript is not a formal dissemination of information by FDA and does not represent agency position or policy.

SUPPLEMENTARY MATERIAL
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets. 2019.00298/full#supplementary-material Supplemental Table 1 | Antimicrobial resistance genes detected by WGS analyses.
Supplemental Table 2 | Virulence-associated genes detected from WGS data analyzed using the PATRIC-VF database.
Supplemental Figure 1 | Cladogram of SNP analyses of WGS data. Isolate numbers in red indicate strains that had significantly greater persistence relative to invasion numbers and those in blue have significantly lower persistence numbers.