Nanopore adaptive sampling effectively enriches bacterial plasmids

Bacterial plasmids play a major role in the spread of antibiotic resistance genes. However, their characterization via DNA sequencing suffers from the low abundance of plasmid DNA in those samples. Although sample preparation methods can enrich the proportion of plasmid DNA before sequencing, these methods are expensive and laborious, and they might introduce a bias by enriching only for specific plasmid DNA sequences. Nanopore adaptive sampling could overcome these issues by rejecting uninteresting DNA molecules during the sequencing process. In this study, we assess the application of adaptive sampling for the enrichment of low-abundant plasmids in known bacterial isolates using two different adaptive sampling tools. We show that a significant enrichment can be achieved even on expired flow cells. By applying adaptive sampling, we also improve the quality of de novo plasmid assemblies and reduce the sequencing time. However, our experiments also highlight issues with adaptive sampling if target and non-target sequences span similar regions.

IMPORTANCE

Antimicrobial resistance causes millions of deaths every year. Mobile genetic elements like bacterial plasmids are key drivers for the dissemination of antimicrobial resistance genes. This makes the characterization of plasmids via DNA sequencing an important tool for clinical microbiologists. Since plasmids are often underrepresented in bacterial samples, plasmid sequencing can be challenging and laborious. To accelerate the sequencing process, we evaluate nanopore adaptive sampling as an in silico method for the enrichment of low-abundant plasmids. Our results show the potential of this cost-efficient method for future plasmid research but also indicate issues that arise from using reference sequences.

Endothelial Indoleamine-2,3-Dioxygenase-1 is not Critically Involved in Regulating Antitumor Immunity in the Central Nervous System

The vascular niche of malignant gliomas is a key compartment that shapes the immunosuppressive brain tumor microenvironment (TME). The blood-brain-barrier (BBB) consisting of specialized endothelial cells (ECs) and perivascular cells forms a tight anatomical and functional barrier critically controlling transmigration and effector function of immune cells. During neuroinflammation and tumor progression, the metabolism of the essential amino acid tryptophan (Trp) to metabolites such as kynurenine has long been identified as an important metabolic pathway suppressing immune responses. Previous studies have demonstrated that indoleamine-2,3-dioxygenase-1 (IDO1), a key rate-limiting enzyme in tryptophan catabolism, is expressed within the TME of high-grade gliomas. Here, we investigate the role of endothelial IDO1 (eIDO1) expression for brain tumor immunity. Single-cell RNA sequencing data revealed that in human glioma tissue, IDO1 is predominantly expressed by activated ECs showing a JAK/STAT signaling pathway-related CXCL11⁺ gene expression signature. In a syngeneic experimental glioma model, eIDO1 is induced by low-dose tumor irradiation. However, cell type-specific ablation of eIDO1 in experimental gliomas did not alter frequency and phenotype of tumor-infiltrating T cells nor tumor growth. Taken together these data argue against a dominant role of eIDO1 for brain tumor immunity.

How to survive pig farming: Mechanism of SCCmec element deletion and metabolic stress adaptation in livestock-associated MRSA

Previous research on methicillin susceptible Staphylococcus aureus (MSSA) belonging to livestock-associated (LA-) sequence type (ST) 398, isolated from pigs and their local surroundings, indicated that differences between these MSSA and their methicillin resistant predecessors (MRSA) are often limited to the absence of the staphylococcal cassette chromosome mec (SCCmec) and few single nucleotide polymorphisms. So far, our understanding on how LA-MRSA endure the environmental conditions associated with pig-farming as well as the putative impact of this particular environment on the mobilisation of SCCmec elements is limited. Thus, we performed in-depth genomic and transcriptomic analyses using the LA-MRSA ST398 strain IMT38951 and its methicillin susceptible descendant. We identified a mosaic-structured SCCmec region including a putative replicative SCCmecVc which is absent from the MSSA chromosome through homologous recombination. Based on our data, such events occur between short repetitive sequences identified within and adjacent to two distinct alleles of the large cassette recombinase genes C (ccrC). We further evaluated the global transcriptomic response of MRSA ST398 to particular pig-farm associated conditions, i.e., contact with host proteins (porcine serum) and a high ammonia concentration. Differential expression of global regulators involved in stress response control were identified, i.e., ammonia-induced alternative sigma factor B-depending activation of genes for the alkaline shock protein 23, the heat shock response and the accessory gene regulator (agr)-controlled transcription of virulence factors. Exposure to serum transiently induced the transcription of distinct virulence factor encoding genes. Transcription of genes reported for mediating the loss of methicillin resistance, especially ccrC, was not significantly different compared to the unchallenged controls. We concluded that, from an evolutionary perspective, bacteria may save energy by incidentally dismissing a fully replicative SCCmec element in contrast to the induction of ccr genes on a population scale. Since the genomic SCCmec integration site is a hot-spot of recombination, occasional losses of elements of 16 kb size may restore capacities for the uptake of foreign genetic material. Subsequent spread of resistance, on the other hand, might depend on the autonomous replication machinery of the deleted SCCmec elements that probably enhance chances for reintegration of SCCmec into susceptible genomes by mere multiplication.

