Francis Muchaamba
Athmanya K. Eshwar
Ueli von Ah
Marc J. A. Stevens
Taurai Tasara- 1Institute for Food Safety and Hygiene, Vetsuisse Faculty, University of Zurich, Zurich, Switzerland
- 2Agroscope, Bern, Switzerland
Listeria monocytogenes associated prosthetic joint infections (PJI) are a rare but increasing clinical problem of listeriosis. We characterized two isolates of the same L. monocytogenes strain isolated within five years of each other from a recurrent human prosthetic joint infection. The two isolates although clonally identical were phenotypically distinct confirming that the original infection strain had evolved within the human host PJI environment giving rise to a phenotypically distinct variant. The recurrent PJI isolate displayed various phenotypic differences compared to the parental original PJI isolate including diminished growth and carbon source metabolism, as well as altered morphology and increased stress sensitivity. The PJI isolates were both diminished in virulence due to an identical truncation mutation in the major virulence regulator PrfA. Genome wide sequence comparison provided conclusive evidence that the two isolates were identical clonal descendants of the same L. monocytogenes strain that had evolved through acquisition of various single nucleotide polymorphisms (SNPs) as well as insertion and deletion events (InDels) during a persistent human PJI. Acquired genetic changes included a specific mutation causing premature stop codon (PMSC) and truncation of RNAse J1 protein. Based on analysis of this naturally truncated as well as other complete RNAse J1 deletion mutants we show that the long-term survival of this specific L. monocytogenes strain within the prosthetic joint might in part be explained by the rnjA PMSC mutation that diminishes virulence and activation of the host immune system in a zebrafish embryo localized infection model. Overall our analysis of this special natural case provides insights into random mutation events and molecular mechanisms that might be associated with the adaptation and short-term evolution of this specific L. monocytogenes strain within a persistent human PJI environment.
Introduction
Listeriosis is a serious foodborne disease caused by Listeria monocytogenes, predominantly affecting people with weakened immune systems including neonates, pregnant women, elderly and immunocompromised individuals leading to serious illness with high mortality rates (European Food Safety Authority [EFSA], 2017; Centers for Disease Control [CDC], 2018; Radoshevich and Cossart, 2018). Although infections usually manifest as meningitis, bacteremia, and feto-maternal complications, listeriosis can in rare cases present as focal infections in different organs including bone and joint (Allerberger and Wagner, 2010; Del Pozo et al., 2013; Chavada et al., 2014; European Food Safety Authority [EFSA], 2017; Radoshevich and Cossart, 2018). Septic arthritis is one of infrequent listeriosis manifestation that usually involves orthopedic implant devices and it is more common in older or immunosuppressed individuals (Cone et al., 2001; Chougle and Narayanaswamy, 2004; Kesteman et al., 2007; Charlier et al., 2012; Del Pozo et al., 2013; Chavada et al., 2014; Bush et al., 2015; Bader et al., 2016).
Listeria monocytogenes associated prosthetic joint infections (PJI) are on the rise probably in part due to the upsurge in prosthetic joint replacements in listeriosis high-risk groups (Charlier et al., 2012; Bush et al., 2015; Bader et al., 2016; Lenguerrand et al., 2018). Such infections are associated with substantial morbidity, impaired joint function, and at times limb amputation (Bush et al., 2015; Lenguerrand et al., 2018). Several antibiotic treatment regimens are reported with ampicillin or amoxicillin, either alone or in combination with gentamicin, being the most commonly used antibiotics (Charlier et al., 2012; Bush et al., 2015). Treatment duration is variable ranging from two weeks to as long as 18 months with some patients being placed on life-long oral antibiosis (Charlier et al., 2012; Bush et al., 2015). Moreover, treatment failures are also quite frequent among Listeria PJI cases receiving antibiotic therapy alone without surgical removal of the infected prosthesis (Charlier et al., 2012; Bush et al., 2015; Bader et al., 2016).
Osteoarticular listeriosis may become more common as the size of the population having specific risk factors related to this condition increases (Bader et al., 2016). Thus, further knowledge of L. monocytogenes involvement in this disease is necessary. We examined two L. monocytogenes strains that were isolated five years apart from a case of recurrent hip PJI in an 84-year-old patient who also received antibiotic treatments for the PJI and chronic obstructive pulmonary disease (COPD) during that period. Strains N843_10 isolated in 2010 and its recurrent strain N843_15 isolated in 2015 were phenotypic variants of the same L. monocytogenes strain. These two strains were characterized using phenome and genome-based approaches to investigate the physiological and molecular adaptive changes that might have evolved during long-term exposure to the human hip prosthetic joint environment resulting in these two phenotypic variants of the same L. monocytogenes clone.
Materials and Methods
Bacterial Strains, Genetic Manipulations, and Culture Conditions
Table 1 lists the strains and plasmids used in this study. Strains N843_10 and N843_15 were isolated by the National Center for Enteropathogenic Bacteria and Listeria (NENT) Switzerland. The ΔrnjA mutants were created through in-frame deletion of the rnjA (lmo1027) gene. A deleted copy for this gene retaining the reading frame, first six and last ten codons, as well as 500 bp of the upstream and downstream flanking sequences was synthesized based on the N843_10 genome sequence (GenScript Biotech, Netherlands). The ΔrnjA DNA fragment was cloned into the pKSV7 plasmid via the EcoR1 and Sa1I sites and used to replace the rnjA chromosomal copies in N843_10 and N2306 strains by homologous recombination as previously described (Smith and Youngman, 1992; Camilli et al., 1993; Schmid et al., 2009). N843_10 and N2306 rnjA locus deletion mutants were confirmed through PCR analysis and DNA sequencing. The Green fluorescent protein (GFP) and mOrange2 fluorescent protein (mO2FP) expressing strain variants were created through site specific PSA-integrase mediated single copy integration of the pPL3-eGFP (Shen and Higgins, 2005) and pIMK-mO2FP plasmids (unpublished and kindly provided by the Loessner lab, Laboratory for Food Microbiology, ETH Zurich) into the tRNA-Arg locus (Lauer et al., 2002; Monk et al., 2008). Bacteria were stored at −80°C in brain heart infusion medium (BHI, Oxoid, United Kingdom) supplemented with 20% glycerol. Strains were initially grown overnight on blood agar or BHI agar plates at 37°C to obtain single colonies, and then cultured twice in 10 ml BHI broth (37°C, 150 rpm) for 16 h generating stationary phase cultures that were routinely used as a starting point for experiments unless otherwise stated.
Table 1. Strains and plasmids used in this study
Quantification of Cell Growth
Stationary cultures prepared from each strain as described above were diluted (1:100) in 10 ml BHI (107 CFU/ml) and incubated at 37°C and 150 rpm. Growth was monitored by viable cell counting and OD600 measurement at defined time points. The experiments performed in duplicate were repeated on three separate occasions. The program DMFit (Baranyi and Roberts, 1994) was used to estimate growth parameters (lag phases and growth rates) from the OD600 data.
Multilocus Sequence Typing (MLST) and Serotyping
MLST based on seven housekeeping gene fragments was performed as described on the Institut Pasteur website1 (last accessed 27 August 2019). The isolates were serotyped with Listeria O and H antisera (Listeria Antisera set, Denka Seiken co., Ltd, Japan), according to the manufacturer’s recommendations and PCR-serogroups were determined using a previously described multiplex PCR assay (Doumith et al., 2004).
Microscopy
100 μl of overnight cultures grown in BHI broth, were fixed onto objective slides, Gram stained and processed for microscopic analysis. Slides were examined and photographed with a Leica DM4000B digital microscope through a 100x/1.3 oil-immersion objective. For electron microscopy, cells were harvested from overnight BHI cultures by centrifugation (6 000 g at 25°C for 5 min), fixed for 2 h at room temperature in 2.5% glutaraldehyde (Electron Microscopy Sciences) buffered in 0.1 M Sodium phosphate buffer (Sigma, Buchs, Switzerland), washed 3 times (0.1 M sodium phosphate buffer), fixed, and stained with 1% osmium tetroxide (Sigma, Buchs, Switzerland). Samples were dehydrated in ascending concentrations of ethanol followed by dehydration with propylene oxide (Sigma, Buchs, Switzerland) and infiltration in 30% and 50% Epon (Epoxy embedding medium, Sigma, Buchs, Switzerland). From each cell pellet, 0.9 mm toluidine blue stained semi thin sections were produced. Representative areas were trimmed and subsequently 90 nm, lead citrate (Merck, Germany) and uranyl acetate (Serva Electrophoresis, Baden-Würtenberg, Germany) contrasted ultrathin sections were produced and viewed under a transmission electron microscope (TEM Phillips CM10) at the Institute for Veterinary Pathology, Zurich, (IVPZ) at the University of Zurich.
Flow Cytometric Analysis
For this analysis, mOrange2 labeled strains (Table 1) were used. 750 μl of the overnight cultures of each strain were mixed with 250 μl of 4% paraformaldehyde and incubated at room temperature for an hour to fix the bacteria. Analysis of the strains was done at the Flow Cytometry Center (University of Zurich) using an Amnis ImageStream X Mk II imaging flow cytometer (Luminex, United States) fitted with a 488 nm excitation laser. Channel 3 that detects fluorescence between 560 to 595 nm was applied for the assay. Ten thousand events were counted for each sample and data was analyzed using the Image Data Exploration and Analysis software (IDEAS) version 6.2.
Phenotypic Microarray Analysis
Strains N843_10 and N843_15 were analyzed for carbon source utilization (PM01 and PM02) as well as resistance to osmotic and pH stress (PM09 and PM10) on Biolog Phenotype Microarrays (PM) (Bochner, 20092). The full list of the tested compounds can be obtained from https://www.biolog.com/wp-content/uploads/2020/04/00A-042-Rev-C-Phenotype-MicroArrays-1-10-Plate-Maps.pdf (accessed July 20, 2020). The experiments were performed in duplicate at 37°C following standard Biolog Inc., protocols with a few modifications as previously described (Muchaamba et al., 2019).
Motility Assay
Five microliter of overnight BHI culture (1 × 109 CFU/ml) of the different strains was spot inoculated on the surface of soft BHI agar (0.25% agar) plates containing 0.05% Triphenyltetrazolium Chloride (TTC) and incubated at 25°C for 48 h. Motility was determined by measuring the diameter of the red zone created by the spreading colony.
Antibiotic Sensitivity
Before each experiment, bacteria were plated by streaking on blood agar plates and grown overnight at 37°C. Tests for antibiotic susceptibility using commercial E test strips against a panel of 8 antimicrobials: amoxicillin and clavulanic acid (Amoxy Clav), cephalothin, tetracycline, ciprofloxacin, penicillin G, sulfamethoxazole, azithromycin and gentamicin, were done in accordance with the manufacturers recommendations (Biomerieux, Lyon, France; Liofilchem, Italy). Briefly, bacteria solutions of 0.5 McFarland standard density were plated to form a lawn on the surface of a Muller Hinton plus blood agar plate. An appropriate antibiotic E test strip was then placed on respective plates and results were assessed after 48 h of incubation at 37°C. The minimum inhibitory concentration (MIC) was defined as the highest concentration at which the colonies touched the E test strip.
Biofilm Assays
Biofilm experiments were performed in Tryptone soy broth (TSB) medium at 37°C. Overnight aliquots of TSB secondary cultures prepared as described above were added (1:40) to fresh TSB medium then 100 μl of this mixture was added to 96-well microtiter plates. Biofilms were grown for 96 h at 37°C. Unbound cells were removed by washing three times with 150 μl sterile deionized water then inversion and tapping of the 96 well plate on absorbent paper. Microplates were dried at 37°C for 30 mins, and adherent cells stained with 150 μl aqueous crystal violet for 20 min. Excess stain was removed by 5 washings with deionized water. The bound stained cells were quantified by addition of 96% ethanol to dissolve the crystal violet and the dissolved stain was measured at an optical density of 595 nm using a Synergy HT OD reader (Biotek Instruments, GmbH, Switzerland). Each biomass was standardized relative to L. monocytogenes N11-1850.
Cell Invasion, Hemolysis and Phosphatidylinositol-Specific Phospholipase C (PI-PLC) Assays
Cell invasion assays using the human enterocyte-like Caco-2 (ATCC® HTB-37TM) cell line and hemolysis using human red blood cells were performed as previously described (Muchaamba et al., 2019). To compare PI-PLC activities the strains were grown on OCLA (Oxoid Chromogenic Listeria Agar) plates. Overnight BHI cultures of each strain were serially diluted in phosphate buffered saline (PBS; 105, 104 and 103 CFU/ml) and spotted (10 μl) on OCLA plates that were incubated at 37°C and visually examined after 48 h for the zone of opacity. L. monocytogenes LL195 and L. innocua JF5051 strains were included as positive and negative controls, respectively. Experiments were conducted in triplicate on three separate occasions.
Zebrafish Microinjection and Localized Infection Assays
Zebrafish husbandry and assays were performed using the Danio rerio wik zebrafish line strains as previously described (Eshwar et al., 2017). Bacteria for microinjection experiments were harvested from secondary stationary phase stage BHI cultures prepared as described above by centrifugation, washed once and diluted to 5 x 108 CFU/ml in DPBS. Two-day post fertilization embryos were injected with approximately 500 CFU in 1–2 nl volume of a bacterial suspension in DPBS into the blood circulation via the caudal vein. Injected CFU numbers were controlled through viable cell counting performed on the microinjection DPBS droplet and five individual embryos immediately disintegrated after microinjection. Post-infection embryos were placed into 24-well plates (one embryo per well) in 1 ml E3 medium per well, incubated at 28°C and observed under a stereomicroscope twice a day up to 72 h post infection (hpi) for signs of disease including developmental delay (especially swim bladder and decreased locomotory rate), necrosis at the inoculation site, subsequent necrosis in other parts of the body, and eventual death. The number of dead larvae was determined visually based on the absence of a heartbeat. The Danio rerio wik zebrafish Fli1: GFP and Danio rerio wik mpeg1: GFP zebrafish line strains with GFP-labeled endothelial cells and macrophages, respectively, were used for the localized infection experiments. Fluorescent-labeled bacteria were prepared and harvested as described above. Two-day post fertilization embryos were injected with approximately 1200 CFU in 1–2 nl volume of a bacterial suspension in DPBS into the otic vesicle to simulate a localized infection. The number of CFU injected and post-infection embryo husbandry was controlled as described above. The embryos were monitored for signs of disease and survival under a stereomicroscope at 8, 24, 48, and 72 hpi. The number of dead larvae was determined visually based on the absence of a heartbeat. At each time point, 10 embryos were fixed and prepared for confocal microscopy. 3D-image stacks of whole mount samples were prepared using a confocal laser-scanning microscope (CLSM, Leica TCS SP8, Leica Microsystems, Heerbrugg, Switzerland). GFP and mOrange were sequentially excited with the 488 nm and 561 nm laser lines, respectively, with emission signals collected within the respective range of wave lengths. 3D image stacks were collected sequentially (to prevent green–red channel crosstalk) according to Nyquist criteria and deconvolved using HuygensPro via the Huygens Remote Manager v2.1.2 (SVI, Netherlands). Images were further analyzed with Imaris 7.6.1 (Bitplane, Zurich, Switzerland). Furthermore, two-day post fertilization zebrafish Fli1: GFP embryos were injected with approximately 2400 CFU in 2 nl volume of a 1:1 mixture of GFP labeled N843_10 and mOrange labeled N845_15 bacterial suspension in DPBS into the otic vesicle to simulate a coinfection. The number of CFU injected and post-infection embryo husbandry was controlled as described above. Plating on selective media to enumerate the individual coinfected strains was done at 0, 8, 24, and 48 hpi.
Genome Analysis
Genomes of L. monocytogenes N843_10 (CP046361), N843_15 (CP046362), N12-1273 (QYFZ00000000), EGDe (NC003210), LL195 (HF558398), Lm3136 (CP013723), Lm3163 (CP013722), N1546 (CP013724), N2306 (CP011004), N16-0044 (CP035187), EGD (HG421741), 10403S (NC_017544), Clip 804259 (AE017262), and F2365 (NC002973) were used in this study. Genomic DNA isolated from the strains N843_10 and N843_15 using the GenElute Bacterial Genomic DNA Kit (Sigma, Buchs, Switzerland) was sequenced using the Pacific Biosciences single-molecule real-time sequencing technology (SMRT) and assembled de novo using the SMRT Analysis 2.3.0 software (ChunLab, Seoul National University). Rapid Annotation Subsystem Technology (RAST) and Seed Viewer3 were used for genome annotation and comparisons. Gepard was used for the dot-plot analysis of the genomes (Krumsiek et al., 2007). MAUVE was used to align the genomes and to derive the coordinates for the positions of the single nucleotide polymorphisms (SNPs), insertions and deletions (InDels) (Darling et al., 2010). Genes of interest were extracted and compared between the genomes using CLC genomics Workbench (Qiagen, Prismet, Denmark) and the BLASTn and BLASTp programs (blast.ncbi.nlm.nih.gov/Blast.cgi). Core genome MLST (cgMLST) analysis was performed using the software package Seqsphere+ 6.0 (Ridom GmbH, Münster, Germany). Complete genomes sequences of N843_10 and N843_15 were compared with a local database containing 567 L. monocytogenes genomes (ILS, unpublished data). The MLST sequence type was determined according to the scheme of the Institut Pasteur4; cgMLST according to scheme at cgMLST.orgm5 (Ruppitsch et al., 2015). Minimal spanning trees were constructed using Seqsphere+ with the option “missing values pairwise ignored.” Relatedness of the two strains (N843_10 and N843_15) was further assessed by SNP comparisons including a selection of other unrelated reference strains. SNPs were identified using parsnp within the harvest suite (Treangen et al., 2014) using standard settings and nucleotide fasta files as input. Each strain was used as a reference strain and compared to the other strains. The output files were converted to variant calling files using harvesttools and a SNP matrix was constructed by taking the sum of the variants compared to the reference strain. The SNP matrix was visualized in a heatmap using clustvis (Metsalu and Vilo, 2015). DuctApe software was used to compare genomes in correlation with PM data, using it to detect genes encoding enzymes that could be involved in the metabolic pathways responsible for the phenotypes observed in carbon source utilization on PM01 and PM02 (Galardini et al., 2014). Genes described in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were considered in this approach. Mutations in selected genes uncovered in WGS were further confirmed through Sanger sequencing at Microsynth AG (Balgach, Switzerland; Supplementary Table S1).
RNA Isolation and Reverse Transcription Quantitative PCR (RT-qPCR)
Secondary bacterial cultures prepared as described above were diluted (1:100) in 50 ml and grown to the late exponential phase (OD600 of 1.4; 109 CFU/ml). Bacteria (1 ml) were harvested in RNA protect Bacteria reagent (Qiagen GmbH, Hilden, Germany) and total RNA was isolated using the RNeasy Plus Mini Kit (Qiagen GmbH, Hilden, Germany) as previously described (Eshwar et al., 2017). RNA yield and quality were assessed using the Quanti Fluor RNA System (Promega, Madison, United States) and the BioAnalyzer (Agilent Technologies, United States). RNA (400 ng and RIN ≥ 7.8) was converted to cDNA using the Quantitect Reverse Transcription Kit (Qiagen GmbH, Hilden, Germany). Samples (2.5 ng cDNAs) were amplified using primers listed in Supplementary Table S1 and the SYBR green I kit (Roche Molecular Diagnostics GmbH, Mannheim, Germany) in the Light Cycler LC480 instrument (Roche Molecular Diagnostics, Rotkreuz, Switzerland). Relative mRNA quantification was performed using the Light Cycler 480 Relative Quantification Software (Roche Molecular Diagnostics). mRNA amounts were normalized using 16S rRNA and expressed relative to a calibrator mRNA sample derived from a stationary phase L. monocytogenes EGDe culture.
Statistical Analyses
All experiments presented were performed independently in duplicate at least three times unless stated otherwise. GraphPad Prism (Version 8.3.0 (328), GraphPad Software, San Diego, CA, United States) was used for the statistical analysis of data. One-way ANOVA with post-hoc Tukey HSD tests was used to assess the statistical significance of differences between the strains. P-values < 0.05 were considered to be statistically significant. PM data was analyzed using the programs ductApe and opm (version 1.3.64) (Galardini et al., 2014; Göker et al., 2016). Area under the curve was the reference parameter for opm analysis, whereas for DuctApe the parameter, activity index (AV) was used as previously described (Galardini et al., 2014; Göker et al., 2016). For each compound tested, the final result was expressed as the mean of two replicates. The bacterium was not able to respire under conditions were AV value was equal to zero, whilst it was able to respire under conditions were the AV values were higher than zero.
Results
N834_15 and N843_10 Are Phenotypic Variants of the Same L. monocytogenes Strain Isolated Five Years Apart From a Recurrent Human Hip PJI
We isolated two L. monocytogenes strains from recurrent human hip PJI episodes recorded five years apart in 2010 (N843_10) and 2015 (N843_15), respectively. Typing of these two isolates based on PCR and seven loci multi-locus sequence typing (MLST) assigned them both to PCR serogroup IIa, MLST sequence type (ST) ST412 and clonal complex (CC) CC412. As such the original (N843_10) and relapse (N843_15) PJI isolates were clonal descendants and derivatives of the same L. monocytogenes strain responsible for the recurrent PJI episodes recorded five years apart. On phenotypic analysis these two isolates were however, phenotypically distinct since compared to the initial PJI isolate N843_10, the recurrent infection isolate N843_15 grew slower on both solid and liquid media as well as aggregated and sedimented in BHI broth (Figures 1A–D). In a swarming motility assays the isolate N843_15 was also significantly less motile than N843_10 (Figure 1E). Both light and electron microscopy examination of liquid and pelleted stationary phase cultures, respectively, revealed that while all N843_15 cells observed occurred in long filamented chains those of N843_10 retained typical L. monocytogenes cell morphology comprising single rods to short chains (Figures 2A,B). Consistent with filamentation in this strain, flow cytometric analysis also revealed increased area and chain length for N843_15 cells compared to N843_10 cells (Figures 2C,D). Electron microscopic examination further revealed improper cell septum formation and incomplete separation suggestive of cell division defects in N843_15 (Figures 3A–C). Cells of this isolate were highly irregular arranged in long chains of varying lengths without a well-defined cell structure, whereas the peptidoglycan layer was disordered with loosely attached fragments projecting from the surface giving it a solar flare appearance in contrast to a well-organized peptidoglycan observed for N843_10. Further physiological comparison of the two isolates using Biolog phenotypic microarray revealed reduced carbon source utilization range as well as increased osmotic and pH stress sensitivity in N843_15, the relapse PJI isolate compared to the initial isolate N843_10 (Table 2 and Figure 4). N843_15 had lost the ability to utilize 13 C-sources including relevant intracellular C-sources such as glycerol, maltose, and cellobiose (Table 2, Supplementary Table S2 and Supplementary Figure S1). Although overall more sensitive to various stress conditions, N843_15 was more resistant to alkaline stress (pH 9.5) than N843_10 in presence of β-Phenylethylamine (Table 2). Overall these analyses thus showed that N843_10 and N843_15 were phenotypic variants of the same L. monocytogenes strain, which suggested that N843_10, the initial PJI strain had evolved within the human hip PJI environment giving rise at some point to strain N843_15 isolated five years later during the relapse infection.
Figure 1. N843_15 exhibits slower growth, atypical aggregation behavior and reduced motility compared to N843_10. (A) Growth of N843_10 and N843_15 on blood agar incubated 24 h at 37°C. (B) N843_15 aggregated and sedimented upon growth in BHI broth. (C and D) Comparison of N843_15, N843_10 and N12-1273 growth kinetics in BHI broth at 37°C. (E) N843_10, N843_15 and N12-1273 swarming motility zones determined on semi-solid BHI agar at 25°C. Results presented are the means and standard deviations from three independent experiments. ∗ Indicates statistically significant difference between the strains (P < 0.05).
Figure 2. N843_10 and N843_15 cells display morphological differences. (A) N843_15 is filamented whereas N843_10 occurs in single rods and short chains under the light microscope. (B) Electron microscopy (EM) images showing that N843_15 filamentation phenotype is due to chains of unseparated cells whilst N843_10 showed single rods and short chains typical of L. monocytogenes morphology. (C and D) Flow cytometry analysis revealed significantly increased area and cell chain length in N843_15 compared to N843_10 and the clonal control strain N12-1273 (∗P < 0.05).
Figure 3. Electron microscopy (EM) reveals variations between N843_10 and N843_15 cell morphology. (A). N843_10 displays typical L. monocytogenes morphology of rod-shaped cells whereas the filamented N843_15 cells varied in shape (B) N843_10 cells displayed normal cell wall morphology and separation during cell division whilst N843_15 cells displayed altered cell wall morphology and lacked clear cell division septum demarcation. (C) N843_10 has a normal cell wall and is undergoing normal septum formation whereas N843_15 has altered cell walls with loosely attached fragments of peptidoglycan protruding from its surface giving it a solar flare appearance.
Figure 4. Overall growth or metabolic activity ring comparison of N843_15 and N843_10 with respect to C-source utilization (PM01 and PM02) and stress (osmolytes and pH; PM09 and PM10) resistance based on phenotypic microarray analysis. The gray inner circles indicate the strains’ order whilst the external circle indicates the PM categories. The metabolic activity referred to as activity index (AV) calculated for N843_15 under each assay condition per well is reported as color stripes going from red (AV = 0, no metabolism) to green (AV = 4, highest metabolic activity). Delta activity: the difference in the metabolic activity (AV) of N843_10 and N843_15 is reported when equal to or higher than 2 AV; gray is no difference; purple indicates a higher metabolic activity of N843_10 whilst orange color indicates that N843_10 has a lower metabolic activity than N843_15 under the assay conditions in that well.
Table 2. Comparison of carbon source utilization ability, pH and osmotic stress tolerance.
N843_15 the Recurrent PJI Isolate Displays Increased Antibiotic Sensitivity and Reduced Biofilm Production Than the Original Infection Parent Isolate N843_10
Assuming changes in either antibiotic susceptibility or biofilm production ability could have facilitated the long-term persistence and survival of the two isolates in the PJI contributing to emergence of the evolved variant N843_15, the antibiotic susceptibilities and biofilm production of the two PJI isolates were compared. Both strains were clinically sensitive to different therapeutic antibiotics associated with antibiotic therapy in the patient as well as gentamicin and sulfamethoxazole. The recurrent infection isolate N843_15 although more sensitive than the parent strain on most of the tested antibiotics was, however, more resistant to sulfamethoxazole (Table 3). N843_15 produced significantly less biofilm than the parent isolate N843_10 when assessed in TSB at 37°C (Figure 5A). The overall biofilm production of the parent isolate N843_10 was similar to other L. monocytogenes strains tested such as LL195 (CC1) and N12-1273 (CC412) but significantly less in comparison to a high biofilm producer reference strain N11-1850 (CC217) (Figure 5A). Our findings thus indicated that increased antibiotic resistance or enhanced biofilm production were unlikely to have played a role in the emergence of the evolved recurrent infection N843_15 isolate from the initial PJI isolate N843_10.
Figure 5. Biofilm production, hemolysis and Caco-2 cell invasion. (A) N843_10 produced more biofilm than N843_15. (B) L. monocytogenes N843_10 and N843_15 were poorly hemolytic compared to LL195 (CC1) and N12-1273 (CC412) strains (C) N843_15 showed lower cell invasion capacity than N843_10 but both strains were significantly less invasive than the reference strain LL195. Presented data shows the mean (bars) and standard deviation (error bars) of three independent biological experiments. ∗Indicates statistically significant differences where P < 0.05 based on one-way ANOVA and Tukey post-hoc test pairwise comparison of all the strains.
Table 3. MICs of different antibiotics determined for the L. monocytogenes N843_10 and N843_15 strains using E-tests.
PJI Isolates N843_10 and N843_15 Are Less Virulent Compared to Other L. monocytogenes Strains
A virulence comparison of the two PJI isolates showed that they had similar hemolytic activity levels on human red blood cells but the relapse infection isolate N843_15 was impaired in human Caco 2 cell invasion compared to the parent PJI isolate N843_10 (Figures 5B,C). On the other hand, both isolates showed significantly lower hemolytic and cell invasion when compared to other L. monocytogenes strains including another MLST clonal complex CC412 sporadic human listeriosis isolate N12-1273, as well as a CC1 previous listeriosis outbreak strain LL195 (Figures 5B,C). Comparing virulence using a zebrafish embryo-based infection model showed that both PJI isolates were avirulent whereas other strains used as positive controls such as L. monocytogenes EGDe, N12-1273, N2306 and LL195 were virulent causing mortality (Figure 6A). In a simulated localized zebrafish embryo infection model both PJI isolates similar to L. innocua the negative control were unable to induce disease signs over 48 h. Interestingly, while the injected N843_10 cells were cleared by the immune system those of the relapse isolate N843_15 remained uncleared nor did they decrease in quantity at the infection site (Figures 6B, 7). On the other hand, the CC412 clonal positive control strain N12-1273 was able to spread from the injection site and induce disease signs. Furthermore, the cells of this strain were not cleared by the immune system 48 hours post infection (hpi) but those of the L. innocua negative control had been cleared already at 24 hpi (Figure 7).
Figure 6. Strains vary in virulence and persistence within zebrafish: (A) In zebrafish embryos infected using the blood stream injection route, both N843_15 and N843_10 were unable to induce any clinical signs or mortality 72 h post infection. The reference strains used N2306 (CC4) and LL195 (CC1) had induced 100%, whilst EGDe (CC9) and N12-1273 (CC412) had induced 100% and 50% mortality of the zebrafish embryos within 24 and 72 h of infection, respectively. (B) GFP zebrafish embryo line coinfected with N843_15 mO2FP and N843_10 GFP. The parent strain N843_10 decreased in CFU numbers over time consistent with it being cleared while N843_15 levels remained unchanged over 48 h of infection. Results show the mean and standard deviation from 3 independent experiments. ∗Indicates significant differences in the level of the two strains at 24 and 48 hpi (P < 0.05 based on repeated measures two-way ANOVA).
Figure 7. Localized infection simulation using Danio rerio wik zebrafish Fli1: GFP. N843_15 was not cleared by the immune system neither did it cause diseases over 48 h of infection whereas N843_10 infection was cleared. N12-1273 used as a clonal positive control caused disease and was not cleared 48 hpi whilst the L innocua negative control used was already cleared at 24 hpi. Arrow denotes the infection site (otic vesicle) and the bacteria are in red (Magnification 20× confocal. Scale bar – 35 μm).
Genome-Wide Comparison Reveals Genetic Changes in the Evolved PJI Isolate L. monocytogenes N843_15
To assess if the different phenotypes between N843_10 and N843_15 were associated with genetic changes during L. monocytogenes N843 evolution within the human prosthetic joint environment we used whole genome sequences (WGS) analysis and compared the two N843 isolates. An in silico classical (seven loci) and core genome based MLST analysis both confirmed the two isolates to belong to MLST sequence type ST412 and clonal complex CC412. In addition, the N843_10 and N843_15 genomes clustered indistinguishably in a WGS-based phylogenetic tree and aligned colinearly without major chromosomal differences in a DNA sequence dot blot (Supplementary Figures S2A,B). Overall the two isolates showed only seven cgMLST allelic differences that are below the 10-allele difference recommended as the strain clonality cut off for L. monocytogenes. In contrast there were 21-24 and 1086-1655 cgMLST allelic profile differences when the two PJI isolates were compared to other epidemiologically unrelated ST412 strains (N12-1273, N18-2578 and N18-2708) and selected L. monocytogenes reference strains (EGDe, 10403S, ScottA and LL195) from other MLST STs, respectively (Figure 8, Supplementary Figure S3). Overall there were 26 genome wide single nucleotide polymorphisms (SNPs) detected in N843_15 compared to N843_10, which included 9 SNPs that cause non-synonymous amino acid changes in proteins with known functions (Table 4, Supplementary Table S3). While the other SNPs were located in intergenic regions and hypothetical proteins their potential functional consequences remain unknown. In addition, there were 44 insertion and deletion events (InDels) detected in the N843_15 genome including 21 InDels that are predicted to cause premature stop codons (PMSC) and truncation of proteins such as RNase J1 and HtrA (Table 5, Supplementary Table S3). A single base deletion was detected in the prfA gene of both N843_10 and N843_15 compared to L. monocytogenes EGDe prfA, which is predicted to cause a PMSC and an eleven amino acid PrfA carboxyl-terminal truncation (Supplementary Figure S4).
Figure 8. Minimum-spanning tree illustrating the phylogenetic relationship based on cgMLST allelic profiles of the 2 PJI isolates (N843_10 and N843_15) and seven selected reference strains including three CC412 isolates (N12-1273, N18-2578 and N18-2708). Each circle represents an allelic profile based on sequence analysis of 1,701 cgMLST target genes. The numbers on the connecting lines illustrate the numbers of target genes with differing alleles. The different groups of strains are distinguished by the colors of the circles. Clonal strains (cutoff 10 allele difference) are shaded in light brown.
Table 4. Single nucleotide polymorphism detected between N843_15 and N843_10 genomes and their predicted consequences1.
Table 5. InDels detected in the N843_15 genome compared to N843_101.
Assessing the Impact of the PrfA Truncation Mutation Predicted in N843_10 and N843_15
The phenotypic impact of the predicted PrfA truncation mutation in both PJI isolates was assessed through qRT-PCR analysis of prfA target genes showing that although containing more prfA mRNA, both N843_10 and N843_15 contained significantly low hly and plcA mRNA levels when compared to a clonal but epidemiologically unrelated positive control strain N12-1273 with an intact full length PrfA (Figure 9). The relapse PJI isolate N843_15 interestingly also contained higher prfA but lower hly and plcA mRNA compared to the parent strain isolate N843_10. Reduced plcA expression in N843_15 was also corroborated through qualitative PI-PLC activity analysis, which showed reduced activity in N843_15 than N843_10 (Supplementary Figure S5). PI-PLC activity of both PJI strains was, however, significantly lower than L. monocytogenes N12-1273 (CC412) and N2306 (CC4) (Supplementary Figures S5A,B). Overall these results were consistent with impaired PrfA activity in N843_10 and N843_15 strains confirming the WGS analysis predicted PrfA truncation mutation.
Figure 9. Impact of the PrfA truncation mutation predicted in N843_10 and N843_15 on virulence gene mRNA levels. Quantification of prfA, plcA, and hly mRNAs using qRT-PCR in the study strains that were cultured in BH broth at 37°C to the late exponential phase. Relative quantities (RQ) of prfA, plcA, and hly mRNA levels were normalized to 16S rRNA and are expressed relative to those of a L. monocytogenes EGDe based mRNA calibrator sample. Presented data shows the mean (bars) and standard deviation (error bars) of three independent biological experiments. ∗Indicates statistically significant differences where P < 0.05 based on one-way ANOVA and Tukey post-hoc test pairwise comparison of all the strains.
Phenotypic Characterization of RNase J1 Mutants in L. monocytogenes N843_10 and N2306 Strains
In order to assess phenotypic consequences of RNase J1 function loss due to a truncation mutation predicted in the relapse strain N843_15, rnjA gene deletion mutants in L. monocytogenes N843_10 (CC412) and N2306 (CC4) were created. These strains represented genetic backgrounds of the parent isolate (N843_10) and a different genetic lineage and MLST clonal complex (N2306), respectively. Recapitulating some of the phenotypic defects displayed in N843_15, both N843_10 and N2306 ΔrnjA mutants had reduced swarming motility, low plcA and hly mRNAs amounts and reduced PI-PLC activity, but showed increased cellular filamentation, antibiotic sensitivity and sulfamethoxazole resistance compared to their parental WT strains (Supplementary Figures S6–S8, Supplementary Table S4). Although the extent of filamentation in these ΔrnjA mutants was less pronounced than that of N843_15. In virulence evaluation using the simulated localized zebrafish embryo infection model, similar to N843_15, cells of N843_10 ΔrnjA were also not cleared by the immune system and did not cause detectable disease signs over 48 h of infection (Figure 10). As with the localized N843_15 infection macrophages were also initially attracted to the site of N843_10 ΔrnjA injection but returned to a normal distribution as was also observed in PBS injected negative controls. In the case of N2306, which unlike the N843_10 isolate has an intact prfA gene, both the WT and ΔrnjA strains were not cleared by the immune system and were able to spread from the infection site, with the WT strain showing a higher spread intensity (Figure 10). The WT strain N2306 caused 100% zebrafish mortality at 24 hpi whereas N2306 ΔrnjA had only caused 10% mortality, and only reached 100% at 72 hpi (Supplementary Figure S9).
Figure 10. In the zebrafish mpeg1: GFP line infection with N843_15 and N843_10 ΔrnjA, macrophages were initially attracted to the injection site (otic vesicle) but returned to the normal distribution similar to PBS injected fish within 24 h without clearing the bacterial from the site of infection. Macrophages are green and the bacteria are in red (Magnification 20x confocal. Scale bar – 35 μm).
Discusion
We examined in this study two L. monocytogenes isolates that were recovered within five years of each other from two episodes of a recurrent human PJI. Genomic analyses showed that the two L. monocytogenes isolates; N843_10 and N843_15 isolated in 2010 and 2015, respectively, were the same strain. Although displaying various minor genetic changes the two isolates clustered indistinguishably in a WGS-based phylogenetic tree and had less than 10 cgMLST allelic differences, which is currently considered a cut off for clonality among listerial strains (Ruppitsch et al., 2015). These two isolates displayed distinctive phenotypic traits suggesting that the relapse infection isolate N843_15 was an evolved variant of the original infection strain N843_10, which had evolved during the 5 years of residence within the human hip PJI environment. Besides impaired growth behavior and cellular morphology, the relapse infection isolate N843_15 also showed reduced motility, virulence and stress resistance as well as decreased carbon source utilization range, low biofilm production and altered antibiotic susceptibility compared to the N843_10 isolate. Comparing the genomes between the two strains revealed several genetic changes in the relapse infection isolate N843_15 including mutations in specific genes that might be linked to some of the phenotypic variations exhibited between this isolate and the original infection parent isolate N843_10. Presumably, these genetic changes arose from random mutation events that might be related to the evolution and adaptation of L. monocytogenes N843 strain to the niche of the infected human prosthetic joint environment. Others have also previously described phenotypically distinguishable clonally identical L. monocytogenes strains differing in phenotypic traits such as hemolysis, but these were isolated simultaneously from the same septic joint fluid (Charlier et al., 2012).
One of the genetic mutations in the evolved variant N843_15 was a single base deletion (106A) in the rnjA gene causing truncation of RNase J1(Table 5). RNase J1 is an RNA nuclease that plays an important role in bacterial RNA processing and degradation (Even et al., 2005; de la Sierra-Gallay et al., 2008; Figaro et al., 2013). The altered cellular morphology phenotypic defects observed in N843_15 such as the disordered peptidoglycan and cell chaining mirrors those previously described in Bacillus subtilis RNase J1 knock-down cells (Hunt et al., 2006; Mader et al., 2008; Durand et al., 2012; Figaro et al., 2013). Furthermore, we could recapitulate some phenotypic defects exhibited by N843_15 such as cellular filamentation as well as reductions in motility, antibiotic resistance and virulence through rnjA gene deletion in L. monocytogenes N843_10 and N2306 strains, which represent genetic backgrounds of the original PJI isolate as well as a different L. monocytogenes genetic lineage, respectively. The N843_10 and N2306 ΔrnjA mutant phenotypes were however less pronounced than those of the N843_15 isolate, which was presumably because the N843_15 phenotypic defects are probably exacerbated due to contributions of other mutations besides the RNase J1 truncation (Table 5). Possible contributing mutations include a single base insertion (769A) induced truncation of the D-alanine carboxypeptidase (PBDB2) protein Lmo2812 previously found to be important in peptidoglycan synthesis and cell wall turnover events in L. monocytogenes (Principe et al., 2009).
Reduction of flagella-based swarming motility observed in N843_15 could be due to altered cellular morphology since incompletely separated chained cells formed by this strain might be too large to be moved through flagella. Besides that, the peptidoglycan structural disordering observed in this isolate could also impair flagella assembly. Although it remains to be investigated, the flagellum specific muraminidase FlgJ in N843_15 also has a SNP induced V15C amino acid change, which might have possible functional consequences that could compromise flagella assembly and motility processes. As previously suggested by others, L. monocytogenes might also strategically down regulate the expression of motility related genes as a way of avoiding the activation of the host immune system, reduced motility in N843_15 might thus also be a deliberate adaptive change to long-term residence within the human host PJI environment (Toledo-Arana et al., 2009).
There was loss in the ability to utilize some carbon sources and increased sensitivity to both acidic and osmotic stress detected in the relapse isolate N843_15 compared to the parent isolate N843_10. Notably such losses included the ability to utilize intracellular and food relevant C-sources such as glycerol and pectin, respectively, as well as increased sensitivity to sodium chloride, sodium lactate, and sodium benzoate stresses. It is tempting to speculate that sacrificing of such phenotypes might be related to the evolution of the L. monocytogenes N843 strain as it adapted to the human host PJ capsule environment. At genome level, the changes in C-source utilization ability in N843_15 might in part be associated with amino acid changing and truncation inducing mutations of carbon metabolism related proteins such as the maltose phosphorylase (loss of D-maltose utilization ability) and galactitol operon BglG anti-terminator, BvrA. Meanwhile increased osmotic and acid stress sensitivities could have arisen from the peptidoglycan structural alteration, as well as ClpQ, HtrA and ArcD truncation mutations found in the N843_15 genome. Proper peptidoglycan architecture is important for bacterial survival in diverse environments and as previously shown in Lactobacillus and Ochrobactrum species its alteration increases osmotic and acid stress sensitivity (Piuri et al., 2005; Principe et al., 2009; Alonzo et al., 2011; Bergholz et al., 2012; Vollmer, 2015). Proteases HtrA and ClpQ promote general stress responses through degradation of misfolded proteins, whereas ArcD promotes acid stress tolerance via the ADI system by facilitating the exchange of intracellular ornithine for extracellular arginine in L. monocytogenes (Stack et al., 2005; Bowman et al., 2008; Kocaman and Sarımehmetoǧlu, 2016). In addition, proteins with DD carboxypeptidase activity are also involved in cell wall turnover as they function to cleave peptidoglycan cross-links (Piuri et al., 2005; Korsak et al., 2010). Similar with previous observations for DD-carboxypeptidase mutants in Ochrobactrum species and L. monocytogenes, the N843_15 isolate that also bears a D-alanine carboxypeptidase (PBDB2-lmo2812) truncation mutation also displayed increased sensitivity to osmotic stress compared to the parent isolate N843_10 (Principe et al., 2009; Bergholz et al., 2012). Although not yet proven, the overall reduction in ability to utilize intracellular carbon sources, as well as low acid stress tolerance and hemolytic ability might be reflective of N843_15 adaptation to an extracellular life within the human prosthetic joint environment.
N843_10 and N843_15 isolates were both clinically sensitive based on Clinical and Laboratory Standards Institute and European Committee on Antimicrobial Susceptibility Testing standards against antibiotics tested including those associated with the PJI therapy (CLSI, 2010; European Committee on Antimicrobial Susceptibility Testing, 2017). This observation rules out the possibility that the evolved mutant isolate N843_15 had emerged through increased antibiotic resistance and subsequent selection during PJI treatment. Both strains could, however, have survived due to antibiotic therapy failure since bactericidal antibiotic levels might not have been achieved in the PJ environment during therapy. On the other hand, N843_15 was more sensitive to most of the tested antibiotics than the parent isolate N843_10 probably due to increased antibiotic permeability due to its peptidoglycan structural alteration or as a result of its various specific acquired mutations. A 50S ribosomal protein L15 subunit truncation mutation found in this isolate might explain the increased sensitivity to the 50S targeting antibiotics tetracycline and azithromycin. The HtrA and RNase J1 truncation mutations could increase sensitivity to penicillin and other antibiotics in N843_15 since previous studies showed loss of these proteins increases antibiotic sensitivity in L. monocytogenes and B. subtilis, respectively (Stack et al., 2005; Figaro et al., 2013). The contribution of RNase J1 loss to increased antibiotic sensitivity was also corroborated since both N843_10 and N2306 ΔrnjA mutants also showed increased antibiotic sensitivity. Similar with observations upon RNase J1 depletion in B. subtilis, L. monocytogenes N843_15 was also more resistant to the dihydropteroate synthetase inhibitor sulfamethoxazole than the parent N843_10 isolate (Hunt et al., 2006). Increased sulfamethoxazole resistance was also observed for N843_10 ΔrnjA but not the N2306 ΔrnjA mutant, which indicates some genetic lineage and strain-specific differences regarding the contribution of RNase J1 loss to this phenotype in L. monocytogenes. Similarly, strain specific variation in sensitivity phenotypes has also been described among RNase J1 mutants created in different B. subtilis strains against the dihydrofolate reductase inhibitor trimethoprim, which like sulfamethoxazole also targets a stage in the folate pathway (Hunt et al., 2006; Figaro et al., 2013).
Prosthetic joint infections pathogenesis varies from that of a native joint since it relies on biofilm formation and L. monocytogenes has a high affinity for implants (Francolini and Donelli, 2010; Bader et al., 2016). Biofilm environments protect the bacterium from exposure to antibiotics and possibly the immune system (Kleemann et al., 2009). Although it remains possible that the relapse isolate N843_15 could also have persisted as a biofilm in the PJI our analysis here showed that it produces significantly less biofilm than the parent isolate N843_10. A possible explanation for this could be the HtrA truncation mutation in this isolate since reduced biofilm production has been described in htrA null mutants of L. monocytogenes (Wilson et al., 2006). We, however, deem it unlikely that N843_15 survived overtime in the PJI due to an increased biofilm production ability compared to the original N843_10 isolate, although it is important to note that the biofilm comparison was conducted under laboratory conditions using attachment surfaces that differ from those encountered within the human prosthetic joint.
Virulence analysis showed that the relapse infection isolate N843_15 has reduced cell invasion capacity than the original infection isolate N843_10, which might probably be associated with altered display and anchoring of cell surface associated virulence factors including internalins due to the filamentous morphology and disordered peptidoglycan in this strain. L. monocytogenes virulence depends on a variety of virulence factors that are transcriptionally regulated through PrfA (Chatterjee et al., 2006; Desvaux and Hebraud, 2006; Cossart and Lebreton, 2014; Radoshevich and Cossart, 2018). One hypothesis is that both N843 isolates might have persisted in the prosthetic joint due to low virulence that limited deleterious effects to the host. In support of this notion both N843 isolates displayed low virulence than other reference L. monocytogenes strains tested based on hemolysis, cell invasion and zebrafish pathogenicity assessments. The low virulence in the N843 isolates was associated with a prfA truncation mutation similar to previous observations in other L. monocytogenes strains (Maury et al., 2017). In both N843 PJI isolates this was due to a prfA deletion mutation that truncates the PrfA protein within the dimerization domain. Site directed mutagenesis of amino acid residues within this PrfA region was previously shown to impair DNA binding and virulence gene expression causing low virulence and LLO secretion (Eiting et al., 2005; Scortti et al., 2007; de las Heras et al., 2011; Good et al., 2016). Meanwhile the slightly elevated prfA mRNA levels detected in N843_10 and N843_15 compared to a control CC412 strain N12_1273 might be indicative of increased production of truncated PrfA to compensate for reduced activity. Such increased expression might be driven through other transcriptional regulators such as alternative sigma factor Sig B (Lebreton and Cossart, 2017). Interestingly, prfA mRNA levels were also higher in N843_15 compared to N843_10. We, however, do not have an explanation for this, but it might be linked to other mutations that N843_15 has in addition to the PrfA truncation.
In a simulation of localized infection using different zebrafish lines, isolate N843_15 and the ΔrnjA mutant of N843_10 were not cleared by the immune system over 48 hrs and neither did they cause any overt disease signs. Notably although macrophages were initially attracted to the bacterial injection site, they returned to normal distribution without clearing cells of these strains from the infection site. Such a response pattern was also observed for the PBS injected negative controls. Although mechanisms behind this phenotype are not yet clear, one possible explanation might the altered expression of motility and virulence associated genes. Furthermore, the detected mutations might have an impact on the expression and secretion of Pathogen-associated molecular patterns (PAMPs) such as flagellin and LLO, thereby inducing a different response from the immune system compared to other normal wild type L. monocytogenes strains (Rose et al., 2001; Way et al., 2004; Schuppler and Loessner, 2010; Wallecha et al., 2013).
A genome comparison between N843_10 and N843_15 revealed various genetic variations some of which might explain the phenotypic differences detected between the two N843 isolates. It is, however, important to note that going forward specific gene targeting mutagenesis and complementation approaches will be necessary to validate some of the observations reported here. Meanwhile, the N843_15 genome changes detected relative to the N843_10 parent strain demonstrates vestigiality suggestive of short-term adaptation through inactivation and retainment of the unnecessary genes. In long-term adaptation, an eventual complete loss of the unnecessary genes would be expected. In contrast to findings by others, our study shows relatively higher level of genetic variability between the two strains which differs from high level genome conservation observed from long-term evolution of L. monocytogenes in other settings such as food processing environments (Orsi et al., 2008; Lomonaco et al., 2015; Moura et al., 2016; Harranda et al., 2020). A likely contributor to the relatively high mutation rate could be that the strains were replicating at a relatively higher rate within the PJI environment at 37°C, than they would in food production plant environments upon which previous observations are based (Orsi et al., 2008; Lomonaco et al., 2015; Harranda et al., 2020). We are, however, aware that such a comparison of mutation rates might be confounded due to differences in strain population sizes and genetic backgrounds as well as environmental niches involved between these other studies and our current study that only involved two isolates within the PJI environment.
In conclusion, we have characterized two clonally identical L. monocytogenes strains N843_10 and N843_15 that were isolated 5 years apart from a persistent rare manifestation of listeriosis. We show that the long-term residence of this bacterium within a human host PJI environment led to a wide range of phenotypic changes in virulence, metabolic flexibility, and stress resilience that were associated with various genetic changes during short-term evolution of L. monocytogenes within a human PJI environment. Overall our observations besides highlighting the phenotypic and genotypic variations between these two strains might have provided some insights into molecular mechanisms associated with adaptation and evolution of L. monocytogenes within the environment of a human host PJI. It is, however, important to note that our case might be a special case involving strains that carry a natural truncation mutation in the main virulence regulator PrfA. As such, more case studies would be required in future in order to generalize our findings regarding mechanism of adaptation of other L. monocytogenes strains during persistence within the human host PJI environment.
Data Availability Statement
The datasets generated for this study can be found in the NCBI GeneBank – Accession numbers CP046361 and CP046362.
Ethics Statement
This study was performed in accordance with the principles and recommendations of the “Ordinance on laboratory animal husbandry, the production of genetically modified animals and the methods of animal experimentation; Animal Experimentation Ordinance” (SR 455.163, April 12, 2010), Swiss Federal Food Safety and Veterinary Office (FSVO/BLV). The maximum age reached by the embryos during experimentation was 5 days post fertilization (dpf) for which no license is required from the cantonal veterinary office in Switzerland, since such embryos will not have yet reached the free-feeding stage. Husbandry and breeding of the adult zebrafish were performed under the supervision of Prof. Stephan Neuhauss, Institute for Molecular Life Sciences, University of Zurich, Zurich, Switzerland. All animal protocols used were in compliance with internationally recognized standards as well as with Swiss legal ethical guidelines for the use of fish in biomedical research. All the experiments were approved by the local authorities (Veterinäramt Zürich Tierhaltungsnummer 150).
Author Contributions
TT and FM designed the study. TT supervised the study. FM and AE performed the experiments. UA assisted in the PM experiments. FM, AE, UA, MS, and TT analyzed the data and wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This project was funded by the University of Zurich and done in collaboration with Agroscope, Bern, Switzerland.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank the NENT for providing the strains and their metadata, the lab of Martin Loessner (Laboratory for Food Microbiology, ETH Zürich) for providing the pIMK-mO2FP plasmid and the Listeria innocua mO2FP strain, the TEM Unit of the Institute for Veterinary Pathology, Zurich, (IVPZ) and the Flow cytometry center both of the University of Zurich for availing their facilities and all their technical assistance with electron microscopy and flow cytometry, respectively. We would also thank Anne Guisolan and Dieter Weik (both of Agroscope, Bern, Switzerland) for their assistance with the PM experiments.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2020.01726/full#supplementary-material
Footnotes
- ^ http://bigsdb.web.pasteur.fr/listeria/listeria.html
- ^ https://biolog.com
- ^ https://rast.nmpdr.org/
- ^ https://bigsdb.pasteur.fr/listeria/
- ^ https://www.cgmlst.org/ncs/schema/690488/
References
Allerberger, F., and Wagner, M. (2010). Listeriosis: a resurgent foodborne infection. Clin. Microbiol. Infect. 16, 16–23. doi: 10.1111/j.1469-0691.2009.03109.x
Alonzo, F., McMullen, D. P., and Freitag, N. E. (2011). Actin polymerization drives septation of Listeria monocytogenes namA hydrolase mutants, demonstrating host correction of a bacterial defect. Infect. Immun. 79, 1458–1470. doi: 10.1128/iai.01140-10
Althaus, D., Lehner, A., Brisse, S., Maury, M., Tasara, T., and Stephan, R. (2014). Characterization of Listeria monocytogenes strains isolated during 2011-2013 from human infections in Switzerland. Foodborne Pathog. Dis. 11, 753–758. doi: 10.1089/fpd.2014.1747
Bader, G., Al-Tarawneh, M., and Myers, J. (2016). Review of prosthetic joint infection from Listeria monocytogenes. Surg. Infect. 17, 739–744. doi: 10.1089/sur.2016.067
Baranyi, J., and Roberts, T. A. (1994). A dynamic approach to predicting bacterial-growth in food. Int. J. Food Microbiol. 23, 277–294. doi: 10.1016/0168-1605(94)90157-0
Bergholz, T. M., Bowen, B., Wiedmann, M., and Boor, K. J. (2012). Listeria monocytogenes shows temperature-dependent and -independent responses to salt stress, including responses that induce cross-protection against other stresses. Appl. Environ. Microbiol. 78, 2602–2612. doi: 10.1128/aem.07658-11
Bille, J. (1990). “Epidemiology of human listeriosis in Europe with special reference to the Swiss outbreak,” in Foodborne Listeriosis, eds A. J. Miller, J. L. Smith, and G. A. Somkuti, (New York, NY: Elsevier), 71–74.
Bochner, B. R. (2009). Global phenotypic characterization of bacteria. FEMS Microbiol. Rev. 33, 191–205. doi: 10.1111/j.1574-6976.2008.00149.x
Bowman, J. P., Bittencourt, C. R., and Ross, T. (2008). Differential gene expression of Listeria monocytogenes during high hydrostatic pressure processing. Microbiology 154, 462–475. doi: 10.1099/mic.0.2007/010314-0
Bush, L. M., Alrifai, A., and Perez, M. T. (2015). Listeria monocytogenes prosthetic joint infections a review a propos a case report. Infect. Dis. Clin. Pract. 23, 66–69. doi: 10.1097/IPC.0000000000000218
Camilli, A., Tilney, L. G., and Portnoy, D. A. (1993). Dual roles of plcA in Listeria monocytogenes pathogenesis. Mol. Microbiol. 8, 143–157. doi: 10.1111/j.1365-2958.1993.tb01211.x
Centers for Disease Control [CDC] (2018). Centre for disease control and prevention. Preliminary incidence and trends of infections with pathogens transmitted commonly through food. Foodborne Diseases Active Surveillance Network, 10 U.S. Sites, 2006-2017. MMWR Morb. Mortal. Wkly. Rep. 67, 324–328. doi: 10.15585/mmwr.mm6711a3
Charlier, C., Leclercq, A., Cazenave, B., Desplaces, N., Travier, L., Cantinelli, T., et al. (2012). Listeria monocytogenes – associated joint and bone infections: a study of 43 consecutive cases. Clin. Infect. Dis. 54, 240–248. doi: 10.1093/cid/cir803
Chatterjee, S. S., Hossain, H., Otten, S., Kuenne, C., Kuchmina, K., Machata, S., et al. (2006). Intracellular gene expression profile of Listeria monocytogenes. Infect. Immun. 74, 1323–1338. doi: 10.1128/IAI.74.2.1323-1338.2006
Chavada, R., Keighley, C., Quadri, S., Asghari, R., Hofmeyr, A., and Foo, H. (2014). Uncommon manifestations of Listeria monocytogenes infection. BMC Infect. Dis. 14:641. doi: 10.1186/s12879-014-0641-x
Chougle, A., and Narayanaswamy, V. (2004). Delayed presentation of prosthetic joint infection due to Listeria monocytogenes. Int. J. Clin. Pract. 58, 420–421. doi: 10.1111/j.1368-5031.2004.0040.x
CLSI, (2010). Methods for Antimicrobial Dilution and Disk Susceptibility Testing for Infrequently Isolated or Fastidious Bacteria: Approved Guideline, 2nd Edn, Vol. 30. Wayne, PA: CLSI.
Cone, L. A., Fitzmorris, A. O., and Hirschberg, J. M. (2001). Is Listeria monocytogenes an important pathogen for prosthetic joints? J. Clin. Rheumatol. 7, 34–37. doi: 10.1097/00124743-200102000-00008
Cossart, P., and Lebreton, A. (2014). A trip in the “new microbiology” with the bacterial pathogen Listeria monocytogenes. FEBS Lett. 588, 2437–2445. doi: 10.1016/j.febslet.2014.05.051
Darling, A. E., Mau, B., and Perna, N. T. (2010). ProgressiveMauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One 5:e11147. doi: 10.1371/journal.pone.0011147
de la Sierra-Gallay, I., Zig, L., Jamalli, A., and Putzer, H. (2008). Structural insights into the dual activity of RNase. J. Nat. Struct. Mol. Biol. 15, 206–212. doi: 10.1038/nsmb.1376
de las Heras, A., Cain, R. J., Bielecka, M. K., and Vázquez-Boland, J. A. (2011). Regulation of Listeria virulence: PrfA master and commander. Curr. Opin. Microbiol. 14, 118–127. doi: 10.1016/j.mib.2011.01.005
Del Pozo, J., Garcia de la Garza, R., Diaz de Rada, P., Enrique, O., and Jose, R. Y. (2013). Listeria monocytogenes septic arthritis in a patient treated with mycophenolate mofetil for polyarteritis nodosa: a case report and review of the literature. Int. J. Infect. Dis. 17, 132–133. doi: 10.1016/j.ijid.2012.11.005
Desvaux, M., and Hebraud, M. (2006). The protein secretion systems in Listeria: inside out bacterial virulence. FEMS Microbiol. Rev. 30, 774–805. doi: 10.1111/j.1574-6976.2006.00035.x
Doumith, M., Buchrieser, C., Glaser, P., Jacquet, C., and Martin, P. (2004). Differentiation of the major Listeria monocytogenes serovars by multiplex PCR. J. Clin. Microbiol. 42, 3819–3822. doi: 10.1128/JCM.42.8.3819-3822.2004
Durand, S., Gilet, L., Bessières, P., Nicolas, P., and Condon, C. (2012). Three essential ribonucleases—RNase Y, J1, and III—control the abundance of a majority of Bacillus subtilis mRNAs. PLoS Genet. 8:e1002520. doi: 10.1371/journal.pgen.1002520
Ebner, R., Stephan, R., Althaus, D., Brisse, S., Maury, M., and Tasara, T. (2015). Phenotypic and genotypic characteristics of Listeria monocytogenes strains isolated during 2011-2014 from different food matrices in Switzerland. Food Control 57, 321–326. doi: 10.1016/j.foodcont.2015.04.030
Eiting, M., Hageluken, G., Schubert, W. D., and Heinz, D. W. (2005). The mutation G145S in PrfA, a key virulence regulator of Listeria monocytogenes, increases DNA-binding affinity by stabilizing the HTH motif. Mol. Microbiol. 56, 433–446. doi: 10.1111/j.1365-2958.2005.04561.x
Eshwar, A. K., Guldimann, C., Oevermann, A., and Tasara, T. (2017). Cold-shock domain family proteins (Csps) are involved in regulation of virulence, cellular aggregation, and flagella-based motility in Listeria monocytogenes. Front. Cell. Infect. Microbiol. 7:453. doi: 10.3389/fcimb.2017.00453
European Committee on Antimicrobial Susceptibility Testing (2017). Breakpoint Tables for Interpretation of MICs and Zone Diameters. Version 7.0.
European Food Safety Authority [EFSA] (2017). The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2016. EFSA J. 15:5077. doi: 10.2903/j.efsa.2017.5077
Even, S., Pellegrini, O., Zig, L., Labas, V., Vinh, J., Bréchemmier-Baey, D., et al. (2005). Ribonucleases J1 and J2: two novel endoribonucleases in B. subtilis with functional homology to E. coli RNase E. Nucleic Acids Res. 33, 2141–2152. doi: 10.1093/nar/gki505
Figaro, S., Durand, S., Gilet, L., Cayet, N., Sachse, M., and Condon, C. (2013). Bacillus subtilis mutants with knockouts of the genes encoding Ribonucleases RNase Y and RNase J1 are viable, with major defects in cell morphology, sporulation, and competence. J. Bacteriol. 195, 2340–2348. doi: 10.1128/jb.00164-13
Francolini, I., and Donelli, G. (2010). Prevention and control of biofilm-based medical-device-related infections. FEMS Immunol. Med. Microbiol. 59, 227–238. doi: 10.1111/j.1574-695X.2010.00665.x
Galardini, M., Mengoni, A., Biondi, E. G., Semeraro, R., Florio, A., Bazzicalupo, M., et al. (2014). DuctApe: a suite for the analysis and correlation of genomic and OmniLog phenotype microarray data. Genomics 103, 1–10. doi: 10.1016/j.ygeno.2013.11.005
Glaser, P., Frangeul, L., Buchrieser, C., Rusniok, C., Amend, A., Baquero, F., et al. (2001). Comparative genomics of Listeria species. Science 294, 849–852. doi: 10.1126/science.1063447
Göker, M., Hofner, B., Montero Calasanz, M. D. C., Sikorski, J., and Vaas, L. A. I. (2016). opm: An R Package for Analysing Phenotype Microarray and Growth Curve Data. Phenotype Microarray Data. Available at: http://www.goeker.org/opm/opm_doc/doc/opm-tutorial.pdf (accessed October 11, 2019).
Good, J. A., Andersson, C., Hansen, S., Wall, J., Krishnan, K. S., Begum, A., et al. (2016). Attenuating Listeria monocytogenes virulence by targeting the regulatory protein PrfA. Cell Chem. Biol. 23, 404–414. doi: 10.1016/j.chembiol.2016.02.013
Guldimann, C., Bärtschi, M., Frey, J., Zurbriggen, A., Seuberlich, T., and Oevermann, A. (2015). Increased spread and replication efficiency of Listeria monocytogenes in organotypic brain-slices is related to multilocus variable number of tandem repeat analysis (MLVA) complex. BMC Microbiol. 15:134. doi: 10.1186/s12866-015-0454-0
Harranda, A. S., Jagadeesanb, B., Baertb, L., Wiedmanna, M., and Orsia, R. H. (2020). Evolution of Listeria monocytogenes in a food-processing plant involves limited single nucleotide substitutions, but considerable diversification by gain and loss of prophages. Appl. Environ. Microbiol. 86:e02493-19.
Hunt, A., Rawlins, J. P., Thomaides, H. B., and Errington, J. (2006). Functional analysis of 11 putative essential genes in Bacillus subtilis. Microbiology 152, 2895–2907. doi: 10.1099/mic.0.29152-0
Kesteman, T., Yombi, J. C., Gigi, J., and Durez, P. (2007). Listeria infections associated with infliximab: case reports. Clin. Rheumatol. 26, 2173–2175. doi: 10.1007/s10067-007-0660-8
Kleemann, P., Domann, E., Chakraborty, T., Bernstein, I., and Lohoff, M. (2009). Chronic prosthetic joint infection caused by Listeria monocytogenes. J. Med. Microbiol. 58, 138–141. doi: 10.1099/jmm.0.004234-0
Kocaman, N., and Sarımehmetoǧlu, B. (2016). Stress responses of Listeria monocytogenes. Ankara Üniv. Vet. Fak. Derg. 63, 421–427. doi: 10.1501/vetfak_0000002764
Korsak, D., Markiewicz, Z., Gutkind, G. O., and Ayala, J. A. (2010). Identification of the full set of Listeria monocytogenes penicillin-binding proteins and characterization of PBPD2 (Lmo2812). BMC Microbiol. 10:239. doi: 10.1186/1471-2180-10-239
Krumsiek, J., Arnold, R., and Rattei, T. (2007). Gepard: a rapid and sensitive tool for creating dotplots on genome scale. Bioinformatics 23, 1026–1028. doi: 10.1093/bioinformatics/btm039
Lauer, P., Chow, M. Y., Loessner, M. J., Portnoy, D. A., and Calendar, R. (2002). Construction, characterization, and use of two Listeria monocytogenes site-specific phage integration vectors. J. Bacteriol. 184, 4177–4186. doi: 10.1128/jb.184.15.4177-4186.2002
Lebreton, A., and Cossart, P. (2017). RNA- and protein-mediated control of Listeria monocytogenes virulence gene expression. RNA Biol. 14, 460–470. doi: 10.1080/15476286.2016.1189069
Lenguerrand, E., Whitehouse, R. M., Beswick, A. D., Kunutsor, S. K., Burston, B., Porter, M., et al. (2018). Risk factors associated with revision for prosthetic joint infection after hip replacement: a prospective observational cohort study. Lancet Infect. Dis. 18, 1004–1014. doi: 10.1016/s1473-3099(18)30345-1
Lomonaco, S., Nucera, D., and Filipello, V. (2015). The evolution and epidemiology of Listeria monocytogenes in Europe and the United States Infect. Genet. Evol. 35, 172–183. doi: 10.1016/j.meegid.2015.08.008
Mader, U., Zig, L., Kretschmer, J., Homuth, G., and Putzer, H. (2008). mRNA processing by RNases J1 and J2 affects Bacillus subtilis gene expression on a global scale. Mol. Microbiol. 70, 183–196. doi: 10.1111/j.1365-2958.2008.06400.x
Maury, M. M., Chenal-Francisque, V., Bracq-Dieye, H., Han, L., Leclercq, A., Vales, G., et al. (2017). Spontaneous loss of virulence in natural populations of Listeria monocytogenes. Infect. Immun. 85:e00541-17.
Metsalu, T., and Vilo, J. (2015). ClustVis: a web tool for visualizing clustering of multivariate data using principal component analysis and heatmap. Nucleic Acids Res. 43, W566–W570. doi: 10.1093/nar/gkv468
Monk, I. R., Gahan, C. G. M., and Hill, C. (2008). Tools for functional postgenomic analysis of Listeria monocytogenes. Appl. Environ. Microbiol. 74, 3921–3934. doi: 10.1128/aem.00314-08
Moura, A., Criscuolo, A., Pouseele, H., Maury, M. M., Leclercq, A., Tarr, C., et al. (2016). Whole genome- based population biology and epidemiological surveillance of Listeria monocytogenes. Nat. Microbiol. 2:16185. doi: 10.1038/nmicrobiol.2016.185
Muchaamba, F., Eshwar, A. K., Stevens, M. J. A., von Ah, U., and Tasara, T. (2019). Variable carbon source utilization, stress resistance, and virulence profiles among Listeria monocytogenes strains responsible for Listeriosis outbreaks in Switzerland. Front. Microbiol. 10:957. doi: 10.3389/fmicb.2019.00957
Orsi, R. H., Borowsky, M. L., Lauer, P., Young, S. K., Nusbaum, C., Galagan, J. E., et al. (2008). Short-term genome evolution of Listeria monocytogenes in a non-controlled environment. BMC Genomics 9:539. doi: 10.1186/1471-2164-9-539
Piuri, M., Sanchez-Rivas, C., and Ruzal, S. M. (2005). Cell wall modifications during osmotic stress in Lactobacillus casei. J. Appl. Microbiol. 98, 84–95. doi: 10.1111/j.1365-2672.2004.02428.x
Principe, A., Jofre, E., Alvarez, F., and Mori, G. (2009). Role of a serine-type D-alanyl-D-alanine carboxypeptidase on the survival of Ochrobactrum sp. 11a under ionic and hyperosmotic stress. FEMS Microbiol. Lett. 295, 261–273. doi: 10.1111/j.1574-6968.2009.01604.x
Radoshevich, L., and Cossart, P. (2018). Listeria monocytogenes: towards a complete picture of its physiology and pathogenesis. Nat. Rev. Microbiol. 16, 32–46. doi: 10.1038/nrmicro.2017.126
Rose, F., Zeller, S. A., Chakraborty, T., Domann, E., Machleidt, T., Kronke, M., et al. (2001). Human endothelial cell activation and mediator release in response to Listeria monocytogenes virulence factors. Infect. Immun. 69, 897–905. doi: 10.1128/IAI.69.2.897-905.2001
Ruppitsch, W., Pietzka, A., Prior, K., Bletz, S., Fernandez, H. L., Allerberger, F., et al. (2015). Defining and evaluating a core genome multilocus sequence typing scheme for whole-genome sequence-based typing of Listeria monocytogenes. J. Clin. Microbiol. 53, 2869–2876. doi: 10.1128/jcm.01193-15
Schmid, B., Klumpp, J., Raimann, E., Loessner, M. J., Stephan, R., and Tasara, T. (2009). Role of cold shock proteins in growth of Listeria monocytogenes under cold and osmotic stress conditions. Appl. Environ. Microbiol. 75, 1621–1627. doi: 10.1128/aem.02154-08
Schuppler, M., and Loessner, M. J. (2010). The opportunistic pathogen Listeria monocytogenes: pathogenicity and interaction with the mucosal immune system. Int. J. Inflamm. 2010:704321.
Scortti, M., Monzo, H. J., Lacharme-Lora, L., Lewis, D. A., and Vazquez-Boland, J. A. (2007). The PrfA virulence regulon. Microbes Infect. 9, 1196–1207. doi: 10.1016/j.micinf.2007.05.007
Shen, A., and Higgins, D. E. (2005). The 5’ untranslated region-mediated enhancement of intracellular listeriolysin O production is required for Listeria monocytogenes pathogenicity. Mol. Microbiol. 57, 1460–1473. doi: 10.1111/j.1365-2958.2005.04780.x
Smith, K., and Youngman, P. (1992). Use of a new integrational vector to investigate compartment-specific expression of the Bacillus subtilis spoIIM gene. Biochimie 74, 705–711. doi: 10.1016/0300-9084(92)90143-3
Stack, H. M., Sleator, R. D., Bowers, M., Hill, C., and Gahan, C. G. M. (2005). Role for HtrA in stress induction and virulence potential in Listeria monocytogenes. Appl. Environ. Microbiol. 71, 4241–4247. doi: 10.1128/AEM.71.8.4241-4247.2005
Stephan, R., Althaus, D., Kiefer, S., Lehner, A., Hatz, C., Schmutz, C., et al. (2015). Foodborne transmission of Listeria monocytogenes via ready-to-eat salad: a nationwide outbreak in Switzerland, 2013–2014. Food Control 57, 14–17. doi: 10.1016/j.foodcont.2015.03.034
Toledo-Arana, A., Dussurget, O., Nikitas, G., Nina Sesto, N., Guet-Revillet, H., Balestrino, D., et al. (2009). The Listeria transcriptional landscape from saprophytism to virulence. Nature 459, 950–956. doi: 10.1038/nature08080
Treangen, T. J., Ondov, B. D., Koren, S., and Phillippy, A. M. (2014). The harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 15, 1–15. doi: 10.1186/s13059-014-0524-x
Vollmer, W. (2015). “Peptidoglycan,” in Molecular Medical Microbiology, 2nd Edn, eds Y. Tang, M. Sussman, D. Liu, I. Poxton, and J. Schwartzman, (Boston, MA: Academic Press), 105–124.
Wallecha, A., Wood, L., Pan, Z. K., Maciag, P. C., Shahabi, V., and Paterson, Y. (2013). Listeria monocytogenes-derived listeriolysin O has pathogen-associated molecular pattern-like properties independent of its hemolytic ability. Clin. Vaccine Immunol. 20, 77–84. doi: 10.1128/cvi.00488-12
Way, S. S., Thompson, L. J., Lopes, J. E., Hajjar, A. M., Kollmann, T. R., Freitag, N. E., et al. (2004). Characterization of flagellin expression and its role in Listeria monocytogenes infection and immunity. Cell. Microbiol. 6, 235–242. doi: 10.1046/j.1462-5822.2004.00360.x
Keywords: Listeria monocytogenes, prosthetic joint infection, phenotype, genome, evolution
Citation: Muchaamba F, Eshwar AK, von Ah U, Stevens MJA and Tasara T (2020) Evolution of Listeria monocytogenes During a Persistent Human Prosthetic Hip Joint Infection. Front. Microbiol. 11:1726. doi: 10.3389/fmicb.2020.01726
Received: 30 January 2020; Accepted: 01 July 2020;
Published: 28 July 2020.
Edited by:
Ludmila Chistoserdova, University of Washington, United StatesReviewed by:
Sandip Paul, Indian Institute of Chemical Biology (CSIR), IndiaYi Chen, United States Food and Drug Administration, United States
Copyright © 2020 Muchaamba, Eshwar, von Ah, Stevens and Tasara. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Taurai Tasara, tasarat@fsafety.uzh.ch; taurai.tasara@uzh.ch