Molecular and phylogenetic analyses of Salmonella Gallinarum trace the origin and diversification of recent outbreaks of fowl typhoid in poultry farms

Salmonella enterica Serovar Gallinarum Biovars Pullorum and Gallinarum in Poultry: Review of Pathogenesis, Antibiotic Resistance, Diagnosis and Control in the Genomic Era

Salmonella enterica subsp. enterica serovar Gallinarum (SG) has two distinct biovars, Pullorum and Gallinarum. They are bacterial pathogens that exhibit host specificity for poultry and aquatic birds, causing severe systemic diseases known as fowl typhoid (FT) and Pullorum disease (PD), respectively. The virulence mechanisms of biovars Gallinarum and Pullorum are multifactorial, involving a variety of genes and pathways that contribute to their pathogenicity. In addition, these serovars have developed resistance to various antimicrobial agents, leading to the emergence of multidrug-resistant strains. Due to their economic and public health significance, rapid and accurate diagnosis is crucial for effective control and prevention of these diseases. Conventional methods, such as bacterial culture and serological tests, have been used for screening and diagnosis. However, molecular-based methods are becoming increasingly important due to their rapidity, high sensitivity, and specificity, opening new horizons for the development of innovative approaches to control FT and PD. The aim of this review is to highlight the current state of knowledge on biovars Gallinarum and Pullorum, emphasizing the importance of continued research into their pathogenesis, drug resistance and diagnosis to better understand and control these pathogens in poultry farms.

Salmonella gallinarum strains from outbreaks of fowl typhoid fever in Southern Africa closely related to SG9R vaccines

Introduction

Salmonella enterica subspecies enterica serovar Gallinarum biovar Gallinarum (SG) is associated with fowl typhoid fever, and the attenuated rough strain SG9R is widely used as a vaccine in many regions. Reversion to virulence of vaccine strains was suspected as the cause during recent fowl typhoid fever outbreaks in poultry in South Africa and Eswatini.

Methods

To compare nine field isolates with global wild-type SG9 strains and the two commercial SG9R vaccines in use, Nobilis^® SG9R and Cevac^®-SG, we used whole-genome comparison with single-nucleotide polymorphism (SNP) detection.

Results

SNP phylogenic analysis showed that all the southern African field isolates were more closely related to the vaccine strains than wild-type SG9 strains. Furthermore, SNPs in the pyruvate dehydrogenase (aceE) and/or lipopolysaccharide 1,2-glucosyltransferase (rfaJ) genes, which are known markers of attenuation, were found in four of the field isolates along with intact spv, SPI-1, and SPI-2 gene clusters, providing conclusive evidence that these four isolates were originally vaccine strains that reverted to virulence. Five other field isolates lacked the SG9R attenuation markers, but variant analysis identified an SNP in the yihX gene, insertions in the ybjX and hydH genes, and deletions in the ftsK and sadA genes that were shared between the field isolates and vaccine strains but absent in wild-type SG9, indicating that these field isolates were also likely revertant vaccines.

Discussion

Overall, this study highlights different mechanisms of reversion of two commercial vaccines, where virulence caused by field isolates closely related to the Nobilis^® SG9R vaccine was associated with the restoration of intact virulence gene clusters, and those derived from the Cevac^®-SG vaccine were characterized by point mutations resulting in restored aceE and rfaJ genes. A possible new marker of attenuation was identified as a point mutation in the yihX gene, as well as four new candidate genes that could potentially be used to distinguish current vaccine strains from wild-type strains using PCR assays.

Genomic Characterization and Genetic Profiles of Salmonella Gallinarum Strains Isolated from Layers with Fowl Typhoid in Colombia

