Bacterial Genomics Reveal the Complex Epidemiology of an Emerging Pathogen in Arctic and Boreal Ungulates

Genomic characterization and virulence gene profiling of Erysipelothrix rhusiopathiae isolated from widespread muskox mortalities in the Canadian Arctic Archipelago

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.

The geographic distribution, and the biotic and abiotic predictors of select zoonotic pathogen detections in Canadian polar bears

Increasing Arctic temperatures are facilitating the northward expansion of more southerly hosts, vectors, and pathogens, exposing naïve populations to pathogens not typical at northern latitudes. To understand such rapidly changing host-pathogen dynamics, we need sensitive and robust surveillance tools. Here, we use a novel multiplexed magnetic-capture and droplet digital PCR (ddPCR) tool to assess a sentinel Arctic species, the polar bear (Ursus maritimus; n = 68), for the presence of five zoonotic pathogens (Erysipelothrix rhusiopathiae, Francisella tularensis, Mycobacterium tuberculosis complex, Toxoplasma gondii and Trichinella spp.), and observe associations between pathogen presence and biotic and abiotic predictors. We made two novel detections: the first detection of a Mycobacterium tuberculosis complex member in Arctic wildlife and the first of E. rhusiopathiae in a polar bear. We found a prevalence of 37% for E. rhusiopathiae, 16% for F. tularensis, 29% for Mycobacterium tuberculosis complex, 18% for T. gondii, and 75% for Trichinella spp. We also identify associations with bear age (Trichinella spp.), harvest season (F. tularensis and MTBC), and human settlements (E. rhusiopathiae, F. tularensis, MTBC, and Trichinella spp.). We demonstrate that monitoring a sentinel species, the polar bear, could be a powerful tool in disease surveillance and highlight the need to better characterize pathogen distributions and diversity in the Arctic.

Erysipelothrix amsterdamensis sp. nov., associated with mortalities among endangered seabirds

Infectious diseases threaten endangered species, particularly in small isolated populations. Seabird populations on the remote Amsterdam Island in the Indian Ocean have been in decline for the past three decades, with avian cholera caused by Pasteurella multocida proposed as the primary driver. However, Erysipelothrix species have also been sporadically detected from albatrosses on Amsterdam Island and may be contributing to some of the observed mortality. In this study, we genomically characterized 16 Erysipelothrix species isolates obtained from three Indian yellow-nosed albatross (Thalassarche carteri) chick carcasses in 2019. Histological analyses suggest that they died of bacterial septicaemia. Two isolates were sequenced using both Illumina short-read and MinION long-read approaches, which - following hybrid assembly - resulted in closed circular genomes. Mapping of Illumina reads from the remaining isolates to one of these new reference genomes revealed that all 16 isolates were closely related, with a maximum of 13 nucleotide differences distinguishing any pair of isolates. The nucleotide diversity of isolates obtained from the same or different carcasses was similar, suggesting all three chicks were likely infected from a common source. These genomes were compared with a global collection of genomes from Erysipelothrix rhusiopathiae and other species from the same genus. The isolates from albatrosses were phylogenetically distinct, sharing a most recent common ancestor with E. rhusiopathiae. Based on phylogenomic analysis and standard thresholds for average nucleotide identity and digital DNA-DNA hybridization, these isolates represent a novel Erysipelothrix species, for which we propose the name Erysipelothrix amsterdamensis sp. nov. The type strain is A18Y020dT (=CIP 112216T=DSM 115297T). The implications of this bacterium for albatross conservation will require further study.

A new multiplexed magnetic capture-Droplet digital PCR tool for monitoring wildlife population health and pathogen surveillance

Anthropogenic stressors are exacerbating the emergence and spread of pathogens worldwide. In regions like the Arctic, where ecosystems are particularly susceptible, marked changes are predicted in regional diversity, intensity, and patterns of infectious diseases. To understand such rapidly changing host-pathogen dynamics and mitigate the impacts of novel pathogens, we need sensitive disease surveillance tools. We developed and validated a novel multiplexed, magnetic capture, and ddPCR tool for the surveillance of multiple pathogens in polar bears, a sentinel species that is considered susceptible to climate change and other stressors with a pan-Arctic distribution. Through sequence-specific magnetic capture, we concentrated five target template sequences from three zoonotic bacteria (Erysipelothrix rhusiopathiae, Francisella tularensis, and Mycobacterium tuberculosis complex) and two parasitic (Toxoplasma gondii and Trichinella spp.) pathogens from large quantities (<100 g) of host tissue. We then designed and validated two multiplexed probe-based ddPCR assays for the amplification and detection of the low-concentration target DNA. Validations used 48 polar bear tissues (muscle and liver). We detected 14, 1, 3, 4, and 22 tissue positives for E. rhusiopathiae, F. tularensis, M. tuberculosis complex, T. gondii, and Trichinella spp., respectively. These multiplexed assays offer a rapid, specific tool for quantifying and monitoring the changing geographical and host distributions of pathogens relevant to human and animal health.

