- Research
- Open access
- Published:
Genomic characterization and virulence gene profiling of Erysipelothrix rhusiopathiae isolated from widespread muskox mortalities in the Canadian Arctic Archipelago
BMC Genomics volume 25, Article number: 691 (2024)
Abstract
Background
Muskoxen are important ecosystem components and provide food, economic opportunities, and cultural well-being for Indigenous communities in the Canadian Arctic. Between 2010 and 2021, Erysipelothrix rhusiopathiae was isolated from carcasses of muskoxen, caribou, a seal, and an Arctic fox during multiple large scale mortality events in the Canadian Arctic Archipelago. A single strain (‘Arctic clone’) of E. rhusiopathiae was associated with the mortalities on Banks, Victoria and Prince Patrick Islands, Northwest Territories and Nunavut, Canada (2010–2017). The objectives of this study were to (i) characterize the genomes of E. rhusiopathiae isolates obtained from more recent muskox mortalities in the Canadian Arctic in 2019 and 2021; (ii) identify and compare common virulence traits associated with the core genome and mobile genetic elements (i.e. pathogenicity islands and prophages) among Arctic clone versus other E. rhusiopathiae genomes; and iii) use pan-genome wide association studies (GWAS) to determine unique genetic contents of the Arctic clone that may encode virulence traits and that could be used for diagnostic purposes.
Results
Phylogenetic analyses revealed that the newly sequenced E. rhusiopathiae isolates from Ellesmere Island, Nunavut (2021) also belong to the Arctic clone. Of 17 virulence genes analysed among 28 Arctic clone isolates, four genes – adhesin, rhusiopathiae surface protein-A (rspA), choline binding protein-B (cbpB) and CDP-glycerol glycerophosphotransferase (tagF) – had amino acid sequence variants unique to this clone when compared to 31 other E. rhusiopathiae genomes. These genes encode proteins that facilitate E. rhusiopathiae to attach to the host endothelial cells and form biofilms. GWAS analyses using Scoary found several unique genes to be overrepresented in the Arctic clone.
Conclusions
The Arctic clone of E. rhusiopathiae was associated with multiple muskox mortalities spanning over a decade and multiple Arctic islands with distances over 1000 km, highlighting the extent of its spatiotemporal spread. This clone possesses unique gene content, as well as amino acid variants in multiple virulence genes that are distinct from the other closely related E. rhusiopathiae isolates. This study establishes an essential foundation on which to investigate whether these differences are correlated with the apparent virulence of this specific clone through in vitro and in vivo studies.
Background
Muskoxen (Ovibos moschatus) are large-bodied herbivores that are well adapted to the Arctic environment. They provide food, economic opportunities, and cultural well-being for Indigenous people in the Canadian Arctic [1]. Accelerated Arctic warming and emerging infectious diseases [2,3,4,5], together with a limited gene pool [6] and other stressors are increasingly threatening the future of this species [7,8,9]. Since 2010, the bacterium Erysipelothrix rhusiopathiae has emerged as a significant cause of mortality in muskoxen in the Canadian Arctic [10]. E. rhusiopathiae is a Gram-positive, opportunistic, facultative intracellular bacterium best known from domestic livestock, where it is generally transmitted via contaminated feed and water [11]. This bacterium is most commonly associated with domestic pigs [12, 13] in which three clinical manifestations are described: a severe acute septicaemic form, a milder subacute urticarial form and a chronic form presenting with arthritis and/or endocarditis [14]. Muskoxen have presented with the severe acute septicaemic form, where animals in good body condition appear to have died rapidly [10]. In wildlife more broadly, E. rhusiopathiae is associated with sporadic individual mortalities and occasional disease outbreaks have been reported in birds [15], rodents [16], and marine mammals [17, 18]. Although E. rhusiopathiae has been circulating in Canadian muskox populations since at least the 1970s [19], it was first identified as a disease-causing agent in 2010. Between 2010 and 2013 a single clonal lineage, the ‘Arctic clone’, part of the phylogenetic clade 3 of E. rhusiopathiae [20], was cultured from muskox carcasses during mass mortality events on Banks and Victoria Islands, Northwest Territories and Nunavut, Canada [10, 21] (Fig. 1). In 2015, this same strain was cultured from a dead seal on Victoria Island, and in summer 2017, it was cultured from carcasses of muskoxen (n = 4), Peary caribou (Rangifer tarandus pearyi) (n = 6) and an Arctic fox (Vulpes lagopus) (n = 1) on Prince Patrick Island, Northwest Territories (NWT), Canada [22]. E. rhusiopathiae was isolated from other sporadic mortalities of muskoxen, caribou and moose in temperate, sub- and mid Arctic Canada and Alaska during this same timeframe; however, whole genome sequencing of the isolates from these sources revealed various genotypes, representative of the wide population diversity of this species [21, 22]. E. rhusiopathiae has been isolated from more recent muskox mortality events in the eastern Canadian Arctic, Nunavut (NU) and western Canadian mainland, Northwest Territories (NWT) and are sequenced and characterized as part of the present study.
Muskox populations affected by the Arctic clone of E. rhusiopathiae experienced substantial population declines during and following the outbreaks [7, 22, 23]. Understanding the drivers of the emergence and apparently high morbidity and mortality rates associated with the Arctic clone is one of the first steps in developing mitigation and management strategies. Considering the epidemiological triad, plausible explanatory factors for the observed mortalities associated with this pathogen include increased host susceptibility due to genetic bottlenecks [6, 24] or increasing stressors [25]; changing environmental factors facilitating bacterial survival and transmission [26, 27]; or increased pathogen virulence. Based on the unusually high mortality and rapidity of deaths in muskoxen infected with the Arctic clone, we hypothesize that this genotype contains unique virulence traits compared to other closely related isolates.
Virulence factors of bacteria influence their pathogenicity, such as the ability to evade the host immune system, leading to infection and potentially disease progression. Virulence factors may be encoded by genes in the chromosomal DNA or by plasmids, and are often carried by different mobile genetic elements such as pathogenicity islands (PAIs), prophages, or transposons [28]. E. rhusiopathiae is a facultative intracellular pathogen that replicates in phagocytic cells, and its pathogenicity is related to cellular attachment, and intracellular invasion and survival [29, 30]. Previously identified PAIs and prophages in E. rhusiopathiae have been shown to harbour transposases, as well as several genes responsible for antimicrobial resistance [31, 32], and could thus play a role in the horizontal transfer of virulence genes.
Although virulence factors of porcine E. rhusiopathiae are generally well characterized, data on those of E. rhusiopathiae isolated from wild animals are lacking. Initial screening for 59 putative virulence genes among E. rhusiopathiae genomes isolated from widespread geographical locations, including the Arctic clone, in 2016 [21], demonstrated that the majority of these genes were found to be present in most genomes (i.e., are core genes). At that time, only gene presence/absence was investigated. However, changes in the DNA sequences encoding these genes could lead to alterations in the amino acid sequences, potentially changing protein structure and function [33]. To understand the unique pathogenic potential associated with the Arctic clone, identifying variations in its virulence genes is important.
The aim of this study was to genomically characterize E. rhusiopathiae isolated from widespread mortalities of muskoxen and other Arctic wildlife, with a focus on genomic differences that might be related to virulence. Our specific objectives were to: (1) determine whether recently recovered E. rhusiopathiae isolates from eastern NU archipelago and western NWT mainland belong to the Arctic clone; (2) compare sequences of virulence genes and assess for specific mobile genetic elements (PAIs and prophages) of interest among Arctic clone and non-Arctic clone isolates; and (3) determine if there is genetic content unique to the Arctic clone.
Methods
Sample collection and processing
Ethical approval for sample collection was obtained from the University of Calgary Animal Care Committee (Protocol number AC18-0093). Wildlife Research Permits were obtained from the Governments of Nunavut (WL2016-058, WL2021-053, WL2019-051, WL2021-048) and the Northwest Territories (WL 500,664). A total of 57 samples from three geographic locations were screened for the presence of E. rhusiopathiae (Fig. 1). (1) Fifty-two samples were collected opportunistically during the investigation of an unusual mortality event of muskoxen on the Fosheim Peninsula, Ellesmere Island, NU, Canada in August 2021. Collections comprised 29 fecal and tissue samples from 10 muskox carcasses (including lung, liver, kidney, heart, spleen, bone marrow, intestine), as well as 22 fecal samples from two apparently healthy muskox herds in the area. (2) Bone marrow samples were collected from the skeletal remains of two dead muskoxen found on Somerset Island, NU, Canada in July 2019 and 2021. (3) One bone marrow sample was collected from a recently dead adult muskox in August 2021, and three tissue samples (lung, liver, kidney) from a dead muskox calf found in September 2021 from the Sahtu region, NWT, Canada. All samples were frozen at -20 ℃ upon collection and shipped within three months from the collection sites. Sample processing and analyses were performed at the Faculty of Veterinary Medicine, University of Calgary, Alberta, Canada.
Isolation and identification of E. rhusiopathiae
For recovery of E. rhusiopathiae, samples were thawed at room temperature (est. 22 ℃). One-gram subsamples of tissue or fecal samples were mechanically homogenized using gentleMACS™ Octo Dissociator (Miltenyi Biotec, Inc., San Diego, CA, USA) in 10 ml brain heart infusion (BHI) broth (BD, Sparks, MD, USA) supplemented with 5% horse serum and incubated for 24 − 48 hours (h) at 37 ℃ with 5% CO2. Subsequently, 300 µl of enriched culture was further incubated for 48 h in 2.7 ml modified BHI (mBHI) supplemented with 5% horse serum, kanamycin (40 µg/ml) (VWR International, Mississauga, Ontario, Canada), neomycin (50 µg/ml) (VWR International) and vancomycin (25 µg/ml) (VWR International), as previously described for the selective culture of this bacterium [20, 34]. To assess the colony morphology, culture positive samples were sub-cultured to selective media containing sodium azide (MilliporeSigma, Oakville, Ontario, Canada) and crystal violet agar (Ward’s Science, Rochester, New York, USA) (SACV) and incubated for 48 h at 37 ℃ with 5% CO2. A single colony from SACV agar that showed very small, transparent, and circular colonies was sub-streaked on 5% sheep blood agar (TSA with sheep blood) (ThermoScientific, USA) to obtain pure colonies for DNA extraction [35]. For molecular confirmation of species identity, one loopful of purified colonies from blood agar was suspended in 200 µl sterile phosphate buffer saline (PBS) for DNA extraction using DNeasy blood and tissue kit (Qiagen, Mississauga, ON, Canada) following manufacturer’s instructions. The confirmation of isolates was performed by qPCR using primers and a species-specific probe targeting the 3’ non-coding region of the rRNA gene cluster as described by Pal et al. [36]. E. tonsillarum was simultaneously ruled out (in multiplex) by including a second species-specific probe described in the same publication.
Library preparation, sequencing, and de novo assembly
Whole genome sequencing of E. rhusiopathiae isolates from Ellesmere Island, Somerset Island and the Sahtu region was outsourced to the Centre de diagnostic vétérinaire de l’Université de Montréal, Saint Hyacinthe, Quebec, Canada. At their facility, DNA extracted from isolates was quantified using Qubit™ dsDNA HS Assay Kit and a Qubit™ fluorometer (Thermofisher Scientific, MA, USA). Sequencing libraries were prepared with Nextera XT DNA Library Preparation Kit (Illumina, CA, USA) following standard protocols. Library quality was assessed using Agilent High Sensitivity DNA Kit in a Bioanalyzer (Agilent, CA, USA). Libraries were sequenced on a v3 600-cycle cartridge using a MiSeq instrument with PhiX spiked in at 1% (Illumina, CA, USA). The reads obtained were assessed for quality using FASTQC v.0.11.9 in Galaxy and the reads were trimmed for quality and adaptors using Trimmomatic v.0.36 [37] within the CLIMB computing platform for microbial genomics [38]. The reads of low-quality were removed using a sliding window (4:20) and MINLEN (50) parameters. The de novo assembly of newly-sequenced genomes was performed with SPAdes tool v.3.15.4 [39] also within CLIMB. Final quality assessment of scaffolds was performed using QUAST v.4.6.3 [40] (i.e., to determine the number of contigs, GC%, assembly length, and N50). Assemblies thus obtained were further assessed using Kmerfinder tool v.3.2 to confirm the species assignation [41]. The good quality assemblies with length (~ 1.8 mb) and GC content (~ 36.6%) consistent with the genome of E. rhusiopathiae and no evidence of contamination based on Kmerfinder were used for further analyses (n = 9); one exception was made for a tenth genome (593-EU04-K), where 2% of the assembly was classified as Pseudomonas aeruginosa. In addition, 45 high quality (coverage > 30X) whole genome sequences of previously characterized E. rhusiopathiae isolates, along with reference strains Fujisawa (GenBank accession no. AP012027), SY1027 (GenBank accession no. CP005079), WH13013 (GenBank accession no. CP017116), and ATCC19414 (GenBank accession no. NZ_ACLK00000000) were included in the core-genome alignment (Supplementary Table S1, Additional File 1). The dataset comprised E. rhusiopathiae isolates from a variety of hosts, including from muskoxen, caribou, moose, seal, swine, poultry, and marine mammals as previously described [20, 21]. Most of the genomes included were from Clade 3, although representation was included from other clades for use as an outgroup.
Core genome alignment
A core genome alignment, consisting of the subset of genetic regions conserved across the 57 E. rhusiopathiae de novo assemblies (clade 1 genomes excluded), was carried out using the Parsnp tool within the Harvest suite [42] through CLIMB. Parsnp identifies the orthologous sequences conserved in all the aligned genomes, from which it estimates the phylogenetic relationship using a maximum likelihood approach. The resulting phylogeny was visualized and edited with iTOL v6 [43].
Identification of amino acid variability among E. rhusiopathiae virulence genes
A total of 17 putative virulence genes were selected based on previous experimental studies and those that have been linked with certain pathogenic mechanisms of E. rhusiopathiae (Table 1). Nucleotide sequences of virulence genes from E. rhusiopathiae strain Fujisawa available on GenBank were downloaded and used to create custom BLAST databases in Geneious v.10.2. [44]. We searched for the presence of these 17 virulence genes in all de novo assemblies (n = 59) by BLASTn searches, where a BLAST hit of 95% pairwise identity was considered positive for a particular virulence gene. Nucleotide sequences were extracted from the positive hits and translations to amino acid sequences were performed within Geneious using translation Table 1 for bacteria as previously described [45]. Amino acid rather than nucleotide sequences were compared in order to focus on differences with a higher likelihood of phenotypic/functional relevance. Multiple sequence alignment to identify the amino acid variability was performed using CLUSTAL-2.1 implemented within Geneious, and variants were recorded in comparison with Fujisawa reference sequence. Phylogenetic trees for each gene were estimated from the multiple sequence alignment using a Neighbor-Joining method implemented in Geneious Tree Builder, using the Jukes-Cantor model.
Genome annotation and identification of PAIs and prophages
De novo assemblies were annotated using Prokka v.1.13 [55] within CLIMB. The resulting Genbank files were used as input into IslandViewer-4 software [56]. Prophages were searched for among the genomes using PHASTER tool (Phage Search Tool- Enhanced Release) [57]. This was done by uploading de novo assemblies in FASTA format to the online version of this tool (https://phaster.ca/). All tools were implemented using default settings.
Identification of unique gene content using Pan-GWAS analysis
A pan-genome wide association study (GWAS) approach was taken to identify gene content unique to the Arctic clone in comparison with other E. rhusiopathiae genomes. The Prokka-annotated genomes in general feature format (.gff) were used to analyse the pan-genome using Roary, using default settings [58]. The gene presence/absence file generated by Roary, along with a trait file, were used to establish gene clusters associated with the Arctic clone using Scoary [59]. The trait file was created in Microsoft excel by classifying all the isolates into binary categories, i.e., Arctic clone isolates were indicated as “1” and others as “0”. The p values for multiple testing and managing error-I were adjusted using the Benjamin-Hochberg method within Scoary. This method orders original p-values from smallest to largest and then the hypotheses are accepted or rejected to minimize the detection of false positives [60]. The list of significant genes for the Arctic clone generated by Scoary comprised both characterized (annotated) and hypothetical genes. We searched for gene function of annotated proteins using UniProt database [61].
Data availability
All newly sequenced genomes in this work were submitted to European Nucleotide Archive under project accession number PRJEB73478.
Results
Culture of E. rhusiopathiae from dead muskoxen from Fosheim Peninsula, Ellesmere Island, and other Arctic locations
Of the 35 samples cultured for E. rhusiopathiae from 14 muskox carcasses, 30 (85.7%) were culture positive (Supplementary Table S2, Additional File 1). Samples from one carcass each from Ellesmere Island and Somerset Island did not culture positive. Samples that showed growth were confirmed as E. rhusiopathiae using qPCR. In general, dead muskoxen were in good body condition. Several had bleeding from the nose and anus, and/or had serosanguinous fluid in the thorax and peritoneal cavities. None of the 22 fecal samples collected from two herds of apparently healthy muskoxen on Ellesmere Island cultured positive for E. rhusiopathiae.
Phylogenetic relationships of newly isolated E. rhusiopathiae compared with previously characterized isolates
Ten whole genome sequences newly-generated in this study were included in our genomic analyses. i) Nine whole genome sequences of E. rhusiopathiae isolated from Ellesmere (6) and Somerset (1) Islands and the Sahtu region (2) passed quality thresholds and were used for downstream analyses (Supplementary Table S3, Additional File 1). The contigs had an average N50 of 106,167 base pairs and an average of 104 contigs per assembly. ii) One additional genome of 593-EU04-K, isolated from an Ellesmere Island muskox, which had a lower N50 (79,868) and a higher number of contigs (257), was also included to assess at least one isolate per culture-positive muskox carcass.
The phylogenetic tree comprising 57 E. rhusiopathiae genomes, including 43 previously-characterised genomes, showed that E. rhusiopathiae isolates from Ellesmere Island fell within clade 3, and clustered together with E. rhusiopathiae from Banks, Victoria, and Prince Patrick Islands; these were isolated from muskox die-offs with symptoms of acute death and previously identified as being part of the Arctic clone (Fig. 2). Conversely, isolates from the Sahtu region (594-1-BM and 599-SAHMX01-Lu) and Somerset Island (591-AW1701-BM) muskoxen, which clustered within clade 3, were genetically distinct from the Arctic clone and from each other. Isolate AKM3-iv, from a muskox in Alaska [21], remained the closest isolate to the Arctic clone, with the next closest isolates from Prince Charles Island caribou, where mortalities were assumed to be caused by a severe weather-related starvation event [62, 63].
Amino acid variance between the Arctic clone and other closely related isolates
Of the 17 virulence genes investigated, 16 were present in all genomes, four of which were highly conserved, i.e., with a maximum of 9 AA variations among all 59 genomes: GAPDH, patatin like phospholipase B, patatin like phospholipase, and phospholipase D. Internalin was present in only 17/59 E. rhusiopathiae isolates examined. Four genes had specific amino acid variants that were present in all or nearly all Arctic clone isolates and absent in most of the other isolates investigated. These were adhesin, choline binding protein-B (cbpB), rhusiopathiae surface protein-A (rspA), and CDP-glycerol glycerophosphotransferase (tagF) (Fig. 3). Adhesin had 15 amino acid variations overall, two of which (V97I and D122G) were unique to the Arctic clone and two other clade 3 isolates (SKnoID1 (British Columbia, CA) and AKM3-iv (Alaska, US). Similarly, choline binding protein-B had 24 variations, two of which (R95L and N495D) were unique to the Arctic clone except genome MX10-1fe (Victoria Island, NWT, CA), which had the N495D variant but not R95L. Rhusiopathiae surface protein-A showed a higher number of variations, with variants at 73 (3.7%) of 1992 amino acid positions. With the exception of isolate 459-ER34-BM (Prince Patrick Island, NWT, CA), the Arctic clone sequences for this protein clustered separately, along with that of AKM3-iv (Alaska, US). Within CDP-glycerol glycerophosphotransferase sequences, a single amino acid variation at position 223 (D223N) was unique to the Arctic clone and 6 other isolates [AKM3-iv (Alaska, US), 600-16626-4bBM, 600-16266-5BM, 600-16626-14BM (Prince Charles Island, NU, CA), F440-C (Alberta, CA) and Z197-ii (Alberta, CA)]. However, F440-C had an additional variation at the position 104 (R104Q) which led to this isolate segregating from others in the tree (Fig. 3D).
Only hylA and spaA had insertions variably present among the genomes. For instance, the Arctic clone and a limited number of non-Arctic clone isolates exhibited the presence of a tandem repeat in the hylA gene consisting of four amino acids (ETKP) repeated four times starting from position 984, whereas other isolates had fewer than four copies of the repeat motif (Supplementary Fig. S1, Additional File 1).
Identification of pathogenicity islands and prophages
We screened for the presence of virulence genes within PAIs using the IslandViewer-4 tool (Supplementary Table S4, Additional File 1). Over 183 PAIs of variable length and location were predicted across all 59 E. rhusiopathiae genomes investigated. In general, genomes isolated from swine showed a higher number of PAIs, whereas the Arctic clone isolates each had 2–3 PAIs. All of the PAIs harboured genes associated with horizontal transfer such as transposase, integrase, recombinase, and prophage genes. PAIs from porcine isolates (B3129, 6028, 8929) harboured genes encoding antibiotic resistance such as aminoglycoside 6-adenylyltransferase, streptomycin 3’-adenylyltransferase, tetracycline resistance (tetM) from transposon TnFO1, bleomycin resistance protein, and kanamycin nucleotidyltransferase. Isolates other than those of swine origin in clade 3 contained putative virulence genes identified as collagen adhesin, CFA/I fimbrial subunit D, multidrug export protein MepA, transcriptional regulatory protein SrrA, and Internalin-J. Finally, two putative toxin-encoding genes were found within PAI of several genomes, including within the Arctic clone: NAD (+)—arginine ADP- ribosyl transferase EFV, and toxin B.
The PHASTER online tool [57] identified prophages present in all 59 E. rhusiopathiae genomes. A total of 137 prophages accounting for 3 intact, 126 incomplete and 8 questionable prophages were found. Incomplete and questionable prophages were present in most of the E. rhusiopathiae genomes and ranged from 5.9 to 62.4 kb in length. All the Arctic clone genomes carried 2–3 prophages each. Meanwhile, non-Arctic clone genomes carried between 1 and 4 prophages. A specific prophage 11.2 kb long was found in the Arctic clone, as well as in one other non-Arctic clone isolate (594-1BM). However, none of the prophages harboured known virulence genes.
Identification of gene content unique to the Arctic clone through pan-GWAS
Scoary identified a total of 101 genes that were significantly (p < 0.05) associated with the Arctic clone (Supplementary Table S5, Additional File 2). Of those, 69 genes were annotated as hypothetical proteins with unknown function. The top 15 annotated genes with highest sensitivity and specificity are shown in Table 2.
Predicted functions of unique genes found in the Arctic clone
Based on analyses using Uniprot [61], the top hit genes from Scoary had predicted functions related to DNA binding and repair. More specifically, four genes are involved in DNA binding: hhalM_1 encodes a DNA methyltransferase enzyme that is a component of restriction modification systems [64]; xerD encodes a phage integrase SAM-like protein that has role in DNA binding [65]; noc_1 encodes a nucleoid occlusion protein; lexA encodes a DNA-binding helix turn helix protein responsible for DNA binding and serine-type endopeptidase activity. The DNA repair gene radD is involved in alteration of one or more nucleotide sites in DNA through a DNA repair helicase protein. The N-acetylhexosamine 1-kinase encoded by nahK_2 regulates phosphorylation responsible for cleaving DNA and integrating phage-related fragments. The gene cpsY encodes an exopolysaccharide phosphotransferase enzyme that has a role in exopolysaccharide synthesis. Group_1036 encodes a putative ABC transporter permease enzyme. Two genes are involved in teichoic acid biosynthesis process: the gene tarF encodes teichoic acid poly (glycerol phosphate) polymerase and tarL gene encodes Ribitol-5-phosphate cytidylyltransferase 2. The gene arnB encodes a UDP-4-amino-4-deoxy-L-arabinoseoxoglutarate that helps in intracellular protein transport. Finally, two toxin-associated genes were among the top Scoary hits with predicted functions: Group_140 encodes NAD (+)—arginine ADP- ribosyltransferase EFV (Uniprot ID Q838U8), and toxB gene encodes toxin B (Uniprot ID P18177).
Discussion
New understanding of the geographical distribution of the Arctic clone
Phylogenetic analysis of E. rhusiopathiae isolates from the Fosheim Peninsula, Ellesmere Island, analyzed alongside a global collection of genomes, confirmed that the Arctic clone is present in this geographical location. The disease presentation was similar to what had been previously described in other locations, i.e. summer mortalities, animals of all age and sex classes and in good condition, septicemia and acute death, and bleeding from natural orifices [10]. Ellesmere Island is located over 1000 km to the northeast of Banks and Victoria Islands, where the Arctic clone was first identified in association with muskox mortalities starting in 2010. Wildlife health surveillance in the Arctic is extremely limited due to vast distances and sparse human populations. It is thus not clear how long the Arctic clone has been present on the Fosheim peninsula, how widespread it is, and what its impacts are on muskoxen or other species. Further research looking into its distribution in the eastern Arctic Archipelago, environmental persistence, host range and transmission, as well as establishing a molecular clock to understand patterns of spread, are underway. The finding of diverse strains of E. rhusiopathiae on Somerset Island and in the Sahtu region is consistent with previous observations that a wide variety of strains of this bacterium are found in wildlife in other locations [21, 22], and emphasizes the novelty of finding this unique clone across such a wide geographic area.
Virulence gene sequence variants unique to the Arctic clone
The variations we observed in the amino acid sequences – as predicted from the whole genome sequence data – have the potential to alter protein production and activity (e.g., through changes in shape or residual charge) [66], and is thus one of the potential mechanisms that could alter the virulence of particular strains. We found that the Arctic clone exhibited unique amino acid sequences in 4 out of the 17 virulence genes investigated in this study. Of these four genes, two (adhesin, and choline binding protein B) are related to adhesion of E. rhusiopathiae to host endothelial cells [49, 50]. In addition, rhusiopathiae surface protein A is involved in adherence to abiotic and biotic surfaces through biofilm formation [52]. For pathogens to withstand host immunological mechanisms and to deploy available virulence genes, adhesion is a universal prerequisite [67]. Furthermore, CDP-glycerol glycerophosphotransferase has been suggested to encode a component of the cell wall, contributing to the structural integrity and pathogenic potential of the bacterium [48].
In addition to these unique virulence gene variants, we identified different copy numbers of a 4-amino acid tandem repeat motif in the hyaluronidase (hylA) gene in the Arctic clone, as well as some non-Arctic clone isolates, whereas other genomes showed 2–3 copies. Variability in tandem repeat units of different genes has been previously associated with virulence and stress-related strategies amongst strains of Salmonella spp. [68]. Hyaluronidase is an important virulence factor produced by bacteria, helminths and arthropods, and is hypothesized to be a spreading factor that disseminates these pathogens into tissues [69]. Although some studies have suggested a role of hylA in E. rhusiopathiae pathogenesis [70], its importance in virulence has been heavily disputed for this species [69, 71]. In addition to the unknown relevance of these differential repeats, there is also some uncertainty in the de novo assembly of repetitive regions of DNA using short sequence reads, such as those generated using Illumina. These limitations include potential assembly errors, incomplete assemblies, and challenges in accurately resolving complex repeat structures due to short read lengths, which complicate the understanding and analysis of these genomic regions [72]. The addition of long-read data, e.g., using Nanopore, would increase the confidence in the assembly of such genomic regions. Overall, further studies to assess potential functional implications of these predicted unique amino acid differences in the Arctic clone could include using structural predictions, as well as in vitro and in vivo studies.
While the Arctic clone genomes harbored 2–3 prophages, none of the genes they carried had putative roles in virulence. These mobile elements therefore seem less likely to be linked to any unique functional differences of this strain. Some of the PAI detected carried potential virulence genes that Scoary detected as being over-represented in the Arctic clone compared with other genomes, as discussed below.
Genes associated with the Arctic clone
Scoary analyses found several unique genes to be overrepresented in the Arctic clone compared to other closely related genomes. A high proportion of genes (69/101 = 68%) were hypothetical, i.e., with unknown function. Of the top 15 genes associated with the Arctic clone having a predicted function, it is perhaps not surprising that most are involved in DNA recombination and repair. E. rhusiopathiae is highly recombinogenic [20] and these genes are necessary for integration of foreign genetic material [73]. Further supporting the importance of recombination in this bacterium, several of these unique genes were among those found on pathogenicity islands, including tyrosine recombinase XerD, nucleoid occlusion protein Noc_1, and putative DNA repair helicase RadD. Moreover, of note was the detection of two genes annotated as putative toxins – toxin B and NAD (+)—arginine ADP- ribosyltransferase EFV – carried on PAIs of E. rhusiopathiae. No toxin has ever been associated with this bacterium previously [11]. The gene annotated as toxin B by Prokka was 1827 nucleotides in length. While the Uniprot reference (P18177) was to an exotoxin produced by Clostridium difficile that causes cytotoxicity (7692 bp) [74], there was low homology between the sequences. This gene has been previously annotated as N-acetylmuramoyl-L-alanine amidase family protein in other E. rhusiopathiae genomes (e.g., RNM30934.1). The gene annotated as NAD (+)—arginine ADP- ribosyl transferase EFV by Prokka was 1581 nucleotides in length. The Uniprot reference for this predicted protein (Q838U8) is associated with Enterococcus faecalis (1461 bp) [75]. This gene has been previously annotated as a minor capsid protein of E. rhusiopathiae (e.g., WP_115354732.1), and is likely a component of a bacteriophage. BLAST searches of both genes did not find any significant homology in other bacterial species. Investigations into the function of these two genes, particularly whether they act as toxins, would be of high interest. One annotated gene, hhaIM_1, showed perfect sensitivity and specificity, and thus represents a good candidate for the development of a PCR assay for the rapid confirmation of the Arctic clone; this is currently under development by our team.
Conclusion
E. rhusiopathiae has emerged as an important and enigmatic pathogen of muskoxen and other Arctic species. It is associated with widespread mortality events, high mortality rates and concomitant to population declines of muskoxen in the Canadian Arctic archipelago and other populations [10, 23]. Our work demonstrates that a single, apparently highly pathogenic strain, the Arctic clone, first detected in the western Arctic Archipelago in 2010, is now present on the Fosheim Peninsula of Ellesmere Island, to the far northeast. This Arctic clone possesses several unique genes, including putative toxin-encoding genes, and distinct amino acid variations in several conserved virulence genes, all of which may be playing a role in the apparent high pathogenicity of this strain. Understanding the functional significance of these differences, together with unraveling the factors related to the pathogen, a genetically depauperate host [6], and a rapidly changing environment [76] in the observed outbreaks, are critical for understanding the impacts of this pathogen on muskoxen, other species, and the broader Arctic ecosystem. The comprehensive comparative genome analysis presented here will guide further studies (e.g., in vitro or challenge studies) to elucidate the genetic and functional basis of the pathogenesis of E. rhusiopathiae.
Data availability
Sequence data generated for this manuscript are available at European Nucleotide Archive under project accession number PRJEB73478. Other publicly available sequences included in this study are described in Supplementary Table S1, Additional File 1.
References
Tomaselli M, Gerlach SC, Kutz SJ, Checkley SL. Community of Iqaluktutiaq. Iqaluktutiaq voices: local perspectives about the importance of Muskoxen, Contemporary and Traditional Use and practices. Arctic. 2018;71:1–14.
Tomaselli M, Ytrehus B, Opriessnig T, Duignan P, Dalton C, van der Meer F et al. Contagious ecthyma dermatitis as a portal of entry for Erysipelothrix rhusiopathiae in muskoxen (Ovibos moschatus) of the Canadian Arctic. J Wildl Dis. 2022;58.
Tomaselli M, Elkin B, Kutz S, Harms NJ, Nymo HI, Davison T, et al. A transdisciplinary approach to Brucella in muskoxen of the western Canadian Arctic 1989–2016. EcoHealth. 2019;16:488–501.
Kutz SJ, Jenkins EJ, Veitch AM, Ducrocq J, Polley L, Elkin B, et al. The Arctic as a model for anticipating, preventing, and mitigating climate change impacts on host-parasite interactions. Vet Parasitol. 2009;163:217–28.
Kafle P, Peller P, Massolo A, Hoberg E, Leclerc L-M, Tomaselli M, et al. Range expansion of muskox lungworms track rapid arctic warming: implications for geographic colonization under climate forcing. Sci Rep. 2020;10:17323.
Prewer E, Kutz S, Leclerc LM, Kyle CJ. Already at the bottom? Demographic declines are unlikely further to undermine genetic diversity of a large Arctic ungulate: muskox, Ovibos moschatus (Artiodactyla: Bovidae). Biol J Linn Soc. 2020;129:459–69.
Kutz S, Rowell J, Adamczewski J, Gunn A, Cuyler C, Aleuy OA et al. Muskox health ecology symposium 2016: Gathering to share knowledge on Umingmak in a time of rapid change. ARCTIC. 2017;70:225–36.
Cuyler C, Rowell J, Adamczewski J, Anderson M, Blake J, Bretten T, et al. Muskox status, recent variation, and uncertain future. Ambio. 2020;49:805–19.
Di Francesco J. Indigenous Knowledge and Biomarkers of Physiological Stress Inform Muskox Conservation in a Rapidly Changing Arctic. 2020.
Kutz S, Bollinger T, Branigan M, Checkley S, Davison T, Dumond M, et al. Erysipelothrix rhusiopathiae associated with recent widespread muskox mortalities in the Canadian Arctic. Can Vet J. 2015;56:560–3.
Wang Q, Chang BJ, Riley TV. Erysipelothrix rhusiopathiae. Vet Microbiol. 2010;140:405–17.
Bender JS, Shen HG, Irwin CK, Schwartz KJ, Opriessnig T. Characterization of Erysipelothrix species isolates from clinically affected pigs, environmental samples, and vaccine strains from six recent swine erysipelas outbreaks in the United States. Clin Vaccine Immunol. 2010;17:1605–11.
Stephenson EH, Berman DT. Isolation of Erysipelothrix rhusiopathiae from tonsils of apparently normal swine by two methods. Am J Vet Res. 1978;39:187–8.
Opriessnig T, Coutinho TA. Erysipelas. Diseases of Swine. Wiley; 2019. pp. 835–43.
Jensen WI, Cotter SE. An outbreak of erysipelas in eared grebes (Podiceps nigricollis). J Wildl Dis. 1976;12:583–6.
Wayson NE. An epizootic among meadow mice in California, caused by the bacillus of mouse septicemia or of swine erysipelas. Public Health Rep. 1927;42:1489–93. (1896–1970).
Danil K, Colegrove KM, Delaney MA, Mena A, Stedman N, Wurster E. Systemic Erysipelas Outbreak among Free-Ranging Bottlenose Dolphins, San Diego, California, USA, 2022. Emerg Infect Dis. 2023;29:2561–3.
IJsseldijk LL, Begeman L, Duim B, Gröne A, Kik MJL, Klijnstra MD, et al. Harbor Porpoise Deaths Associated with Erysipelothrix rhusiopathiae, the Netherlands, 2021. Emerg Infect Dis. 2023;29:835–8.
Mavrot F, Orsel K, Hutchins W, Adams LG, Beckmen K, Blake JE, et al. Novel insights into serodiagnosis and epidemiology of Erysipelothrix rhusiopathiae, a newly recognized pathogen in muskoxen (Ovibos moschatus). PLoS ONE. 2020;15:e0231724.
Forde T, Biek R, Zadoks R, Workentine ML, De Buck J, Kutz S, et al. Genomic analysis of the multi-host pathogen Erysipelothrix rhusiopathiae reveals extensive recombination as well as the existence of three generalist clades with wide geographic distribution. BMC Genomics. 2016;17:461.
Forde TL, Orsel K, Zadoks RN, Biek R, Adams LG, Checkley SL et al. Bacterial genomics reveal the complex epidemiology of an emerging pathogen in Arctic and boreal ungulates. Front Microbiol. 2016;7.
Mavrot F, Forde T, Tomaselli M, Anderson M, Leclerc L-M, Fenton H, et al. Studying rare occurrences: long-term phylogenetic monitoring of the emerging pathogen Erysipelothrix rhusiopathiae during wildlife mortality events. Online; 2020.
Davison T, Baryluk S. Aerial survey of Muskoxen (Ovibos moschatus) and Peary Caribou (Rangifer tarandus pearyi) on Banks Island, July and August 2019. Environment and Natural Resources, Government of the Northwest Territories; 2021.
Pečnerová P, Lord E, Garcia-Erill G, Hanghøj K, Rasmussen MS, Meisner J, et al. Population genomics of the muskox’ resilience in the near absence of genetic variation. Mol Ecol. 2023;33:e17205.
Di Francesco J, Hanke A, Milton T, Leclerc L-M, Association KA, Gerlach C, et al. Documenting Indigenous Knowledge to identify and understand the stressors of Muskoxen (Ovibos moschatus) in Nunavut, Canada. Arctic. 2022;74:418–36.
Hellberg RS, Chu E. Effects of climate change on the persistence and dispersal of foodborne bacterial pathogens in the outdoor environment: a review. Crit Rev Microbiol. 2016;42:548–72.
Kock RA, Orynbayev M, Robinson S, Zuther S, Singh NJ, Beauvais W, et al. Saigas on the brink: multidisciplinary analysis of the factors influencing mass mortality events. Sci Adv. 2018;4:eaao2314.
Juhas M, van der Meer JR, Gaillard M, Harding RM, Hood DW, Crook DW. Genomic islands: tools of bacterial horizontal gene transfer and evolution. FEMS Microbiol Rev. 2009;33:376–93.
Ogawa Y, Ooka T, Shi F, Ogura Y, Nakayama K, Hayashi T, et al. The genome of Erysipelothrix rhusiopathiae, the causative agent of swine erysipelas, reveals new insights into the evolution of Firmicutes and the organism’s intracellular adaptations. J Bacteriol. 2011;193:2959–71.
Shimoji Y. Pathogenicity of Erysipelothrix rhusiopathiae: virulence factors and protective immunity. Microbes Infect. 2000;2:965–72.
Yang L, Zhu Y, Peng Z, Ding Y, Jie K, Wang Z, et al. Comparative genome analysis of a pathogenic Erysipelothrix rhusiopathiae isolate WH13013 from pig reveals potential genes involve in bacterial adaptions and pathogenesis. Veterinary Sci. 2020;7:74.
Gu J, Li Y, Xu C, Xie X, Li P, Ma G, et al. Genome sequence of multidrug-resistant Erysipelothrix rhusiopathiae ZJ carrying several acquired antimicrobial resistance genes. J Global Antimicrob Resist. 2020;21:13–5.
Sotomayor-Vivas C, Hernández-Lemus E, Dorantes-Gilardi R. Linking protein structural and functional change to mutation using amino acid networks. PLoS ONE. 2022;17:e0261829.
Bender JS, Kinyon JM, Kariyawasam S, Halbur PG, Opriessnig T. Comparison of conventional direct and enrichment culture methods for Erysipelothrix spp. from experimentally and naturally infected swine. J Vet Diagn Invest. 2009;21:863–8.
Packer RA. The use of sodium azide (NaN (3)) and crystal violet in a selective medium for Streptococci and Erysipelothrix rhusiopathiae. J Bacteriol. 1943;46:343–9.
Pal N, Bender JS, Opriessnig T. Rapid detection and differentiation of Erysipelothrix spp. by a novel multiplex real-time PCR assay. J Appl Microbiol. 2010;108:1083–93.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Connor TR, Loman NJ, Thompson S, Smith A, Southgate J, Poplawski R et al. CLIMB (the Cloud Infrastructure for Microbial Bioinformatics): an online resource for the medical microbiology community. Microb Genomics. 2016;2.
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:455–77.
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: quality assessment tool for genome assemblies. Bioinformatics. 2013;29:1072–5.
Larsen MV, Cosentino S, Lukjancenko O, Saputra D, Rasmussen S, Hasman H, et al. Benchmarking of methods for genomic taxonomy. J Clin Microbiol. 2014;52:1529–39.
Treangen TJ, Ondov BD, Koren S, Phillippy AM. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 2014;15:524.
Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49:W293–6.
Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647–9.
Forde TL, Kollanandi Ratheesh N, Harvey WT, Thomson JR, Williamson S, Biek R, et al. Genomic and Immunogenic Protein Diversity of Erysipelothrix rhusiopathiae isolated from pigs in Great Britain: implications for Vaccine Protection. Front Microbiol. 2020;11:418.
Wang Y, Li J, Zhang A, Zhu W, Zhang Q, Xu Z, et al. iTRAQ-based quantitative proteomic analysis reveals potential virulence factors of Erysipelothrix rhusiopathiae. J Proteom. 2017;160:28–37.
Zhu Weifeng C, Chengzhi W, Ya L, Jingtao W, Chao K, Chao, et al. Characterization of roles of SpaA in Erysipelothrix rhusiopathiae adhesion to porcine endothelial cells. Microb Pathog. 2017;113:176–80.
Shimoji Y, Ogawa Y, Tsukio M, Shiraiwa K, Nishikawa S, Eguchi M. Genome-wide identification of virulence genes in Erysipelothrix rhusiopathiae: use of a mutant deficient in a tagF homolog as a safe oral vaccine against swine erysipelas. Infect Immun. 2019;87.
Zhu W, Cai C, Huang J, Liu L, Xu Z, Sun X, et al. Characterization of pathogenic roles of two Erysipelothrix rhusiopathiae surface proteins. Microb Pathog. 2018;114:166–8.
Zhu W, Cai C, Li J, Zhang Q, Huang J, Jin M. Characterization of protective antigen CbpB as an adhesin and a plasminogen-binding protein of Erysipelothrix rhusiopathiae. Res Vet Sci. 2019;124:352–6.
Ding Y, Zhu D, Zhang J, Yang L, Wang X, Chen H, et al. Virulence determinants, antimicrobial susceptibility, and molecular profiles of Erysipelothrix rhusiopathiae strains isolated from China. Emerg Microbes Infections. 2015;4:1–5.
Shimoji Y, Ogawa Y, Osaki M, Kabeya H, Maruyama S, Mikami T, et al. Adhesive surface proteins of Erysipelothrix rhusiopathiae bind to polystyrene, fibronectin, and type I and IV collagens. J Bacteriol. 2003;185:2739–48.
Nakato H, Shinomiya K, Mikawa H. Adhesion of Erysipelothrix rhusiopathiae to cultured rat aortic endothelial cells. Pathol - Res Pract. 1987;182:255–60.
Zhu Weifeng Z, Qiang L, Jingtao W, Yanmin C, Chengzhi L, Liang, et al. Glyceraldehyde-3-phosphate dehydrogenase acts as an adhesin in Erysipelothrix rhusiopathiae adhesion to porcine endothelial cells and as a receptor in recruitment of host fibronectin and plasminogen. Vet Res. 2017;48:16.
Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.
Bertelli C, Laird MR, Williams KP, Lau BY, Hoad G, Winsor GL, et al. IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic Acids Res. 2017;45:W30–5.
Arndt D, Marcu A, Liang Y, Wishart DS. PHAST, PHASTER and PHASTEST: tools for finding prophage in bacterial genomes. Brief Bioinform. 2018;20:1560–7.
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.
Brynildsrud O, Bohlin J, Scheffer L, Eldholm V. Rapid scoring of genes in microbial pan-genome-wide association studies with Scoary. Genome Biol. 2016;17:238.
Chen S-Y, Feng Z, Yi X. A general introduction to adjustment for multiple comparisons. J Thorac Disease. 2017;9:1725–9.
Bateman A, Martin M-J, Orchard S, Magrane M, Agivetova R, Ahmad S, et al. UniProt: the universal protein knowledgebase in 2021. Nucleic Acids Res. 2021;49:D480–9.
Kutz S. Diagnostic services Unit Case Report (Accession 16–626). University of Calgary; 2016.
Dolant C, Montpetit B, Langlois A, Brucker L, Zolina O, Johnson CA, et al. Assessment of the barren ground caribou die-off during winter 2015–2016 using passive microwave observations. Geophys Res Lett. 2018;45:4908–16.
Casselli T, Tourand Y, Scheidegger A, Arnold WK, Proulx A, Stevenson B et al. DNA methylation by restriction modification systems affects the global transcriptome profile in Borrelia burgdorferi. J Bacteriol. 2018;200.
Subramanya Hosahalli S, Arciszewka K, Rachel L, Louise B, David B, Dale S. Crystal structure of the site-specific recombinase,XerD. e EMBO J. 1997;16:5178–87.
Biro J. Amino acid size, charge, hydropathy indices and matrices for protein structure analysis. Theoretical Biology Med Modelling. 2006;3:15.
Krachler AM, Orth K. Targeting the bacteria–host interface. Virulence. 2013;4:284–94.
Das G, Das S, Dutta S, Ghosh I. In silico identification and characterization of stress and virulence associated repeats in Salmonella. Genomics. 2018;110:23–34.
Shimoji Y, Asato H, Seizaki T, Mori Y, Yokomizo Y. Hyaluronidase is not essential for the lethality of Erysipelothrix rhusiopathiae infection in mice. J Vet Med Sci. 2002;64:173–6.
Mann S. Cell wall composition of Listeria monocytogenes and Erysipelothrix rhusiopathiae. Zentralbl Bakteriol Orig. 1969;209:510–22.
Nørrung V. Studies on Erysipelothrix Insidiosa S. Rhusiopathiae. Acta Vet Scand. 1970;11:577–85.
Treangen TJ, Salzberg SL. Repetitive DNA and next-generation sequencing: computational challenges and solutions. Nat Rev Genet. 2011;13:36–46.
Kowalczykowski SC. An overview of the molecular mechanisms of recombinational DNA repair. Cold Spring Harb Perspect Biol. 2015;7:a016410.
Lyras D, O’Connor JR, Howarth PM, Sambol SP, Carter GP, Phumoonna T, et al. Toxin B is essential for virulence of Clostridium difficile. Nature. 2009;458:1176–9.
Paulsen IT, Banerjei L, Myers GSA, Nelson KE, Seshadri R, Read TD, et al. Role of mobile DNA in the evolution of Vancomycin-resistant Enterococcus faecalis. Science. 2003;299:2071–4.
IPCC. Climate Change 2022: Impacts, Adaptation, and Vulnerability. Contribution of Working Group II to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, UK and New York, USA: Cambridge University Press; 2022.
Acknowledgements
We thank Tabitha Mullin and Matt Fredlund with Government of Nunavut, Don Lavalee with Environment Canada Eureka Weather Station, Ronan Donovan and BBC Natural History Unit for supporting logistics of 2021 field work; Kevin Chan and Naima Jutha of Government of Northwest Territories for support during sample collection; Angela Schneider for laboratory assistance; and Sreejith Radhakrishnan for assistance in depositing sequence data.
Funding
This work was supported by funding to SJK: Polar Knowledge Canada (NST-2122-0049), NSERC Canada Research Chair in Arctic One Health (CRC-2020-00315), NSERC Northern Supplement (RGPNS/02796 − 2021), Canada North Outfitting, and Weber Arctic. Logistical and in-kind support was provided by the Government of Nunavut and the Government of Northwest Territories.
Author information
Authors and Affiliations
Contributions
LVS, TLF, YDN and SJK contributed to study conceptualization. LVS, TLF and SJK were responsible for data curation. LVS led the formal analysis and figure development, with support from TLF. LVS conducted the lab analyses, and wrote the original draft. AR-C and FM contributed to field investigations. TLF, YDN and SJK contributed to supervision. SJK led the acquisition of funding. LVS, TLF, FM, YDN and SJK contributed to manuscript review and editing. All authors read and approved the final manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
All work was done under University of Calgary Animal Care Protocol number AC18-0093. Wildlife Research Permits were obtained from the Governments of Nunavut (WL2016-058, WL2021-053, WL2019-051, WL2021-048) and the Northwest Territories (WL 500664).
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
12864_2024_10592_MOESM1_ESM.docx
Supplementary Material 1. Additional File 1, Table S1. Erysipelothrix rhusiopathiae genomes included in this study. Table S2. Culture results for Erysipelothrix rhusiopathiae from tissue samples collected from muskox mortalities. Table S3. Quality metrics of 10 newly sequenced E. rhusiopathiae genome assemblies, as determined by QUAST. Table S4. Prophages and pathogenicity islands (PAIs) detected in 59 Erysipelothrix rhusiopathiae genomes. Fig. S1. Multiple sequence alignment of hyaluronidase amino acid sequence from 57 E. rhusiopathiae genomes.
12864_2024_10592_MOESM2_ESM.csv
Supplementary Material 2. Additional File 2. Table S5. Scoary results for the significant genes associated with the Arctic clone.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Seru, L.V., Forde, T.L., Roberto-Charron, A. et al. Genomic characterization and virulence gene profiling of Erysipelothrix rhusiopathiae isolated from widespread muskox mortalities in the Canadian Arctic Archipelago. BMC Genomics 25, 691 (2024). https://doi.org/10.1186/s12864-024-10592-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12864-024-10592-9