Salmonella Gallinarum (SG) is the causative agent of fowl typhoid (FT), a disease that is harmful to the poultry industry. Despite sanitation and prophylactic measures, this pathogen is associated with frequent disease outbreaks in developing countries, causing high morbidity and mortality. We characterized the complete genome sequence of Colombian SG strains and then performed a comparative genome analysis with other SG strains found in different regions worldwide. Eight field strains of SG plus a 9R-derived vaccine were subjected to whole-genome sequencing (WGS) and bioinformatics analysis, and the results were used for subsequent molecular typing; virulome, resistome, and mobilome characterization; and a comparative genome study. We identified 26 chromosome-located resistance genes that mostly encode efflux pumps, and point mutations were found in gyrase genes (gyrA and gyrB), with the gyrB mutation S464T frequently found in the Colombian strains. Moreover, we detected 135 virulence genes, mainly in 15 different Salmonella pathogenicity islands (SPIs). We generated an SPI profile for SG, including C63PI, CS54, ssaD, SPI-1, SPI-2, SPI-3, SPI-4, SPI-5, SPI-6, SPI-9, SPI-10, SPI-11, SPI-12, SPI-13, and SPI-14. Regarding mobile genetic elements, we found the plasmids Col(pHAD28) and IncFII(S) in most of the strains and 13 different prophage sequences, indicating a frequently obtained profile that included the complete phage Gifsy_2 and incomplete phage sequences resembling Escher_500465_2, Shigel_SfIV, Entero_mEp237, and Salmon_SJ46. This study presents, for the first time, the genomic content of Colombian SG strains and a profile of the genetic elements frequently found in SG, which can be further studied to clarify the pathogenicity and evolutionary characteristics of this serotype.

Serotype Screening of Salmonella enterica Subspecies I by Intergenic Sequence Ribotyping (ISR): Critical Updates

(1) Background: Foodborne illness from Salmonella enterica subspecies I is most associated with approximately 32 out of 1600 serotypes. While whole genome sequencing and other nucleic acid-based methods are preferred for serotyping, they require expertise in bioinformatics and often submission to an external agency. Intergenic Sequence Ribotyping (ISR) assigns serotype to Salmonella in coordination with information freely available at the National Center for Biotechnology Information. ISR requires updating because it was developed from 26 genomes while there are now currently 1804 genomes and 1685 plasmids. (2) Methods: Serotypes available for sequencing were analyzed by ISR to confirm primer efficacy and to identify any issues in application. Differences between the 2012 and 2022 ISR database were tabulated, nomenclature edited, and instances of multiple serotypes aligning to a single ISR were examined. (3) Results: The 2022 ISR database has 268 sequences and 40 of these were assigned new NCBI accession numbers that were not previously available. Extending boundaries of sequences resolved hdfR cross-alignment and reduced multiplicity of alignment for 37 ISRs. Comparison of gene cyaA sequences and some cell surface epitopes provided evidence that homologous recombination was potentially impacting results for this subset. There were 99 sequences that still had no match with an NCBI submission. (4) The 2022 ISR database is available for use as a serotype screening method for Salmonella enterica subspecies I. Finding that 36.9% of the sequences in the ISR database still have no match within the NCBI Salmonella enterica database suggests that there is more genomic heterogeneity yet to characterize.

Lacticaseibacillus rhamnosus alleviates intestinal inflammation and promotes microbiota-mediated protection against Salmonella fatal infections

The fatal impairment of the intestinal mucosal barrier of chicks caused by Salmonella significantly resulting economic losses in the modern poultry industry. Probiotics are recognized for beneficially influencing host immune responses, promoting maintenance of intestinal epithelial integrity, antagonistic activity against pathogenic microorganisms and health-promoting properties. Some basic studies attest to probiotic capabilities and show that Lacticaseibacillus rhamnosus could protect intestinal mucosa from injury in animals infected with Salmonella Typhimurium. However, the mechanisms underlying its protective effects in chicks are still not fully understood. Here, we used the chick infection model combined with histological, immunological, and molecular approaches to address this question. The results indicated that L. rhamnosus significantly reduced the diarrhea rate and increased the daily weight gain and survival rate of chicks infected with S. Typhimurium. Furthermore, we found that L. rhamnosus markedly improved the immunity of gut mucosa by reducing apoptotic cells, hence effectively inhibiting intestinal inflammation. Notably, pre-treatment chicks with L. rhamnosus balanced the expression of interleukin-1β and interleukin-18, moderated endotoxin and D-lactic acid levels, and expanded tight junction protein levels (Zonula occluden-1 and Claudin-1), enhanced the function of the intestinal mucosal epithelial cells. Additionally, investigations using full-length 16S rRNA sequencing also demonstrated that L. rhamnosus greatly weakened the adhesion of Salmonella, the mainly manifestation is the improvement of the diversity of intestinal microbiota in infected chicks. Collectively, these results showed the application of L. rhamnosus against Salmonella fatal infection by enhancing barrier integrity and the stability of the gut microbiota and reducing inflammation in new hatch chicks, offering new antibiotic alternatives for farming animals.