An integrative and multi-indicator approach for wildlife health applied to an endangered caribou herd

Assessing wildlife health in remote regions requires a multi-faceted approach, which commonly involves convenient samplings and the need of identifying and targeting relevant and informative indicators. We applied a novel wildlife health framework and critically assessed the value of different indicators for understanding the health status and trends of an endangered tundra caribou population. Samples and data from the Dolphin and Union caribou herd were obtained between 2015 and 2021, from community-based surveillance programs and from captured animals. We documented and categorized indicators into health determinants (infectious diseases and trace elements), processes (cortisol, pathology), and health outcomes (pregnancy and body condition). During a recent period of steep population decline, our results indicated a relatively good body condition and pregnancy rates, and decreasing levels of stress, along with a low adult cow survival. We detected multiple factors as potential contributors to the reduced survival, including Brucella suis biovar 4, Erysipelothrix rhusiopathiae and lower hair trace minerals. These results remark the need of targeted studies to improve detection and investigations on caribou mortalities. We also identified differences in health indicators between captured and hunter sampled caribou, highlighting the importance of accounting for sampling biases. This integrative approach that drew on multiple data sources has provided unprecedented knowledge on the health in this herd and highlights the value of documenting individual animal health to understand causes of wildlife declines.

Cook Inlet beluga whale Delphinapterus leucas with valvular endocarditis, encephalitis, rhabdomyolysis, heavy parasitism and fungal dermatitis

This is a case report of a Cook Inlet beluga whale Delphinapterus leucas found dead stranded on September 28, 2020 in Turnagain Arm, Alaska. This subadult male had valvular endocarditis, encephalitis, rhabdomyolysis, myoglobinuric nephropathy, severe parasitism and fungal dermatitis. Erysipelothrix rhusiopathiae was detected in the heart lesion, eye and external swabs. The level of infection and parasitism in this individual is markedly higher than what has been found in other Cook Inlet belugas, suggesting immunosuppression.

Population Structure and Genomic Characteristics of Australian Erysipelothrix rhusiopathiae Reveals Unobserved Diversity in the Australian Pig Industry

Erysipelothrix rhusiopathiae is a bacterial pathogen that is the causative agent of erysipelas in a variety of animals, including swine, emus, turkeys, muskox, caribou, moose, and humans. This study aims to investigate the population structure and genomic features of Australian isolates of E. rhusiopathiae in the Australian pig industry and compare them to the broader scope of isolates worldwide. A total of 178 isolates (154 Australian, seven vaccine isolates, six international isolates, and 11 of unknown origin) in this study were screened against an MLST scheme and publicly available reference isolates, identifying 59 new alleles, with isolates separating into two main single locus variant groups. Investigation with BLASTn revealed the presence of the spaA gene in 171 (96%) of the isolates, with three main groups of SpaA protein sequences observed amongst the isolates. Novel SpaA protein sequences, categorised here as group 3 sequences, consisted of two sequence types forming separate clades to groups 1 and 2, with amino acid variants at positions 195 (D/A), 303 (G/E) and 323(P/L). In addition to the newly identified groups, five new variant positions were identified, 124 (S/N), 307 (Q/R), 323 (P/L), 379 (M/I), and 400 (V/I). Resistance screening identified genes related to lincomycin, streptomycin, erythromycin, and tetracycline resistance. Of the 29 isolates carrying these resistance genes, 82% belonged to SpaA group 2-N101S (n = 22) or 2-N101S-I257L (n = 2). In addition, 79% (n = 23) of these 29 isolates belonged to MLST group ST 5. Our results illustrate that Australia appears to have a unique diversity of E. rhusiopathiae isolates in pig production industries within the wider global context of isolates.

Association of Environmental Factors with Seasonal Intensity of Erysipelothrix rhusiopathiae Seropositivity among Arctic Caribou

Several caribou (Rangifer tarandus) populations have been declining concurrently with increases in infectious diseases in the Arctic. Erysipelothrix rhusiopathiae, a zoonotic bacterium, was first described in 2015 as a notable cause of illness and death among several Arctic wildlife species. We investigated epidemiologic and environmental factors associated with the seroprevalence of E. rhusiopathiae in the Arctic and found that seropositivity was highest during warmer months, peaking in September, and was highest among adult males. Summer seroprevalence increases tracked with the oestrid index from the previous year, icing and snowing events, and precipitation from the same year but decreased with growing degree days in the same year. Seroprevalence of E. rhusiopathiae varied more during the later years of the study. Our findings provide key insights into the influence of environmental factors on disease prevalence that can be instrumental for anticipating and mitigating diseases associated with climate change among Arctic wildlife and human populations.

A SCOPING REVIEW OF THE RANGIFER TARANDUS INFECTIOUS DISEASE LITERATURE: GAP BETWEEN INFORMATION AND APPLICATION

