- 1School of Biosciences, University of Nottingham, Sutton Bonington Campus, Sutton Bonington, United Kingdom
- 2Department of Genetics and Genome Biology, University of Leicester, Leicester, United Kingdom
- 3Sequenceanalysis.Co.uk, Innovation Centre, Norwich, United Kingdom
- 4School of Biological Sciences, University of East Anglia, Norwich, United Kingdom
Salmonella Typhimurium carrying the multidrug resistance (MDR) plasmid pMG101 was isolated from three burns patients in Boston United States in 1973. pMG101 was transferrable into other Salmonella spp. and Escherichia coli hosts and carried what was a novel and unusual combination of AMR genes and silver resistance. Previously published short-read DNA sequence of pMG101 showed that it was a 183.5Kb IncHI plasmid, where a Tn7-mediated transposition of pco/sil resistance genes into the chromosome of the E. coli K-12 J53 host strain had occurred. We noticed differences in streptomycin resistance and plasmid size between two stocks of E. coli K-12 J53 pMG101 we possessed, which had been obtained from two different laboratories (pMG101-A and pMG101-B). Long-read sequencing (PacBio) of the two strains unexpectedly revealed plasmid and chromosomal rearrangements in both. pMG101-A is a non-transmissible 383Kb closed-circular plasmid consisting of an IncHI2 plasmid sequence fused to an IncFI/FIIA plasmid. pMG101-B is a mobile closed-circular 154 Kb IncFI/FIIA plasmid. Sequence identity of pMG101-B with the fused IncFI/IncFIIA region of pMG101-A was >99%. Assembled host sequence reads of pMG101-B showed Tn7-mediated transposition of pco/sil into the E. coli J53 chromosome between yhiM and yhiN. Long read sequence data in combination with laboratory experiments have demonstrated large scale changes in pMG101. Loss of conjugation function and movement of resistance genes into the chromosome suggest that even under long-term laboratory storage, mobile genetic elements such as transposons and insertion sequences can drive the evolution of plasmids and host. This study emphasises the importance of utilising long read sequencing technologies of plasmids and host strains at the earliest opportunity.
Introduction
Clinical resistance to silver nitrate (AgNO3) in Salmonella Typhimurium was first reported by McHugh and colleagues in 1975 (Mchugh et al., 1975). The S. Typhimurium strain had been responsible for three deaths within the burns ward of Massachusetts General Hospital in 1973. Each of the three patients who died was receiving topical applications of 0.5% AgNO3 solution during treatment of severe burns. Blood- and wound-cultures all consistently produced AgNO3-resistant S. Typhimurium, which was also resistant to ampicillin, tetracycline, chloramphenicol, sulphonamides, streptomycin, and mercuric chloride (HgCl2). Conjugation studies demonstrated transmission of metal (AgNO3/HgCl2) and antibiotic resistances to other Salmonella spp. and to Escherichia coli K-12 J53. The plasmid responsible was designated pMG101 and is historically important as Mchugh et al. (1975) gave the first description of transmissible AgNO3-resistance in the clinic. Although the original S. Typhimurium isolate has not been the focal point of further published research (and may no longer exist in available culture collections), the E. coli J53/pMG101 transconjugant strain has been studied further. Pulse-field gel electrophoresis (PFGE) of pMG101 DNA indicated an approximate size of 180Kb, and a novel silver resistance operon was identified (Gupta et al., 1999). Assignment of pMG101 to the IncHI group and comparison of AgNO3-resistance encoding genes of the sil operon with other plasmids has given insights into the distribution of similar elements in IncH plasmids (Gupta et al., 2001). A whole genome sequence of a stock of E. coli J53/pMG101 (NCTC 50110) obtained from NCTC (National Collection of Type Cultures, United Kingdom) using Illumina MiSeq technology (Accession No. ASRI00000000) was previously published (Randall et al., 2015). The calculated size of pMG101 from the published short read sequence assembly was estimated as 183.5Kb, which agreed with the S1-endonuclease PFGE size estimate for pMG101 of 180Kb (Gupta et al., 1999). The sequence data showed evidence of Tn7-mediated transposition of the silver (sil) and copper (pco) resistance operons as a single unit (Tn7/sil/pco) into the E. coli J53 chromosome (Randall et al., 2015). The antimicrobial resistance gene (ARG) and metal resistance gene (MRG) content of pMG101 was however, congruent with the observed phenotypic resistances originally reported by Mchugh et al. (1975). A class D β-lactamase (bla-OXA-1), tetracycline resistance (tetA), chloramphenicol acetyltransferase (catA), sulphonamide resistance (folP), and the mercury resistance (mer) operon were all reported as being present in pMG101 (Randall et al., 2015).
IncHI plasmids constitute a diverse collection of large (75–400Kb), self-transmissible, low-copy number mobile genetic elements commonly found in Gram negative Enterobacteriaceae (Rozwandowicz et al., 2018). Molecular typing of IncHI plasmids provides discrimination between subgroups IncHI1, IncHI2, and IncHI3. IncHI2 plasmids are prevalent in Enterobacteriaceae associated with causing disease in humans and animals. As such, plasmids from this incompatibility group constitute a threat to public health and animal welfare. Of significance, Salmonella spp., E. coli, Klebsiella pneumoniae, and Enterobacter cloacae have all been identified as carriers of IncHI2 plasmids (Wong et al., 2016; Matamoros et al., 2017; Zhao et al., 2018; Chaudhry et al., 2020). Gene operons encoding metal resistances, such as tellurium (terZABCDE/Y3Y2XY1W), mercury (merRTPCADE), copper (pcoABCDRSE), and silver (silESRCBAP) are also common-place in this plasmid family (Fang et al., 2016). Co-occurrence of MRGs with multiple ARGs is widely reported (Pal et al., 2015; Li et al., 2017), and a recent study has highlighted the diverse array of ARGs encoded in IncHI2 plasmids isolated from poultry, ducks, and pigs in China (Fang et al., 2016). Extended spectrum β-lactamases (ESBL), plasmid-mediated quinolone resistance determinants (PMQRs), and Tn7/pco/sil transposition units are frequently seen in IncHI plasmids (Fang et al., 2016; Wang et al., 2017; Billman-Jacobe et al., 2018; Branger et al., 2018). Further expansion of the genetic content of IncHI plasmids can be driven via the acquisition of IncFI, IncFII, and IncN plasmids co-residing in the same bacterial cell. Fusion between two separate incompatibility groups (e.g., IncHI2/IncFII) may provide unique opportunities for plasmid host-range expansion allowing dissemination into novel hosts.
Given the historical interest of pMG101 as an archetype of AgNO3-resistance within the clinic and its relevance to current trends in AMR research, this mobile genetic element represents an important aspect of IncHI plasmid biology. Technological developments in long-read DNA sequencing provide a useful tool for resolving features of bacterial chromosomes and plasmids. Resolution of genomic repeat regions, multi-copy insertion sequences, and transposable elements is quite often beyond the reach of short-read sequencing platforms (Rhoads and Au, 2015), but long-read sequencing can overcome such issues, and in order to build upon the previously published sequence data of pMG101 (Randall et al., 2015) we decided to reanalyze this important mobile genetic element. Due to the advances in long-read and hybrid assembly technologies, it is possible to resolve some of these repeat regions. The technology may also be able to allow identification of genomic rearrangements and changes in the strain genome after longer term storage. The changes after long term storage would be expected to be small but support the theory that strains will have the capacity to undergo recombination after years in culture. Genome sequencing of microorganisms offers a snapshot to the genome and plasmidome as part of a continuous landscape (Hobman et al., 2007).
Using two separate stocks of E. coli J53/pMG101 obtained from different laboratories, PacBio sequencing of both strains and plasmids was used to provide complete closed-circular chromosomal and plasmid DNA sequences. Data obtained provided confirmation of previously published work regarding pMG101 such as ARG content and Tn7/pco/sil transposition into the E. coli chromosome. However, previously undescribed features of pMG101 including two distinct genotypes were identified from the sequence data, and further confirmed via S1-endonuclease PFGE. Here, we present analysis of each isoform – pMG101-A is a 383Kb IncHI2/IncFI/IncFII plasmid, whereas pMG101-B is a 154Kb IncFI/IncFII plasmid. Evidence obtained in the resequencing of these two plasmids, has allowed us to suggest a model for the transposition of the Tn7/pco/sil module into the chromosome of E. coli K-12 J53.
Given the medical relevance of silver as an antimicrobial able to enhance the effects of antibiotics, surveillance of resistance operons such as sil in a Tn7 transposable element across conjugative plasmids and novel hosts in lab and environmental hosts is of significance (Zhou et al., 2015).
Materials and Methods
Bacterial Strains and Growth Conditions
Two separate stocks of E. coli J53 harbouring pMG101 were obtained from Anne Summers, University of Georgia, United States (pMG101-A) and Simon Silver, University of Chicago, United States (pMG101-B). Frozen (−80°C) glycerol stocks of both strains were used to inoculate LB agar plates prior to incubation overnight at 37°C. Escherichia coli J53 carrying a chromosomal insertion of mrfp1 and kanr was used as a recipient in conjugation experiments (J53 RFP). For liquid growth, a single well-isolated colony was used to inoculate LB broth prior to overnight incubation at 37°C with shaking, unless stated otherwise.
Pulse Field Gel Electrophoresis
Pulse-field gel electrophoresis was used to determine the sizes of the two plasmids pMG101-A and pMG101-B (Hooton et al., 2011). A single colony of E. coli J53 pMG101-A and E. coli J53 pMG101-B were inoculated in to separate 10ml Miller’s lysogeny broth -LB (Merck, United Kingdom). Broth cultures were incubated at 25°C for 48h and shaking at 180rpm. Following incubation cultures were centrifuged at 4,400×g for 5min, the supernatant was decanted, pellets were re-suspended in 1ml LB Broth, and 55μl of cell suspension was transferred to a 1.5ml micro-centrifuge tube. For PFGE plug preparation, 5μl of 20mg/ml proteinase K (Sigma Aldrich, United Kingdom) was added followed by mixing 50μl of molten 1.2% PFGE grade agarose (Biorad, United States) in TE [10mM Tris-HCl (pH 7.5), 1mM EDTA] buffer. The suspension was mixed thoroughly and then set in PFGE plug moulds (Biorad, United States). The plugs were then incubated overnight at 55°C shaking at 300rpm in 1ml lysis buffer [50mM Tris-HCl (pH 8), 50mM EDTA, 1% N-lauroyl sarcosine, and 100μg/ml proteinase K]. Plugs were then washed three times for 1h in 1ml wash buffer [20mM Tris-HCl (pH 8), 50mM EDTA] at 55°C shaking at 300rpm. A 3mm slice of each plug was then loaded to a 1% PFGE grade agarose gel (100ml 1 X TAE buffer) along with PFGE lambda ladder (New England Biolabs, United States). The lanes were sealed with 1% PFGE grade agarose. A CHEF-DR II system was used for the run with a 10–30s switch time over 18h at 6V/cm, and 1 X TAE running buffer was circulated at 14°C. The gel was then stained using ethidium bromide (1μg/ml) and imaged using a Biorad ChemiDoc MP System (BioRad, United States).
Conjugation of pMG101-A and pMG101-B
A single colony of E. coli J53 pMG101-A, E. coli J53 pMG101-B and E. coli J53 RFP were inoculated in to separate 20ml Miller’s LB and incubated at 25°C with shaking until an OD600nm 1.0 was reached. About 1ml of each culture was then pelleted at 16,000g for 1min, the supernatant was removed before being resuspended and washed in 1ml maximum recovery diluent – MRD (Oxoid, United Kingdom). The wash step was repeated once more. Around 50μl of each culture containing pMG101 was added to 50μlJ53 RFP, mixed thoroughly by pipetting, and deposited on to non-selective LB agar plates and incubated at 25°C for 48h. The lawn of growth was then scraped off each plate and resuspended in 1ml MRD and vortexed thoroughly to end mating. The cells were then pelleted at 16,000g for 1min and supernatant removed and resuspended in 1ml MRD. The wash step was repeated once. The OD600nm was then measured and adjusted to 1.0 with MRD. About 10-fold dilutions of the cell suspension were then made and spread plated in triplicate onto LB agar supplemented with 25μg/ml HgCl2 and incubated at 25°C for 48h. Successful transconjugants, that produced red pigment, and were able to grow in the presence of HgCl2, were counted. Transconjugants were then selected at random, and patch plated to LB agar supplemented with 25μg/ml HgCl2 and non-selective LB agar to confirm conjugation.
DNA Extraction
Genomic DNA was extracted from overnight E. coli J53 LB broth cultures using a GenElute™ Bacterial Genomic DNA Kit (Sigma Aldrich, United Kingdom) as per manufacturer’s instructions. DNA integrity was checked by electrophoresing 5μl of each preparation in a 1% agarose gel with Quick-Load® 1Kb DNA ladder (New England Biolabs, United Kingdom) as a molecular weight marker. Genomic DNA extracts were quantified using an Invitrogen™ Qubit™ 4 Fluorometer (Fisher Scientific, United Kingdom) to determine overall yield.
Pacific Biosciences Library Preparation and Sequencing
Genomic DNA samples were purified with 1x cleaned AMPure beads (Agencourt, United Kingdom), and the quantity and quality was assessed using a Nanodrop spectrophotometer and Qubit assay. Additionally, the Fragment Analyser (Agilent, United Kingdom; using a high sensitivity genomic kit) was used to determine the average size of the DNA and the extent of degradation. This procedure was also used at the steps indicated below to determine average fragment size of the DNA. One microgram of DNA for each sample was sheared to approximately 10Kb using a Covaris G tube (Beckman Coulter, United Kingdom), as per manufacturer’s instruction. The average size was checked using a Fragment Analyser (Agilent, United Kingdom). Samples were bead purified and subjected to a DNA damage repair reaction at 37°C for 1h followed by an end repair reaction at 25°C for 5min. Bead cleaning was performed as described above. Samples were individually ligated to barcoded blunt-end adapters overnight at 25°C. After a heat inactivation step at 65°C for 10min, the samples were pooled and treated with exonucleases 3 and 7 to remove any non-circular DNA. The library was size selected on Sage BluePippin (Beverly MA, United States) using a 0.75% cassette between 5 and 50Kb to remove any shorter material. The final library concentration was measured by Qubit assay, and DNA sizes were determined by Fragment analyser. Two SMRT cells were run for this DNA using diffusion loading and sequenced on the PacBio RS system. Libraries were sequenced using 600-min movies at the Centre for Genomic Research, University of Liverpool, Crown Street, Liverpool, L69 3BX, United Kingdom.
Sequence Assembly and Bioinformatics
The resulting reads had adapters removed prior to the trimmed reads being assembled using the “canu” assembler tool, and subsequently polished using the “arrow” algorithm (Chin et al., 2013; Koren et al., 2017). Polished assemblies were quality checked for contiguity using “QUAST” with “busco” being employed to determine completeness of assemblies using Benchmarking Universal Single-Copy Orthologs (Gurevich et al., 2013; Simão et al., 2015). In this case, BUSCO V2 was used with the BUSCO proteobacteria database used as the reference. The assembled genomes for both E. coli K-12 J53 pMG101 variants were subsequently annotated using Prokka (Seemann, 2014) and screened for ARGs and MRGs using ResFinder (Bortolaia et al., 2020) and BacMet 2.0 database (Pal et al., 2014). Plasmid incompatibility groups were identified using PlasmidFinder (Carattoli et al., 2014).
Chromosome and plasmid sequences have been deposited at NCBI under Bioprojects PRJNA 701544 and PRJNA701546. Individual Genbank accession numbers host strain and plasmid sequences are E. coli J53/pMG101-A (JAFFIC010000001.1), pMG101-A (JAFFIC010000003.1), E. coli J53/pMG101-B/Tn7 insertion (CP070962), and pMG101-B (CP070963).
Antimicrobial Susceptibility Testing
The two E. coli J53 pMG101 variants were tested for antibiotic susceptibility according to the Clinical and Laboratory Standards Institute standard protocol (CLSI, 2019). The following antibiotics were tested: ampicillin 10μg (AMP), amoxicillin-clavulanic acid 20 and 10μg (AMC), cefoxitin 30μg (FOX), ceftazidime 30μg (CAZ), cefpodoxime 10μg (CPD), aztreonam 30μg (ATM), imipenem 10μg (IPM), streptomycin 10μg (S10), tetracycline 30μg (T), ciprofloxacin 5μg (CIP), nalidixic Acid 30μg (NA), trimethoprim-sulfamethoxazole 1.25μg and 23.75μg (SXT), chloramphenicol 30μg (C), nitrofurantoin 300μg (F), and azithromycin 15μg (AZM). All discs were obtained from Pro-Lab Diagnostics (Birkenhead, United Kingdom). Antimicrobial metal resistance testing for HgCl2 (184μM) and AgNO3 (400μM) was carried out according to standard techniques (Yang et al., 2020).
Results
Pulsed Field Gel Electrophoresis
Imaging of the plasmids extracted from the two J53 strains by PFGE showed a significant difference in size between the two large, low copy number plasmids. pMG101-A was identified as approximately 380Kb, more than double the size of pMG101-B at approximately 170Kb. Due to the size and the copy number of the plasmids, the band intensity was significantly reduced although was still clearly visible (Supplementary Figure S1, Supplementary data).
Sequencing of IncHI2/IncFI/II Plasmid pMG101-A and IncFI/FII Plasmid pMG101-B
Previously, we had sequenced the ~380Kb pMG101-A plasmid using Illumina MiSeq. Assembly of short reads produced an approximately 300Kb plasmid sequence, albeit fragmented in 27 contigs ranging from 120 to 34 and 220bp in length (Unpublished data). To improve assembly of pMG101-A, we utilised long-read PacBio sequencing. Assembly of 10Kb PacBio reads yielded a single closed-circular 383,246bp dsDNA molecule representative of pMG101-A. The size is in agreement with S1-endonuclease PFGE analysis of cultures of the E. coli J53 strain harbouring pMG101-A (Figure 1). Analysis of the nucleotide composition revealed a GC content of 47.8%. Identification of plasmid replicons in pMG101-A was achieved using the online tool PlasmidFinder (Carattoli et al., 2014). Five separate replicons were identified in pMG101-A representing incompatibility groups IncH and IncF (Tables 1, 2). Several features haracteristic of IncH plasmids are apparent in pMG101-A (Figure 1A) including the tellurium resistance operon (terZABCDE/terY3Y2XY1W), smr loci 0018/0199, IncH conjugal transfer (tra1/tra2) regions, and plasmid replicons of the IncHI2/IncHI2A sequence types (Carattoli et al., 2014). A single 4,634,108bp contig for the host strain E. coli K-12 J53 harbouring pMG101-A was obtained.
Figure 1. (A) Circular genetic map of IncH/IncF plasmid pMG101-A – the IncHI2/HI2A replicon (yellow), IncFIA/FIB/FII component (blue), Tn7/pco/sil mobile unit (grey), ter and mer metal resistance operons (red), Tn3/strAB antimicrobial resistance gene (ARGs; black), tetBCD/catA1 (light blue), class I integron (light green), and two copies (dark green) of ISVSa5 (IS4-like transposase) are all highlighted. (B) Circular genetic map of IncF plasmid pMG101-B – the IncFIA/FIB/FII backbone of pMG101-B (blue), class I integron (light green), tetBCD/catA1 (light blue), mer operon (red), and a single copy (dark green) of ISVSa5 (IS4-like transposase).
As was found with pMG101-A, Illumina sequencing of plasmid pMG101-B yielded a fragmented assembly (Unpublished data) containing gaps, unresolvable small contigs etc. The pco/sil genes were noticeably absent from the plasmid assembly. We used PacBio sequencing of the J53 pMG101-B strain to understand the biology of the plasmid and to confirm chromosomal integration of the pco/sil operon. A single closed-circular 154,497bp dsDNA molecule representative of pMG101-B was assembled from PacBio reads as described above (Figure 1B). A single 4,683,867bp contig representative of the host E. coli K12 J53/Tn7 insertion harbouring pMG101-B was also assembled. Analysis of the nucleotide content of pMG101-B indicates a higher GC content (53.28%) compared to pMG101-A (47.8% GC). As can be seen from Table 1, pMG101-B is an IncF plasmid containing replicon-types IncFIA, IncFIB, and IncFII.
pMG101-A and pMG101-B ARG and MRG Content
Resistance to toxic metals and antibiotics is encoded by MRGs and ARGs present in both pMG101-A and pMG101-B. Table 3 shows resistances, nucleotide positions, and associated mobile genetic elements (e.g., transposons/insertion sequences). Antibiotic sensitivity tests showed that resistance phenotypes of E. coli J53 carrying pMG101-A and pMG101-B were congruent with ARGs identified in the PacBio sequence assembly. Both J53 pMG101 strains were resistant to ampicillin, aminoglycosides, tetracycline, chloramphenicol, nalidixic acid, and sulphonamides, as well as AgNO3 and HgCl2. The discernible difference between the two variant pMG101 plasmids was strAB-encoded streptomycin resistance in pMG101-A, which was absent in pMG101-B.
The shared ARG content of pMG101-A and pMG101-B consists of a classical Tn21-family mercury resistance operon and associated class I integron harbouring antibiotic resistance gene cassettes (Figure 2). The mer operon (merRTPCADE) and class I integron structure is 20,676bp in length and forms a composite transposable element of the Tn3 family (Tn21 subgroup). Both sides of the transposon are flanked by 38bp imperfect inverted repeats (IRs) with the left-IR positioned immediately upstream of merR and the right-IR immediately downstream of the Tn21 transposon. Moving inwards from the mer operon, a further set of embedded transposable elements is seen. A Mu-like transposon (tniA/tniB) and an IS26-like element (istA/istB) constitute the region between the mer operon and the class I integron. The 3'-conserved segment of the class I integron contains a typical arrangement of sul1/qacE that encodes resistance to sulphonamides and quaternary ammonium compounds, respectively. Resistance to aminoglycoside antibiotics is encoded by the aadA1 gene cassette and is located immediately downstream of sul1. A class D β-lactamase gene, blaOXA-1, a progenitor of the well-established and studied oxacillinases is the next antibiotic resistance gene cassette in the array of the class I integron. This suggests blaOXA-1 was the most recently captured cassette in the integron cassette array. Finally, the diagnostic feature of the class I integron, the intI1 integrase, is located at the 5'-region. Immediately downstream of intI1 is the previously described Tn21 family transposon unit. Three genes encoding the modulation function (tnpM), a resolvase (tnpR), and the transposase (tnpA) form the Tn21 transposable element. This region is followed by the right-IR, which defines the end of the composite transposon. Immediately, downstream of the Tn21 unit is a chloramphenicol acetyltransferase (catA1) gene that confers resistance to phenicol antibiotics. IS1-like transposon genes (insA/insB) and a tetracycline resistance transposable unit (tetRBCD) are the final ARGs in this multidrug and metal resistance region.
Figure 2. Genetic map of Tn21 present in both pMG101-A and pMG101-B – the backbone consisting of the mer operon (red) and transposase genes (green). The insertion of the Tn402-like element containing Class I integron with cassette array encoding the following resistance genes: blaOXA-1, aadA1, qacE∆1, and sul1. The Tn21 mobile genetic element is flanked by two 38bp imperfect inverted repeats.
IS-Mediated Fusion of IncHI/IncF pMG101-A
As discussed previously, the 383,246bp pMG101-A is composed of IncH1/H12A/FIA/FII/FIB plasmid sequences. Based on sequence analysis and mapping of insertion sequences and transposable elements, it is possible to hypothesise a chain of events leading to fusion of the various Inc. plasmid regions. Downstream of conjugal transfer region 2 of the IncH1/H12A element is an ISVSa5 (IS4-like transposon) that defines the beginning of the IncF region. Approximately, 157Kb further downstream is a further ISVSa5 transposon that delineates the end of the IncF plasmid sequence. This ISVSa5 element is located immediately downstream of tetD. Due to the presence of these elements, at the junctions of the IncH components of the plasmid, it is possible that the integration event has been driven by ISVSa5-mediated recombination. An interesting feature of the IncF region is its similarity to a chromosomal genomic island (GI-DT12) previously observed in S. Typhimurium T000240. Chromosomal insertion of this large genomic island is reported as being mediated by IS1-like transposable elements (Izumiya et al., 2011). The major differences between GI-DT12 and the homologous regions found in both pMG101 variants are the presence of IncF replicons in the plasmid nucleotide sequences. The two copies of ISVSa5 that flank the IncF region of pMG101-A may have arisen following a duplication event that occurred during recombination. Analysis of pMG101-B shows that there is only a single copy of ISVSa5 present in this plasmid, which is positioned downstream of tetD.
Factors Influencing Conjugal Transfer of pMG101-A
Multiple attempts to transfer pMG101-A via conjugation into different E. coli recipient strains at 25 and 37°C were undetected, whereas pMG101-B transferred into E. coli recipients at 25°C. The pMG101-B plasmid was transferred into the recipient strain E. coli K-12 J53 RFP kanr at a frequency of 1.25×10−3 conjugants/ml. Analysis of the conjugal transfer regions of pMG101-A highlighted the presence of an ISEc27 element that had integrated into orf009 of the transfer region 2 (tra2). It has recently been shown that orf009 (herein designated rsp) encodes an essential component of the conjugation apparatus – the RSP/R27-secreted protein (Huttener et al., 2019), suggestive that ISEc27 has inactivated the transfer ability of pMG101-A. ISEc27 is found as multiple copies in both E. coli K12 J53 chromosomes with a total of nine in the strain harbouring pMG101-A (plus one insertion in pMG101-A rsp), and seven copies in E. coli J53 pMG101-B. No insertions of ISEc27 were identified in pMG101-B. Expansion of mobile genetic elements such as ISEc27 in long-term stocks derived from strains such as E. coli K12 J53 should be monitored for such events.
Tn7-Mediated sil/pco Chromosomal Integration
As previously reported the NCTC stock of E. coli J53 pMG101 strain carries a Tn7/sil/pco module chromosomally, conferring resistance to copper- and silver-containing compounds (Randall et al., 2015). In most instances, mobility of Tn7-elements is driven via Tn7-mediated transposition into the canonical attTn7 site located at the 3'-end of glmS in the E. coli chromosome. Analysis of the Tn7/sil/pco region in pMG101-A shows that the complete module (tnsABCDE/silESRCFBA/pcoEABCDRS) is located in the plasmid DNA assembly – specifically within the IncHI region of the plasmid. IncHI plasmids and their relatives are recognised as providing an environmental reservoir of MRGs, including pco/sil within the family Enterobacteriaceae (Gupta et al., 2001; Fang et al., 2016). For pMG101-B, the Tn7/sil/pco module has transposed into the chromosome of the E. coli J53 host. The transposition event has mediated the integration of Tn7/sil/pco in between yhiN (NAD[P]/FAD-dependent oxidoreductase) and yhiM (inner membrane protein; Figure 3). The overall effect of this is that yhiN/yhiM is now split by the Tn7/sil/pco module such that the two genes are separated by a genomic island of approximately 33.7Kb in length. The absence of Tn7/sil/pco in the plasmid sequence of pMG101-B is also coupled to the complete absence of all the recognisable features of the IncH component of pMG101-A. Analysis of glmS in E. coli J53/pMG101-B shows that the site is unoccupied with respect to the presence of Tn7-like elements. The use of alternate insertion sites in the E. coli chromosome for the Tn7/sil/pco that are not the canonical attTn7 has also been observed in pig E. coli isolates (Chalmers et al., 2018).
Figure 3. Tn7-mediated transposition of metal resistance operons pco and sil into Escherichia coli J53 chromosome – the transposition event has resulted in insertion of the mobile unit between yhiN and yhiM.
Discussion
Observable differences in antibiotic resistances and PFGE size profiles of plasmid pMG101 between what was expected and previously reported data (Randall et al., 2015) drove us to investigate the underlying mechanisms driving these changes. Previously, published Illumina MiSeq data for pMG101 hinted at some unusual biology but the fragmented assembly (Randall et al., 2015) and difficulties faced when resolving repeat sequences (Cahill et al., 2010) suggested a long-read sequencing approach may improve our current understanding of pMG101. Therefore, we applied PacBio long-read sequencing in order to resolve the important features of pMG101, gain complete closure of the circular genetic element, and determine any alterations to the host E. coli K12 J53 chromosome. Working from two separately obtained stocks of pMG101, it was possible to show evidence of plasmid recombination between IncH and IncF plasmids (pMG101-A and pMG101-B), Tn7-mediated transposition of MRGs into a non-canonical attTn7 site in the chromosome of E. coli J53 (pMG101-B), and loss of mobilisation (pMG101-A). The transposition event of pMG101-B Tn7/sil/pco is consistent with other reports that the Tn7 element can insert into the chromosome of other E. coli strains at non-canonical attTn7 sites (Chalmers et al., 2018). Transposition of the Tn7 element resulted in the host strain gaining chromosomal resistance to silver and copper in the form of silESRCFBAGP and pcoEABCDRS, which form with the Tn7 transposition machinery – a recognised, likely ancient, and widely distributed, mobile genetic element (Hao et al., 2015; Hobman and Crossman, 2015; Staehlin et al., 2016). The use of this attTn7 suggests that the chromosomal insertion occurred via the TnsABC+D rather than the TnsABC+E system. Whilst both of the mechanisms are able to induce transposition to the chromosome, TnsABC+D promotes insertion by “cut and paste” mechanism to the specific attTn7 sites (5'-CCCGC-3'), which is duplicated and normally 23 bases from the C-terminal end of the glmS gene (Waddell and Craig, 1989; Bainton et al., 1993; Mitra et al., 2010). In the case of pMG101-B, the same attTn7 site is observed between the yhiN and M genes in the E. coli J53 chromosome. Tn7 and Tn7-like elements undergo chromosomal insertion by inserting downstream (3') adjacent to glmS within the host chromosome, as it is considered a safe site that does not interfere with host cell function and is therefore more likely to propagate to the daughter cells (Peters and Craig, 2001). The mechanism for insertion has been exploited and used to produce high copy number expression vectors such as pRM2, a pUC18 backbone with a Tn7 insertion site also presenting resistance to ampicillin (Mckown et al., 1988).
From the assembled data, it is seen here that from the original location in pMG101-B, Tn7/pco/sil has inserted using this mechanism between yhiN and yhiM. Chromosomal insertion of Tn7 is documented to be less frequent than transposition to conjugal plasmids. In this case, pMG101-A has lost mobilisation ability due to orf009, a part of the tra genes needed for conjugative plasmid transfer, having an integrative and conjugative element (ICE) inserting in the middle of it and separating it into two parts. Over time, it is possible that the loss of conjugative ability of the plasmid may have increased the likelihood that the TnsABC+D mechanism for chromosomal insertion of the Tn7 element as a means to make sure it is passed onto daughter cells (Craig, 1996). This hypothesis is in agreement with the reported function of TnsC, which may evaluate insertion target sites for attTn7 sites acting as a regulator for the TnsD and E proteins. TnsD then binds to a consensus sequence allowing for insertion of the Tn7 (Hauer and Shapiro, 1984). While a similar chain of events has allowed pMG101-B/Tn7 to integrate at an alternative chromosomal locus, limited sequence homology between attTn7/glmS and the insertion site of yhiNM is observed.
Insertion of mobile genetic elements to host DNA may be advantageous, but equally may have less favourable implications. Transmissibility of IncHI plasmids from host to recipient cells is determined by the presence of specific transfer regions – Tra1 and Tra2. In this instance, the ISEc27 element that has undergone expansion in numbers in the E. coli K12 J53 (pMG101-A) host strain is inserted on the opposite strand of Tra2 rsp, is transcribed in the opposite direction to all other genes encoded in the gene operon and may well also be exerting polar effects on genes downstream that are recognised as being essential for conjugal transfer (trhFWUN). For IncHI plasmids, R27 isolated from S. Typhi and prototypical founder for the IncHI1 plasmid group remains one of the best studied plasmids of this type (Sherburne et al., 2000; Lawley et al., 2003; Huttener et al., 2019; Luque et al., 2019). In S. Typhimurium SL1344, the 155.4kDa RSP was initially identified in cell-free supernatants and found to be expressed in a temperature-dependent manner (optimal expression 25°C). RSP is exported to the extracellular surface via the plasmid-encoded type IV secretion system in a trhC-dependent manner. Externally RSP has been identified as primarily interacting with flagella, which consequently reduces cell-motility, an initial step towards bacterial conjugation. However, parallel studies with non-flagellated cells also indicate other external cellular components may interact with RSP to reduce motility (Huttener et al., 2019). Disruption of essential transfer genes such as rsp would therefore likely have a negative effect on conjugation ability as observed for E. coli K12 J53 (pMG101-A).
From the two E. coli J53 pMG101 cell stocks we analysed and the similarity of one stock (of pMG101-B) to the NCTC50110 stock sequenced by Randall et al. (2015), it is quite clear that there are some significant variations between the two plasmids and host strains. The fusion of a whole plasmid replicon larger than the first reported analysis of pMG101 in one stock (Gupta et al., 2001) and loss of transmissibility resulting in chromosomal insertion of the Tn7/pco/sil mobile element in the other suggest that bacterial evolution during long term storage and use in laboratories has occurred. This has broader, but unsurprising, implications for our understanding of isolates of plasmid containing strains, suggesting that in the real world these are likely to be transient forms, which can also evolve in the laboratory. Whilst sequence analysis demonstrates what the plasmid structure is, what the sequence of events that led up to these significant differences in pMG101 was, is less easy to reconstruct in a timeline. Over 20years passed between the original reports of transfer of pMG101 into E. coli K-12 and reports of the size of the plasmid by PFGE. It is clear that at some point in the history of pMG101, two plasmids fused, one of which carried a Tn7/sil/pco mobile genetic element and the other of which carried SGI11. It is also clear that the Tn7/sil/pco element transposed into a non-standard attTn7 site in the E. coli chromosome. As there is now cheap and readily accessible long-read sequencing available for bacterial strains, we argue that researchers should take advantage of this technology for early sequencing of minimally manipulated isolates and use sequencing as a tool to regularly check strain integrity.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.
Author Contributions
SH, AM, DS, and JH conceived and designed the experiments. SH, AP, and C-GH performed experimental work. SH, AP, LC, C-GH, AM, and JH: bioinformatic analysis. SH, AP, KA, C-GH, LC, DS, AM, and JH: wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by UK NERC grants (EVAL-farms NE/N019881/1 – AM, DS, and JH and NERC NBAF1115 – SH and AM) and funding from UK BBSRC to KA, C-GH, and AP. Bioinformatic analysis was carried out on MRC-CLIMB infrastructure (MR/L015080/1).
Conflict of Interest
LC is the director of SequenceAnalysis.co.uk and was paid to perform this work.
The remaining 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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Simon Silver and Anne Summers for providing the pMG101 plasmids used in this study, and for their continued support and encouragement.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.723322/full#supplementary-material
Supplementary Figure S1 | PFGE analysis of S1-endonuclease digestion of pMG101-A (left image) and pMG101-B (right image).
References
Bainton, R. J., Kubo, K. M., Feng, J. N., and Craig, N. L. (1993). Tn7 transposition: target DNA recognition is mediated by multiple Tn7-encoded proteins in a purified in vitro system. Cell 72, 931–943. doi: 10.1016/0092-8674(93)90581-A
Billman-Jacobe, H., Liu, Y., Haites, R., Weaver, T., Robinson, L., Marenda, M., et al. (2018). pSTM6-275, a conjugative IncHI2 plasmid of salmonella enterica That confers antibiotic and heavy-metal resistance under changing physiological conditions. Antimicrob. Agents Chemother. 62, e02357–e02317. doi: 10.1128/AAC.02357-17
Bortolaia, V., Kaas, R. F., Ruppe, E., Roberts, M. C., Schwarz, S., Cattoir, V., et al. (2020). ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 75, 3491–3500. doi: 10.1093/jac/dkaa345
Branger, C., Ledda, A., Billard-Pomares, T., Doublet, B., Fouteau, S., Barbe, V., et al. (2018). Extended-spectrum beta-lactamase-encoding genes are spreading on a wide range of Escherichia coli plasmids existing prior to the use of third-generation cephalosporins. Microb. Genom. 4:e000203. doi: 10.1099/mgen.0.000203
Cahill, M. J., Koser, C. U., Ross, N. E., and Archer, J. A. (2010). Read length and repeat resolution: exploring prokaryote genomes using next-generation sequencing technologies. PLoS One 5:e11518. doi: 10.1371/journal.pone.0011518
Carattoli, A., Zankari, E., Garcia-Fernandez, A., Voldby Larsen, M., Lund, O., Villa, L., et al. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents Chemother. 58, 3895–3903. doi: 10.1128/AAC.02412-14
Chalmers, G., Rozas, K. M., Amachawadi, R. G., Scott, H. M., Norman, K. N., Nagaraja, T. G., et al. (2018). Distribution of the pco gene cluster and associated genetic determinants among swine Escherichia coli from a controlled feeding trial. Gene 9:504. doi: 10.3390/genes9100504
Chaudhry, T. H., Aslam, B., Arshad, M. I., Alvi, R. F., Muzammil, S., Yasmeen, N., et al. (2020). Emergence of Bla NDM-1 harboring Klebsiella pneumoniae ST29 and ST11 in veterinary settings and waste of Pakistan. Infect. Drug Resist. 13, 3033–3043. doi: 10.2147/IDR.S248091
Chin, C.-S., Alexander, D. H., Marks, P., Klammer, A. A., Drake, J., Heiner, C., et al. (2013). Nonhybrid, finished microbial genome assemblies from long-read SMRT sequencing data. Nat. Methods 10, 563–569. doi: 10.1038/nmeth.2474
CLSI (2019). Clinical and Laboratory Standards Institute CLSI performance standards for antimicrobial disk and dilution susceptibility tests for bacteria isolated from animals; second informational supplement; VET01-S2; Clinical and Laboratory Standards Institute. Available at: https://clsi.org/standards/products/veterinary-medicine/documents/vet01/ (Accessed July 02, 2019).
Craig, N. L. (1996). Transposon Tn7. Curr. Top. Microbiol. Immunol. 204, 27–48. doi: 10.1007/978-3-642-79795-8_2
Fang, L., Li, X., Li, L., Li, S., Liao, X., Sun, J., et al. (2016). Co-spread of metal and antibiotic resistance within ST3-IncHI2 plasmids from E. coli isolates of food-producing animals. Sci. Rep. 6:25312. doi: 10.1038/srep25312
Gupta, A., Matsui, K., Lo, J. F., and Silver, S. (1999). Molecular basis for resistance to silver cations in salmonella. Nat. Med. 5, 183–188. doi: 10.1038/5545
Gupta, A., Phung, L. T., Taylor, D. E., and Silver, S. (2001). Diversity of silver resistance genes in IncH incompatibility group plasmids. Microbiology 147, 3393–3402. doi: 10.1099/00221287-147-12-3393
Gurevich, A., Saveliev, V., Vyahhi, N., and Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinformatics 29, 1072–1075. doi: 10.1093/bioinformatics/btt086
Hao, X., Luthje, F. L., Qin, Y., Mcdevitt, S. F., Lutay, N., Hobman, J. L., et al. (2015). Survival in amoeba--a major selection pressure on the presence of bacterial copper and zinc resistance determinants? Identification of a “copper pathogenicity island.” Appl. Microbiol. Biotechnol. 99, 5817–5824. doi: 10.1007/s00253-015-6749-0
Hauer, B., and Shapiro, J. A. (1984). Control of Tn7 transposition. Mol. Gen. Genet. 194, 149–158. doi: 10.1007/BF00383510
Hobman, J. L., and Crossman, L. C. (2015). Bacterial antimicrobial metal ion resistance. J. Med. Microbiol. 64, 471–497. doi: 10.1099/jmm.0.023036-0
Hobman, J. L., Penn, C. W., and Pallen, M. J. (2007). Laboratory strains of Escherichia coli: model citizens or deceitful delinquents growing old disgracefully? Mol. Microbiol. 64, 881–885. doi: 10.1111/j.1365-2958.2007.05710.x
Hooton, S. P., Atterbury, R. J., and Connerton, I. F. (2011). Application of a bacteriophage cocktail to reduce salmonella Typhimurium U288 contamination on pig skin. Int. J. Food Microbiol. 151, 157–163. doi: 10.1016/j.ijfoodmicro.2011.08.015
Huttener, M., Prieto, A., Aznar, S., Bernabeu, M., Glaria, E., Valledor, A. F., et al. (2019). Expression of a novel class of bacterial Ig-like proteins is required for IncHI plasmid conjugation. PLoS Genet. 15:e1008399. doi: 10.1371/journal.pgen.1008399
Izumiya, H., Sekizuka, T., Nakaya, H., Taguchi, M., Oguchi, A., Ichikawa, N., et al. (2011). Whole-genome analysis of salmonella enterica serovar Typhimurium T000240 reveals the acquisition of a genomic island involved in multidrug resistance via IS1 derivatives on the chromosome. Antimicrob. Agents Chemother. 55, 623–630. doi: 10.1128/AAC.01215-10
Koren, S., Walenz, B. P., Berlin, K., Miller, J. R., Bergman, N. H., and Phillippy, A. M. (2017). Canu: scalable and accurate long-read assembly via adaptive k-mer weighting and repeat separation. Genome Res. 27, 722–736. doi: 10.1101/gr.215087.116
Lawley, T. D., Klimke, W. A., Gubbins, M. J., and Frost, L. S. (2003). F factor conjugation is a true type IV secretion system. FEMS Microbiol. Lett. 224, 1–15. doi: 10.1016/S0378-1097(03)00430-0
Li, L. G., Xia, Y., and Zhang, T. (2017). Co-occurrence of antibiotic and metal resistance genes revealed in complete genome collection. ISME J. 11, 651–662. doi: 10.1038/ismej.2016.155
Luque, A., Paytubi, S., Sanchez-Montejo, J., Gibert, M., Balsalobre, C., and Madrid, C. (2019). Crosstalk between bacterial conjugation and motility is mediated by plasmid-borne regulators. Environ. Microbiol. Rep. 11, 708–717. doi: 10.1111/1758-2229.12784
Matamoros, S., Van Hattem, J. M., Arcilla, M. S., Willemse, N., Melles, D. C., Penders, J., et al. (2017). Global phylogenetic analysis of Escherichia coli and plasmids carrying the mcr-1 gene indicates bacterial diversity but plasmid restriction. Sci. Rep. 7:15364. doi: 10.1038/s41598-017-15539-7
Mchugh, G. L., Moellering, R. C., Hopkins, C. C., and Swartz, M. N. (1975). Salmonella typhimurium resistant to silver nitrate, chloramphenicol, and ampicillin. Lancet 1, 235–240. doi: 10.1016/s0140-6736(75)91138-1
Mckown, R. L., Orle, K. A., Chen, T., and Craig, N. L. (1988). Sequence requirements of Escherichia coli attTn7, a specific site of transposon Tn7 insertion. J. Bacteriol. 170, 352–358. doi: 10.1128/jb.170.1.352-358.1988
Mitra, R., Mckenzie, G. J., Yi, L., Lee, C. A., and Craig, N. L. (2010). Characterization of the TnsD-attTn7 complex that promotes site-specific insertion of Tn7. Mob. DNA 1:18. doi: 10.1186/1759-8753-1-18
Pal, C., Bengtsson-Palme, J., Kristiansson, E., and Larsson, D. G. (2015). Co-occurrence of resistance genes to antibiotics, biocides and metals reveals novel insights into their co-selection potential. BMC Genomics 16:964. doi: 10.1186/s12864-015-2153-5
Pal, C., Bengtsson-Palme, J., Rensing, C., Kristiansson, E., and Larsson, D. G. (2014). BacMet: antibacterial biocide and metal resistance genes database. Nucleic Acids Res. 42, D737–D743. doi: 10.1093/nar/gkt1252
Peters, J. E., and Craig, N. L. (2001). Tn7: smarter than we thought. Nat. Rev. Mol. Cell Biol. 2, 806–814. doi: 10.1038/35099006
Randall, C. P., Gupta, A., Jackson, N., Busse, D., and O'neill, A. J. (2015). Silver resistance in Gram-negative bacteria: a dissection of endogenous and exogenous mechanisms. J. Antimicrob. Chemother. 70, 1037–1046. doi: 10.1093/jac/dku523
Rhoads, A., and Au, K. F. (2015). PacBio sequencing and its applications. Genomics Proteomics Bioinformatics 13, 278–289. doi: 10.1016/j.gpb.2015.08.002
Rozwandowicz, M., Brouwer, M. S. M., Fischer, J., Wagenaar, J. A., Gonzalez-Zorn, B., Guerra, B., et al. (2018). Plasmids carrying antimicrobial resistance genes in Enterobacteriaceae. J. Antimicrob. Chemother. 73, 1121–1137. doi: 10.1093/jac/dkx488
Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics 30, 2068–2069. doi: 10.1093/bioinformatics/btu153
Sherburne, C. K., Lawley, T. D., Gilmour, M. W., Blattner, F. R., Burland, V., Grotbeck, E., et al. (2000). The complete DNA sequence and analysis of R27, a large IncHI plasmid from salmonella typhi that is temperature sensitive for transfer. Nucleic Acids Res. 28, 2177–2186. doi: 10.1093/nar/28.10.2177
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V., and Zdobnov, E. M. (2015). BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212. doi: 10.1093/bioinformatics/btv351
Staehlin, B. M., Gibbons, J. G., Rokas, A., O'halloran, T. V., and Slot, J. C. (2016). Evolution of a heavy metal homeostasis/resistance island reflects increasing copper stress in Enterobacteria. Genome Biol. Evol. 8, 811–826. doi: 10.1093/gbe/evw031
Waddell, C. S., and Craig, N. L. (1989). Tn7 transposition: recognition of the attTn7 target sequence. Proc. Natl. Acad. Sci. U. S. A. 86, 3958–3962.
Wang, J., Li, Y., Xu, X., Liang, B., Wu, F., Yang, X., et al. (2017). Antimicrobial resistance of Salmonella enterica Serovar Typhimurium in Shanghai, China. Front. Microbiol. 8:510. doi: 10.3389/fmicb.2017.00510
Wong, M. H., Chan, E. W., Xie, L., Li, R., and Chen, S. (2016). IncHI2 plasmids are the key vectors responsible for oqxAB transmission among Salmonella species. Antimicrob. Agents Chemother. 60, 6911–6915. doi: 10.1128/AAC.01555-16
Yang, H., Wei, S. H., Hobman, J. L., and Dodd, C. E. R. (2020). Antibiotic and metal resistance in Escherichia coli isolated from pig slaughterhouses in the United Kingdom. Antibiotics 9:746. doi: 10.3390/antibiotics9110746
Zhao, H., Chen, W., Xu, X., Zhou, X., and Shi, C. (2018). Transmissible ST3-IncHI2 plasmids are predominant carriers of diverse complex IS26-class 1 integron arrangements in multidrug-resistant Salmonella. Front. Microbiol. 9:2492. doi: 10.3389/fmicb.2018.02492
Keywords: plasmid, silver resistance, pMG101, PacBio, recombination
Citation: Hooton SPT, Pritchard ACW, Asiani K, Gray-Hammerton CJ, Stekel DJ, Crossman LC, Millard AD and Hobman JL (2021) Laboratory Stock Variants of the Archetype Silver Resistance Plasmid pMG101 Demonstrate Plasmid Fusion, Loss of Transmissibility, and Transposition of Tn7/pco/sil Into the Host Chromosome. Front. Microbiol. 12:723322. doi: 10.3389/fmicb.2021.723322
Edited by:
Elisabeth Grohmann, Beuth Hochschule für Technik Berlin, GermanyReviewed by:
Pablo Vinuesa, National Autonomous University of Mexico, MexicoPriscilla Brunetto, Université de Fribourg, Switzerland
Copyright © 2021 Hooton, Pritchard, Asiani, Gray-Hammerton, Stekel, Crossman, Millard and Hobman. 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: Jon L. Hobman, jon.hobman@nottingham.ac.uk
†Present address: Karishma Asiani, Diaceutics PLC, Belfast, United Kingdom; Charlotte J. Gray-Hammerton,Department of Zoology, University of Oxford, Oxford, UK; Lisa C. CrossmanNRP Innovation Centre, Norwich, United Kingdom