Using unique ORFan genes as strain-specific identifiers for Escherichia coli
BMC Microbiology volume 22, Article number: 135 (2022)
Bacterial identification at the strain level is a much-needed, but arduous and challenging task. This study aimed to develop a method for identifying and differentiating individual strains among multiple strains of the same bacterial species. The set used for testing the method consisted of 17 Escherichia coli strains picked from a collection of strains isolated in Germany, Spain, the United Kingdom and Vietnam from humans, cattle, swine, wild boars, and chickens. We targeted unique or rare ORFan genes to address the problem of selective and specific strain identification. These ORFan genes, exclusive to each strain, served as templates for developing strain-specific primers.
Most of the experimental strains (14 out of 17) possessed unique ORFan genes that were used to develop strain-specific primers. The remaining three strains were identified by combining a PCR for a rare gene with a selection step for isolating the experimental strains. Multiplex PCR allowed the successful identification of the strains both in vitro in spiked faecal material in addition to in vivo after experimental infections of pigs and recovery of bacteria from faecal material. In addition, primers for qPCR were also developed and quantitative readout from faecal samples after experimental infection was also possible.
The method described in this manuscript using strain-specific unique genes to identify single strains in a mixture of strains proved itself efficient and reliable in detecting and following individual strains both in vitro and in vivo, representing a fast and inexpensive alternative to more costly methods.
The tracing of microbes in complex biological systems is indispensable to answer many scientific questions in applied, clinical, and environmental microbiology. Bacterial identification at the strain level is highly challenging since closely related bacteria may have similar morphologic and physiologic profiles, and strains belonging to the same bacterial species and even to the same bacterial family, are indistinguishable solely by morphological methods . Techniques to identify strains have evolved and advanced in the past decades, but all available methods have limitations and flaws, and none of them is 100% reliable or accurate . Differentiation of bacterial strains based on their genome information has become the preferred approach due to its excellent resolution, high reliability, and easy availability . Genotyping methods can be grouped into three major categories . The first is based on DNA fragment patterns, in which amplification and/or enzymatic digestion of bacterial DNA is followed by electrophoretic resolution of differently sized fragments, the pattern of which serves as specific species or strain identifier. For amplicon-based classification, the DNA sequence of a reference genome has to be at least partially known. Second, DNA-hybridisation systems deploy nucleic acid-based probes, labelled fragments of known sequence complementary to their corresponding targets, which are detected after probe binding. Third, DNA sequence-based genotyping, the most powerful tool currently being used to classify bacteria, utilizes strain-specific variations, such as single nucleotide polymorphisms (SNP)  as well as deletion or addition of genetic material [5, 6]. These methods all rely on specific differences between individual strains, which have to be at least qualitatively, but preferably quantitatively, detectable within the context of a complex microbial environment, in which one or more strains of the target species may already be present.
Escherichia coli is a commensal member of the vertebrate gut, but certain strains, grouped into E. coli pathovars, have acquired virulence genes, mainly located on mobile elements, on multiple occasions resulting in a high degree of genomic flux [7, 8]. E. coli has an open pan-genome, i.e., the number of genes it contains (consisting of core and accessory genes) increases with the number of additional genomes sequenced. Many accessory genes are only present in one strain (also known as singletons) or in a few strains [9, 10]. Such genes with no known relatives or homologues in species belonging to other lineages are known as orphan or ORFan genes . Improved genotyping and sequencing techniques during the last two decades have led to the discovery of large numbers of ORFan genes present within bacterial genomes [11, 12].
New studies analysing the genomes of pathogenic and non-pathogenic bacteria have shown that all genomes of a single species have their own share of specific and unique ORFan genes, which are common to some strains but are not found consistently in all members of the species: this part of the genome is also considered to belong to the “variable or accessory genome” [13,14,15]. Several ORFan genes are lineage-specific, while others can even be strain-specific and seem to contribute to a particular strain’s characteristics, including potential pathogenicity [16, 17]. These features qualify specific ORFan genes to serve as molecular tracer for bacteria at the strain level. In this study, we describe an ORFan gene-based identification method characterized by the generation and development of PCR primers specific at the bacterial strain level. An extensive E. coli library containing 1198 whole-genome-sequenced strains collected by the consortium served as database to identify suitable ORFan genes. Specific regions within these ORFan genes were tested against larger sequence databases to increase the marker’s probability of being restricted to a single strain. Specific PCR and quantitative PCR (qPCR) primers derived from the search results were successfully applied in vitro and in vivo in an animal infection experiment.
Materials and methods
One thousand one hundred and ninety-eight E. coli strains were provided by all HECTOR project partners and whole-genome-sequenced using short-read technology (van der Putten, B., Tiwari, S.K., HECTOR consortium, Semmler, T. and Schultsz, C., unpublished data https://doi.org/10.1101/2022.02.08.479532). The sequencing was executed using Illumina MiSeq (2 × 150 bp, 2 × 250 bp, 2 × 300 bp paired reads) and the Illumina HiSeq 4000 system (2 × 150 bp paired reads). As quality control, adapter sequences and low-quality bases within raw reads were trimmed. For genome assembly and annotation, adapter-trimmed reads were assembled with SPAdes v3.13.1 using read correction . Scaffolds smaller than 500 bp were discarded. QUAST v5.0.0 was used to assess assembly quality using default parameters . Data on the assemblies is given in Supplementary Table S1. Seventeen strains, carrying extended-spectrum beta-lactamase (ESBL) genes, were selected from this collection for an animal study to assess their colonization properties in livestock animals in parallel with an infection approach using a mixture of strains (cocktail) (Table 1).
To facilitate experimental strain detection in complex, non-sterile matrices by adding a selective culturing step, all 17 strains were artificially selected for rifampicin resistance. To this end, a 50 ml overnight culture of each strain was centrifuged at 4000 x g for 10 min, the pellet resuspended in 1 ml of LB media, and then plated on LB-agar plates containing 100 μg/ml rifampicin (Carl Roth, Karlsruhe, Germany). After overnight incubation at 37 °C, the plates were inspected for colony growth. The rpoB genes of the rifampicin-resistant strains were sequenced and one mutant was selected per experimental strain.
ORFan gene identification
The DNA isolated from each strain was sequenced using Illumina short read technology (see Supplemental Table S1). The draft genomes were annotated by Prokka v1.13  using a genus-specific blast for E. coli. The pan-genome was constructed at 95% amino acid identity by using Roary v3.12.0 . Genes found in 99% of the strains within our collection were considered to represent core genes and the remaining genes classified as accessory genes. Paralogs were split into different orthologous groups. The strain-specific genes were identified by in-house scripts based on the binary matrix of gene presence or absence obtained from Roary. The nucleotide sequences of these strain-specific genes were extracted for each strain, and their specificity was further confirmed by using BLAST . Each strain-specific gene was scanned against the entire gene pool of the HECTOR strain collection. Genes found in other strains or in more than one copy in the same strain at 90% identity and 90% coverage were discarded. A gene with a single copy present only in one strain was considered as strain-specific in this study. Its corresponding sequence was extracted based on the respective locus-tag. These genes were further examined individually via online BLAST analysis  against the GenBank  database specific for E. coli in order to assess exclusivity. As databases, the “Standard database” and the “Nucleotide collection (nr/nt)” were selected, and the organism was specified as “Escherichia coli (taxid: 562)”. As program selection, the search was optimized for “highly similar sequences (megablast)” and when no results were given, a new search was performed with the “somewhat similar sequence (blastn)” option selected. At the moment of the analysis (between February and June 2019), the GenBank database was composed of 4 representative E. coli genomes, more than 8000 complete E. coli genomes and more than 31,000,000 drafted E. coli genomes. If the top results from the BLAST  analysis of the GenBank  database belonged to known plasmids and fully matched the sequences being compared, they were considered plasmid-derived, and excluded from further analysis. If no hits were found from the search, a more exhaustive search against the non-redundant nucleotide database from GenBank was carried out. Those genes with few (< 10) or no hits were considered ORFans in this analysis.
PCR and qPCR primer design
The ORFan genes detected in silico were used to design primers. If no “low-hit” ORFans were available for a specific strain, those regions within a selected ORFan gene, identified in another strain, but showing a higher sequence variation when blasted against the GenBank database were used to generate the primers. All PCR and qPCR primers were manually designed for each strain according to the following criteria: primer length and melting temperature, avoidance of dimers, hairpins, and self-complementarity. Specificity verification was performed using the tool Primer-BLAST . The PCR primers were designed so that four multiplex PCR reactions could be performed in order to qualitatively (i.e., presence/absence from a given sample) trace all 17 experimental strains with a minimum number of reactions. For this, both melting temperatures of the primers and sizes of the products within a multiplex were selected so that the temperature should not differ more than 1 °C between primer pairs, and product sizes should differ by at least 100 bp from each other, if possible. The oligonucleotide primers were synthesized by Eurofins Genomics (Ebersberg, Germany). Their sequences are listed in Table 2.
Primers for qPCR were optimized for use with the Luna qPCR mastermix (New England Biolabs Inc., Ipswich, MA, USA). Amplicon sizes ranged between 100 and 250 bp, with a GC content from 40 to 60%, a melting temperature not greater than 61 °C with less than 1 °C difference between primers of the same pair, and a primer length of 19 to 25 nucleotides. If possible, ORFan gene and ORFan gene region used to design the qPCR primers were the same as the ones used for PCR primer design. The oligonucleotide primers were manufactured by Eurofins Genomics (Ebersberg, Germany). Their sequences are listed in Table 3.
Four multiplexes were needed in order to detect all 17 experimental strains. As much as possible, multiplexes targeted strains originating from the same host. In the case of the “Mix multiplex”, due either to fragment size or melting temperature three strains were included that didn’t fit in the other multiplexes. Multiplex PCR conditions were optimised following the recommendations published by Zangenberg et al.  or the PCR mastermix manufacturer. PCR was performed in a total volume of 25 μl containing 2 μl of purified DNA, 12.5 μl of OneTaq 2x Master Mix with Standard Buffer (New England Biolabs Inc., Ipswich, MA, USA), 9.5 μl of nuclease-free water, and 0.5 μl each of 10 μM forward and reverse primer. The reactions were performed in a Biometra T3 Thermocycler System (Analytik Jena, Jena, Germany) using the conditions specified in Table 4.
Gel electrophoresis was performed by using 1.0% all-purpose, high-purity agarose (VWR International, Radnor, PA, USA) gels with 0.25X SERVA DNA stain clear G (SERVA Electrophoresis GmbH, Heidelberg, Germany) in 1X Tris-borate-EDTA buffer (VWR International, Radnor, PA, USA) in a Perfect Blue Gel System (VWR International, Radnor, PA, USA) at 150 V for 1 h. Two microliters of amplified DNA were mixed with 4 μl of gel loading dye (New England Biolabs Inc., Ipswich, MA, USA) for analysis. For reference, a Quick-Load Purple 100 bp DNA Ladder (New England Biolabs Inc., Ipswich, MA, USA) was used (bands every 100 bp up to 1000 bp, plus two additional bands at 1200 bp and 1500 bp).
Quantitative PCR reaction setup
Quantitative PCR conditions were optimised following the qPCR mastermix manufacturer’s specifications. Assays for qPCR were performed on a CFX96 Touch Real-Time PCR Detection System (Bio-Rad Laboratories, Hercules, CA, USA). Reactions contained a total volume of 20 μl, in which 2 μl of purified DNA were used together with 10 μl of Luna Universal qPCR Master Mix (New England Biolabs Inc., Ipswich, MA, USA), 7 μl of nuclease-free water, and 0.5 μl each of 10 μM forward and reverse primer. Each reaction was performed in triplicate. The cycling conditions included an initial denaturation step of 1 min at 95 °C followed by 40 cycles of 95 °C for 15 s and 60 °C for 30 s. No-template controls (2 μl of nuclease-free water instead of DNA extract) and an internal calibrator control for each strain (2 μl of each strain’s purified DNA at a concentration of 10 ng/μl with a known Ct value, ranging between 12 and 14 Ct, used to account for possible variations between plate runs) were performed with each batch of samples tested. The uidA gene encoding a β-glucuronidase specific for E. coli was included in the qPCR assay as housekeeping gene .
Specificity and efficiency testing
Primer pairs were individually tested with their respective target strain first by simplex PCR and afterwards together in multiplexes. For this, each strain was streaked onto a Gassner agar plate (Sifin, Berlin, Germany) containing ceftiofur at 4 μg/ml (ceftiofur hydrochloride, VETRANAL®, St. Louis, MO, USA) and allowed to grow overnight at 37 °C. After verifying a pure culture, a single colony of each strain was then picked and used to inoculate 10 ml of fresh LB medium. Liquid cultures were incubated overnight at 37 °C. After 18 h, 3 ml of the overnight culture were used to isolate genomic DNA (peqGOLD Bacterial DNA Mini Kit, Peqlab, Erlangen, Germany). A Nanodrop microvolume spectrophotometer (NanoDrop One, Thermo Fisher) was used to monitor DNA quality and concentration. Serial dilutions (100 ng/μl to 0.1 ng/μl) were made from the DNA solutions and used as templates to test primer specificity. Culture mixtures of different strains were also used to ensure no cross-detection occurred. Faecal spiking was performed to verify that the primers were sufficiently specific to detect individual strains in the background of an intestinal microbiome. For this, pure cultures of each strain were allowed to grow overnight at 37 °C in liquid culture. A hundred microliters of the overnight culture were used to inoculate 10 millilitres of fresh media. When the culture reached an OD600 of 0.5, each individual culture was adjusted to a bacterial concentration of 5.88 × 108 cells, based on previously collected data, and 10 microliters of the concentration-adjusted culture were used to spike 1 g of a porcine faecal sample. Single-strain spiking and spiking with mixtures containing 2–10 of the experimental bacterial strains were performed. Afterwards, DNA was extracted using the Quick-DNA Faecal/Soil Microbe Miniprep Kit (Zymo Research, Irvine, CA, USA).
The efficiency of the qPCR primers was calculated following the recommendations published by Svec et al. . Ten-fold dilutions ranging from 10 ng down to 0.0001 ng of DNA were tested in triplicate. The mean average of the triplicates was plotted on a logarithmic scale along with the corresponding template concentrations. A linear regression curve was applied to the data points, and the slope of the trend line was calculated. Finally, efficiency was calculated by using the equation: E = − 1 + 10(− 1/slope). All primer pairs tested showed efficiency values between 90 and 95%. Non-template controls did not show any amplification, and internal calibrator values always stayed within the determined Ct value range of 12–14 cycles depending on the respective calibrator used.
Tracing of bacterial strains after cocktail infection of piglets
For in vivo analysis, pigs were inoculated with a bacterial cocktail containing 17 different strains (Table 1). The animal experiment was approved by the competent authority (State Office for Agriculture, Food Safety and Fisheries of Mecklenburg-Western Pomerania, Rostock, Germany, reference no. 7221.3–1-034/19). Eight German landrace pigs, 42–45 days old and healthy as per veterinary guidelines, were purchased from a conventional pig farm (bhzp Garlitz, Langenheide, Germany) and housed in an environmentally controlled animal facility at the Friedrich-Loeffler-Institut (FLI) on the Isle of Riems, Greifswald. The animals adapted to the environmental conditions for 3 weeks prior to experimental infection. Meanwhile, faecal samples were collected to determine the resistance status of the coliform bacterial population in the intestinal tract of the pigs. Some samples tested positive for ceftiofur-resistance, but all samples tested negative for ceftiofur/rifampicin double-resistant bacteria. All inoculation strains were grown individually on Gassner plates containing 4 μg/ml ceftiofur and 50 μl/ml rifampicin and stored at 4 °C. The bacterial cocktail was prepared before inoculation, by mixing liquid cultures of all 17 strains at equal numbers (5.88 × 108 cells per strain) in order to reach a total of 1010 bacteria per inoculation dose. Mixtures were gently centrifuged, the media removed, and the bacterial pellets resuspended in 10 ml of physiological saline containing 10% sodium bicarbonate to buffer stomach acid. After re-suspension, cooled individual doses were immediately transported to the animal facility. For inoculation, all animals were lightly sedated intramuscularly with azaperon (Stresnil®, Elanco, Greenfield, IN, USA) using 0.5 ml / 20 kg of body weight. The inoculation of the strain cocktail was performed intra-gastrically using a gastric tube (B. Braun, Melsungen, Germany). The animals recovered quickly and were fed immediately after the procedure. Post-inoculation, clinical observation of the animals was performed once per day during the entire experiment. Rectal swabs, in addition to faecal samples from the pen, were collected daily from day 1 to 14 post-infection, and every 2 days from day 15 until day 56 at the end of the experiment.
Rectal swabs were suspended in 1 ml of LB medium and allowed to rest at 37 °C for 30 min. The swab wash-offs were serially diluted from 10− 1 to 10− 4 and used for plate spotting on Gassner agar plates containing either no antibiotics or ceftiofur (4 μg/ml) and rifampicin (50 μg/ml). For spotting, 10 μl droplets of the 10− 1 to 10− 4 dilutions were gently spread on each plate in duplicate and plates were left open inside the bench for 1–2 min to allow excess media to be absorbed by the agar. After overnight aerobic incubation at 37 °C, colonies were counted in each droplet. In addition, 100 μl of the suspension from the rectal swabs were plated on Gassner agar plates containing ceftiofur (4 μg/ml) and rifampicin (50 μg/ml). After overnight incubation at 37 °C, plates were washed off using 2 ml of LB, and the suspension was used to isolate DNA with a commercial kit (peqGOLD Bacterial DNA Mini Kit, Peqlab, Erlangen, Germany).
To collect colonic content and tissue, four animals each were euthanized by intravenous administration of Pentobarbital (Release® 500 mg/ml, WDT) on days 43 p.i. and 56 p.i., respectively. At post-mortem examination, the intestinal tract of each animal was removed and the colon section was separated by a double ligation. After opening, approximately 50 ml of the content was collected and a large piece of intestinal tissue (approximately 2–3 cm) sampled and gently washed to remove any remaining content.
One gram of content was weighed, diluted in 9 ml of LB medium and allowed to rest at 37 °C for 30 min. One gram of tissue was weighed, finely chopped, suspended in 9 ml of LB medium and allowed to rest at 37 °C for 30 min. One hundred microliters of each suspension were serially diluted from 10− 1 to 10− 4 and used for plate spotting (see detailed description above) on Gassner agar plates containing either no antibiotics or ceftiofur (4 μg/ml) and rifampicin (50 μg/ml). After overnight incubation at 37 °C, colonies were counted.
One millilitre of the initial dilution (10-1) was used to isolate DNA with a commercial kit (peqGOLD Bacterial DNA Mini Kit, Peqlab, Erlangen, Germany) prior to enrichment.
Sample enrichment was performed by adding rifampicin to the initial content suspensions, which was further incubated overnight at 37 °C. The next day, 100 μl of the overnight suspension were plated on Gassner agar plates containing ceftiofur (4 μg/ml) and rifampicin (50 μg/ml). After overnight aerobic incubation at 37 °C, plates were washed off using 2 ml of LB. One millilitre of the suspension was used for generating glycerol stocks and the remaining millilitre was used to isolate DNA with a commercial kit (peqGOLD Bacterial DNA Mini Kit, Peqlab, Erlangen, Germany).
Statistical analyses were performed with GraphPad Prism Software (GraphPad Prism version 9.0.2 for Windows, GraphPad Software, San Diego, California USA, www.graphpad.com). The one-way ANOVA on ranks (Kruskal-Wallis) test for multiple comparisons was used to determine the significance of differences between ∆Ct (∆Ct = Ct (gene of interest) – Ct (housekeeping gene)) values from strains. Values with p ≤ 0.05 were considered significant.
ORFan gene identification
A total of 299 ORFan genes were identified in the whole genome sequences of the 17 experimental strains (Table 5).
The number of strain-specific genes per strain was highly variable, ranging from a maximum of 158 in one strain to zero in two strains. Approximately one-third (86) of the 299 strain-specific genes were classified as plasmid-borne and, therefore, not deemed suitable for strain identification. They were not further analysed as to their ORFan status. Nearly all of the strain-specific genes, that were not potentially plasmid-encoded (213), turned out to be ORFan genes (206). For the three remaining strains that showed no ORFan genes, a comparison within the 17 strains selected for the bacterial cocktail was made. Genes that were unique among the cocktail strains, and showed less than 5 hits with other E. coli in the GenBank search, were then used as strain identifiers. These primers were tested with DNA isolated from faecal samples from clinically healthy (non-inoculated) pigs with no positive matches, and were deemed specific enough to be used during the animal experiment.
PCR primer specificity and multiplex functionality
A primer pair was designed for each experimental strain (Table 2). Only primer pairs that exclusively amplified DNA of the corresponding strain were considered target-specific. Among the 34 initially designed primers, seven pairs failed (41.2%), either by not-amplifying their respective target or by yielding non-specific bands when tested against other strains. Consequently, seven new primers pairs were designed and two pairs again failed (28.5%) for the same reasons. The third pair of primers yielded specific signals for the last two strains. Combining the primer pairs into multiplexes affected neither their specificity nor their ability to identify the matching strain (Fig. 1). No non-specific bands were seen after multiplexing the primers; however, it was observed that the bands of the smallest fragments in the multiplex (smaller than 150 bp) had reduced intensity in those multiplexes containing a larger number of primer pairs (e.g., Cattle multiplex with six pairs of primers and Pig multiplex with five pairs of primers – see Table 4).
During multiplex testing with DNA isolated from pure cultures, the detection limit of the reactions ranged between 1 to 10 ng/μl DNA, which is equivalent to 200 to 2000 genome copies. However, DNA isolated directly from faecal samples spiked with 108 bacteria did not yield positive results, indicating that inhibitory substances or the high endogenous background of bacterial and eukaryotic DNA might be problematic in the aforementioned set up. To overcome this obstacle, a selection step was added for experimental strain detection from faecal samples. Since all experimental strains possessed ESBL genes and were rifampicin resistant, a plating step on Gassner agar containing ceftiofur (4 μg/ml) and rifampicin (50 μg/ml) was added to limit strain isolation to the experimental strains and to eliminate the presence of inhibitory faecal substances. This eliminated growth of other endogenous ESBL bacteria that were already present in the animals pre-inoculation, and facilitated specific detection of experimental strains in a multiplex PCR approach using boiled lysates or DNA prepared from the pooled bacteria.
Monitoring of faecal samples
The suitability of the ORFan approach to qualitatively monitor shedding of different E. coli strains by pigs when inoculated with 17 strains simultaneously was assessed by conducting colony counting of resistant E. coli and detection of strains by classical PCR. All experimental strains were confirmed to be present in the initial inoculation cocktail prepared to be given to the animals (Fig. 2). After inoculation of the animals with the bacterial cocktail, faecal samples were collected for 56 days. Results after the first 24 h post-inoculation were highly variable between animals. After 48 h, 12 of the 17 strains were detected in at least one animal, with a minimum of four strains and a maximum of 12 strains detected in the eight animals. The remaining five strains were not detected at all in faecal matter during the experiment.
Experimental strain counts, i.e., the total number of colonies, which grew without differentiating the single strains, remained high up to day 3–4 post-inoculation, after which they slowly declined. A significant number of strains from the inoculum mixture were detectable at day 3 (Fig. 3), with three animals being positive for four strains, one animal for six strains, one animal for nine strains, two animals for eleven strains, and one animal positive for twelve strains. Inoculated bacteria were shed and identified by PCR up to day 21 p.i., when their counts on selective plates had declined to single-digit numbers per gram of faeces or were absent, requiring an enrichment protocol for further detection of experimental strains. By day 29 p.i., only 5–6 different strains were shed, even after enrichment, and by day 53 p.i., only four experimental strains were detected.
qPCR system results
The qPCR system was used to prove the general suitability of the ORFan approach to compare different E. coli strains quantitatively within complex intestinal porcine microbiomes. To this end, strains that gave positive signals in the PCR Multiplex setup were further analysed via qPCR for their presence in the dissection sample contents. Delta-Ct (ΔCt) values obtained from the qPCRs showed significant differences in quantities between the strains. The strain classified as P3 (39533) showed the lowest ΔCt values, indicating a high presence in the colonic content. This strain’s presence in the samples was significantly different from all other strains tested via qPCR (Fig. 4; Table 6). The strains classified as C2 (R45) and M1 (21225_2#178) showed the highest ΔCt values on average, denoting both strains’ lower presence in the intestinal content of the inoculated animals. Also, M1, together with P2 (IMT28138), showed the lowest degree of significance when compared against the other strains tested (Fig. 4; Table 6).
Current strain identification methods based on DNA fragment patterns, including pulsed-field gel electrophoresis (PFGE), restriction fragment length polymorphisms (RFLP), repetitive sequencing-based PCR (REP-PCR), Enterobacterial Repetitive Intergenic Consensus PCR (ERIC-PCR), and multiple-locus variable-number tandem repeat analysis (MLVA) , cannot be used to characterize mixtures of strains because it is not possible to assign individual bands to their cognate isolate unambiguously. DNA sequencing-based methods cannot resolve individual strains in a mixture except if sequences with differentiating SNPs have been identified for all strains. The analysis would then require prior testing of the animals to ensure the absence of these discriminatory SNPs in the endogenous bacterial population. DNA hybridization-based methods, such as cDNA and oligonucleotide microarrays, permit the detection of individual genes or gene fragments and, consequently, strains, but they also require previous knowledge of strain-specific sequences and have to be individually adapted for each new combination of to-be-detected strains. Other strategies, commonly used in infection models, like introducing artificial selection markers, such as antimicrobial resistance genes or genes coding for fluorescent proteins, are of limited feasibility since they allow only a limited number of strains to be introduced simultaneously .
The ORFan gene targeting approach utilized herein has the advantage of allowing the introduction of multiple strains simultaneously into an experimental setup, including an in vivo animal trial, as demonstrated. As opposed to existing methods for strain identification, which are either extremely time-consuming, expensive or need specific equipment, the ORFan identification system can be implemented relatively fast and is accessible to everyone with standard laboratory equipment and moderate knowledge of bioinformatics. This method can also be flexibly and easily scaled up, because it only requires the identification of a specific ORFan gene for any novel strain to be introduced.
The ORFan gene approach allows the combination with other visualization techniques to expand strain-specific detection to other types of biological sample. Recently, a HiPR-FISH (high phylogenetic resolution - fluorescence in situ hybridization)-based identification technique was described . In the aforementioned study, HiPR-FISH was employed to identify (i) over 1000 E. coli isolates, using artificially introduced barcoded sequences to generate fluorescent probes for individual strain identification and (ii) bacterial genera present in the murine gut microbiome or in human plaque biofilms using 16S rRNA sequences to identify different genera . Using ORFan genes instead of 16S rDNA as targets for such strain-specific fluorescent hybridization probes could generate an alternative and complementary data set on individual strains. Fluorescent detection of ORFan gene presence via hybridization of fixed and permeabilized cells would allow rapid quantification of large samples via fluorescence-activated cell sorting. Another possible application could be specific strain detection in tissue sections from infected animals, which would even permit comparative spatial-temporal resolution of multiple strains in a host coupled with respective niche identification.
The ORFan gene identification workflow pinpointed a large enough number of ORFan genes for most strains to allow the successful selection of specific primers. Previous E. coli phylogenetic studies unveiled that ORFan genes could compose approximately 1% of the bacteria’s total core genome, with an increase of approximately 26 genes per new genome sequenced . The average number of 18 specific genes per strain was lower for the set of strains selected here for the animal test, and the respective total number of strain-specific genes was lower for 15 of the 17 strains. An additional complication is that these genes could also be located on plasmids and, thus, subjected to horizontal transmission, disqualifying them as stable markers. In our study, approximately one-third of the strain-specific genes were classified as plasmid-derived, but all affected strains, except one, harboured alternate ORFans for which a specific detection system could be designed. Primer implementation was an iterative process, requiring up to three rounds of sequence selection until all 17 strains could be explicitly identified in multiplex PCR assays. Each round required up to 2 weeks of time spent between ORFan gene selection and successful testing of the specific primer pairs. For the three strains, for which no ORFan genes were identified, an additional search of their accessory genomes was performed to identify unique genes within the sequence context of the 17 experimental strains. ORFan genes found using this approach were also subjected to the workflow presented in the “Materials and Methods” section to assess exclusivity. As expected, most of the newly identified genes did have more than 10 hits in the GenBank search. For this reason, only nucleotide stretches within a selected ORFan gene that showed a higher degree of sequence variation between the experimental strain’s sequence and the sequences available in GenBank, were used to generate primers. Furthermore, designated cocktail strains were artificially rendered rifampicin-resistant to distinguish them from cephalosporin-resistant E. coli in the endogenous animal microbiota. The combination of these measures was successful for expanding the limits of the ORFan approach, at least for the list of bacteria under study and in the group of animals used.
An aliquot from the strain cocktail used to infect the animals was immediately stored at −80 °C after preparation. DNA was extracted, and the four multiplexes performed to corroborate the experimental strains’ presence. All 17 strains were detected, confirming their presence in the inoculation cocktail. Their respective band intensities were also similar, indicating that each strain had been added in a similar quantity to the cocktail. Among the 17 strains used to infect the pigs, five strains were not detected in the faecal samples throughout the experiment and were also not detected in the intestinal content at necropsy. Eight strains were only detectable intermittently during the entire experiment, with all of them displaying higher prevalence during the first 14 days p.i.. The remaining four strains were consistently detected throughout the entire experiment. The presence or absence of the 12 successfully detected strains was closely monitored during the entire experiment via the multiplex PCRs, demonstrating that the detection system can be used to follow dynamic changes in a strain’s presence.
All qPCR primers tested showed high efficiency with values ranging from 90 to 95%. Several tests with various samples corroborated that the primers were indeed sufficiently specific to allow the direct use of faecal DNA to detect experimental bacteria from the mixed background of total bacteria. Strain presence in gut content from the animals’ colon, which had already been demonstrated in the multiplex PCRs, was also seen with qPCR primers. The qPCR data indicated colonization differences between experimental strains. Based on the qPCR results, strains with high colonization capacity were also easily re-isolated from the faecal samples, demonstrating that the qPCR results are accurate at monitoring specific strain abundancies in the colonic samples.
Next-generation metagenome sequencing provides vital information on microbial populations and genetic diversity at all taxonomic levels. A fast and easy but robust and reliable method for individual strain identification based on information derived from whole-genome sequencing has yet to be described. The usage of ORFan genes, specific for individual strains, could be such a valuable tool, as it allows to develop precise PCR markers for tracing the strains in complex mixed-culture experiments. The method has the potential to be applied in multiple ways to foster our understanding of, for example, the population dynamics of closely related strains of a pathogen. A probiotic or any other strain of interest can be evaluated as to its colonization ability, general strain fitness, or zoonotic risk and might guide the development of intervention strategies . This method could also potentially be used as a fast and powerful tool for back-tracking the identification of specific pathogens in the event of an outbreak. At present, whole-genome-sequencing is an essential part of outbreak investigations. After acquiring the sequence data of the suspected outbreak agent, the strain detection method presented here could be used for rapid identification of a specific outbreak strain following identification of ORFan genes unique to the specific outbreak agent and the design of strain-specific primers. Without the need to amplify a whole set of virulence markers characteristic of an outbreak strain or to isolate the pathogen from multiple samples, ORFan genes may be used for specific-strain identification via PCR, allowing the rapid pre-screening of many different samples to narrow down the potential sources that could have served as origin or as potential transmission route of the pathogen in the outbreak. In the ensuing second step of outbreak analysis, PCR-positive samples would be subjected to classical approaches involving strain isolation and characterization for unambiguous identification of the culprit strain. This would allow to perform large-scale surveys of many different samples to identify potential outbreak sources and clusters. Similar approaches have been used in the past, such as the one described by Bielaszewska et al. , where an outbreak strain was successfully identified by detecting a specific virulence gene profile. A similar approach was described by Kiel et al. , where two pipelines were simultaneously run, comparing Shiga-toxin expressing E. coli (STEC) genomes versus control genomes and an STEC core proteome versus control proteomes. Lineage- and serotype-specific genes were identified this way and used for monitoring specific STEC strains.
The method described in this manuscript using single unique genes to identify specific strains proved easy to implement and very reliable in identifying and following individual strains and their dynamics in an in vivo animal model of experimental infection, thus representing a fast, inexpensive and reliable alternative to more costly and laborious identification methods.
Availability of data and materials
The datasets generated and/or analysed during the current study are available in the NCBI repository, under the title “The impact of Host restriction of Escherichia coli on Transmission dynamics and spread of antimicrobial Resistance”, BioProject number PRJNA739205.
Delta threshold cycle
Basic Local Alignment Search Tool
Enterobacterial Repetitive Intergenic Consensus PCR
Fluorescence in situ hybridization
Multiple-locus variable-number tandem repeat analysis
Open reading frames
Polymerase chain reaction
Pulsed-field gel electrophoresis
Quantitative polymerase chain reaction
Repetitive sequencing-based PCR
Restriction fragment length polymorphisms
Shiga-toxin expressing E. coli
Witkowska E, Korsak D, Kowalska A, Janeczek A, Kamińska A. Strain-level typing and identification of bacteria – a novel approach for SERS active plasmonic nanostructures. Anal Bioanal Chem. 2018;410:5019–31.
Janda JM, Abbott SL. Bacterial identification for publication: when is enough enough? J Clin Microbiol. 2002;40:1887-91.
Li W, Raoult D, Fournier P-E. Bacterial strain typing in the genomic era. FEMS Microbiol Rev. 2009;33:892–916.
Schork NJ, Fallin D, Lanchbury JS. Single nucleotide polymorphisms and the future of genetic epidemiology. Clin Genet. 2000;58:250–64.
Lupski JR, Weinstock GM. Short, interspersed repetitive DNA sequences in prokaryotic genomes. J Bacteriol. 1992;174:4525–9.
van Belkum A, Scherer S, van Alphen L, Verbrugh H. Short-sequence DNA repeats in prokaryotic genomes. Microbiol Mol Biol Rev. 1998;62:275–93.
Denamur E, Clermont O, Bonacorsi S, Gordon D. The population genetics of pathogenic Escherichia coli. Nat Rev Microbiol. 2021;19:37–54.
Touchon M, Hoede C, Tenaillon O, Barbe V, Baeriswyl S, Bidet P, et al. Organised genome dynamics in the Escherichia coli species results in highly diverse adaptive paths. PLoS Genet. 2009;5:e1000344.
Land M, Hauser L, Jun S, Nookaew I, Leuze MR, Ahn T, et al. Insights from 20 years of bacterial genome sequencing. Funct Integr Genomics. 2015;15:141–61.
Gordienko EN, Kazanov MD, Gelfand MS. Evolution of pan-genomes of Escherichia coli, Shigella spp., and Salmonella enterica. J Bacteriol. 2013;195:2786–92.
Wilson GA, Bertrand N, Patel Y, Hughes JB, Feil EJ, Field D. Orphans as taxonomically restricted and ecologically important genes. Microbiology. 2005;151:2499–501.
Yu G, Stoltzfus A. Population diversity of ORFan genes in Escherichia coli. Genome Biol Evol. 2012;4:1176–87.
Medini D, Donati C, Tettelin H, Masignani V, Rappuoli R. The microbial pan-genome. Curr Opin Genet Dev. 2005;15:589–94.
Medini D, Serruto D, Parkhill J, Relman DA, Donati C, Moxon R, et al. Microbiology in the post-genomic era. Nat Rev Microbiol. 2008;6:419–30.
Lukjancenko O, Wassenaar TM, Ussery DW. Comparison of 61 sequenced Escherichia coli genomes. Microb Ecol. 2010;60:708–20.
Pallen MJ, Wren BW. Bacterial pathogenomics. Nature. 2007;449:835–42.
Raskin DM, Seshadri R, Pukatzki SU, Mekalanos JJ. Bacterial genomics and pathogen evolution. Cell. 2006;124:703–14.
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19(5):455-77. https://doi.org/10.1089/cmb.2012.0021.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29(8):1072-5. https://doi.org/10.1093/bioinformatics/btt086.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.
Page AJ, Cummins CA, Hunt M, Wong VK, Reuter S, Holden MTG, et al. Roary: rapid large-scale prokaryote pan-genome analysis. Bioinformatics. 2015;31:3691–3.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:1–11.
Zangenberg G, Saiki RK, Reynolds R. Multiplex PCR: Optimization guidelines. In: Innis M, Gelfand D, Sninsky J, editors. PCR applications. 1st ed. San Diego: Academic Press; 1999. p. 73–94.
Farnleitner AH, Kreuzinger N, Kavka GG, Grillenberger S, Rath J, Mach RL. Simultaneous detection and differentiation of Escherichia coli populations from environmental freshwaters by means of sequence variations in a fragment of the β-d-glucuronidase gene. Appl Environ Microbiol. 2000;66:1340–6.
Svec D, Tichopad A, Novosadova V, Pfaffl MW, Kubista M. How good is a PCR efficiency estimate: recommendations for precise and robust qPCR efficiency assessments. Biomol Detect Quantif. 2015;3:9–16.
Hamm K, Barth SA, Stalb S, Geue L, Liebler-Tenorio E, Teifke JP, et al. Experimental infection of calves with Escherichia coli O104:H4 outbreak strain. Sci Rep. 2016;6:32812.
Shi H, Shi Q, Grodner B, Lenz JS, Zipfel W, Brito I, et al. Highly multiplexed spatial mapping of microbial communities. Nature. 2020;588:676–81.
Touchon M, Perrin A, Moura de Sousa JA, Vangchhia B, Burn S, O’Brien CL, et al. Phylogenetic background and habitat drive the genetic diversification of Escherichia coli. PLoS Genet. 2020;16:e1008866.
Bielaszewska M, Mellmann A, Zhang W, Köck R, Fruth A, Bauwens A, et al. Characterisation of the Escherichia coli strain associated with an outbreak of haemolytic uraemic syndrome in Germany, 2011: a microbiological study. Lancet Infect Dis. 2011;11:671–6.
Kiel M, Sagory-Zalkind P, Miganeh C, Stork C, Leimbach A, Sekse C, et al. Identification of novel biomarkers for priority serotypes of Shiga toxin-producing Escherichia coli and the development of multiplex PCR for their detection. Front Microbiol. 2018;9:1321.
HECTOR consortium members:
Marta Ferrandis-Vila1, Sumeet K. Tiwari2, Boas van der Putten3,4, Nguyen V. Trung5, Rik Oldenkamp3, Martin Bootsma6, Sebastien Matamoros4, Hoa T. Ngo5,7, Julio Alvarez8,9, Jennifer M. Ritchie10, Amanda Fivian-Hughes10, Angelika Fruth11, Joy Leng12, Roberto M. La Ragione12, Maria Ugarte-Ruiz8,9, Astrid Bethe13,14, Stefan Schwarz13,14, Torsten Semmler2, Constance Schultsz3,4, Christian Menge1 & Christian Berens1.
1Friedrich-Loeffler-Institut – Federal Research Institute for Animal Health, Institute of Molecular Pathogenesis, Jena (Germany); 2Robert Koch Institute, Microbial Genomics, Berlin (Germany); 3Amsterdam UMC, Department of Global Health, Amsterdam (the Netherlands); 4Department of Medical Microbiology, Amsterdam (the Netherlands); 5Oxford University Clinical Research Unit (Vietnam); 6Mathematical Institute and UMCU, Utrecht University, Utrecht (the Netherlands); 7Nuffield Department of Medicine, University of Oxford, Oxford OX3 7BN (United Kingdom); 8VISAVET, Health Surveillance Center, Complutense University of Madrid (Spain); 9Animal Health Department, Veterinary Faculty, Complutense University of Madrid (Spain); 10University of Surrey, School of Biosciences and Medicine, Guildford, Surrey (United Kingdom); 11Robert Koch Institute, Enteropathogenic Bacteria and Legionella, Wernigerode (Germany); 12University of Surrey, School of Veterinary Medicine, Guildford, Surrey (United Kingdom); 13Institute for Microbiology and Epizootics, Freie Universität Berlin, Berlin (Germany); 14Veterinary Centre for Resistance Research (TZR), Freie Universität Berlin, Berlin (Germany).
Open Access funding enabled and organized by Projekt DEAL. The HECTOR research project was supported under the framework of the JPIAMR - Joint Programming Initiative on Antimicrobial Resistance – through the 3rd joint call, thanks to the generous funding by the Netherlands Organization for Health Research and Development (ZonMw, grant number 547001012), the Federal Ministry of Education and Research (BMBF/DLR grant numbers 01KI1703A, 01KI1703B and 01KI703C), the State Research Agency (AEI) of the Ministry of Science, Innovation and Universities (MINECO, grant number PCIN-2016-096), and the Medical Research Council (MRC, grant number MR/R002762/1). The funders were not involved in (i) the design of the study, (ii) the collection, analysis, and interpretation of data and (iii) in writing the manuscript.
Ethics approval and consent to participate
The experimental animal study involving pigs was approved by the competent authority (State Office for Agriculture, Food Safety and Fisheries of Mecklenburg-Western Pomerania, Rostock, Germany) under the reference no. 7221.3–1-034/19–1.
All experiments involving animals were performed according to the guidelines for the care and use of experimental animals established by the State Office for Agriculture, Food Safety and Fisheries of Mecklenburg-Western Pomerania.
The study was carried out in compliance with the ARRIVE guidelines.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Ferrandis-Vila, M., Tiwari, S.K., Mamerow, S. et al. Using unique ORFan genes as strain-specific identifiers for Escherichia coli. BMC Microbiol 22, 135 (2022). https://doi.org/10.1186/s12866-022-02508-y
- ORFan gene
- Strain identification
- Multiplex PCR
- E. coli
- Animal experiment