The role and impact of infectious diseases in wildlife population dynamics are increasingly recognized, yet disease information is variably incorporated into wildlife management frameworks. This discrepancy is particularly relevant for Rangifer tarandus (caribou or reindeer), a keystone circumarctic species experiencing widespread population declines. The primary objective of this review was to characterize the available peer-reviewed literature on infectious diseases of Rangifer by using a scoping review methodology. Three databases of peer-reviewed literature-Web of Science, BIOSIS previews, and Scopus-were searched and 695 articles met the criteria for initial review. After screening for relevance and language, 349 articles, published between 1967 and 2020, remained. More than half of the excluded articles (181/346; 52%) were left out because they were not published in English; the majority of these excluded articles (120) were in Russian. From the 349 included articles, 137 (39%) pertained to wild (as opposed to semidomesticated or captive) Rangifer populations. Articles on infectious disease in wild Rangifer were published in 40 different journals across various disciplines; the most common journals were disease and parasitology oriented, accounting for 55% of included articles. Most studies were descriptive (87%), followed by experimental (9%). Of the pathogen taxa investigated, helminths were the most common, comprising 35% of articles. Rangifer subspecies were not equally represented in the literature, with barren-ground caribou (R. t. groenlandicus; n=40) and woodland caribou (R. t. caribou; n=39) having the greatest abundance and diversity of infectious disease information available. Few studies explicitly examined individual or population-level impacts of disease, or related disease to vital population rates, and only 27 articles explicitly related results to management or conservation. Findings from this review highlight an unbalanced distribution of studies across Rangifer ecotypes, a preference for dissemination in disease-specialized publication venues, and an opportunity for investigating population-level impacts that may be more readily integrated into caribou conservation frameworks.

Population genomics of Bacillus anthracis from an anthrax hyperendemic area reveals transmission processes across spatial scales and unexpected within-host diversity

Genomic sequencing has revolutionized our understanding of bacterial disease epidemiology, but remains underutilized for zoonotic pathogens in remote endemic settings. Anthrax, caused by the spore-forming bacterium Bacillus anthracis, remains a threat to human and animal health and rural livelihoods in low- and middle-income countries. While the global genomic diversity of B. anthracis has been well-characterized, there is limited information on how its populations are genetically structured at the scale at which transmission occurs, critical for understanding the pathogen's evolution and transmission dynamics. Using a uniquely rich dataset, we quantified genome-wide SNPs among 73 B. anthracis isolates derived from 33 livestock carcasses sampled over 1 year throughout the Ngorongoro Conservation Area, Tanzania, a region hyperendemic for anthrax. Genome-wide SNPs distinguished 22 unique B. anthracis genotypes (i.e. SNP profiles) within the study area. However, phylogeographical structure was lacking, as identical SNP profiles were found throughout the study area, likely the result of the long and variable periods of spore dormancy and long-distance livestock movements. Significantly, divergent genotypes were obtained from spatio-temporally linked cases and even individual carcasses. The high number of SNPs distinguishing isolates from the same host is unlikely to have arisen during infection, as supported by our simulation models. This points to an unexpectedly wide transmission bottleneck for B. anthracis, with an inoculum comprising multiple variants being the norm. Our work highlights that inferring transmission patterns of B. anthracis from genomic data will require analytical approaches that account for extended and variable environmental persistence, as well as co-infection.

Qiviut cortisol is associated with metrics of health and other intrinsic and extrinsic factors in wild muskoxen (Ovibos moschatus)

Glucocorticoid (GC) levels are increasingly and widely used as biomarkers of hypothalamic-pituitary-adrenal (HPA) axis activity to study the effects of environmental changes and other perturbations on wildlife individuals and populations. However, identifying the intrinsic and extrinsic factors that influence GC levels is a key step in endocrinology studies to ensure accurate interpretation of GC responses. In muskoxen, qiviut (fine woolly undercoat hair) cortisol concentration is an integrative biomarker of HPA axis activity over the course of the hair's growth. We gathered data from 219 wild muskoxen harvested in the Canadian Arctic between October 2015 and May 2019. We examined the relationship between qiviut cortisol and various intrinsic (sex, age, body condition and incisor breakage) and extrinsic biotic factors (lungworm and gastrointestinal parasite infections and exposure to bacteria), as well as broader non-specific landscape and temporal features (geographical location, season and year). A Bayesian approach, which allows for the joint estimation of missing values in the data and model parameters estimates, was applied for the statistical analyses. The main findings include the following: (i) higher qiviut cortisol levels in males than in females; (ii) inter-annual variations; (iii) higher qiviut cortisol levels in a declining population compared to a stable population; (iv) a negative association between qiviut cortisol and marrow fat percentage; (v) a relationship between qiviut cortisol and the infection intensity of the lungworm Umingmakstrongylus pallikuukensis, which varied depending on the geographical location; and (vi) no association between qiviut cortisol and other pathogen exposure/infection intensity metrics. This study confirmed and further identified important sources of variability in qiviut cortisol levels, while providing important insights on the relationship between GC levels and pathogen exposure/infection intensity. Results support the use of qiviut cortisol as a tool to monitor temporal changes in HPA axis activity at a population level and to inform management and conservation actions.