Transmission pathways of campylobacter spp. at broiler farms and their environment in Brandenburg, Germany

Broiler meat is widely known as an important source of foodborne Campylobacter jejuni and Campylobacter coli infections in humans. In this study, we thoroughly investigated transmission pathways that may contribute to possible Campylobacter contamination inside and outside broiler houses. For this purpose we carried out a comprehensive longitudinal sampling approach, using a semi-quantitative cultivation method to identify and quantify transmissions and reservoirs of Campylobacter spp.. Three german broiler farms in Brandenburg and their surrounding areas were intensively sampled, from April 2018 until September 2020. Consecutive fattening cycles and intervening downtimes after cleaning and disinfection were systematically sampled in summer and winter. To display the potential phylogeny of barn and environmental isolates, whole genome sequencing (WGS) and bioinformatic analyses were performed. Results obtained in this study showed very high Campylobacter prevalence in 51/76 pooled feces (67.1%) and 49/76 boot swabs (64.5%). Average counts between 6.4 to 8.36 log10MPN/g were detected in pooled feces. In addition, levels of 4.7 and 4.1 log10MPN/g were detected in boot swabs and litter, respectively. Samples from the barn interior showed mean Campyloacter values in swabs from drinkers 2.6 log10MPN/g, walls 2.0 log10MPN/g, troughs 1.7 log10MPN/g, boards 1.6 log10MPN/g, ventilations 0.9 log10MPN/g and 0.7 log10MPN/g for air samples. However, Campylobacter was detected only in 7/456 (1.5%) of the environmental samples (water bodies, puddles or water-filled wheel tracks; average of 0.6 log10MPN/g). Furthermore, WGS showed recurring Campylobacter genotypes over several consecutive fattening periods, indicating that Campylobacter genotypes persist in the environment during downtime periods. However, after cleaning and disinfection of the barns, we were unable to identify potential sources in the broiler houses. Interestingly, alternating Campylobacter genotypes were observed after each fattening period, also indicating sources of contamination from the wider environment outside the farm. Therefore, the results of this study suggest that a potential risk of Campylobacter transmission may originate from present environmental sources (litter and water reservoirs). However, the sources of Campylobacter transmission may vary depending on the operation and farm environmental conditions.

Phenotypic and genetic analyses of two Campylobacter fetus isolates from a patient with relapsed prosthetic valve endocarditis

Campylobacter fetus can cause intestinal and systemic disease in humans and are well established veterinary and economic pathogens. We report the complete genomic sequences of two C. fetus subsp. fetus (Cff) isolates recovered in 2017 (CITCf01) and 2018 (CITCf02) from a case of recurrent prosthetic valve endocarditis. Both were capable of growth aerobically. Their genomes were found to be highly conserved and syntenic with 99.97% average nucleotide identity (ANI) while differences in their respective sap loci defined the temporal separation of their genomes. Based on core genome phylogeny and ANI of 83 Cff genomes belonging to the previously described human-associated Cff lineage, CITCf01 and CITCf02 grouped in a clade of eleven sequence type (ST)3 Cff (including the Cff type strain NCTC 10842T). CITCf01 and CITCf02 were marked for their lack of unique genomic features when compared to isolates within the subspecies and the type strain in particular. We identified point mutations in oxidative stress response genes, among others, that may contribute to aerobiosis. We report a case of Cff causing relapsed prosthetic valve endocarditis and we highlight the sap island as a polymorphic site within the genetically stable ST3 lineage, central to pathogenicity.

