Prevalence of virulence- and antibiotic resistance-associated genotypes and phenotypes in Staphylococcus aureus strains from the food sector compared to clinical and cow mastitis isolates
Andrea Jurado
Lucía Fernández
Ana Rodríguez
Pilar García- 1Department of Technology and Biotechnology, Instituto de Productos Lácteos de Asturias (IPLA-CSIC), Asturias, Spain
- 2DairySafe Group. Instituto de Investigación Sanitaria del Principado de Asturias (ISPA), Oviedo, Spain
Background: Infections by the pathogen Staphylococcus aureus currently represent one of the most serious threats to human health worldwide, especially due to the production of enterotoxins and the ability to form biofilms. These structures and the acquisition of antibiotic resistance limit the action of antibiotics and disinfectants used to combat this microorganism in the industry and the clinic.
Methods: This work reports a comparative phenotypic and genotypic study of 18 S. aureus strains from different origins: clinical samples, milk from mastitic cows and food industry surfaces, most of which were isolated in Northern Spain.
Results: Genetically, the strains were very diverse but, in most cases, a closer proximity was observed for those from the same source. Notably, the average number of virulence genes was not significantly different in strains from the food sector. Of the 18 strains, 10 coded for at least one enterotoxin, and four of them carried 6 or 7 enterotoxin genes. The latter were all veterinary or clinical isolates. Most strains carried prophages, plasmids and/or pathogenicity islands. Regarding antibiotic resistance, although phenotypically all strains showed resistance to at least one antibiotic, resistance genes were only identified in 44.5% of strains, being mastitis isolates those with the lowest prevalence. Virulence-related phenotypic properties such as haemolytic activity, staphyloxanthin production, biofilm-forming capacity and spreading ability were widely distributed amongst the isolates.
Conclusions: Our results indicate that production of virulence factors, antibiotic resistance and biofilm formation can be found in S. aureus isolates from diverse environments, including the food industry, although some of these traits are more prevalent in strains isolated from infections in cows or humans. This emphasizes on the importance of monitoring the spread of these determinants not only in samples from the clinical environment, but also along the food chain, a strategy that falls under the prism of a one-health approach.
1 Introduction
Foodborne diseases remain one of the most serious health problems worldwide. According to the latest available EFSA report, 4,005 foodborne outbreaks (FBOs) and 32,543 human cases, including 2,495 hospitalizations and 31 deaths, were reported by the EU Member States in 2021 (EFSA and ECDC (European Food Safety Authority and European Centre for Disease Prevention and Control), 2022). These figures represent a notable increase compared to 2020 (between 29-62%). The number of FBOs involving bacterial toxins has also increased by 152 cases in 2021 (EFSA and ECDC (European Food Safety Authority and European Centre for Disease Prevention and Control), 2022). Amongst toxin-derived outbreaks, those produced by coagulase-positive staphylococci, such as Staphylococcus aureus, led to the highest hospitalization rates, and are amongst the most common foodborne diseases in many parts of the world (Hennekinne et al., 2012). This microorganism commonly inhabits the bodily surfaces of humans and other animals, most frequently without causing disease. Indeed, about 70-90% of the population is considered to be, at least transiently, asymptomatic carriers of S. aureus (Hasse-Cieślińska, 2007). With this in mind, it is not surprising that the main routes of food contamination are workers, during food processing, and farm animals (Argudín et al., 2010). Additionally, dairy animals, such as cows, goats or sheep, can also suffer from staphylococcal infections, especially mastitis, which causes large economic losses in the livestock sector (Ajose et al., 2022). Some of the foods that pose a greater risk of contamination with this bacterium include meat and meat products, poultry and egg products, milk and dairy products, bakery products (cream-filled pastries and cakes) and ready-to-eat foods (Kadariya et al., 2014). Contamination with S. aureus frequently takes place during food processing. As a result, the isolation of this microorganism from food contact surfaces in the dairy, meat and seafood industries is highly documented (Gutiérrez et al., 2012), and is considered an indicator of inadequate cooking and/or food processing (Castro et al., 2018).
Apart from food poisoning, S. aureus can also cause a wide range of infections in humans, especially when their immune system is compromised, being currently considered one of the most notorious nosocomial pathogens (Kadariya et al., 2014; WHO (World Health Organization), 2021). This ability to cause disease is due to the production of diverse toxins (e.g. enterotoxins, leukotoxins, hemolysins, exfoliative toxin, etc.) and the possession of immune evasion mechanisms (e.g. antiphagocytic capsule, sequestration of host antibodies, antigen masking by protein A, intracellular survival, blocking of leukocyte chemotaxis, etc.) (Hennekinne et al., 2012; Malak et al., 2020).
S. aureus enterotoxins are the main cause of food poisoning, as they are able to persist in food even after denaturing treatments that kill bacteria, and are resistant to the action of human enzymes such as pepsin and trypsin (Fisher et al., 2018). As superantigens, their effect leads to the activation of a high number of T cells, causing cytokine release and, sometimes, systemic shock (Fisher et al., 2018). The currently described enterotoxins can be classified into staphylococcal enterotoxins with emetic activity (SEA-SEE, SEG-SEI, SER-SET) and staphylococcal-like enterotoxins, which do not have a demonstrated emetic activity in a primate model or have not been tested yet (SElL, SElQ, SElJ, SElK, SElM to SElP, SElU, SElU2 and SElV). They can also be divided into classical (SEA to SEE) and new (SEG to SElU2) enterotoxins, the latter being encoded by the enterotoxin gene cluster (Argudín et al., 2010; Johler et al., 2015).
Genes encoding virulence factors, including the aforementioned enterotoxins, and antibiotic resistance are often located in mobile genetic elements, such as plasmids or bacteriophages, increasing their distribution potential (Ayliffe et al., 1996; Hennekinne et al., 2012). These phages capable of exchanging genes between bacterial cells belong mostly to the so-called temperate phages, and have as a main feature the ability of integrating into the bacterial chromosome as prophages. Based on their role in horizontal gene transfer, temperate phages pose a serious problem regarding the treatment of these strains; hence the importance of their detection and characterization (Goerke et al., 2009). Antibiotic resistance development can also be associated with the acquisition of chromosomal cassettes, transposons located in conjugative plasmids or chromosomal mutations. Worryingly, some S. aureus strains have evolved resistance to many of the currently available antibiotics, including methicillin and even vancomycin (McGuinness et al., 2017; Lee et al., 2018). As a result, the WHO has categorized this microbe as a high priority pathogen (WHO (World Health Organization), 2021). Indeed, healthcare-acquired methicillin-resistant S. aureus (HA-MRSA) strains constitute a serious problem in the hospital environment, causing serious diseases, such as bacteremia, ventilator-associated pneumonia and surgical wound infections, that often result in therapeutic failure (Kanafani and Fowler, 2006; Pofahl et al., 2011; Hurley, 2018). In addition, the spread of these MRSA strains has extended the problem to the out-of-hospital community (CA-MRSA) (de Kraker et al., 2011). This bacterium is also well known for its ability to persist in hostile environments, including body tissues, by forming biofilms, in which bacterial cells are surrounded by an extracellular matrix consisting of polysaccharides, proteins and/or DNA (Flemming and Wingender, 2010). This complex structure confers biofilms with high resistance to antibiotics, as well as to the cleaning and disinfection systems commonly used in the food industry (Bridier et al., 2015; Olsen, 2015). Staphylococcal biofilms have been associated with endovascular, bone and joint infections, as they can grow and persist on catheters and prostheses. Moreover, biofilms are one of the best ways for this bacterium to survive on food industry surfaces and, ultimately, contaminate foods (Castro et al., 2018).
The determination of genetic markers associated with strains from specific environments or geographical areas can be used to study the dispersion path of a strain or its mobile elements and identify its source. This would allow the implementation of more refined measures to prevent its expansion. Also, the levels of certain virulence factors can help establish the degree of severity of the pathology, as in the case of leukocidin LukM/F, which is associated with zoonotic infections and can be detected in the milk of mastitic cows (Vrieling et al., 2016).
In the present work, we carried out an in depth analysis of several strains from Northern Spain from different sources: the food (meat and dairy) industry and mastitic cows, many of which had been isolated as part of the same study. Additionally, these strains were compared to well-characterized reference strains of clinical origin. Regarding the food and mastitis strains, our most significant finding was that despite having the same geographic origin and, in some cases, having been isolated within a short time period, they are actually quite genetically diverse.
2 Materials and methods
2.1 Bacterial strains, bacteriophages and culture conditions
Eighteen different S. aureus strains from different origins were used in this study (Table 1). All the bacteria were routinely cultured in tryptic soy broth (TSB; Scharlau, Barcelona, Spain) at 37°C with shaking or on TSB plates containing 2% (wt/vol) bacteriological agar (TSA).
Table 1 Typing and origin of the different S. aureus strains used in this study.
2.2 Genome sequencing, assembly and annotation
The sequences of strains SH1000, MW2, 132, RN4220, Newman, JE2 and V329 were available in public databases with accession numbers JANFOB010000001, NC_003923, NZ_ACOT01000046, NZ_WWFP01000001, NC_009641, NZ_CP020619 and JAGTJH000000000, respectively.
The following S. aureus strains were grown on TSA plates at 37 °C: Sa9, Sa5, Sa7, IPLA19, IPLA1, IPLA3, IPLA5, IPLA11, IPLA13, IPLA15, IPLA16. A single colony from each strain was then streaked out on a fresh plate and incubated overnight. All cells grown on this plate were harvested and resuspended in 1 ml PBS to obtain approximately 5 x 109 cells. To estimate the number of cells, dilutions of the resulting suspension were made to obtain an OD600 of 0.1, which corresponds to 1 x 107 cells. The final number was calculated by correcting by the dilution factor. Cells were pelleted and washed with PBS and, after centrifuging again at maximum speed for 3 minutes, resuspended in 0.5 ml of 1 x DNA/RNA Shield buffer (Zymo Research, Irvine, USA) until gDNA isolation. Samples were treated with 0.2 mg/ml lysostaphin and 0.1 mg/ml RNase A (ITW Reagents, Spain) for 25 min at 37 °C. These samples were further incubated for 5 min at 65 °C with 0.1 mg/ml proteinase K (VWR Chemicals, Ohio, USA) and 0.5% v/v SDS (Sigma-Aldrich, Missouri, USA). Purification of gDNA was carried out using solid-phase reversible immobilization (SPRI) beads (Beckman Coulter, Brea, USA), and the genomic libraries were prepared according to the manufacturer’s protocol using Nextera XT Library Prep Kit (Illumina, San Diego, USA). The resulting libraries were then sequenced on an Illumina HiSeq platform using a 250 bp paired end protocol. Additionally, strains Sa9, IPLA15 and IPLA16 were sequenced with an Oxford Nanopore sequencing platform to obtain long-read sequencing data. Genome sequencing was provided by MicrobesNG (http://www.microbesng.com). Preparation of the DNA libraries was carried out with Oxford Nanopore SQK-LSK109 kit with Native Barcoding EXP-NBD104/114 (ONT, United Kingdom) using 400-500 ng of high-molecular-weight DNA. Next, the barcoded samples were pooled together into a single sequencing library before being loaded in a FLO-MIN106 (R.9.4.1) or FLO-MIN111 (R10.3) flow cell in a GridION (ONT, United Kingdom).
Unless otherwise noted, genome analysis was carried out using default parameters for all software. Quality of the reads was checked with FASTQC v. 0.11.3 (Andrews, 2010) and trimming was performed with Trimmomatic v. 0.39 by using a sliding window quality cutoff of Q15 (Bolger et al., 2014). De novo assembly of the genomes was achieved using SPAdes v. 3.14.1 (Bankevich et al., 2012). In the case of strains Sa9, IPLA15 and IPLA16, hybrid assemblies were generated using Unicycler version 0.4.0 (Wick et al., 2017). Quality of the assemblies was assessed with QUAST v. 5.0.2 (Gurevich et al., 2013). Sequences shorter than 500 bp were removed prior to genome annotation with the NCBI Prokaryotic Genome Annotation Pipeline v. 5.1 (Tatusova et al., 2016) for strains Sa9, IPLA15 and IPLA16 and Prokka (Seemann, 2014) for the rest.
2.3 Analysis of mobile genetic elements
The annotated genomes of all the strains were screened for the presence of genomic islands with the tool IslandViewer 4 (IslandViewer 4 (https://www.pathogenomics.sfu.ca/islandviewer) (Bertelli et al., 2017). Identification and annotation of prophage sequences and plasmids was performed by using PHASTER (PHASTER, http://phaster.ca/) (Arndt et al., 2016) and PlasmidFinder 2.1 (https://cge.cbs.dtu.dk/services/PlasmidFinder-2.0/) (Carattoli et al., 2014), respectively. For the analysis of the SCCmec elements, the tool SCCmecFinder1.2 from the same webserver was used with the following parameters: 90% minimum identity and 60% minimum coverage (CGE Server (https://cge.dtu.dk/services/SCCmecFinder) (International Working Group on the Classification of Staphylococcal Cassette Chromosome Elements (IWG-SCC), 2009). This latter analysis was only carried out in mecA positive strains.
2.4 Phylogenetic analysis
A single nucleotide polymorphism (SNP) tree was constructed using CSI Phylogeny 1.4 (CGE Server (https://cge.food.dtu.dk/services/CSIPhylogeny/) (Kaas et al., 2014) and the assembled FASTA files of all S. aureus strains sequenced in this study together with the following reference strains from GenBank: S. aureus SH1000 (JANFOB010000001.1), MW2 (NC_003923.1), JE2 (NZ_CP020619.1), 132 (NZ_ACOT01000046.1), RN4220 (NZ_WWFP01000001.1), Newman (NC_009641.1) and V329 (JAGTJH000000000). This analysis was carried out using default settings and excluding heterozygous SNPs. The following criteria for high-quality SNP calling and filtering were chosen: (i) a minimum depth of 10× at SNP positions; (ii) a minimum relative depth of 10% at SNP positions; (iii) a minimum distance of 10 bp between SNPs; (iv) a minimum SNP quality of 30; (v) a minimum read mapping quality of 25; and (vi) a minimum Z score of 1.96.
2.5 Typing and identification of virulence and antimicrobial resistance genes
The presence of antimicrobial resistance and virulence genes was evaluated by using two resources, namely, ResFinder4.1 (https://cge.cbs.dtu.dk/services/ResFinder/) (Bortolaia et al., 2020) and VFDB (VFDB: Virulence Factor Database (http://www.mgc.ac.cn/cgi-bin/VFs/v5/main.cgi) (Liu et al., 2019), respectively, using the default parameters. This last tool was also used to carry out a comparative analysis of our strain collection with S. aureus strains from different origins and geographical locations (Table S3). The presence of phenol-soluble-modulin genes (PSMs-α, PSM-β, hld) and mutations in selected genes of interest was analyzed with BLAST https://blast.ncbi.nlm.nih.gov/Blast.cgi) (Camacho et al., 2009). To determine the sequence type (ST) of the isolated strains, in silico MLST analysis was performed using the MLST 2.0 database (MLST 2.0 (https://cge.cbs.dtu.dk/services/MLST/) (Larsen et al., 2012). The S. aureus MLST scheme consisted of 7 housekeeping genes: arcC (carbamate kinase), aroE (shikimate dehydrogenase), glpF (glycerol kinase), gmk (guanylate kinase), pta (phosphate acetyltransferase), tpi (triosephosphate isomerase), and yqi (acetyl coenzyme A acetyltransferase). The parameter 5X for minimum depth for an allele was selected and the assembled genomes were used as input.
Spa types were predicted using spaTyper v1.0 webserver from the Center of Genomic Epidemiology (https://cge.cbs.dtu.dk/services/spatyper) (Bartels et al., 2014). In this typing technique, 21 to 27 polymorphic VNTR in the 3-coding region of staphylococcal protein A (spa) were compared in order to assign a unique repeat code corresponding to its spa type.
2.6 Biofilm formation by different bacterial strains
Biofilm assays were carried out as described previously, with some modifications (Herrera et al., 2007). Briefly, overnight cultures of different S. aureus strains were diluted 100-fold in fresh TSBg (TSB supplemented with 0.25% w/v D-(+)-glucose) and 200 μl-aliquots were poured into each well of a 96-well plate (Thermo Scientific, NUNC, Madrid, Spain). Biofilms were grown for 24 h at 37 °C and, subsequently stained with 0.1% crystal violet. To do that, the planktonic phase was removed and the adhered cells were washed twice with PBS (137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4 and 2 mM KH2PO4; pH 7.4). Afterwards, 200 μl of crystal violet were added and, following incubation for 15 minutes at room temperature, excess of the dye was removed by washing the biofilms with water. Finally, the crystal violet retained in the biofilm was solubilized with 33% (v/v) acetic acid, and the total biomass was quantified by measuring absorbance at 595 nm (A595) with a Bio-Rad Benchmark plus microplate spectrophotometer (Bio-Rad Laboratories, Hercules, CA, USA).
2.7 Hemolysin, staphyloxanthin and protease production
To determine staphyloxanthin production, 5 μl from an overnight culture of each strain were placed on a TSA plate, which was then incubated for 20-24 hours at 37 °C. The next day, staphyloxanthin production was determined by visual inspection as the presence of an orange-yellow color. On the other hand, for the characterization of hemolysin production, 5 μl from an overnight culture of each strain were laid on a blood agar plate and after 24 hours of incubation at 37 °C, the action of alpha hemolysin was detected by visual analysis as the presence or absence of a halo surrounding the colonies. The plates were observed again after another 24 hours of incubation at 4 °C in order to observe the effect of the beta hemolysin (Smyth et al., 1975). The detection of protease activity was carried out on plates containing nutrient broth (Becton, Dickinson and Company, Sparks, USA), supplemented with 1.5% agar and 1.5% skim milk (OXOID, Hampshire, England). These plates were inoculated with 5 μl from S. aureus overnight cultures and subsequently incubated for 24 h at 37 °C. After incubation, the proteolytic activity was considered positive when a clear zone was observed around the bacterial colonies.
2.8 Colony spreading ability of different strains
Fresh culture plates were prepared with autoclaved TSB medium containing 0.25% agar. Once dried, a volume of 2 μl from an overnight culture of each strain was inoculated on the center of each plate and incubated overnight at 37°C. After incubation, the diameter of each bacterial colony was measured.
2.9 Antimicrobial susceptibility test
The antimicrobial susceptibility profile of the S. aureus isolates was determined by the Kirby–Bauer disk diffusion method (Biemer, 1973). Briefly, overnight cultures of the different strains were diluted 10-fold in PBS and the suspension was streaked on 20-ml TSA plates and allowed to dry. Antibiotic disks were then placed on each plate with a considerable separation between them and incubated at 37°C for 24 h. After the incubation time, the zones of inhibition were measured with a ruler and interpreted in order to compare the antimicrobial susceptibility of the strains. The following antimicrobial agents were used: ampicillin (10 µg), methicillin (10 µg), tetracycline (30 µg), chloramphenicol (30 µg), gentamicin (10 µg), tobramycin (10 µg), erythromycin (15 µg), ciprofloxacin (5 µg), vancomicin (5 µg), streptomycin (10 µg). Antibiotic disks were purchased from Oxoid (Wade Road, Basingstoke, Hants, UK).
2.10 Statistical analysis
All experiments were performed in triplicate. Different groups of samples were compared by using one-way analysis of variance (ANOVA) followed by a post hoc Tukey HDS test. P-values <0.05 were considered statistically significant.
2.11 Genome accession numbers
The genome assemblies of 11 S. aureus strains were deposited in the DDBJ/ENA/GenBank database under the following accession numbers: CP134618-CP134619 (IPLA15), CP134617 (IPLA16), CP134620 (Sa9), JAVRXU000000000 (IPLA1), JAVRXT000000000 (IPLA3), JAVRXS000000000 (IPLA5), JAVRXR000000000 (IPLA11), JAVRXQ000000000 (IPLA13), JAVRXV000000000 (IPLA19), JAVRXX000000000 (Sa5) and JAVRXW000000000 (Sa7).
3 Results
3.1 Genomic analysis of S. aureus strains
The aim of this study was the genomic and phenotypic comparison of 18 S. aureus strains from different origins, including 5 from cows with mastitis (V329, Sa9, Sa5, Sa7, IPLA19), 4 from the meat industry (IPLA11, IPLA13, IPLA15, IPLA16), 3 from the dairy industry (IPLA1, IPLA3, IPLA5), 4 from the clinic (JE2, MW2, 132, Newman) and 2 laboratory strains (SH1000 and RN4220) derived from the clinical strain S. aureus NCTC 8325 (Table 1). In order to do that, we sequenced the genomes of 11 strains, namely Sa9, Sa5, Sa7, IPLA19, IPLA1, IPLA3, IPLA5, IPLA11, IPLA13, IPLA15, IPLA16, whose genomes were not available in public databases. Supplementary Table 1 summarizes the results regarding genome assembly and annotation. The genome size ranged from 2,708,529 to 2,805,322 bp, and the G+C content was very similar between the strains, with 32.70% and 32.87% being the minimum and the highest values, respectively.
Once all the genome sequences were available, we carried out a comparative genomic analysis of the 18 strains. All of them showed high diversity regarding their multilocus sequence typing (MLST) profile (Table 1). A similar trend was observed for the spa types, which were also very diverse. Only one type (t008) was found in several clinical strains (JE2, Newman, 132) (Table 1). Regarding capsule type, 39% (7/18) of the strains had capsule type 8, while 61% (11/18) were type 5 (Table 1). It is noteworthy that all the strains from the meat industry were type 8, whereas those from the dairy industry were type 5. Amongst the strains of clinical origin, most were type 5 except for MW2, which had capsule type 8, and those from mastitic cows presented in similar numbers capsule type 5 (V329, Sa5, Sa7) and 8 (Sa9, IPLA19). The most prevalent agr type was type I (50% of the strains), followed by type II (39%) and type III (11%) (Table 1). Again, MW2 (agr type III) was the only clinical strain that did not have agr type I. Also, strains from the dairy industry all were type I, whereas most strains from mastitis and the meat industry were type II except for Sa5 (type I) and IPLA15 (type III). The phylogenetic tree obtained based on SNP analysis revealed close proximity between all the clinical strains, including the laboratory strains, with the exception of MW2, which clustered with the meat industry isolates (Figure 1). Some mastitis isolates (Sa5, Sa9 and IPLA19) grouped fairly close to those from the dairy industry (IPLA5, IPLA1, IPLA3) and in close proximity to the clinical isolates. In contrast, the other two strains from milk of cows with mastitis (Sa7 and V329) formed a separate cluster.
Figure 1 Phylogenetic analysis of S. aureus isolates. The genome-wide SNP-based maximum likelihood phylogenetic tree of S. aureus isolates obtained from different origins (Clinical sector: Green, Food industry: Purple (meat industry), Blue (Dairy Industry), mastitis: Red). The genome of strain NCTC 8325 was used as the reference to map and screen the SNPs. The phylogenetic tree was constructed using the CSI Phylogeny v1.4 and visualized using FigTree.
3.2 Identification of virulence and antibiotic resistance genes
The presence of different types of virulence factors in 17 of the strains was examined and the results are summarized in Supplementary Table 2. Strain SH1000 was not included in this analysis as it is practically identical to RN4220.The virulence factors with the highest prevalence in all isolates regardless of origin were toxin-coding genes (19.11%), followed by adhesion-related genes (15.35%), while those involved in evasion showed the lowest prevalence (3.82%) (Figure 2A). In terms of origin, strains of clinical origin exhibited the highest average prevalence of virulence factors of all types (adherence: 16%, enzymes: 14%, evasion factors: 4.2%, secretion: 11.6%, toxins: 20.4%) (Figure 2A). Mastitis isolates possessed a higher mean number of genes involved in secretion (10.60%) and toxin production (20.2%) than food isolates, but a lower prevalence of those related to adherence (14.60%), enzymes (12.2%) and evasion (3.4%) (Figure 2A). Nonetheless, it must be noted that none of these differences were statistically significant, indicating that food industry strains, exhibited an average number of virulence determinants similar to that of strains from infected cows or humans.
Figure 2 Distribution and pairwise associations of virulence genes in S. aureus isolates. (A) Distribution of genes encoding different types of virulence factors in S. aureus isolated from the clinical sector (white), food industry (black), and mastitis (grey). Data represents the average and standard deviation. (B) Pairwise associations of selected virulence genes, computed using phi coefficients (cut-off 0.7), with a colour gradient representing the type of association. The intensity of the colour and the size of the circle indicate the strength of the association.
Multiple virulence genes were present in all strains, including adhesion factors (atl, ebp, eap, efb, fnbA, icaA, icaC, icaR and spa), enzymes (sspB, sspC, sspA, hysA, geh, lip, coa, nuc and aur), evasion factors (adsA, sbi, and isd family genes), and toxins (hly/hla, hlgA, hlgB, and hlgC). However, other virulence genes were only observed in a few isolates. These included the enzyme-encoding gene sak (JE2, MW2, 132 and Newman from clinical origin; Sa7 from mastitis and IPLA13 from the meat industry) and the evasion factor gene chp (JE2 and Newman from clinical origin and IPLA11 from the meat industry). Strains of meat origin showed the lowest possession of type VII secretion genes (esxB, esxC, esxD, esaA and esaE), which were only present in strain IPLA15. On the other hand, the exfoliative toxin genes were not detected in any strain. It is worth highlighting that the Panton-Valentine leukocidin genes lukF0PV and lukS0PV were only present in two strains of clinical origin (JE2 and MW2), and the gene coding for the toxic shock syndrome toxin (tsst) was only present in the isolate strain Sa9, which together with IPLA19 are the only strains carrying the lukM gene. Another interesting finding was that only clinical strains were positive for phenol soluble modulins type α (PSM-α), whereas PSM-β modulins type I and II were present in all strains in the collection. Regarding enterotoxin genes, 55.5% (10/18) of the strains carried at least one gene, with sea (27.7%), sec (27.7%) and selk (27.7%) being the most prevalent (Table 2). Four strains possessed 6 or 7 enterotoxin genes, of which three were isolated from mastitic cows and one was a clinical strain. By contrast, seven strains did not carry any enterotoxin genes, four of which were isolated from food industry samples. Some enterotoxin genes (seb, sed, see, sei, sej, selp and selr) were absent from all strains, while enterotoxin-like genes were more prevalent in mastitis isolates (Table 2). With the aim of determining if these trends could be extrapolated to other studies, we examined the presence of different enterotoxins in the genomes of 75 strains from different origins (Supplementary Table 3). Just like in our study, no strain carried enterotoxins sed, see, sej, selp or selr, and none of the mastitis isolates produced sea or seb. However, in this case, enterotoxin-like genes were not more prevalent in mastitis isolates and were also very frequent in clinical and even in food strains. Another significant difference is that the strains with the highest number of enterotoxin genes were actually food isolates.
Table 2 Presence of enterotoxin genes in S. aureus isolates.
Analysis of the correlation between the different types of virulence genes in the eighteen strains showed that the adhesion gene fnbB exhibited a negative association with genes coding for serine proteases (splA, splB, splC) and type VII secretion systems, except essA. In turn, genes coding for serine proteases showed a strong positive association among themselves and with the leukocidin lukD gene. Most type VII secretion system genes (esaD, esaE, esxB, esxC, esxD) showed a pattern of positive association with each other. In contrast, essA exhibited a negative association with the other secretion genes and a positive correlation with the presence of toxin hld. Regarding enterotoxin genes, many gene pairs (i.e selm-seln, yent1-yent2, selo-selu, seg-selo, selu-set20) and groups of three (i.e set3-set4-set5, set22-set23-set25, set30-set35-set36) showed a strong positive association between them and a similar pattern of association with genes of the other categories. Regarding leukocidin genes, lukF0PV and lukS0PV showed a strong positive association with each other and with the selq, set21 toxins genes, whereas lukM is positively associated with the toxin genes selU and set20 (Figure 2B).
Regarding antibiotic resistance determinants, we found that 47% (8/17) of all strains presented genes that might lead to resistance to at least one class of antibiotics. However, only two of them (132 and IPLA3) showed a multiresistance profile (resistance to three or more classes of antibiotics). In addition, 53% (9/17) did not carry any detectable antibiotic resistance genes in their genome, being the strains from mastitis the ones that showed the lowest prevalence of such genes (Table 3). Genes related to beta-lactam resistance were a common denominator in all strains with resistance genes, being in five strains (IPLA15, IPLA11, IPLA5, IPLA3, Sa7) due to the blaZ gene and in three (132, JE2, MW2) due to mecA, genes that cause decreased susceptibility to penicillin and methicillin, respectively. Strain IPLA3 possessed both genes. Presence of the staphylococcal cassette chromosome was observed in all mecA-positive strains. On the other hand, potential resistance to quinolones was also identified in two strains of clinical origin (JE2 and 132), which had mutations in the genes gyrA (p.S84L) and grlA (p.S80Y), to aminoglycosides in three strains from different origins (132, Sa7, IPLA3) mediated by the genes ant(4’)-Ia, str, aph(3’)-Ia, and finally to tetracycline and streptomycin only in IPLA3 strain due to the presence of the genes tetM, tetK and str, respectively. It is noteworthy that the bleO gene was found in strain 132, which confers resistance to bleomycin, an antibiotic that causes DNA strand breaks (Table 3).
Table 3 Antimicrobial resistance genes of S. aureus isolates.
3.3 High abundance of pathogenicity islands and other mobile elements in S. aureus strains
The prevalence of intact phages among the different groups of strains according to their origin was very homogeneous. Only 28% (5/18) of the strains did not carry any intact prophages in their genomes, being some strains of clinical origin the ones that showed the highest number (Newman: 4, 132: 3). The % GC of the intact prophages detected ranged between 31.99% and 35.09% (Table 4). All strains had at least one pathogenicity island in their genomes, the number being variable among the different strains. Of note, the dairy industry strain IPLA3 possessed the highest number (13), followed by the clinical strain MW2, which had 8 (Table 4). The presence of plasmids was also detected in several strains, although 7 did not carry any predicted replicon according to PlasmidFinder (Table 4). Some plasmids could be identified by homology with the available databases, and their sizes ranged from 4 to 36 kb. It is worth noting that some of the resistance determinants described above were located in the identified plasmids.
Table 4 Genomic analysis of mobile genetic elements.
3.4 Phenotypic analysis of virulence properties
In addition to genomic characterization, we also studied several phenotypes of interest in the different isolates, namely biofilm formation, hemolysin activity, staphyloxanthin production, colony spreading, proteolytic activity and antibiotic resistance. Regarding biofilm formation, all strains displayed some degree of adherence to the surface but only strains V329, SH1000 and 132 were strong biofilm formers (Figure 3A). These three strains were a mastitis isolate, a laboratory strain of clinical origin, and a clinical isolate, respectively. Two of them presented exclusive genes like bap in V329 and sasG in Newman. They also lacked the icaD gene, belonging to the staphylococcal icaADBC operon involved in the production of PIA/PNAG (polysaccharide intercellular adhesin). In the case of strain 132 the quorum sensing gene agrC was missing. Also noteworthy was the lack of icaB in the IPLA11 strain. Additionally, the genes clfA, clfB, related to the clumping mechanism, were found in some strains (i.e SH1000, JE2, MW2, RN4220, Newman, Sa9, IPLA5) (Supplementary Table 2).
Figure 3 Biofilm and staphyloxantin production of S. aureus isolates. (A) Biofilm production of 17 S. aureus strains from different origin (Clinical sector, Food industry, Mastitis) by conventional crystal violet staining technique. Bars represent mean ± standard deviation of three biological replicates of each strain. Bars within each strain having distinct lower case letter indicate that the biofilm production is statistically different (p < 0.05). Pink colour represents the most biofilm-producing strains. (B) Staphyloxantin production of 18 S. aureus strains from different origin.
A greater diversity was found in terms of staphyloxanthin production. This carotenoid pigment has antioxidant properties and helps the microbe evade the action of reactive oxygen species produced by the immune system (Clauditz et al., 2006). Two strains (V329, RN4220) did not produce any pigment, one of them (Sa9) only synthesized an intermediary pigment and the rest produced varying levels of staphyloxanthin. The absence of production in the above mentioned strains could be explained because of the inactive SigB caused by a deletion in the rsbU gene in RN4220 and a mutation in crtN in V329, a gene that encodes a dehydrosqualene desaturase which dehydrogenates dehydrosqualene (the first compound in the staphyloxanthin biosynthesis pathway) (Pelz et al., 2005). Also, in strain Sa9 the light yellow production is likely caused by a mutation in gene crtP, whose product is responsible for the oxidation of the terminal methyl group of 4,4′-diaponeurosporene to form 4,4′-diaponeurosporenic acid (Pelz et al., 2005) (Figure 3B).
All but two strains (132, IPLA11) displayed haemolytic activity on blood agar. This is consistent with the presence in all strains of the hla (alpha-hemolysin) and hlb (beta-hemolysin) genes except in strains JE2 and 132, which lacked an intact hlb. In the case of strain 132, the absence of haemolytic activity is likely due to a mutation in agrC, which might prevent expression of the quorum-sensing-dependent gene hla gene (Le and Otto, 2015). We could not pinpoint the mutation(s) behind the negative phenotype of strain IPLA11.The highest proteolytic activity was observed in the clinical strains SH1000 and Newman, and the food industry isolates IPLA5, IPLA13 and IPLA15. In contrast, strains Sa9 and IPLA19 (mastitis) together with IPLA11 and IPLA16 (meat industry) showed very low or no protease production, respectively (Figure 4A). Proteases allow Staphylococcus to degrade host tissue or food components (in the case of industry isolates) to obtain nutrients (Dinges et al., 2000). All strains showed some degree of colony spreading, as can be observed in Figure 4B. This phenotype provides an indirect measure of the production of PSMs, which are known virulence factors. The strains with the largest colony diameter were both from mastitis, namely Sa9 (27.33 mm) and IPLA19 (25 mm). These two strains showed a mutation in gene hld that may lead to the production of a truncated protein, which has been demonstrated to cause overproduction of PSMs. Mutations in hld were also found in 20% of the mastitis isolates included in Supplementary Table 3, as well as in 8% of strains from the food industry. In contrast, the strains with the smallest diameter both belonged to the clinical sector, SH1000 (8.67 mm) and 132 (8.67 mm); this latter strain had a mutation in the quorum-sensing gene agrC.
Figure 4 Protease production and spreading capacity of S. aureus isolates. (A) Protease production of 18 S. aureus strains from different origin. (B) Expansion of 18 S. aureus strains from different origin. Bars represent mean ± standard deviation of three biological replicates of each strain. Bars within each strain having distinct lower case letter indicate that the biofilm production is statistically different (p < 0.05). Pink colour represents the strains with the largest displacement diameter (mm).
3.5 Antibiotic resistance
All strains displayed total (R) o partial resistance (I), although some were only resistant to one (16.7%, 3/18) or two (16.6%, 3/18) antibiotics, while others exhibited a multiresistance profile, being resistant to three (22.2%, 4/18), four (27.8%, 5/18), five (11.1%, 2/18) or even seven (1/18, IPLA3) antibiotics. The antibiotics against which we found the highest number of resistant strains were gentamicin (13 strains) and tobramycin (12 strains), followed by ciprofloxacin to which all strains showed intermediate or total resistance (Table 3). It is worth noting that all mastitis isolates were susceptible to methicillin. IPLA3 was the only tetracycline resistant strain. In most cases, these phenotypes were in agreement with the genotyping analysis (Table 3), although there were some exceptions. For instance, strain IPLA15 was methicillin resistant despite lacking gene mecA, and some tetracycline-resistant strains did not carry any specific genetic determinant.
4 Discussion
S. aureus is a ubiquitous microbe that can cause infections in different animals and contaminate the food environment, potentially leading to food poisoning. As strains spread from one environment to another, they carry their antibiotic resistance and virulence determinants. However, adaptation to the new milieu will involve changes in the bacterial genome due to mutations and genetic exchanges with neighboring bacteria. In this context, it is important to determine how isolates from diverse sources differ from each other and identify genetic signatures that characterize a given environmental niche. This is especially significant from the perspective of food safety. Strains isolated from human patients typically produce virulence factors and exhibit resistance to antibiotics commonly used in therapy. However, if these same traits are shared by isolates spread along the food chain this would pose a serious danger to the general public. Here, several strains from food surfaces and from cows with mastitis were analyzed and compared to reference strains of clinical origin.
In silico genotyping revealed that the studied strains are highly diverse, although some patterns can be observed. The most frequent MLST type is ST8, which appears only in several strains of clinical origin. This type was also one of the most frequent in hospital MRSA and MSSA strains in a genomic sequencing study carried out by (Khan et al. (2022a; 2022b). After that, the most common genotypes are ST97 and ST1. ST97 is the most frequent genotype involved in bovine mastitis, but rare in human infections. However, Manara et al. (2018) identified this MLST type in clinical strains from an Italian paediatric hospital. In turn, ST1 has been associated with MRSA community-acquired infections (Yamamoto et al., 2013). The predominant spa type in this work was t008, as was the case for strains of the same origin in previous studies (Khan et al., 2022a; Khan et al., 2022b). This type is frequently associated with ST8 and SCCmec type IV (Manara et al., 2018), as observed for strains JE2 and 132. Other frequent genotypes, t359 and t529, were found in mastitis isolates. T529 was also prevalent in previous studies concerning S. aureus presence in raw milk such as Oliveira et al. (2022). Regarding capsule type, all strains were CP5 or CP8, often referred to as microencapsulated nonmucoid strains, which are typical in both clinical and animal strains (O'Riordan and Lee, 2004; Salimena et al., 2016). The most prevalent agr type was type I, followed by type II and type III. Type I was also the most abundant in Chinese S. aureus strains of food and human origin (Wang et al., 2021). These source-dependent grouping based on genotyping is also reflected in the phylogenetic tree obtained by SNP analysis.
Clinical strains possess, on average, a larger number of virulence genes than isolates from other origins. In the case of toxins, mastitis strains also have a higher number of determinants than food industry strains. Even though these differences are not statistically significant, they do correlate with the fact that these strains were successful in causing disease. Given the importance of enterotoxins in staphylococcal food poisoning, we looked more closely at the prevalence of enterotoxin-coding genes. Notably, one human and three mastitis isolates have six or seven enterotoxin genes. Mastitis is known to be often related to enterotoxin production (Rall et al., 2014). The high prevalence of enterotoxins in the milk of mastitic cows highlights the importance of controlling this disease in dairy cows and preventing the use of contaminated milk down the dairy chain. It is worth noting that one of the two most frequent enterotoxins involved in staphylococcal foodborne diseases, SEB, is not present in any of the strains studied here (Pinchuk et al., 2010). SEC, the most common enterotoxin involved in bovine mastitis, is present in only two out of five mastitis isolates (Fang et al., 2019). In our work, mastitis isolates had a greater prevalence of non-classical enterotoxin genes than strains from other origins. However, we did not observe this trend in strains from other studies (Supplementary Table 3).
In terms of antibiotic resistance determinants, 6 strains from our collection as well as four of the reference strains have no detectable resistance markers in their genomes. Additionally, none of the strains analyzed in this work display resistance to vancomycin, erythromycin, linezolid or chloramphenicol. In the case of erythromycin, this differs from other studies on strains from food and clinical settings where the prevalence was high (Fernandes et al., 2022; Li et al., 2022). In contrast, the results obtained for linezolid and chloramphenicol are similar to those reported for S. aureus from human and food origin (Wang et al., 2021) and milk (Oliveira et al., 2022) where resistance to these compounds was absent or low. This is important as linezolid is a new bacteriostatic antimicrobial agent with activity against drug-resistant staphylococci, including VISA and VRSA strains (Lowy, 2003). We observed that resistance to beta-lactams is the most prevalent, in some cases due to blaZ and, in others, to mecA, the gene responsible for methicillin resistance. Only isolate (IPLA3) possessed the latter gene, but IPLA15 was also resistant to this antibiotic. This phenomenon is known as borderline resistance, and is explained in the literature as the result of beta-lactamase overproduction and/or modification of PBP genes (Chambers, 1997). Since IPLA15 does have blaZ, the first alternative appears as a plausible explanation. While MRSA clinical strains carry type IV SCCmec, the dairy isolate IPLA3 has type V SCCmec. Type IV and V are relatively new types that replaced the previously prevalent types I and II. Interestingly, strain 132 has plasmid pUB110 (McKenzie et al., 1986; Gennimata et al., 1996) inserted into SCCmec type IVc that harbours kanamycin/tobramycin and bleomycin resistance genes as observed in previous isolates (Manara et al., 2018; Wang et al., 2021). Resistance to aminoglycosides was also frequent, generally due to the presence of ant(4’)-Ia, aph(3’)Ia and aph(6)-Ic (str) genes. In strains lacking these genes, adaptive resistance may be the reason for the resistant phenotypes, as described in previous studies (Chandrakanth et al., 2008). Some clinical strains also carry mutations that confer resistance to quinolones, but none of the food or mastitis isolates do.
The spread of virulence and antibiotic resistance genes can be greatly facilitated by mobile genetic elements. In this study, all strains had at least one pathogenicity island in their genomes, with IPLA13 being the strain with the highest number (13 SaPIs). However, these S. aureus pathogenicity islands (SaPIs) require a helper prophage for their mobilization (Novick and Ram, 2017). For instance, Dearborn and Dokland (2012) demonstrated that of the 4 prophages present in S. aureus Newman, wNM1 and wNM2 can mobilize SaPI1 and SaPIbov1. In this regard, most strains harbour intact temperate phages in their genomes. In addition to the aforementioned function, temperate phages might contribute to resistance gene transfer by transduction, and/or alter the expression of virulence factors in the host bacterium (Jurado et al., 2022). It is also worth mentioning the potential presence of plasmids in all but 7 strains of this study.
In addition to the prevalence of virulence and antibiotic resistance determinants, we also explored some phenotypes of the different strains, namely biofilm formation, staphyloxanthin biosynthesis, hemolysin production, proteolytic activity and colony spreading. Only three of the tested strains are strong biofilm formers, all of which were reference strains. These structures increase the persistence of S. aureus on surfaces, materials and food, complicating their elimination and treatment (Avila-Novoa et al., 2018). The mastitis strain V329 carries the biofilm-associated gene, bap, which allows a non-dependent form of polysaccharide intercellular adhesin (PIA) production for the formation of resistant cell aggregates (Lasa and Penadés, 2006). This gene is not present in any of the food or mastitis isolates from our collection. On the other hand, the strain 132 mutation in agrC could explain its high biofilm formation, given that the agr system participates in biofilm dispersion (Yarwood and Schlievert, 2003; Le and Otto, 2015). We should also mention the presence of genes involved in biofilm attachment (clfA, clfB, atl, ebp, efb, fnbA, eap-map) in most of the strains in our collection, regardless of their origin. However, it is also noteworthy that approximately 50% of the strains lack an intact icaD gene, which might explain their weak biofilm development. The strains were also quite homogeneous regarding other phenotypes. For instance, most strains produce staphyloxanthin, a pigment with antioxidant properties that defends the microorganism from the host immune system, and display both haemolytic and proteolytic activity. Fernandes et al. (2022) also demonstrated haemolytic activity in S. aureus strains from food handlers. However, this same study also found that 14.4% of the strains exhibited proteolytic activity, far from the percentage found in our collection of food strains (71%). In the food industry sector, this activity represents an important adaptation to the environment (Shaw et al., 2004), allowing the strains to use milk as a substrate (Anderson, 1976). Another study carried out with a CA-MRSA mutant strain demonstrated the importance of proteases for growth in peptide-rich environments, serum, in the presence of antimicrobial peptides (AMPs), and in human blood (Kolar et al., 2013).
Although S. aureus is regarded as a non-motile microorganism, we now know that cells can passively expand across the surface of soft agar plates through the production of surfactants (PSMs), a phenomenon called sliding or colony spreading (Pollitt et al., 2015). Apart for this function, PSMs exert many others related to pathogenesis (lysis of leukocytes and erythrocytes, stimulation of inflammatory responses, contribute to biofilm development), increasing the virulence of S. aureus strains (Gordon et al., 2014). Depending on their length, PSMs can be classified into alpha-type (20 to 25 amino acids) and beta-type (43 to 45 amino acids), with α PSMs conferring greater virulence (Wang et al., 2007). In our collection, all strains showed β-type PSMs, but only the clinical strains carried the alpha-type PSMs. This would explain why all strains in the study showed some degree of spreading. However, it is noteworthy that strains Sa9 and IPLA19, both from milk of mastitic cows, showed a high spreading capacity. This may be due to a mutation in the hld gene present in both strains. This PSM is regulated by the quorum-sensing system, and its production has been linked to a partial inhibition of spreading (Omae et al., 2012; Gordon et al., 2014). Mutations of the hld gene were also found in strains from other studies (Karki et al., 2020; Merda et al., 2020; Naushad et al., 2020). Given the role of PSMs in virulence, it would be interesting to track the prevalence of hld truncation in strains involved in cattle and human infections.
5 Concluding remarks
Tracking the prevalence of antibiotic resistance and virulence genes and phenotypes in S. aureus strains is very valuable in order to predict the potential risk of such traits spreading. Most studies of this kind focus on clinical isolates due to their direct impact on human health. However, it would be useful to intensify the characterization of S. aureus strains isolated in the food environment, since food would make it easier for isolates with a high pathogenic potential to reach the wider community. Here, we have focused on food isolates from Northern Spain and our results highlight the genetic diversity of such S. aureus isolates, including those isolated from the same location or within a short period of time. Nonetheless, there is a clear correlation between genetic proximity and source of the isolate. Some of the genotypic and phenotypic trends that we found in isolates coming from the same environment are similar to those identified in other studies. Nonetheless, other characteristics, such as the prevalent types of enterotoxins in mastitis isolates, seem to be specific to this particular collection, at least considering the information available at the moment. Finally, another interesting finding was the relatively low frequency of toxin-encoding genes and antibiotic resistance determinants in the food strains compared to strains from mastitic cows and to reference strains of clinical origin. Overall, our results suggest that possession of genes related to virulence and antibiotic resistance might not confer an advantage to S. aureus in the food environment. Therefore, strains carrying these determinants in foods might likely come from infected humans or animals. Nonetheless, it is important to remain vigilant and perform studies that track the prevalence of these genes and their corresponding phenotypes in food industry isolates and not only in clinical or veterinary strains.
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 below: https://www.ncbi.nlm.nih.gov/genbank/, CP134618; https://www.ncbi.nlm.nih.gov/genbank/, CP134619; https://www.ncbi.nlm.nih.gov/genbank/, CP134617; https://www.ncbi.nlm.nih.gov/genbank/, CP134620; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXU000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXT000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXS000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXR000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXQ000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXV000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXX000000000; https://www.ncbi.nlm.nih.gov/genbank/, JAVRXW000000000.
Author contributions
AJ: Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft. LF: Conceptualization, Data curation, Formal analysis, Methodology, Supervision, Visualization, Funding acquisition, Writing – review & editing. AR: Conceptualization, Funding acquisition, Writing – review & editing. PG: Conceptualization, Funding acquisition, Supervision, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by grants MCIN/AEI/10.13039/501100011033/FEDER, UE, grant number PID2022-140988OB-I00, awarded to PG and LF, PID2019-105311RB-I00 (MICIU/AEI/FEDER, UE, Spain) to PG and AR, and AYUD/2021/52120 (Program of Science, Technology and Innovation 2021-2023, FICYT and FEDER EU, Principado de Asturias, Spain). AJ was awarded a JAE-Intro Scholarship (CSIC, Spain) and is currently a fellow of the program “Ayudas Severo Ochoa” of the Principality of Asturias.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2024.1327131/full#supplementary-material
References
Ajose, D. J., Oluwarinde, B. O., Abolarinwa, T. O., Fri, J., Montso, K. P., Fayemi, O. E., et al. (2022). Combating bovine mastitis in the dairy sector in an era of antimicrobial resistance: ethno-veterinary medicinal option as a viable alternative approach. Front. Vet. Sci. 9. doi: 10.3389/fvets.2022.800322
Anderson, J. C. (1976). Mechanisms of staphylococcal virulence in relation to bovine mastitis. Br. Vet. J. 132 (3), 229–245. doi: 10.1016/s0007-1935(17)34682-1
Andrews, S. (2010) FastQC: a quality control tool for high throughput sequence data. Available at: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (Accessed November 22, 2023).
Argudín, M. Á., Mendoza, M. C., Rodicio, M. R. (2010). Food poisoning and Staphylococcus aureus enterotoxins. Toxins 2 (7), 1751–1773. doi: 10.3390/toxins2071751
Arndt, D., Grant, J. R., Marcu, A., Sajed, T., Pon, A., Liang, Y., et al. (2016). HASTER: a better, faster version of the PHAST phage search tool. Nucleic. Acids. Res. 44 (W1), W16–W21. doi: 10.1093/nar/gkw387
Avila-Novoa, M.-G., Iniguez-Moreno, M., Solis-Velazquez, O. A., Gonzalez-Gomez, J. P., Guerrero-Medina, P. J., Gutierrez-Lomeli, M. (2018). Biofilm formation by staphylococcus aureus isolated from food contact surfaces in the dairy industry of jalisco, Mexico. J. Food. Qual. 2018, 1746139. doi: 10.1155/2018/1746139
Ayliffe, A. J., World Health Organization, Division of Emerging and other Communicable Diseases Surveillance Control (1996) Recommendations for the control of methicillin-resistant Staphylococcus aureus (MRSA) (World Health Organization). Available at: https://apps.who.int/iris/handle/10665/62984 (Accessed November 22, 2023).
Bankevich, A., Nurk, S., Antipov, D., Gurevich, A. A., Dvorkin, M., Kulikov, et al. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.: J. Comput. Mol. Cell Biol. 19 (5), 455–477. doi: 10.1089/cmb.2012.0021
Bartels, M. D., Petersen, A., Worning, P., Nielsen, J. B., Larner-Svensson, H., Johansen, H. K., et al. (2014). Comparing whole-genome sequencing with Sanger sequencing for spa typing of methicillin-resistant Staphylococcus aureus. J. Clin. Microbiol. 52 (12), 4305–4308. doi: 10.1128/JCM.01979-14
Bertelli, C., Laird, M. R., Williams, K. P., Fraser, S., Lau, B. Y., Hoad, G., et al. (2017). IslandViewer 4: expanded prediction of genomic islands for larger-scale datasets. Nucleic. Acids. Res. 45 (W1), W30–W35. doi: 10.1093/nar/gkx343
Biemer, J. J. (1973). Antimicrobial susceptibility testing by the Kirby-Bauer disc diffusion method. Ann. Clin. Lab. Sci. 3 (2), 135–140.
Bolger, A. M., Lohse, M., Usadel, B. (2014). Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinf. (Oxford England). 30 (15), 2114–2120. doi: 10.1093/bioinformatics/btu170
Bortolaia, V., Kaas, R. S., Ruppe, E., Roberts, M. C., Schwarz, S., Cattoir, V., et al. (2020). ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother. 75 (12), 3491–3500. doi: 10.1093/jac/dkaa345
Bridier, A., Sanchez-Vizuete, P., Guilbaud, M., Piard, J. C., Naïtali, M., Briandet, R. (2015). Biofilm-associated persistence of food-borne pathogens. Food. Microbiol. 45 (Pt B), 167–178. doi: 10.1016/j.fm.2014.04.015
Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al. (2009). BLAST+: architecture and applications. BMC. Bioinf. 10, 421. doi: 10.1186/1471-2105-10-421
Carattoli, A., Zankari, E., García-Fernández, A., Voldby Larsen, M., Lund, O., Villa, L. (2014). In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob. Agents. Chemother. 58 (7), 3895–3903. doi: 10.1128/AAC.02412-14
Castro, A., Silva, J., Teixeira, P. (2018). Staphylococcus aureus, a food pathogen: Virulence factors and antibiotic resistance. Foodborne. Dis. 15, 213–238. doi: 10.1016/B978-0-12-811444-5.00008-7
Centers for Disease Control and Prevention (CDC). (1999). Four pediatric deaths from community-acquired methicillin-resistant Staphylococcus aureus — Minnesota and North Dakota 1997-1999. MMWR. Morbid. Mortal. Weekly. Rep. 48 (32), 707–710.
Chambers, H. F. (1997). Methicillin resistance in staphylococci: molecular and biochemical basis and clinical implications. Clin. Microbiol. Rev 10(4), 781–791. doi: 10.1128/CMR.10.4.781
Chandrakanth, R. K., Raju, S., Patil, S. A. (2008). Aminoglycoside-resistance mechanisms in multidrug-resistant Staphylococcus aureus clinical isolates. Curr. Microbiol. 56 (6), 558–562. doi: 10.1007/s00284-008-9123-y
Clauditz, A., Resch, A., Wieland, K. P., Peschel, A., Götz, F. (2006). Staphyloxanthin plays a role in the fitness of Staphylococcus aureus and its ability to cope with oxidative stress. Infect. Immun. 74, 4950–4953. doi: 10.1128/IAI.00204-06
Cucarella, C., Solano, C., Valle, J., Amorena, B., Lasa, I., Penadés, J. R. (2001). ). Bap, a Staphylococcus aureus surface protein involved in biofilm formation. J. Bacteriol. 183 (9), 2888–2896. doi: 10.1128/JB.183.9.2888-2896.2001
Dearborn, A. D., Dokland, T. (2012). Mobilization of pathogenicity islands by Staphylococcus aureus strain Newman bacteriophages. Bacteriophage 2 (2), 70–78. doi: 10.4161/bact.20632
de Kraker, M. E., Wolkewitz, M., Davey, P. G., Koller, W., Berger, J., Nagler, J., et al. (2011). Clinical impact of antimicrobial resistance in European hospitals: excess mortality and length of hospital stay related to methicillin-resistant Staphylococcus aureus bloodstream infections. Antimicrob. Agents. Chemother. 55 (4), 1598–1605. doi: 10.1128/AAC.01157-10
Dinges, M. M., Orwin, P. M., Schlievert, P. M. (2000). Exotoxins of Staphylococcus aureus. Clin. Microbiol. Rev. 13, 16–34. doi: 10.1128/cmr.13.1.16
Duthie, E. S., Lorenz, L. L. (1952). Staphylococcal coagulase; mode of action and antigenicity. J. Gen. Microbiol. 6, 95–107. doi: 10.1099/00221287-6-1-2-95
EFSA and ECDC (European Food Safety Authority and European Centre for Disease Prevention and Control) (2022). The european union one health 2021 zoonoses report. EFSA. J. 2022. 20 (12), 7666, 273. doi: 10.2903/j.efsa.2022.7666
Fang, R., Cui, J., Cui, T., Guo, H., Ono, H. K., Park, C. H., et al. (2019). Staphylococcal enterotoxin C is an important virulence factor for mastitis. Toxins 11 (3), 141. doi: 10.3390/toxins11030141
Fernandes, A., Ramos, C., Monteiro, V., Santos, J., Fernandes, P. (2022). Virulence potential and antibiotic susceptibility of S. aureus strains isolated from food handlers. Microorganisms 10 (11), 2155. doi: 10.3390/microorganisms10112155
Fey, P. D., Endres, J. L., Yajjala, V. K., Widhelm, T. J., Boissy, R. J., Bose, J. L., et al. (2013). A genetic resource for rapid and comprehensive phenotype screening of nonessential Staphylococcus aureus genes. mBio 4 (1), e00537–e00512. doi: 10.1128/mBio.00537-12
Fisher, E. L., Otto, M., Cheung, G. Y. C. (2018). Basis of virulence in enterotoxin-mediated staphylococcal food poisoning. Front. Microbiol. 9. doi: 10.3389/fmicb.2018.00436
Flemming, H. C., Wingender, J. (2010). The biofilm matrix. Nature reviews. Microbiology 8 (9), 623–633. doi: 10.1038/nrmicro2415
García, P., Madera, C., Martínez, B., Rodríguez, A. (2007). Biocontrol of Staphylococcus aureus in curd manufacturing processes using bacteriophages. Int. Dairy. J. 17, 1232–1239. doi: 10.1016/j.idairyj.2007.03.014
Gennimata, D., Davies, J., Tsiftsoglou, A. S. (1996). Bleomycin resistance in Staphylococcus aureus clinical isolates. J. Antimicrob. Chemother. 37 (1), 65–75. doi: 10.1093/jac/37.1.65
Goerke, C., Pantucek, R., Holtfreter, S., Schulte, B., Zink, M., Grumann, D., et al. (2009). Diversity of prophages in dominant Staphylococcus aureus clonal lineages. J. Bacteriol. 191 (11), 3462–3468. doi: 10.1128/JB.01804-08
Gordon, Y. S., Joo, H. S., Chatterjee, S. S., Otto, M. (2014). Phenol-soluble modulins – critical determinants of staphylococcal virulence. FEMS Microbiol. Rev 38(4), 698–719. doi: 10.1111/1574-6976.1205
Gurevich, A., Saveliev, V., Vyahhi, N., Tesler, G. (2013). QUAST: quality assessment tool for genome assemblies. Bioinf. (Oxford England) 29 (8), 1072–1075. doi: 10.1093/bioinformatics/btt086
Gutiérrez, D., Delgado, S., Vázquez-Sánchez, D., Martínez, B., Cabo, M. L., Rodríguez, A., et al. (2012). Incidence of Staphylococcus aureus and analysis of associated bacterial communities on food industry surfaces. Appl. Environ. Microbiol. 78 (24), 8547–8554. doi: 10.1128/AEM.02045-12
Hennekinne, J. A., De Buyser, M. L., Dragacci, S. (2012). Staphylococcus aureus and its food poisoning toxins: characterization and outbreak investigation. FEMS Microbiol. Rev. 36 (4), 815–836. doi: 10.1111/j.1574-6976.2011.00311.x
Hasse-Cieślińska, M. (2007). Staphyloccoccus aureus skin colonization in atopic dermatitis patients. Postepy Dermatologii i Alergologii 24, 107–115.
Herrera, J. J., Cabo, M. L., González, A., Pazos, I., Pastoriza, L. (2007). Adhesion and detachment kinetics of several strains of Staphylococcus aureus subsp. aureus under three different experimental conditions. Food. Microbiol. 24 (6), 585–591. doi: 10.1016/j.fm.2007.01.001
Horsburgh, M. J., Aish, J. L., White, I. J., Shaw, L., Lithgow, J. K., Foster, S. J. (2002). sigmaB modulates virulence determinant expression and stress resistance: characterization of a functional rsbU strain derived from Staphylococcus aureus 8325-4. J. Bacteriol. 184 (19), 5457–5467. doi: 10.1128/JB.184.19.5457-5467.2002
Hurley, J. C. (2018). Unusually High Incidences of Staphylococcus aureus Infection within Studies of Ventilator Associated Pneumonia Prevention Using Topical Antibiotics: Benchmarking the Evidence Base. Microorganisms 6 (1), 2. doi: 10.3390/microorganisms6010002
International Working Group on the Classification of Staphylococcal Cassette Chromosome Elements (IWG-SCC) (2009). Classification of staphylococcal cassette chromosome mec (SCCmec): guidelines for reporting novel SCCmec elements. Antimicrob. Agents. Chemother. 53 (12), 4961–4967. doi: 10.1128/AAC.00579-09
Johler, S., Giannini, P., Jermini, M., Hummerjohann, J., Baumgartner, A., Stephan, R. (2015). Further evidence for staphylococcal food poisoning outbreaks caused by egc-encoded enterotoxins. Toxins 7 (3), 997–1004. doi: 10.3390/toxins7030997
Jurado, A., Fernández, L., Rodríguez, A., García, P. (2022). Understanding the mechanisms that drive phage resistance in staphylococci to prevent phage therapy failure. Viruses 14 (5), 1061. doi: 10.3390/v14051061
Kaas, R. S., Leekitcharoenphon, P., Aarestrup, F. M., Lund, O. (2014). Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PloS. One. 9 (8), e104984. doi: 10.1371/journal.pone.0104984
Kadariya, J., Smith, T. C., Thapaliya, D. (2014). Staphylococcus aureus and staphylococcal food-borne disease: an ongoing challenge in public health. BioMed. Res. Int. 2014, 827965. doi: 10.1155/2014/827965
Kanafani, Z. A., Fowler, V. G., Jr. (2006). Staphylococcus aureus infections: new challenges from an old pathogen. Enfermedades. Infecciosas. y Microbiol. Clin. 24 (3), 182–193. doi: 10.1157/13086552
Karki, A. B., Neyaz, L., Fakhr, M. K. (2020). Comparative genomics of plasmid-bearing staphylococcus aureus strains isolated from various retail meats. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.574923
Khan, S. A., Gudeta, D. D., Aljahdali, N., Chon, J., Snippes Vagnone, P., Nawaz, M. S., et al. (2022a). Draft genome sequences of 27 hospital-associated methicillin-resistant Staphylococcus aureus strains isolated in minnesota. Microbiol. Resour. Announcements. 11 (2), e0118621. doi: 10.1128/mra.01186-21
Khan, S. A., Gudeta, D. D., Chon, J., Snippes Vagnone, P., Nawaz, M. S., Foley, S. L., et al. (2022b). Whole-genome sequences of nine hospital-associated methicillin-susceptible Staphylococcus aureus strains. Microbiol. Resour. Announcements. 11 (2), e0118521. doi: 10.1128/mra.01185-21
Kolar, S. L., Ibarra, J. A., Rivera, F. E., Mootz, J. M., Davenport, J. E., Stevens, S. M., et al. (2013). Extracellular proteases are key mediators of Staphylococcus aureus virulence via the global modulation of virulence-determinant stability. MicrobiologyOpen 2 (1), 18–34. doi: 10.1002/mbo3.55
Kreiswirth, B. N., Löfdahl, S., Betley, M. J., O'Reilly, M., Schlievert, P. M., Bergdoll, M. S., et al. (1983). The toxic shock syndrome exotoxin structural gene is not detectably transmitted by a prophage. Nature 305 (5936), 709–712. doi: 10.1038/305709a0
Larsen, M. V., Cosentino, S., Rasmussen, S., Friis, C., Hasman, H., Marvig, R. L., et al. (2012). Multilocus sequence typing of total-genome-sequenced bacteria. J. Clin. Microbiol. 50 (4), 1355–1361. doi: 10.1128/JCM.06094-11
Lasa, I., Penadés, J. R. (2006). Bap: a family of surface proteins involved in biofilm formation. Res. Microbiol. 157 (2), 99–107. doi: 10.1016/j.resmic.2005.11.003
Le, K. Y., Otto, M. (2015). Quorum-sensing regulation in staphylococci-an overview. Front. Microbiol. 6, 1174. doi: 10.3389/fmicb.2015.01174
Lee, A. S., de Lencastre, H., Garau, J., Kluytmans, J., Malhotra-Kumar, S., Peschel, A., et al. (2018). Methicillin-resistant Staphylococcus aureus. Nature reviews. Dis. Primers. 4, 18033. doi: 10.1038/nrdp.2018.33
Li, X., Zhang, J., Zhang, H., Shi, X., Wang, J., Li, K., et al. (2022). Genomic analysis, antibiotic resistance, and virulence of Staphylococcus aureus from food and food outbreaks: A potential public concern. Int. J. Food Microbiol. 377, 109825. doi: 10.1016/j.ijfoodmicro.2022.109825
Liu, B., Zheng, D., Jin, Q., Chen, L., Yang, J. (2019). VFDB 2019: a comparative pathogenomic platform with an interactive web interface. Nucleic. Acids. Res. 47 (D1), D687–D692. doi: 10.1093/nar/gky1080
Lowy, F. D. (2003). Antimicrobial resistance: the example of Staphylococcus aureus. J. Clin. Invest. 111 (9), 1265–1273. doi: 10.1172/JCI18535
Malak, H. A., Abulreesh, H. H., Organji, S. R., Elbanna, K., Shaaban, M. R., Samreen, , et al. (2020). Immune system evasion mechanisms in Staphylococcus aureus: current understanding. J. Pure. Appl. Microbiol. 14 (4), 2219–2234. doi: 10.22207/JPAM.14.4.01
Manara, S., Pasolli, E., Dolce, D., Ravenni, N., Campana, S., Armanini, F., et al. (2018). Whole-genome epidemiology, characterisation, and phylogenetic reconstruction of Staphylococcus aureus strains in a paediatric hospital. Genome. Med. 10, 82. doi: 10.1186/s13073-018-0593-7
McGuinness, W. A., Malachowa, N., DeLeo, F. R. (2017). Vancomycin resistance in staphylococcus aureus. Yale. J. Biol. Med. 90 (2), 269–281.
McKenzie, T., Hoshino, T., Tanaka, T., Sueoka, N. (1986). The nucleotide sequence of pUB110: some salient features in relation to replication and its regulation. Plasmid 15, 93–103. doi: 10.1016/0147-619x(86)90046-6
Merda, D., Felten, A., Vingadassalon, N., Denayer, S., Titouche, Y., Decastelli, L., et al. (2020). NAuRA: genomic tool to identify staphylococcal enterotoxins in Staphylococcus aureus strains responsible for foodBorne outbreaks. Front. Microbiol. 11. doi: 10.3389/fmicb.2020.01483
Naushad, S., Nobrega, D. B., Naqvi, S. A., Barkema, H. W., De Buck, J. (2020). Genomic analysis of bovine Staphylococcus aureus isolates from milk to elucidate diversity and determine the distributions of antimicrobial and virulence genes and their association with mastitis. mSystems 5 (4), e00063–e00020. doi: 10.1128/mSystems.00063-20
Novick, R. P., Ram, G. (2017). Staphylococcal pathogenicity islands-movers and shakers in the genomic firmament. Curr. Opin. Microbiol. 38, 197–204. doi: 10.1016/j.mib.2017.08.001
O'Riordan, K., Lee, J. C. (2004). Staphylococcus aureus capsular polysaccharides. Clin. Microbiol. Rev. 17 (1), 218–234. doi: 10.1128/CMR.17.1.218-234.2004
Oliveira, R., Pinho, E., Almeida, G., Azevedo, N. F., Almeida, C. (2022). Prevalence and diversity of staphylococcus aureus and staphylococcal enterotoxins in raw milk from Northern Portugal. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.846653
Olsen, I. (2015). Biofilm-specific antibiotic tolerance and resistance. European journal of clinical microbiology & infectious diseases: official publication of the European. Soc. Clin. Microbiol. 34 (5), 877–886. doi: 10.1007/s10096-015-2323-z
Omae, Y., Sekimizu, K., Kaito, C. (2012). Inhibition of colony-spreading activity of Staphylococcus aureus by secretion of δ-hemolysin. J. Biol. Chem. 287 (19), 15570–15579. doi: 10.1074/jbc.M112.357848
Pelz, A., Wieland, K. P., Putzbach, K., Hentschel, P., Albert, K., Götz, F. (2005). Structure and biosynthesis of staphyloxanthin from Staphylococcus aureus. J. Biol. Chem. 280 (37), 32493–32498. doi: 10.1074/jbc.M505070200
Pinchuk, I. V., Beswick, E. J., Reyes, V. E. (2010). Staphylococcal enterotoxins. Toxins 2 (8), 2177–2197. doi: 10.3390/toxins2082177
Pofahl, W. E., Ramsey, K. M., Nobles, D. L., Cochran, M. K., Goettler, C. (2011). Importance of methicillin-resistant Staphylococcus aureus eradication in carriers to prevent postoperative methicillin-resistant Staphylococcus aureus surgical site infection. Am. Surgeon. 77 (1), 27–31. doi: 10.1177/0003134811077001
Pollitt, E. J., Crusz, S. A., Diggle, S. P. (2015). Staphylococcus aureus forms spreading dendrites that have characteristics of active motility. Sci. Rep. 5, 17698. doi: 10.1038/srep17698
Rall, V. L., Miranda, E. S., Castilho, I. G., Camargo, C. H., Langoni, H., Guimarães, F. F., et al. (2014). Diversity of Staphylococcus species and prevalence of enterotoxin genes isolated from milk of healthy cows and cows with subclinical mastitis. J. Dairy. Sci. 97 (2), 829–837. doi: 10.3168/jds.2013-7226
Salimena, A. P., Lange, C. C., Camussone, C., Signorini, M., Calvinho, L. F., Brito, M. A., et al. (2016). Genotypic and phenotypic detection of capsular polysaccharide and biofilm formation in Staphylococcus aureus isolated from bovine milk collected from Brazilian dairy farms. Vet. Res. Commun. 40 (3-4), 97–106. doi: 10.1007/s11259-016-9658-5
Seemann, T. (2014). Prokka: rapid prokaryotic genome annotation. Bioinf. (Oxford England) 30 (14), 2068–2069. doi: 10.1093/bioinformatics/btu153
Shaw, L., Golonka, E., Potempa, J., Foster, S. J. (2004). The role and regulation of the extracellular proteases of Staphylococcus aureus. Microbiol. (Reading England) 150 (Pt 1), 217–228. doi: 10.1099/mic.0.26634-0
Smyth, C. J., Möllby, R., Wadström, T. (1975). Phenomenon of hot-cold hemolysis: chelator-induced lysis of sphingomyelinase-treated erythrocytes. Infect. Immun. 12 (5), 1104–1111. doi: 10.1128/iai.12.5.1104-1111.1975
Tatusova, T., Di Cuccio, M., Badretdin, A., Chetvernin, V., Nawrocki, E. P., Zaslavsky, L., et al. (2016). NCBI prokaryotic genome annotation pipeline. Nucleic. Acids. Res. 44, 6614–6624. doi: 10.1093/nar/gkw569
Vergara-Irigaray, M., Valle, J., Merino, N., Latasa, C., García, B., Ruiz de Los Mozos, I., et al. (2009). Relevant role of fibronectin-binding proteins in Staphylococcus aureus biofilm-associated foreign-body infections. Infect. Immun. 77 (9), 3978–3991. doi: 10.1128/IAI.00616-09
Vrieling, M., Boerhout, E. M., van Wigcheren, G. F., Koymans, K. J., Mols-Vorstermans, T. G., de Haas, C. J., et al. (2016). LukMF' is the major secreted leukocidin of bovine Staphylococcus aureus and is produced in vivo during bovine mastitis. Sci. Rep. 6, 37759. doi: 10.1038/srep37759
Wang, W., Baker, M., Hu, Y., Xu, J., Yang, D., Maciel-Guerra, A., et al. (2021). Whole-genome sequencing and machine learning analysis of Staphylococcus aureus from multiple heterogeneous sources in China reveals common. Genetic traits of antimicrobial resistance. mSystems 6 (3), e0118520. doi: 10.1128/mSystems.01185-20
Wang, R., Braughton, K. R., Kretschmer, D., Bach, T. H., Queck, S. Y., Li, M., et al. (2007). Identification of novel cytolytic peptides as key virulence determinants for community-associated MRSA. Nat. Med. 13 (12), 1510–1514. doi: 10.1038/nm1656
WHO (World Health Organization) (2021) Global antimicrobial resistance and use surveillance system (GLASS) report. Available at: https://www.who.int/publications/i/item/9789240027336 (Accessed November 22, 2023).
Wick, R. R., Judd, L. M., Gorrie, C. L., Holt, K. E. (2017). Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PloS. Comput. Biol. 13, e1005595. doi: 10.1371/journal.pcbi.1005595
Yamamoto, T., Hung, W. C., Takano, T., Nishiyama, A. (2013). Genetic nature and virulence of communityassociated methicillin-resistant Staphylococcus aureus. Biomedicine 3, 2–18. doi: 10.1016/j.biomed.2012.12.001
Keywords: bacteriophages, antibiotic resistance, food industry, biofilms, Staphylococcus aureus
Citation: Jurado A, Fernández L, Rodríguez A and García P (2024) Prevalence of virulence- and antibiotic resistance-associated genotypes and phenotypes in Staphylococcus aureus strains from the food sector compared to clinical and cow mastitis isolates. Front. Cell. Infect. Microbiol. 14:1327131. doi: 10.3389/fcimb.2024.1327131
Received: 24 October 2023; Accepted: 12 January 2024;
Published: 29 January 2024.
Edited by:
Ala E. Tabor, The University of Queensland, AustraliaReviewed by:
Muhammad Noman Naseem, The University of Queensland, AustraliaMichal Bukowski, Jagiellonian University, Poland
Copyright © 2024 Jurado, Fernández, Rodríguez and García. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lucía Fernández, lucia.fernandez@ipla.csic.es