Contagious Ecthyma Dermatitis as a Portal of Entry for Erysipelothrix rhusiopathiae in Muskoxen (Ovibos moschatus) of the Canadian Arctic

Erysipelothrix rhusiopathiae was detected immunohistochemically in contagious ecthyma (orf virus) dermatitis in two muskoxen (Ovibos moschatus), harvested and found dead in 2014 and 2015, respectively, on Victoria Island, Canada. This may help target further research on E. rhusiopathiae epidemiology and mechanisms of infection in muskoxen, recently associated with widespread mortalities in Canada's Arctic.

Zoonotic pathogens in wild muskoxen (Ovibos moschatus) and domestic sheep (Ovis aries) from Greenland

The present study aimed to estimate the prevalence of zoonotic pathogens Giardia duodenalis, Cryptosporidium spp., Toxoplasma gondii and Erysipelothrix in muskoxen (Ovibos moschatus) and sheep (Ovis aries) from Greenland. In 2017 and 2018, faecal samples were collected from wild muskoxen from three distinct populations (Zackenberg, Kangerlussuaq, and Ivittuut) and from domestic sheep from southwest Greenland. Blood samples were collected from muskoxen from Kangerlussuaq and Ivittuut and from sheep. Faecal samples were tested for specific DNA of G. duodenalis and Cryptosporidium spp., and blood samples were tested for antibodies against T. gondii and Erysipelothrix. The estimated prevalence of G. duodenalis was 0% (0/58), 17% (7/41) and 0% (0/55) in muskoxen from Zackenberg, Kangerlussuaq and Ivittuut, respectively, and 37% (16/43) in sheep. The estimated prevalence of Cryptosporidium was 0% (0/58), 2% (1/41), 7% (4/55) in muskoxen from Zackenberg, Kangerlussuaq, Ivittuut, respectively, and 2% (1/43) in sheep. Neither Giardia nor Cryptosporidium were detected in winter samples (0/78). Of the positive samples, Giardia from one muskox sample only was successfully typed as G. duodenalis assemblage A, and Cryptosporidium from two muskoxen was successfully typed as C. parvum, subtype IIdA20G1e. The estimated T. gondii seroprevalence was 2% (1/44) and 0% (0/8) in muskoxen from Kangerlussuaq and Ivittuut, respectively, and 1% (1/155) in sheep. The estimated Erysipelothrix seroprevalence was 2% (1/45) and 13% (1/8) in muskoxen from Kangerlussuaq and Ivittuut, respectively, and 7% (10/150) in sheep. The results of this study add to the scarce knowledge on zoonotic pathogens in the Arctic.

Corynebacterium freneyi Bacterial Septicemia Secondary to Contagious Ecthyma in a Wild Muskox (Ovibos moschatus)

A wild muskox (Ovibos moschatus) with dermatitis typical of contagious ecthyma had secondary bacterial septicemia with Corynebacterium freneyi that included laminitis, hepatitis, and suppurative encephalitis. This case supports the association between orf virus infection and fatal secondary infections, which may have contributed to population declines on Victoria Island, Canada.

Implications of Zoonoses From Hunting and Use of Wildlife in North American Arctic and Boreal Biomes: Pandemic Potential, Monitoring, and Mitigation

The COVID-19 pandemic has re-focused attention on mechanisms that lead to zoonotic disease spillover and spread. Commercial wildlife trade, and associated markets, are recognized mechanisms for zoonotic disease emergence, resulting in a growing global conversation around reducing human disease risks from spillover associated with hunting, trade, and consumption of wild animals. These discussions are especially relevant to people who rely on harvesting wildlife to meet nutritional, and cultural needs, including those in Arctic and boreal regions. Global policies around wildlife use and trade can impact food sovereignty and security, especially of Indigenous Peoples. We reviewed known zoonotic pathogens and current risks of transmission from wildlife (including fish) to humans in North American Arctic and boreal biomes, and evaluated the epidemic and pandemic potential of these zoonoses. We discuss future concerns, and consider monitoring and mitigation measures in these changing socio-ecological systems. While multiple zoonotic pathogens circulate in these systems, risks to humans are mostly limited to individual illness or local community outbreaks. These regions are relatively remote, subject to very cold temperatures, have relatively low wildlife, domestic animal, and pathogen diversity, and in many cases low density, including of humans. Hence, favorable conditions for emergence of novel diseases or major amplification of a spillover event are currently not present. The greatest risk to northern communities from pathogens of pandemic potential is via introduction with humans visiting from other areas. However, Arctic and boreal ecosystems are undergoing rapid changes through climate warming, habitat encroachment, and development; all of which can change host and pathogen relationships, thereby affecting the probability of the emergence of new (and re-emergence of old) zoonoses. Indigenous leadership and engagement in disease monitoring, prevention and response, is vital from the outset, and would increase the success of such efforts, as well as ensure the protection of Indigenous rights as outlined in the United Nations Declaration on the Rights of Indigenous Peoples. Partnering with northern communities and including Indigenous Knowledge Systems would improve the timeliness, and likelihood, of detecting emerging zoonotic risks, and contextualize risk assessments to the unique human-wildlife relationships present in northern biomes.

