Clustered Regularly Interspaced Short Palindromic Repeats Genotyping of Multidrug-Resistant Salmonella Heidelberg Strains Isolated From the Poultry Production Chain Across Brazil

Salmonella enterica subsp. enterica serovar Heidelberg has been associated with a broad host range, such as poultry, dairy calves, swine, wild birds, environment, and humans. The continuous evolution of S. Heidelberg raises a public health concern since there is a global dispersal of lineages harboring a wide resistome and virulome on a global scale. Here, we characterized the resistome, phylogenetic structure and clustered regularly interspaced short palindromic repeats (CRISPR) array composition of 81 S. Heidelberg strains isolated from broiler farms (n = 16), transport and lairage (n = 5), slaughterhouse (n = 22), and retail market (n = 38) of the poultry production chain in Brazil, between 2015 and 2016 using high-resolution approaches including whole-genome sequencing (WGS) and WGS-derived CRISPR genotyping. More than 91% of the S. Heidelberg strains were multidrug-resistant. The total antimicrobial resistance (AMR) gene abundances did not vary significantly across regions and sources suggesting the widespread distribution of antibiotic-resistant strains from farm to market. The highest AMR gene abundance was observed for fosA7, aac(6′)-Iaa, sul2, tet(A), gyrA, and parC for 100% of the isolates, followed by 88.8% for blaCMY–2. The β-lactam resistance was essentially driven by the presence of the plasmid-mediated AmpC (pAmpC) blaCMY–2 gene, given the isolates which did not carry this gene were susceptible to cefoxitin (FOX). Most S. Heidelberg strains were classified within international lineages, which were phylogenetically nested with Salmonella strains from European countries; while CRISPR genotyping analysis revealed that the spacer content was overall highly conserved, but distributed into 13 distinct groups. In summary, our findings underscore the potential role of S. Heidelberg as a key pathogen disseminated from farm to fork in Brazil and reinforce the importance of CRISPR-based genotyping for salmonellae. Hence, we emphasized the need for continuous mitigation programs to monitor the dissemination of this high-priority pathogen.

Salmonella enterica subsp. enterica serovar Heidelberg has been associated with a broad host range, such as poultry, dairy calves, swine, wild birds, environment, and humans. The continuous evolution of S. Heidelberg raises a public health concern since there is a global dispersal of lineages harboring a wide resistome and virulome on a global scale. Here, we characterized the resistome, phylogenetic structure and clustered regularly interspaced short palindromic repeats (CRISPR) array composition of 81 S. Heidelberg strains isolated from broiler farms (n = 16), transport and lairage (n = 5), slaughterhouse (n = 22), and retail market (n = 38) of the poultry production chain in Brazil, between 2015 and 2016 using high-resolution approaches including whole-genome sequencing (WGS) and WGS-derived CRISPR genotyping. More than 91% of the S. Heidelberg strains were multidrug-resistant. The total antimicrobial resistance (AMR) gene abundances did not vary significantly across regions and sources suggesting the widespread distribution of antibiotic-resistant strains from farm to market. The highest AMR gene abundance was observed for fosA7, aac(6 )-Iaa, sul2, tet(A), gyrA, and parC for 100% of the isolates, followed by 88.8% for bla CMY −2 . The β-lactam resistance was essentially driven by the presence of the plasmid-mediated AmpC (pAmpC) bla CMY −2 gene, given the isolates which did not carry this gene were susceptible to cefoxitin (FOX). Most S. Heidelberg strains were classified within international lineages, which were phylogenetically nested with Salmonella strains from European countries; while CRISPR genotyping analysis revealed that the spacer content INTRODUCTION Salmonella enterica subsp. enterica serovar Heidelberg is most often associated with eggs and poultry (Hennessy et al., 2004;Chittick et al., 2006;Foley et al., 2008;Folster et al., 2012;Jackson et al., 2013). However, this scenario has changed since Salmonella Heidelberg has been associated with other sources from a broad host range, such as dairy calves (Centers for Disease Control and Prevention [CDC], 2017), swine (Cabral et al., 2017), wild birds (Liakopoulos et al., 2016a), environmental sources (Antony et al., 2018), human-derived clinical specimens (Abdullah et al., 2021), and outbreaks (Antony et al., 2018), which denotes their importance as a high-priority pathogen.
One of the most important risk factors surrounding foodborne illness is the international food trade that has been circumstantially accompanied by Salmonellae dispersal beyond borders. In this context, there is a global dispersal of multidrugresistant lineages of serovar S. Heidelberg, reaching various countries in North America (Andrysiak et al., 2008;Centers for Disease Control and Prevention, 2014; Public Health Agency of Canada, 2014; Deblais et al., 2018;Cox et al., 2021), South America (Kipper et al., 2021), Europe (Liakopoulos et al., 2016b;Campos et al., 2018), and Asia (Wu et al., 2013). Therefore, the simultaneous increase and extended protraction of S. Heidelberg in many parts of the world have favored their genetic acquisition of virulence and antimicrobial resistance (AMR) genes through horizontal gene transfer (HGT), which has ultimately led to one of the most pressing global concerns.
Owing to their importance as a key poultry producer globally, Brazil quickly became the hotspot of S. Heidelberg and urgent actions were needed from the food safety authorities to mitigate this pathogen in order to reduce the economic losses in the poultry sector. In this context, most investigations to detect Salmonella enterica serovars in the poultry sector still had important methodological gaps, since the food industry focuses especially on Salmonella Typhimurium and Salmonella Enteritidis, which demonstrate the need for a combined approach between classical microbiology and high-resolution methods such as whole-genome sequencing (WGS) and clustered regularly interspaced short palindromic repeats (CRISPR) genotyping (Barrangou and Dudley, 2016;Thompson et al., 2018;Yousfi et al., 2020). Indeed, the use of these high-throughput sequencing analyses exemplifies a useful means, not only for identifying Salmonella serovars but also to trace back the origin of the contamination conferring a substantial aid in decision-making to the poultry sector. In this regard, we demonstrated the usefulness of WGS-based identification in our previous study for genotyping rare Salmonella enterica serovars isolated from food and related sources (Monte et al., 2021). This previous survey demonstrated that the CRISPR arrays were highly conserved, and this genomic inspection provides high-resolution genotyping of Salmonella serovars. Hence, we performed a genomic study by combining WGS and CRISPR genotyping to characterize S. Heidelberg isolates from different sources at broiler farms, slaughterhouses, transport, lairages, and retail markets in Brazil.