Characterization of Campylobacter spp. Strains Isolated From Wild Birds in Turkey

Turkey is an important stopover site for many migrating birds between Europe, Asia and Africa. Campylobacter spp. are frequently found in wildlife, in particular waterfowl, and distinct strains are disseminated within this reservoir. In this study, 183 wild birds of hunting areas in Turkey were collected and thermophilic Campylobacter spp. from cloacal swabs were isolated at a prevalence of 5.2% from song thrushes (6/116) and 93% from Eurasian coots (41/44). After PCR species differentiation and flaA restriction profiles determination, C. jejuni and C. coli strains were further investigated by whole genome sequencing. PCR target amplification of the ceuE gene, commonly used for C. coli species-identification was inefficient and even hampered in one isolate. A close look on the ceuE sequence revealed that various mismatches in the ceuE oligo annealing sites caused less efficient diagnostic detection. All C. coli isolates belonged to the environmental clade II and clade III, for which thirty-six novel MLST types were identified. Further single nucleotide polymorphism (SNP) analysis showed a high genomic divergence between the C. coli isolates. High variability was also implicated for putative plasmid-located genes detected in 51% of the C. coli isolates. Distinct gene variants in clades II and III C. coli were identified by a k-mer analysis. After substracting k-mers in common with C. coli clade I database, 11 and 35 distinct genes were identified in clades II and III isolates, mainly involved in surface structures and modifications as well as signal transduction, suggesting niche adaptation of C. coli strains in wild birds. All strains were susceptible against (fluoro-)quinolones, erythromycin, tetracycline, gentamicin and only one isolate was resistant against streptomycin, suggesting that the sensitive phenotype was due to absence of selective pressure and niche separation in wild birds in Turkey. We conclude that Campylobacter spp. isolates from wildlife and environmental sources are still scarce in the databases and that there is a need for more studies on thermophilic Campylobacter spp. from different places all over the world in order to complement our understanding on dissemination and adaptation to distinct niches of this global food-borne pathogen.

SCORE: Smart Consensus Of RNA Expression-a consensus tool for detecting differentially expressed genes in bacteria

SUMMARY

RNA-sequencing (RNA-Seq) is the current method of choice for studying bacterial transcriptomes. To date, many computational pipelines have been developed to predict differentially expressed genes from RNA-Seq data, but no gold-standard has been widely accepted. We present the Snakemake-based tool Smart Consensus Of RNA Expression (SCORE) which uses a consensus approach founded on a selection of well-established tools for differential gene expression analysis. This allows SCORE to increase the overall prediction accuracy and to merge varying results into a single, human-readable output. SCORE performs all steps for the analysis of bacterial RNA-Seq data, from read preprocessing to the overrepresentation analysis of significantly associated ontologies. Development of consensus approaches like SCORE will help to streamline future RNA-Seq workflows and will fundamentally contribute to the creation of new gold-standards for the analysis of these types of data.

AVAILABILITY AND IMPLEMENTATION

https://github.com/SiWolf/SCORE.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Frequency, Local Dynamics, and Genomic Characteristics of ESBL-Producing Escherichia coli Isolated From Specimens of Hospitalized Horses