Exposure of the Gulf of St. Lawrence grey seal Halichoerus grypus population to potentially zoonotic infectious agents

The population of grey seals Halichoerus grypus in Canadian waters is currently used as a commercial source of meat for human consumption. As with domestic livestock, it is important to understand the occurrence in these seals of infectious agents that may be of public health significance and thus ensure appropriate measures are in place to avoid zoonotic transmission. This study examined the prevalence of antibodies against Brucella spp., Erysipelothrix rhusiopathiae, 6 serovars of Leptospira interrogans, and Toxoplasma gondii in 59 grey seals and determined by polymerase chain reaction (PCR) the presence of these potentially zoonotic agents in specific organs and tissues of seropositive animals. The presence of encysted Trichinella spp. larvae was also investigated by digestion of tongue, diaphragm and other muscle samples, but none were detected. Seroprevalence against Brucella spp. and E. rhusiopathiae was low (5 and 3%, respectively). All 59 seals tested had antibodies against L. interrogans, but no carrier of this bacterium was detected by PCR. Seroprevalence against T. gondii was 53%, and DNA of this protozoan was detected by PCR in 11/30 (37%) seropositive animals. Standard sanitary measures mandatory for commercialization of meat products for human consumption should greatly reduce the potential for exposure to these infectious agents. However, special consideration should be given to freezing seal meat for at least 3 d to ensure destruction of tissue cysts of T. gondii.

Serovars and SpaA Types of Erysipelothrix rhusiopathiae Isolated from Pigs in Japan from 2012 to 2019

Erysipelothrix rhusiopathiae causes swine erysipelas (SE), which results in considerable economic loss on pig farms. During SE outbreaks that occurred sporadically from 2008 to 2011 in Japan, new E. rhusiopathiae strains were isolated with a specific surface protective antigen (Spa)A protein characterized by methionine at position 203 and isoleucine at position 257 (M203/I257 SpaA type). To determine whether strains with the M203/I257 SpaA type are still prevalent in Japan, we collected 79 strains of E. rhusiopathiae from pigs showing various SE symptoms from 2012 to 2019 and classified them based on serovar typing, spaA gene sequence analysis, and lineage typing. We found that the majority of recent E. rhusiopathiae strains (59/79) belonged to the serovar 1a strain, and that the M203/I257 SpaA type (56/59) was predominant continuing from 2008 to 2011. Furthermore, serovar 1a strains with IVb-1 and IVb-2 lineages that had been isolated in specific regions of Japan were no longer local but were found across Japan. The pathogenicity of recent isolates tested in mice was not significantly changed when compared to that of previously isolated strains. Our results suggest that recent SE outbreaks were not due to changes in the SpaA protein or to altered virulence of E. rhusiopathiae but were rather caused by the persistent presence of E. rhusiopathiae with the M203/I257 SpaA type.

Development of a Multiplex PCR-Based Assay for Rapid Serotyping of Erysipelothrix Species

The Gram-positive bacterium Erysipelothrix rhusiopathiae is a zoonotic pathogen that causes erysipelas in a wide range of mammalian and avian species. Historically, E. rhusiopathiae has been differentiated from other Erysipelothrix species by serotyping. Among 28 serovars of Erysipelothrix species, specific serovars, namely, 1a, 1b, and 2 of E. rhusiopathiae, are associated mainly with the disease in pigs, poultry, and humans; however, other serovar strains are often simultaneously isolated from diseased and healthy animals, indicating the importance of isolate serotyping for epidemiology. The traditional serotyping protocol, which uses heat-stable peptidoglycan antigens and type-specific rabbit antisera in an agar-gel precipitation test, is time-consuming and labor-intensive. To develop a rapid serotyping scheme, we analyzed sequences of the 12- to 22-kb chromosomal region, which corresponds to the genetic region responsible for virulence of serovar 1a and 2 strains of E. rhusiopathiae, of the 28 serovars of Erysipelothrix species. We confirmed that the serovar 13 strain lacks the genomic region and that some serovar strains possess very similar or the same genetic structure, prohibiting differentiation of the serovars. We created 4 multiplex PCR sets allowing the simultaneous detection and differentiation of the majority of Erysipelothrix serovars. Together with a previously reported multiplex PCR that can differentiate serovars 1a, 1b, 2, and 5, the multiplex PCR-based assay developed in this study covers all but one (serovar 13) of the reported serovars of Erysipelothrix species and should be a valuable tool for etiological as well as epidemiological studies of Erysipelothrix infections.

Novel insights into serodiagnosis and epidemiology of Erysipelothrix rhusiopathiae, a newly recognized pathogen in muskoxen (Ovibos moschatus)