Emergence, Dissemination and Antimicrobial Resistance of the Main Poultry-Associated Salmonella Serovars in Brazil

Simple Summary

Salmonellosis is a human and animal disease caused by Salmonella, a bacterial genus classified into different species, subspecies, and serological variants (serovars) according to adaptation to one or more different hosts (animals and humans), pathogenicity profiles, and antigenic properties. Some specific Salmonella serovars can spread more easily in the enteric microbiota of avian species, often causing disease in birds and/or being transmitted to humans through food (such as chicken and eggs). Antimicrobial resistance (AMR) has also been reported in poultry-associated Salmonella isolates due to the widespread use of antimicrobials on farms. The availability of comprehensive data on the emergence and spread of Salmonella serovars, as well as their AMR profiles in farms and food products in Brazil (a major producer of poultry in the World), is necessary to understand their relevance in all avian production chains and also occurrence in poultry-derived foods. This article aims to provide an overview of the genus Salmonella and the main serovars that emerged in Brazilian poultry over time (Gallinarum, Typhimurium, Enteritidis, Heidelberg, and Minnesota), reviewing the scientific literature and suggesting more effective prevention and control for the future.

Abstract

Salmonella infects poultry, and it is also a human foodborne pathogen. This bacterial genus is classified into several serovars/lineages, some of them showing high antimicrobial resistance (AMR). The ease of Salmonella transmission in farms, slaughterhouses, and eggs industries has made controlling it a real challenge in the poultry-production chains. This review describes the emergence, dissemination, and AMR of the main Salmonella serovars and lineages detected in Brazilian poultry. It is reported that few serovars emerged and have been more widely disseminated in breeders, broilers, and layers in the last 70 years. Salmonella Gallinarum was the first to spread on the farms, remaining as a concerning poultry pathogen. Salmonella Typhimurium and Enteritidis were also largely detected in poultry and foods (eggs, chicken, turkey), being associated with several human foodborne outbreaks. Salmonella Heidelberg and Minnesota have been more widely spread in recent years, resulting in frequent chicken/turkey meat contamination. A few more serovars (Infantis, Newport, Hadar, Senftenberg, Schwarzengrund, and Mbandaka, among others) were also detected, but less frequently and usually in specific poultry-production regions. AMR has been identified in most isolates, highlighting multi-drug resistance in specific poultry lineages from the serovars Typhimurium, Heidelberg, and Minnesota. Epidemiological studies are necessary to trace and control this pathogen in Brazilian commercial poultry production chains.

Recent Evolution and Genomic Profile of Salmonella enterica Serovar Heidelberg Isolates from Poultry Flocks in Brazil

Salmonella enterica serovar Heidelberg is isolated from poultry-producing regions around the world. In Brazil, S. Heidelberg has been frequently detected in poultry flocks, slaughterhouses, and chicken meat. The goal of the present study was to assess the population structure, recent temporal evolution, and some important genetic characteristics of S. Heidelberg isolated from Brazilian poultry farms. Phylogenetic analysis of 68 S. Heidelberg genomes sequenced here and additional whole-genome data from NCBI demonstrated that all isolates from the Brazilian poultry production chain clustered into a monophyletic group, here called S. Heidelberg Brazilian poultry lineage (SH-BPL). Bayesian analysis defined the time of the most recent common ancestor (tMRCA) as 2004, and the overall population size (N_e) was constant until 2008, when an ∼10-fold N_e increase was observed until circa 2013. SH-BPL presented at least two plasmids with replicons ColpVC (n = 68; 100%), IncX1 (n = 66; 97%), IncA/C2 (n = 65; 95.5%), ColRNAI (n = 43; 63.2%), IncI1 (n = 32; 47%), ColMG828, Col156, IncHI2A, IncHI2, IncQ1, IncX4, IncY, and TrfA (each with n < 4; <4% each). Antibiotic resistance genes were found, with high frequencies of fosA7 (n = 68; 100%), mdf(A) (n = 68; 100%), tet(34) (n = 68; 100%), sul2 (n = 64; 94.1%), and blaCMY-2 (n = 56; 82.3%), along with an overall multidrug resistance (MDR) profile. Ten Salmonella pathogenicity islands (SPI1 to SPI5, SPI9, and SPI11 to SPI14) and 139 virulence genes were also detected. The SH-BPL profile was like those of other previous S. Heidelberg isolates from poultry around the world in the 1990s. In conclusion, the present study demonstrates the recent introduction (2004) and high level of dissemination of an MDR S. Heidelberg lineage in Brazilian poultry operations. IMPORTANCE S. Heidelberg is the most frequent serovar in several broiler farms from the main Brazilian poultry-producing regions. Therefore, avian-source foods (mainly chicken carcasses) commercialized in the country and exported to other continents are contaminated with this foodborne pathogen, generating several national and international economic losses. In addition, isolates of this serovar are usually resistant to antibiotics and can cause human invasive and septicemic infection, representing a public health concern. This study demonstrates the use of whole-genome sequencing (WGS) to obtain epidemiological information for one S. Heidelberg lineage highly spread among Brazilian poultry farms. This information will help to define biosecurity measures to control this important Salmonella serovar in Brazilian and worldwide poultry operations.