Salmonella Heidelberg Strains and Antimicrobial Susceptibility Testing
A total of 79 non-duplicate Salmonella enterica subsp. enterica serovar Heidelberg from our collection that included isolates obtained from broiler farms (n = 16), transport and lairage (n = 5), slaughterhouses (n = 22), and retail markets (n = 38) in Brazil between 2015 and 2016 were used in this study (refer Table 1). We also included two S. Heidelberg strains (SH159 and SSc139) from our previous work (Monte et al., 2019) for comparative purposes, totalizing eighty-one isolates. The Salmonella isolation was performed according to the International Organization for Standardization (Anonymous, 2007(Anonymous, , 2017. These isolates were serotyped on the basis of somatic O, phase 1, and phase 2 of H flagellar antigens by agglutination tests with antisera as specified in the Kauffmann-White-Le Minor scheme (Grimont and Weil, 2007;Guibourdenche et al., 2010).

Genomic Analysis
All S. Heidelberg isolates (n = 81) underwent DNA extraction performed by using a commercial kit (QiAmp tissue, Qiagen,

ST15
Germany) per manufacturer's guidelines. Genomic DNA of eighty-one Salmonella isolates was sequenced at a 300-bp pairedend-read using the Nextera XT library preparation kit at the MiSeq platform (Illumina, San Diego, CA, United States).
Resulted raw sequence reads underwent strict quality control by using default settings in CLC workbench 10.1.1 (Qiagen) as per Monte et al. (2019), while assemblies were annotated with PROKKA version 1.14-dev (Seemann, 2014). A core genome phylogeny was constructed with an alignment of the core genes determined by the software version 3.11.2; the BlastP threshold was set to 95% (Page et al., 2015). A pan-genome genes presenceabsence information from Roary was visualized with Phandango (Hadfield et al., 2018). The single nucleotide polymorphisms were extracted from the alignment using SNP-sites version 2.3.3 (Page et al., 2016). The phylogeny was reconstructed using RAxML version 8.2.12, using a General Time Reversible Model and Gamma distribution for rate heterogeneity (Stamatakis, 2014). The resulting phylogeny was tested against 1,000 bootstrap replications, as determined by implementing the majority rule, autoMR convergence criteria in the RAxML software (Pattengale et al., 2010). The phylogeny was visualized and annotated using iTol version 3 (Letunic and Bork, 2016).
Lastly, the assemblies were analyzed for acquired AMR genes and chromosomal point mutations using default settings of ResFinder 4.1 database available at the Center for Genome Epidemiology. 1 In addition, we used MLST 2.0 to detect multilocus sequence typing (MLST), and the PlasmidFinder software version 2.0.1 was run with database version 2018-11-20 (Carattoli et al., 2014). A minimum identity threshold of 95% was used as a filter for identification.

Clustered Regularly Interspaced Short Palindromic Repeats Genotyping and Phylogenetic Analysis
An automated high-throughput processing pipeline previously described by Nethery and Barrangou (2019) was used to identify the CRISPR loci within each strain. Using CRISPR Visualizer, we extracted and imported CRISPR loci into the web interface for visualization and alignment of all CRISPR spacer and repeat sequences. 2

Salmonella Heidelberg Strains Isolated From Brazil Are Genetically Related to South American, European, and Asian Isolates
While assessing the phylogenetics of S. Heidelberg strains (n = 81) sequenced in this study, we noticed that the core genome, calculated from WGS data, represented 74% of the pan-genome (4,139 out of 5,582 total genes). Little genomic variation was present among the core genome, as only 704 SNP sites were detected. Bootstrap values varied across the phylogeny, likely attributed to the small genomic variation among strains. S. Heidelberg strains did not cluster by year, source, or geographic location across the phylogeny suggesting the widespread distribution [regions (n = 6), sources (n = 4), years of isolation (n = 2)] and persistence of Salmonella strains in Brazil (Figure 1), which validate the previous surveys (Monte et al., 2019). Furthermore, 11 different clusters were identified as shown in Figure 1. Of these, 14 strains appear to be from independent lineages, given that they nested out of the main branches.  We also observed through SNP clustering (PDS000037185.127; n = 765 isolates; NCBI pathogen detection tool), cases of international clustering of S. Heidelberg from our collection (n = 77) with strains isolated from a variety of sources (food, human, and environment) from Brazil (n = 201), Chile (n = 9), the United Kingdom (n = 444), Germany (n = 1), the Netherlands (n = 2), South Korea (n = 1), and China (n = 1), which suggest a common ancestor origin (Supplementary Figure 1). Another two strains (SH265 and SH268) from this study nested in the SNP cluster (PDS000029160.10) with strains isolated from Brazil (n = 11), the United Kingdom (n = 4), and Chile (n = 1) (Supplementary Figure 2).

DISCUSSION
There has been a great interest in surveying the adaptation of Salmonella serovars to the poultry production chain because of their extensive persistence in the past, notably with S. Typhimurium and S. Enteritidis, which have caused significant economic losses to this sector. Furthermore, the prevalence of S. Heidelberg shown in this study is not the only issue, but the fact that highly drug-resistant and/or MDR isolates are being recovered in most steps of the poultry production chain, particularly in Brazil could be considered a public health threat, as there is a risk of it becoming globalized.
Based on AMR results, the β-lactam resistance was essentially driven by the presence of plasmid-mediated AmpC (pAmpC) bla CMY −2 gene, given the isolates which did not carry this gene were susceptible to FOX, while QRDR such as gyrA and parC genes drove quinolone resistance (Table 1). Indeed, the presence of strains harboring bla CMY −2 gene could have implications on a one health interface, since this plasmid is more likely to persist (Teunis et al., 2018). Besides that, all strains harbored chromosomal mutations in gyrA and parC genes promoting high-level resistance against quinolones that could have implications on human health as treatment options become limited. Disturbingly, this result corroborates the findings by van den Berg et al. (2019) that found 98.4% of the S. Heidelberg isolates resistant to fluoroquinolones. On the other hand, all S. Heidelberg strains from our collection, displayed susceptibility to azithromycin, which could be considered a promising agent against Salmonellae infections (Crump et al., 2015;Wen et al., 2017). Azithromycin has been used as an alternative treatment option for enteric fever even when the guidelines on susceptibility testing were not available. Like fluoroquinolones, azithromycin is an antimicrobial agent with efficient intracellular penetration (Crump et al., 2015;Wen et al., 2017).
The total AMR load also included encoding resistance genes for aminoglycoside [aac(3)-VIa, aph(3 )-Ic, aadA1, aadA2], chloramphenicol (cmlA1), macrolides [Inu(G)], trimethoprim (dfrA12), and ammonium quaternary compounds (qacL) (Figure 1). Interestingly, the unique strain (SSc139) that carried the qacL gene was isolated from the retail market, which makes large use of such compounds to disinfect surfaces. Moreover, this strain nested in the same cluster with five strains isolated from broiler farms and slaughterhouse sources, in different states such as São Paulo, Santa Catarina, and Parana, denoting the successful establishment of this lineage (ST15) in the Southern and South regions of Brazil. In this regard, we visualized the co-occurrence of AMR genes in all S. Heidelberg strains.
The presence of AMR genes and transmissible plasmids demonstrated little variation across the strains. The broad distribution and abundance of S. Heidelberg in broiler farms, slaughterhouses, transport, lairages, and retail markets suggest the high adaptability of this serovar in the poultry production chain in Brazil. Similarly, a study by Edirmanasinghe et al. (2017) examining FOX-resistant S. Heidelberg strains isolated FIGURE 2 | Reconstructed phylogeny based on the core genome, distribution of spacers composition, clustered regularly interspaced short palindromic repeats (CRISPR) loci, and CRISPR profiles among Salmonella Heidelberg strains. The location of isolation of each strain is labeled on its respective branch. Color strips depict the year (Y) and source (S) of isolation, respectively. Brazilian States: PR, Paraná; SC, Santa Catarina; SP, São Paulo; MG, Minas Gerais; DF, Distrito Federal; MS, Mato Grosso do Sul.
FIGURE 3 | Salmonella enterica subsp. enterica serovar Heidelberg type I-E CRISPR locus architecture. This system contains two distinct CRISPR arrays-one associated with the cas genes and one disparate locus upstream (5 ) from the cas genes. from human, abattoir poultry, and retail poultry sources in Canada revealed a potential common source among strains, which suggest the simultaneous dispersal of S. Heidelberg strains carrying CMY-2 gene in several sources and different geographical locations. In convergence with our results, another study revealed a high occurrence of S. Heidelberg in imported poultry meat in the Netherlands containing bla CMY −2 gene (van den Berg et al., 2019). Although we could not find colistinresistant strains, it is worthwhile to note the presence of the IncX4 plasmid since it is a promiscuous plasmid with a high capacity of self-transmissibility that is commonly associated with the presence of plasmid-mediated colistin-resistance (mcr-1) in Brazil (Moreno et al., 2019).
Consistent with observations obtained in our study, recent surveys strongly support that S. Heidelberg may originate from a common ancestor, which circulates and persists in the Brazilian poultry production chain since at least 2004 (Kipper et al., 2021). More importantly, these strains are also nested with strains isolated from several sources and countries around the world supporting the hypotheses of intercontinental spread, which demonstrate that probably the common ancestor underwent diversification through genetic changes over time (Supplementary Material).
Protracted dissemination of S. Heidelberg via poultry might be a risk for a globalized food trade era. The introduction and clonal expansion of Salmonella strains across borders remain challenging due to the difficulties of identifying the origins of contamination. In light of this, the continued need for combined approaches between classical microbiology and high-resolution methods such as WGS and CRISPR genotyping truly illustrate to us what is hidden in plain sight.
For the purpose of discussion, other studies provide compelling validation data to support the usefulness of highresolution methods for genotyping rare Salmonella enterica serovars (Monte et al., 2021) and/or to resolve S. Heidelberg isolates involved in foodborne outbreaks (Vincent et al., 2018). Of the latter, while assessing the CRISPR array of 145 S. Heidelberg isolates, Vincent et al. (2018) found 15 different CRISPR profiles endorsing our results.
This study further illustrates the potential of CRISPR for the tracking of variable genotypes in diverse Salmonella strains, as previously determined (DiMarzio et al., 2013;Shariat et al., 2013aShariat et al., ,b, 2015Monte et al., 2021), with noteworthy methodological convenience. Indeed, CRISPR-based analyses have proven relevant for subtyping of Salmonella enterica serovars Typhimurium and Heidelberg strains involved in outbreaks (Shariat et al., 2013b) and occasionally associated with antibiotic resistance (DiMarzio et al., 2013). Our findings underscore the potential role of S. Heidelberg as a key pathogen in the poultry production chain, particularly in Brazil.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.