BACKGROUND

Muskoxen are a key species of Arctic ecosystems and are important for food security and socio-economic well-being of many Indigenous communities in the Arctic and Subarctic. Between 2009 and 2014, the bacterium Erysipelothrix rhusiopathiae was isolated for the first time in this species in association with multiple mortality events in Canada and Alaska, raising questions regarding the spatiotemporal occurrence of the pathogen and its potential impact on muskox populations.

MATERIALS AND METHODS

We adapted a commercial porcine E. rhusiopathiae enzyme-linked immunosorbent assay to test 958 blood samples that were collected from muskoxen from seven regions in Alaska and the Canadian Arctic between 1976 and 2017. The cut-off between negative and positive results was established using mixture-distribution analysis, a data-driven approach. Based on 818 samples for which a serological status could be determined and with complete information, we calculated trends in sample seroprevalences in population time-series and compared them with population trends in the investigated regions.

RESULTS

Overall, 219/818 (27.8%, 95% Confidence Interval: 24.7-31.0) samples were classified as positive for exposure to E. rhusiopathiae. There were large variations between years and regions. Seropositive animals were found among the earliest serum samples tested; 1976 in Alaska and 1991 in Canada. In Alaskan muskoxen, sample seroprevalence increased after 2000 and, in two regions, peak seroprevalences occurred simultaneously with population declines. In one of these regions, concurrent unusual mortalities were observed and E. rhusiopathiae was isolated from muskox carcasses. In Canada, there was an increase in sample seroprevalence in two muskox populations following known mortality events that had been attributed to E. rhusiopathiae.

CONCLUSION

Our results indicate widespread exposure of muskoxen to E. rhusiopathiae in western Canada and Alaska. Although not new to the Arctic, we documented an increased exposure to the pathogen in several regions concurrent with population declines. Understanding causes for the apparent increased occurrence of this pathogen and its association with large scale mortality events for muskoxen is critical to evaluate the implications for wildlife and wildlife-dependent human populations in the Arctic.

Erysipelothrix Spp.: Past, Present, and Future Directions in Vaccine Research

Erysipelothrix spp. comprise a group of small Gram-positive bacteria that can infect a variety of hosts including mammals, fish, birds, reptiles and insects. Among the eight Erysipelothrix species that have been described to date, only Erysipelothrix rhusiopathiae plays a major role in farmed livestock where it is the causative agent of erysipelas. E. rhusiopathiae also has zoonotic potential and can cause erysipeloid in humans with a clear occupational link to meat and fish industries. While there are 28 known Erysipelothrix serovars, over 80% of identified isolates belong to serovars 1 or 2. Vaccines to protect pigs against E. rhusiopathiae first became available in 1883 as a response to an epizootic of swine erysipelas in southern France. The overall vaccine repertoire was notably enlarged between the 1940s and 1960s following major outbreaks of swine erysipelas in the Midwest USA and has changed little since. Traditionally, E. rhusiopathiae serovar 1a or 2 isolates were inactivated (bacterins) or attenuated and these types of vaccines are still used today on a global basis. E. rhusiopathiae vaccines are most commonly used in pigs, poultry, and sheep where the bacterium can cause considerable economic losses. In addition, erysipelas vaccination is also utilized in selected vulnerable susceptible populations, such as marine mammals in aquariums, which are commonly vaccinated at regular intervals. While commercially produced erysipelas vaccines appear to provide good protection against clinical disease, in recent years there has been an increase in perceived vaccine failures in farmed animals, especially in organic outdoor operations. Moreover, clinical erysipelas outbreaks have been reported in animal populations not previously considered at risk. This has raised concerns over a possible lack of vaccine protection across various production species. This review focuses on summarizing the history and the present status of E. rhusiopathiae vaccines, the current knowledge on protection including surface antigens, and also provides an outlook into future directions for vaccine development.

Genomic and Immunogenic Protein Diversity of Erysipelothrix rhusiopathiae Isolated From Pigs in Great Britain: Implications for Vaccine Protection

Erysipelas, caused by the bacterium Erysipelothrix rhusiopathiae, is re-emerging in swine and poultry production systems worldwide. While the global genomic diversity of this species has been characterized, how much of this genomic and functional diversity is maintained at smaller scales is unclear. Specifically, while several key immunogenic surface proteins have been identified for E. rhusiopathiae, little is known about their presence among field strains and their divergence from vaccines, which could result in vaccine failure. Here, a comparative genomics approach was taken to determine the diversity of E. rhusiopathiae strains in pigs in Great Britain over nearly three decades, as well as to assess the field strains’ divergence from the vaccine strain most commonly used in British pigs. In addition, the presence/absence and variability of 13 previously described immunogenic surface proteins was determined, including SpaA which is considered a key immunogen. We found a high diversity of E. rhusiopathiae strains in British pigs, similar to the situation described in European poultry but in contrast to swine production systems in Asia. Of the four clades of E. rhusiopathiae found globally, three were represented among British pig isolates, with Clade 2 being the most common. All British pig isolates had one amino acid difference in the immunoprotective domain of the SpaA protein compared to the vaccine strain. However, we were able to confirm using in silico structural protein analyses that this difference is unlikely to compromise vaccine protection. Of 12 other known immunogenic surface proteins of E. rhusiopathiae examined, 11 were found to be present in all British pig isolates and the vaccine strain, but with highly variable degrees of conservation at the amino acid sequence level, ranging from 0.3 to 27% variant positions. Moreover, the phylogenetic incongruence of these proteins suggests that horizontal transfer of genes encoding for antigens is commonplace for this bacterium. We hypothesize that the sequence variants in these proteins could be responsible for differences in the efficacy of the immune response. Our results provide the necessary basis for testing this hypothesis through in vitro and in vivo studies.