Comparative genome analysis of Salmonella enterica serovar Gallinarum biovars Pullorum and Gallinarum decodes strain specific genes

Salmonella enterica serovar Gallinarum biovar Pullorum (bvP) and biovar Gallinarum (bvG) are the etiological agents of pullorum disease (PD) and fowl typhoid (FT) respectively, which cause huge economic losses to poultry industry especially in developing countries including India. Vaccination and biosecurity measures are currently being employed to control and reduce the S. Gallinarum infections. High endemicity, poor implementation of hygiene and lack of effective vaccines pose challenges in prevention and control of disease in intensively maintained poultry flocks. Comparative genome analysis unravels similarities and dissimilarities thus facilitating identification of genomic features that aids in pathogenesis, niche adaptation and in tracing of evolutionary history. The present investigation was carried out to assess the genotypic differences amongst S.enterica serovar Gallinarum strains including Indian strain S. Gallinarum Sal40 VTCCBAA614. The comparative genome analysis revealed an open pan-genome consisting of 5091 coding sequence (CDS) with 3270 CDS belonging to core-genome, 1254 CDS to dispensable genome and strain specific genes i.e. singletons ranging from 3 to 102 amongst the analyzed strains. Moreover, the investigated strains exhibited diversity in genomic features such as virulence factors, genomic islands, prophage regions, toxin-antitoxin cassettes, and acquired antimicrobial resistance genes. Core genome identified in the study can give important leads in the direction of design of rapid and reliable diagnostics, and vaccine design for effective infection control as well as eradication. Additionally, the identified genetic differences among the S. enterica serovar Gallinarum strains could be used for bacterial typing, structure based inhibitor development by future experimental investigations on the data generated.

Optimized Detoxification of a Live Attenuated Vaccine Strain (SG9R) to Improve Vaccine Strategy against Fowl Typhoid

The live attenuated vaccine strain, SG9R, has been used against fowl typhoid worldwide, but it can revert to the pathogenic smooth strain owing to single nucleotide changes such as nonsense mutations in the rfaJ gene. As SG9R possesses an intact Salmonella plasmid with virulence genes, it exhibits dormant pathogenicity and can cause fowl typhoid in young chicks and stressed or immunocompromised brown egg-laying hens. To tackle these issues, we knocked out the rfaJ gene of SG9R (named Safe-9R) to eliminate the reversion risk and generated detoxified strains of Safe-9R by knocking out lpxL, lpxM, pagP, and phoP/phoQ genes to attenuate the virulence. Among the knockout strains, live ΔlpxL- (Dtx-9RL) and ΔlpxM-9R (Dtx-9RM) strains induced remarkably less expression of inflammatory cytokines in chicken macrophage cells, and oil emulsion (OE) Dtx-9RL did not cause body weight loss in chicks. Live Dtx-9RM exhibited efficacy against field strain challenge in one week without any bacterial re-isolation, while the un-detoxified strains showed the development of severe liver lesions and re-isolation of challenged strains. Thus, SG9R was optimally detoxified by knockout of lpxL and lpxM, and Dtx-9RL and Dtx-9RM might be applicable as OE and live vaccines, respectively, to prevent fowl typhoid irrespective of the age of chickens.

Genomic Characterization of Salmonella Minnesota Clonal Lineages Associated with Poultry Production in Brazil