Previous research identified veterinary clinics as hotspots with respect to accumulation and spread of multidrug resistant extended-spectrum β-lactamase (ESBL)-producing Escherichia coli (EC). Therefore, promoting the prudent use of antibiotics to decrease selective pressure in that particular clinical environment is preferable to enhance biosecurity for animal patients and hospital staff. Accordingly, this study comparatively investigated the impact of two distinct perioperative antibiotic prophylaxis (PAP) regimens (short-term versus prolonged) on ESBL-EC carriage of horses subjected to colic surgery. While all horses received a combination of penicillin/gentamicin (P/G) as PAP, they were assigned to either the “single-shot group” (SSG) or the conventional “5-day group” (5DG). Fecal samples collected on arrival (t₀), on the 3rd (t₁) and on the 10th day after surgery (t₂) were screened for ESBL-EC. All isolates were further investigated using whole genome sequences. In total, 81 of 98 horses met the inclusion criteria for this study. ESBL-EC identified in samples available at t₀, t₁ and t₂ were 4.8% (SSG) and 9.7% (5DG), 37% (SSG) and 47.2% (5DG) as well as 55.6% (SSG) and 56.8% (5DG), respectively. Regardless of the P/G PAP regimen, horses were 9.12 times (95% CI 2.79–29.7) more likely to carry ESBL-EC at t₁ compared to t₀ (p < 0.001) and 15.64 times (95% CI 4.57–53.55) more likely to carry ESBL-EC at t₂ compared to t₀ (p < 0.001). ESBL-EC belonging to sequence type (ST) 10, ST86, ST641, and ST410 were the most prevalent lineages, with bla_CTX–M–1 (60%) being the dominant ESBL gene. A close spatio-temporal relationship between isolates sharing a particular ST was revealed by genome analysis, strongly indicating local spread. Consequently, hospitalization itself has a strong impact on ESBL-EC isolation rates in horses, possibly masking differences between distinct PAP regimens. The results of this study reveal accumulation and spread of multi-drug resistant ESBL-EC among horses subjected to colic surgery with different P/G PAP regimens, challenging the local hygiene management system and work-place safety of veterinary staff. Moreover, the predominance of particular ESBL-EC lineages in clinics providing health care for horses needs further investigation.

Population Biology and Comparative Genomics of Campylobacter Species

The zoonotic pathogen Campylobacter is the leading cause for bacterial foodborne infections in humans. Campylobacters are most commonly transmitted via the consumption of undercooked poultry meat or raw milk products. The decreasing costs of whole genome sequencing enabled large genome-based analyses of the evolution and population structure of this pathogen, as well as the development of novel high-throughput molecular typing methods. Here, we review the evolutionary development and the population diversity of the two most clinically relevant Campylobacter species; C. jejuni and C. coli. The state-of-the-art phylogenetic studies showed clustering of C. jejuni lineages into host specialists and generalists with coexisting lifestyles in chicken and livestock-associated hosts, as well as the separation of C. coli isolates of riparian origin (waterfowl, water) from C. coli isolated from clinical and farm-related samples. We will give an overview of recombination between both species and the potential impact of horizontal gene transfer on host adaptation in Campylobacter. Additionally, this review briefly places the current knowledge of the population structure of other Campylobacter species such as C. lari, C. concisus and C. upsaliensis into perspective. We also provide an overview of how molecular typing methods such as multilocus sequence typing (MLST) and whole genome MLST have been used to detect and trace Campylobacter outbreaks along the food chain.

High-Zinc Supplementation of Weaned Piglets Affects Frequencies of Virulence and Bacteriocin Associated Genes Among Intestinal Escherichia coli Populations

To prevent economic losses due to post-weaning diarrhea (PWD) in industrial pig production, zinc (Zn) feed additives have been widely used, especially since awareness has risen that the regular application of antibiotics promotes buildup of antimicrobial resistance in both commensal and pathogenic bacteria. In a previous study on 179 Escherichia coli collected from piglets sacrificed at the end of a Zn feeding trial, including isolates obtained from animals of a high-zinc fed group (HZG) and a corresponding control group (CG), we found that the isolate collection exhibited three different levels of tolerance toward zinc, i.e., the minimal inhibitory concentration (MIC) detected was 128, followed by 256 and 512 μg/ml ZnCl2. We further provided evidence that enhanced zinc tolerance in porcine intestinal E. coli populations is clearly linked to excessive zinc feeding. Here we provide insights about the genomic make-up and phylogenetic background of these 179 E. coli genomes. Bayesian analysis of the population structure (BAPS) revealed a lack of association between the actual zinc tolerance level and a particular phylogenetic E. coli cluster or even branch for both, isolates belonging to the HZG and CG. In addition, detection rates for genes and operons associated with virulence (VAG) and bacteriocins (BAG) were lower in isolates originating from the HZG (41 vs. 65% and 22 vs. 35%, p < 0.001 and p = 0.002, resp.). Strikingly, E. coli harboring genes defining distinct pathotypes associated with intestinal disease, i.e., enterotoxigenic, enteropathogenic, and Shiga toxin-producing E. coli (ETEC, EPEC, and STEC) constituted 1% of the isolates belonging to the HZG but 14% of those from the CG. Notably, these pathotypes were positively associated with enhanced zinc tolerance (512 μg/ml ZnCl2 MIC, p < 0.001). Taken together, zinc excess seems to influence carriage rates of VAGs and BAGs in porcine intestinal E. coli populations, and high-zinc feeding is negatively correlated with enteral pathotype occurrences, which might explain earlier observations concerning the relative increase of Enterobacterales considering the overall intestinal microbiota of piglets during zinc feeding trials while PWD rates have decreased.