Muskox status, recent variation, and uncertain future

Muskoxen (Ovibos moschatus) are an integral component of Arctic biodiversity. Given low genetic diversity, their ability to respond to future and rapid Arctic change is unknown, although paleontological history demonstrates adaptability within limits. We discuss status and limitations of current monitoring, and summarize circumpolar status and recent variations, delineating all 55 endemic or translocated populations. Acknowledging uncertainties, global abundance is ca 170 000 muskoxen. Not all populations are thriving. Six populations are in decline, and as recently as the turn of the century, one of these was the largest population in the world, equaling ca 41% of today's total abundance. Climate, diseases, and anthropogenic changes are likely the principal drivers of muskox population change and result in multiple stressors that vary temporally and spatially. Impacts to muskoxen are precipitated by habitat loss/degradation, altered vegetation and species associations, pollution, and harvest. Which elements are relevant for a specific population will vary, as will their cumulative interactions. Our summaries highlight the importance of harmonizing existing data, intensifying long-term monitoring efforts including demographics and health assessments, standardizing and implementing monitoring protocols, and increasing stakeholder engagement/contributions.

Conservation Genomics in a Changing Arctic

Although logistically challenging to study, the Arctic is a bellwether for global change and is becoming a model for questions pertinent to the persistence of biodiversity. Disruption of Arctic ecosystems is accelerating, with impacts ranging from mixing of biotic communities to individual behavioral responses. Understanding these changes is crucial for conservation and sustainable economic development. Genomic approaches are providing transformative insights into biotic responses to environmental change, but have seen limited application in the Arctic due to a series of limitations. To meet the promise of genome analyses, we urge rigorous development of biorepositories from high latitudes to provide essential libraries to improve the conservation, monitoring, and management of Arctic ecosystems through genomic approaches.

An unusual case of Erysipelothrix rhusiopathiae prosthetic joint infection from the Canadian Arctic: whole genome sequencing unable to identify a zoonotic source

Background

Erysipelothrix rhusiopathiae is a zoonotic pathogen that causes erysipeloid and is most frequently associated with exposure to domestic swine. Infection of native and prosthetic joints is a rarely reported manifestation.

Case presentation

We describe a case of E. rhusiopathiae prosthetic joint infection in a woman with a history of exposure to wild animals in the Canadian Arctic. Patient management involved a 1-stage surgical revision exchange with an antibiotic impregnated cement spacer and 6 weeks of intravenous penicillin G followed by 6 weeks of oral amoxicillin. Ten previously reported cases of E. rhusiopathiae joint infection are reviewed. Recent increases in mortality due to infection with this organism among host animal populations in the Canadian Arctic have generated concern regarding a potential increase in human infections. However, whole genome sequencing (WGS) of the organism was unable to identify a zoonotic origin for this case.

Conclusions

Consideration should be given to E. rhusiopathiae as a cause of joint infections if the appropriate epidemiologic and host risk factors exist. Expanded use of WGS in other potential animal hosts and environmental sources may provide important epidemiologic information in determining the source of human infections.

Identification of serovar 1a, 1b, 2, and 5 strains of Erysipelothrix rhusiopathiae by a conventional gel-based PCR

Among the four species of the genus Erysipelothrix, Erysipelothrix rhusiopathiae is the main species that causes disease in swine and poultry and has also been isolated from human patients. Recently, E. rhusiopathiae infections in domesticated animals have increased in many countries and are also the cause of emerging wildlife disease in arctic and boreal ecosystems. Historically, E. rhusiopathiae has been differentiated from other Erysipelothrix species by their serovars, which are determined based on cell wall antigens. Serotyping of Erysipelothrix is important, as specific E. rhusiopathiae serovars (1a, 1b, and 2) are associated with disease in pigs, poultry, and humans. However, serotyping is laborious and time-consuming and requires a full set of serovar reference strains and strain-specific antiserum. In this study, to develop a conventional gel-based PCR assay that can detect the main disease-associated serovars of E. rhusiopathiae, the draft genome sequences of E. rhusiopathiae strains of serovars 1a, 1b, 2, and 5, the last of which is often isolated from wild animals, were analyzed. Primers were designed based on the serovar-specific sequences of the strains and tested for field strains isolated from extensive origins. Among two hundred and ninety-seven isolates of various serovar strains of E. rhusiopathiae and other Erysipelothrix species, the PCR assay identified serovar 1a, 1b, 2, and 5 strains of E. rhusiopathiae. This conventional gel-based PCR assay should be useful for serovar surveillance of E. rhusiopathiae isolates in domesticated and wild animals as well as in humans.