Salmonella serotype Minnesota has been increasingly detected in Brazilian poultry farms and food products (chicken meat, eggs) in recent years. In addition, S. Minnesota isolates from poultry are generally resistant to several antibiotics and persistent in farm environments. The present study aimed to assess phylogenomic diversity of S. Minnesota isolates from the poultry production chain in Brazil. In total, 107 worldwide S. Minnesota whole genomes (including 12 from Brazil) were analyzed using a comparative approach. Phylogenetic analysis demonstrated two clades more related to poultry production in Brazil: S. Minnesota poultry lineages I and II (SM-PLI and SM-PLII). Phylodynamic analysis demonstrated that SM-PLI had a common ancestor in 1915, while SM-PLII originated circa 1971. SM-PLII encompassed a higher number of isolates and presented a recent increase in effective population size (mainly from 2009 to 2012). Plasmids IncA/C2 and ColRNA, antimicrobial resistance genes (aph(3')-Ia, blaCMY-2, qnrB19, sul2, and tet(A)) and mainly a virulence genetic cluster (including the yersiniabactin operon) were detected in isolates from SM-PLI and/or SM-PLII. This study demonstrates the dissemination of two distinct S. Minnesota lineages with high resistance to antibiotics and important virulence genetic clusters in Brazilian poultry farms.

Salmonella serotype assignment by sequencing analysis of intergenic regions of ribosomal RNA operons

Salmonella laboratorial detection is usually carried out by bacteriological culture and serological methods. Salmonella isolates are then classified into more than 2,650 serotypes with different somatic (O) and flagellar (H) antigenic combinations. More recently, DNA analysis methods were developed and applied for the identification of Salmonella serotypes, including intergenic spacer regions (ISRs) that separates DNA-encoding ribosomal subunits (rRNA gene) in Salmonella genomes. The present study aimed to evaluate the nucleotide diversity of the ISRs in 2 rRNA operons (rrnB and rrnH) for the assignment of Salmonella serotypes. A total of 63 Salmonella isolates (bacterial cultures) from 21 serotypes were obtained in the period of 2014 to 2017. DNA was extracted, and PCRs were used to detect the genus Salmonella and 4 important serotypes: Enteritidis, Gallinarum, Heidelberg, and Typhimurium. ISRs of the operons rrnB and rrnH were amplified by PCR and then sequenced. All sequence results were submitted to BLASTn search and were aligned in comparison to 72 Salmonella reference nucleotide sequences. The results demonstrated that 60 (95.2%) samples returned a sequence of the same serotype (determined by the traditional serological procedure) after searching in BLASTn and/or in the alignment with the reference sequences for both operons (rrnB and rrnH). These PCR-sequencing procedures had a general agreement index of 0.792 based on the Kappa analysis, 98.7% sensitivity value, 100% specificity, and 98.4% accuracy. Three different phylogenetic trees were generated from the alignments with the sequences of rrnH, rrnB, and rrnH plus rrnB and isolates clustered in specific branches according to the serotypes.

High-resolution diffusion pattern of human infections by Salmonella enterica serovar Napoli in Northern Italy explained through phylogeography

Salmonella enterica serovar Napoli (serovar Napoli) is an emerging cause of human salmonellosis in Northern Italy. No specific reservoirs of serovar Napoli have been identified in Italy, so far. However, the environment, especially surface waters, has been hypothesized as an important source of infection based on the observation that genotypically different clusters of serovar Napoli are detected in different geographical macro-areas. To further support the hypothesis of a spatially-restricted pattern of serovar Napoli diffusion, a spatial segregation of serovar Napoli lineages should be observed also at smaller geographical scale. However, classical genotyping techniques used for Salmonella, such as pulsed-field gel electrophoresis (PFGE), did not possess enough discriminatory power to highlight spatial clustering of serovar Napoli within the macro-areas. To this purpose, we performed phylogeographical analyses based on genome-wide single nucleotide polymorphisms to test whether spatio-temporal evolution patterns of serovar Napoli in Northern Italy could be recognized with high geographical resolution, i.e. at local level. Specifically, we analyzed the local spread of the main PFGE clonal group, responsible for more than 60% of human infections in the study area, that did not show any geographical differentiation by PFGE within Northern Italy, i.e. the macro-area considered in the study. Both discrete and continuous phylogeography highlighted the existence of two main geographically-restricted clades: a Southern clade corresponding to the Po Valley and a Northern clade corresponding to the Pre-Alps area. Furthermore, the phylogeographical analyses suggested that the most probable site of origin of the clone was in an area of the Po Valley at the confluence of the Po and Ticino rivers, one of the most important Italian wetlands. These findings provide further support to the hypothesis that environmental transmission may play an important role in the ecology of serovar Napoli.