Whole genome sequencing reveals extended natural transformation in Campylobacter impacting diagnostics and the pathogens adaptive potential

Campylobacter is the major bacterial agent of human gastroenteritis worldwide and represents a crucial global public health burden. Species differentiation of C. jejuni and C. coli and phylogenetic analysis is challenged by inter-species horizontal gene transfer. Routine real-time PCR on more than 4000 C. jejuni and C. coli field strains identified isolates with ambiguous PCR results for species differentiation, in particular, from the isolation source eggs. K-mer analysis of whole genome sequencing data indicated the presence of C. coli hybrid strains with huge amounts of C. jejuni introgression. Recombination events were distributed over the whole chromosome. MLST typing was impaired, since C. jejuni sequences were also found in six of the seven housekeeping genes. cgMLST suggested that the strains were phylogenetically unrelated. Intriguingly, the strains shared a stress response set of C. jejuni variant genes, with proposed roles in oxidative, osmotic and general stress defence, chromosome maintenance and repair, membrane transport, cell wall and capsular biosynthesis and chemotaxis. The results have practical impact on routine typing and on the understanding of the functional adaption to harsh environments, enabling successful spreading and persistence of Campylobacter.

Comparison of different technologies for the decipherment of the whole genome sequence of Campylobacter jejuni BfR-CA-14430

Background

Campylobacter jejuni is a zoonotic pathogen that infects the human gut through the food chain mainly by consumption of undercooked chicken meat, raw chicken cross-contaminated ready-to-eat food or by raw milk. In the last decades, C. jejuni has increasingly become the most common bacterial cause for food-born infections in high income countries, costing public health systems billions of euros each year. Currently, different whole genome sequencing techniques such as short-read bridge amplification and long-read single molecule real-time sequencing techniques are applied for in-depth analysis of bacterial species, in particular, Illumina MiSeq, PacBio and MinION.

Results

In this study, we analyzed a recently isolated C. jejuni strain from chicken meat by short- and long-read data from Illumina, PacBio and MinION sequencing technologies. For comparability, this strain is used in the German PAC-CAMPY research consortium in several studies, including phenotypic analysis of biofilm formation, natural transformation and in vivo colonization models. The complete assembled genome sequence most likely consists of a chromosome of 1,645,980 bp covering 1665 coding sequences as well as a plasmid sequence with 41,772 bp that encodes for 46 genes. Multilocus sequence typing revealed that the strain belongs to the clonal complex CC-21 (ST-44) which is known to be involved in C. jejuni human infections, including outbreaks. Furthermore, we discovered resistance determinants and a point mutation in the DNA gyrase (gyrA) that render the bacterium resistant against ampicillin, tetracycline and (fluoro-)quinolones.

Conclusion

The comparison of Illumina MiSeq, PacBio and MinION sequencing and analyses with different assembly tools enabled us to reconstruct a complete chromosome as well as a circular plasmid sequence of the C. jejuni strain BfR-CA-14430. Illumina short-read sequencing in combination with either PacBio or MinION can substantially improve the quality of the complete chromosome and epichromosomal elements on the level of mismatches and insertions/deletions, depending on the assembly program used.

Effects of a Four-Week High-Dosage Zinc Oxide Supplemented Diet on Commensal Escherichia coli of Weaned Pigs