A Comparative Study on the Faecal Bacterial Community and Potential Zoonotic Bacteria of Muskoxen (Ovibos moschatus) in Northeast Greenland, Northwest Greenland and Norway

Muskoxen (Ovibos moschatus) are ruminants adapted to a high-fibre diet. There is increasing interest in the role that gut microbes play in the digestion and utilization of these specialized diets but only limited data available on the gut microbiome of high-Arctic animals. In this study, we metabarcoded the 16S rRNA region of faecal samples from muskoxen of Northeast Greenland, Northwest Greenland and Norway, and quantified the effects of physiological and temporal factors on bacterial composition. We found significant effects of body mass, year of sampling and location on the gut bacterial communities of North East Greenland muskoxen. These effects were however dwarfed by the effects of location, emphasizing the importance of the local ecology on the gut bacterial community. Habitat alterations and rising temperatures may therefore have a considerable impact on muskoxen health and reproductive success. Moreover, muskoxen are hunted and consumed in Greenland, Canada and Alaska; therefore, this study also screened for potential zoonoses of food safety interest. A total of 13 potentially zoonotic genera were identified, including the genera Erysipelothrix and Yersinia implicated in recent mass die-offs of the muskoxen themselves.

Multiplex PCR assay for the simultaneous detection and differentiation of clonal lineages of Erysipelothrix rhusiopathiae serovar 1a strains currently circulating in Japan

The species Erysipelothrix rhusiopathiae displays genetic heterogeneity; however, E. rhusiopathiae serovar 1a strains currently circulating in Japan exhibit remarkably low levels of genetic diversity and group into clonal sublineages of Lineage IVb (IVb-1 and IVb-2). In the present study, based on whole genome sequencing data, we designed primers for a multiplex PCR assay to simultaneously detect and differentiate the sublineages of E. rhusiopathiae strains. Among the one hundred and twenty-seven isolates of various serovar strains, including isolates from a wide range of hosts and geographic origins, the PCR assay could successfully detect and differentiate the serovar 1a strains belonging to the sublineages.

Clonal Lineages of Erysipelothrix rhusiopathiae Responsible for Acute Swine Erysipelas in Japan Identified by Using Genome-Wide Single-Nucleotide Polymorphism Analysis

Erysipelothrix rhusiopathiae causes swine erysipelas, an important infectious disease in the swine industry. In Japan, the incidence of acute swine erysipelas due to E. rhusiopathiae serovar 1a has recently increased markedly. To study the genetic relatedness of the strains from the recent cases, we analyzed 34 E. rhusiopathiae serovar 1a swine isolates collected between 1990 and 2011 and further investigated the possible association of the live Koganei 65-0.15 vaccine strain (serovar 1a) with the increase in cases. Pulsed-field gel electrophoresis analysis revealed no marked variation among the isolates; however, sequencing analysis of a hypervariable region in the surface-protective antigen A gene (spaA) revealed that the strains isolated after 2007 exhibited the same spaA genotype and could be differentiated from older strains. Phylogenetic analysis based on genome-wide single-nucleotide polymorphisms (SNPs) revealed that the Japanese strains examined were closely related, showing a relatively small number of SNPs among them. The strains were classified into four major lineages, with Koganei 65-0.15 (lineage III) being phylogenetically separated from the other three lineages. The strains isolated after 2007 and the two older strains constituted one major lineage (lineage IV) with a specific spaA genotype (M203/I257-SpaA), while the recent isolates were further divided into two geographic groups. The remaining older isolates belonged to either lineage I, with the I203/L257-SpaA type, or lineage II, with the I203/I257-SpaA type. These results indicate that the recent increased incidence of acute swine erysipelas in Japan is associated with two sublineages of lineage IV, which have independently evolved in two different geographic regions.IMPORTANCE Using large-scale whole-genome sequence data from Erysipelothrix rhusiopathiae isolates from a wide range of hosts and geographic origins, a recent study clarified the existence of three distinct clades (clades 1, 2, and 3) that are found across multiple continents and host species, representing both livestock and wildlife, and an "intermediate" clade between clade 2 and the dominant clade 3 within the species. In this study, we found that the E. rhusiopathiae Japanese strains examined exhibited remarkably low levels of genetic diversity and confirmed that all of the Japanese and Chinese swine isolates examined in this study belong to clonal lineages within the intermediate clade. We report that spaA genotyping of E. rhusiopathiae strains is a practical alternative to whole-genome sequencing analysis of the E. rhusiopathiae isolates from eastern Asian countries.