Strategies to reduce economic losses associated with post-weaning diarrhea in pig farming include high-level dietary zinc oxide supplementation. However, excessive usage of zinc oxide in the pig production sector was found to be associated with accumulation of multidrug resistant bacteria in these animals, presenting an environmental burden through contaminated manure. Here we report on zinc tolerance among a random selection of intestinal Escherichia coli comprising of different antibiotic resistance phenotypes and sampling sites isolated during a controlled feeding trial from 16 weaned piglets: In total, 179 isolates from "pigs fed with high zinc concentrations" (high zinc group, [HZG]: n = 99) and a corresponding "control group" ([CG]: n = 80) were investigated with regard to zinc tolerance, antimicrobial- and biocide susceptibilities by determining minimum inhibitory concentrations (MICs). In addition, in silico whole genome screening (WGSc) for antibiotic resistance genes (ARGs) as well as biocide- and heavy metal tolerance genes was performed using an in-house BLAST-based pipeline. Overall, porcine E. coli isolates showed three different ZnCl₂ MICs: 128 μg/ml (HZG, 2%; CG, 6%), 256 μg/ml (HZG, 64%; CG, 91%) and 512 μg/ml ZnCl₂ (HZG, 34%, CG, 3%), a unimodal distribution most likely reflecting natural differences in zinc tolerance associated with different genetic lineages. However, a selective impact of the zinc-rich supplemented diet seems to be reasonable, since the linear mixed regression model revealed a statistically significant association between "higher" ZnCl₂ MICs and isolates representing the HZG as well as "lower ZnCl₂ MICs" with isolates of the CG (p = 0.005). None of the zinc chloride MICs was associated with a particular antibiotic-, heavy metal- or biocide- tolerance/resistance phenotype. Isolates expressing the 512 μg/ml MIC were either positive for ARGs conferring resistance to aminoglycosides, tetracycline and sulfamethoxazole-trimethoprim, or harbored no ARGs at all. Moreover, WGSc revealed a ubiquitous presence of zinc homeostasis and - detoxification genes, including zitB, zntA, and pit. In conclusion, we provide evidence that zinc-rich supplementation of pig feed selects for more zinc tolerant E. coli, including isolates harboring ARGs and biocide- and heavy metal tolerance genes - a putative selective advantage considering substances and antibiotics currently used in industrial pork production systems.

The gut microbiome of horses: current research on equine enteral microbiota and future perspectives

Understanding the complex interactions of microbial communities including bacteria, archaea, parasites, viruses and fungi of the gastrointestinal tract (GIT) associated with states of either health or disease is still an expanding research field in both, human and veterinary medicine. GIT disorders and their consequences are among the most important diseases of domesticated Equidae, but current gaps of knowledge hinder adequate progress with respect to disease prevention and microbiome-based interventions. Current literature on enteral microbiomes mirrors a vast data and knowledge imbalance, with only few studies tackling archaea, viruses and eukaryotes compared with those addressing the bacterial components.Until recently, culture-dependent methods were used for the identification and description of compositional changes of enteral microorganisms, limiting the outcome to cultivatable bacteria only. Today, next generation sequencing technologies provide access to the entirety of genes (microbiome) associated with the microorganisms of the equine GIT including the mass of uncultured microbiota, or "microbial dark matter".This review illustrates methods commonly used for enteral microbiome analysis in horses and summarizes key findings reached for bacteria, viruses and fungi so far. Moreover, reasonable possibilities to combine different explorative techniques are described. As a future perspective, knowledge expansion concerning beneficial compositions of microorganisms within the equine GIT creates novel possibilities for early disorder diagnostics as well as innovative therapeutic approaches. In addition, analysis of shotgun metagenomic data enables tracking of certain microorganisms beyond species barriers: transmission events of bacteria including pathogens and opportunists harboring antibiotic resistance factors between different horses but also between humans and horses will reach new levels of depth concerning strain-level distinctions.

SeroBA: rapid high-throughput serotyping of Streptococcus pneumoniae from whole genome sequence data

Streptococcus pneumoniae is responsible for 240 000–460 000 deaths in children under 5 years of age each year. Accurate identification of pneumococcal serotypes is important for tracking the distribution and evolution of serotypes following the introduction of effective vaccines. Recent efforts have been made to infer serotypes directly from genomic data but current software approaches are limited and do not scale well. Here, we introduce a novel method, SeroBA, which uses a k-mer approach. We compare SeroBA against real and simulated data and present results on the concordance and computational performance against a validation dataset, the robustness and scalability when analysing a large dataset, and the impact of varying the depth of coverage on sequence-based serotyping. SeroBA can predict serotypes, by identifying the cps locus, directly from raw whole genome sequencing read data with 98 % concordance using a k-mer-based method, can process 10 000 samples in just over 1 day using a standard server and can call serotypes at a coverage as low as 15–21×. SeroBA is implemented in Python3 and is freely available under an open source GPLv3 licence from: https://github.com/sanger-pathogens/seroba