Microsoft Word - 11546 NSB Bacila 2023.06.26.docx Received: 05 Apr 2023. Received in revised form: 16 Jun 2023. Accepted: 20 Jun 2023. Published online: 26 Jun 2023. From Volume 13, Issue 1, 2021, Notulae Scientia Biologicae journal uses article numbers in place of the traditional method of continuous pagination through the volume. The journal will continue to appear quarterly, as before, with four annual numbers. SHSTSHSTSHSTSHST Horticulture and ForestryHorticulture and ForestryHorticulture and ForestryHorticulture and Forestry Society of TransylvaniaSociety of TransylvaniaSociety of TransylvaniaSociety of Transylvania Băcilă I et al. (2023) Notulae Scientia BiologicaeNotulae Scientia BiologicaeNotulae Scientia BiologicaeNotulae Scientia Biologicae Volume 15, Issue 2, Article number 11546 DOI:10.15835/nsb15211546 Research ArticleResearch ArticleResearch ArticleResearch Article.... NSBNSBNSBNSB Notulae Scientia Notulae Scientia Notulae Scientia Notulae Scientia BiologicaeBiologicaeBiologicaeBiologicae Epidemiology of Epidemiology of Epidemiology of Epidemiology of Enterococcus faeciumEnterococcus faeciumEnterococcus faeciumEnterococcus faecium isolates sampled from different isolates sampled from different isolates sampled from different isolates sampled from different sources in Romania using MLST technique and eBURST algorithmsources in Romania using MLST technique and eBURST algorithmsources in Romania using MLST technique and eBURST algorithmsources in Romania using MLST technique and eBURST algorithm Ioan BĂCILĂ1a, Endre JAKAB2,3,5b, Dana ŞUTEU1*, Octavian POPESCU4,5,6 1National Institute of Research and Development for Biological Sciences, Institute of Biological Research, Department of Experimental Biology; 48 Republicii St., 400015 Cluj-Napoca, Romania; ioan.bacila@icbcluj.ro; dana.suteu@icbcluj.ro; dana.suteu@icbcluj.ro (*corresponding author) 2Babeș-Bolyai University, Faculty of Biology and Geology, Hungarian Department for Biology and Ecology, 5-7 Clinicilor St., 400006 Cluj-Napoca, Romania 3Babeș-Bolyai University, Centre for Systems Biology, Biodiversity and Bioresources, 5-7 Clinicilor St., 400006 Cluj-Napoca, Romania; endre.jakab@ubbcluj.ro; 4Institute of Biology Bucharest of the Romanian Academy, 296 Splaiul Independenței St., Bucharest, 060031, Romania 5Babeș-Bolyai University, Molecular Biology Centre, Interdisciplinary Research Institute on Bio-Nano-Sciences, 42 Treboniu Laurian St., Cluj-Napoca, 400394, Romania 6Babeș-Bolyai University, Emil G. Racoviță Institute, 5-7 Clinicilor St., Cluj-Napoca, 400006, Romania; opopescu.ubbcluj@gmail.com a,bThese authors contributed equally to the work AbstractAbstractAbstractAbstract Enterococcus faecium is emerging as an important cause of multidrug resistance and hospital acquired infections, special attention being paid to the vancomycin resistant species. Therefore, the characterization of pathogenic strains/isolates plays an important role in the epidemiology of infectious diseases. The enterococcal rate was determined from wastewaters in Cluj-Napoca area. As presence of E. faecium was detected, a number of isolates from wastewater, birds and humans were epidemiologically analyzed according to the MLST website. Comparisons were performed against a collection of available isolates, with multiple origins, contained in the MLST database. Out of the Enterococcus isolates collected from wastewater, 11 were identified as E. faecalis (40.74%); 8 as E. casseliflavus (29.62%); 5 as E. faecium (18.50%); 2 as E. gallinarum (7.40%) and one isolate as E. durans. Based on the MLST data and using the eBURST algorithm, the isolates of E. faecium sampled from Romania were split in three groups: one group comprised isolates from human hosts and wastewater (Cj316, 106/6, Cj197, Cj22, 129/6, Cj117, Cj24, 284/7, and 43/7), while the second (G9, G10-2, G7, G3-2, and G9-1) and the third group (G8, G6, and 40/7) originated from bird hosts. The rest of the isolates were not joined in a particular group, assuming the lack of a phylogenetic bond between these isolates. The obtained data suggested the existence of at least two phylogenetic lines of E. faecium in Romania: a line that had mainly human host prevalence, while in the other line the animal hosts dominated. Keywords:Keywords:Keywords:Keywords: eBURST; Enterococcus faecium; epidemiology; MLST https://www.notulaebiologicae.ro/index.php/nsb/index Băcilă I et al. (2023). Not Sci Biol 15(2):11546 2 IntroductionIntroductionIntroductionIntroduction Bacteria of the genus Enterococcus (formerly the ‘fecal’ or Lancefield group D streptococci) are ubiquitous microorganisms. They occur in large numbers in different types of soil, surface waters, vegetables, plant material and foods, especially those of animal origin (Giraffa, 2014), but have a predominant habitat in the gastrointestinal tract of humans and animals. E. faecium and E. faecalis are the predominant Gram-positive cocci in human stools, while E. faecium is the prevalent species in production animals like poultry, cattle, and pig, and E. mundtii and E. casseliflavus are found in plant sources (Klein, 2003). From an ecological point of view, the distribution of Enterococcus species varies throughout Europe: in Spain and the UK E. faecalis and E. faecium are the most commonly isolated species from both clinical and environmental sources, in Sweden E. faecium has a lower incidence and E. hirae a higher isolation rate, whereas in Denmark E. hirae is the dominant species, isolated mainly from slaughtered animals (Kühn et al., 2003). Enterococcus species play an important role in food industry (dairy production, storage of meat and vegetables - Foulquié Moreno et al., 2006) or as probiotics to treat diarrhea and improve immunity (Franz et al., 2011). Nevertheless, some species of Enterococcus, such as E. faecalis and E. faecium, were also reported to be associated with many infections, including urinary tract infection, bacteremia, endocarditis, neonatal infection and infection of the central nervous system (O’Driscoll and Crank, 2015). A much bigger problem is represented by the fact that, with the extensive use of antibiotics, Enterococcus species have developed resistance to many antibiotics (Zhong et al., 2017). Currently, E. faecium is emerging as an important cause of multidrug resistance and hospital acquired infections, special attention being paid to VRE (vancomycin resistant Enterococcus) (Lebreton et al., 2014; Adegoke et al., 2022; Correa-Martínez et al., 2022; Toc et al., 2022). Even if E. faecalis is responsible for about 80% of all enterococcal infections in humans, while E. faecium causes only about 20% of the infections, E. faecium represents the majority of the VRE (Sievert et al., 2013). Ecology and epidemiology studies of Enterococcus have reported that E. faecalis and E. faecium are being regularly isolated from cheese, fish, sausages, minced beef, and pork (Klein, 2003; Foulquié Moreno et al., 2006). The characterization of pathogenic strains plays an important role in the epidemiology of infectious diseases, generating the necessary information for the identification, tracking and intervention against epidemics (Tavanti et al., 2005). MLST (Multilocus Sequence Typing) is one of the techniques used in global epidemiology, aiming to identify strains with a high pathogenicity and, thus, providing an improved picture of the activity of bacteria within environment and human populations. Studies on microbial populations using the MLST technique are generally intended to estimate the genetic diversity (usually counted as the relative contribution of recombination and mutations per allele or per locus), as well as evaluate the relative impact of genetic dispersion and natural selection in the evolutionary history of these pathogens (Stefani and Agodi, 2000; Pérez-Losada et al., 2005). MLST is based on the sequence of housekeeping genes which exhibit in each strain a distinct numerical allelic profile, abbreviated to a unique identifier: the sequence type (ST). The relatedness between two strains can be then inferred by differences between the allelic profiles (Francisco et al., 2009). The possible patterns of evolutionary descent obtained within MLST can be further analyzed by eBURST. eBURST represents an advanced algorithm that involves partitioning the large number of MLST data sets into groups of related, non-overlapping STs, clonal complexes, and then discerns the best-fit pathway within each clonal complex to the founder genotype (Spratt et al., 2004). Thus, following a set of very simple rules, eBURST can be used to find out how the diversification of bacterial clones has occurred and to provide evidence for the emergence of clinically relevant clones (Feil et al., 2004). The aims of this study were i) – to determine the enterococcal ratio in wastewater from the Cluj-Napoca area; ii) - to determine using MLST the genetic relatedness of randomly selected E. faecium isolates from wastewaters, birds and humans, and confront the results against an international MLST database. Băcilă I et al. (2023). Not Sci Biol 15(2):11546 3 Materials and MethodsMaterials and MethodsMaterials and MethodsMaterials and Methods Bacterial isolates sampled from wastewater in Cluj-Napoca area. Cultures and identification Between June 2006 and August 2008, 27 water samples were collected from the Someşeni wastewater treatment plant in Cluj-Napoca, Romania. Residual water samples were isolated in pre-sterilized containers and analyzed within 6 hours of isolation, according to international standards. The water samples divided into 6/100 mL portions were subjected to filtration on a 0.45 μm pore diameter membrane. To ease the filtration process and to reduce the number of colonies, the samples were diluted with H2O UV/UP, the initial dilution being 1:4, followed by a subsequent dilution of 1:9. After the filtering process, the samples were placed directly on the culture medium in Petri dishes. Two types of selective culture media were used. The filters were initially placed on Petri dishes containing M Enterococcus Agar (MEA) (Fluka, Buchs, Switzerland), because the sodium azide, captan and nalidixic acid contained in this formula, are known to inhibit the growth of many species of bacteria and fungi, conferring selectivity to the medium. After casting the plates, 300 μL of 1% TTC solution (2,3,5-triphenyl tetrazolium chloride) were added to the surface of the medium. Within the bacterial cell, TTC was reduced to insoluble formazan, which conferred a dark pink color to the colonies. Cultivation of filters on MEA was performed in the oven, for 48 hours at a temperature of 41 °C. As a safety measure, with respect to the identity of the cultivated microorganisms, cultivation was carried out on a second type of medium, namely Esculin Iron Agar (Fluka, Buchs, Switzerland). This second medium was used to confirm the identity of the colonies based on the ability of enterococci to hydrolyze esculin. The procedure consisted in transferring the M Enterococcus Agar filter membranes onto Esculin Iron Agar and incubating them in the oven at 41 °C for 20 minutes. The hydrolysis of esculin resulted in black or dark brown colonies, confirming the presence of enterococci. For species identification of enterococcal isolates (Table 3) the following methods were used: classical biochemical identification by observing hemolysis pattern on blood agar and two standardized methods: API® 20 Strep from BioMérieux (a kit includes strips that contain up to 20 miniature biochemical tests for manual identification of microorganisms) and VITEK (a fully automated system which performs bacterial identification). DNA extraction of the bacterial isolates sampled from wastewater in Cluj-Napoca area The bacterial isolates were cultivated on LB (Luria-Bertani) liquid media overnight at a temperature of 37 °C. In order to disrupt the cell walls, a digestion with lysozyme (Sigma, St. Louis, USA) and proteinase K (Promega, Madison, USA) was performed, following the protocol described in Băcilă et al. (2007). DNA purification were performed using Wizard Genomic DNA Purification Kit (Promega, Madison, USA), according to the manufacturer’s recommendations. DNA quality was estimated on a 1% agarose gel stained with ethidium bromide, and DNA concentration was quantified using a NanoDrop 1000 Spectrophotometer (Thermo Fisher Scientific, Wilmington, USA). Identification of Enterococcus faecium isolates by PCR The predominant species in cultures are: E. faecalis and E. faecium. To distinguish between the two species of enterococci, one fragment of the ddl gene was amplified within PCR. The amplification reaction was carried out on isolates from wastewater (one sample from Bucharest was added to the 27 samples previously collected from Cluj-Napoca area), as well as on isolates from other sources: human and birds (Table 1). As positive controls, the strains E. faecium ATTC 35667 and E. faecalis ATTC 51299 were used. The PCR program was the following: initial denaturation for 2 min at 94 °C, followed by 30 cycles of 1 min at 94 °C, 1 min at 54 °C, and 1 min at 72 °C; and a final extension for 10 min at 72 °C. The specific primers were: ddlF1: 5’-GCAAGCCTTCTTAGAGA-3’ and ddlF2: 5’-CATCGTGTAAGCTAACTTC-3’ (Dutka-Malen et al., 1995). Băcilă I et al. (2023). Not Sci Biol 15(2):11546 4 MLST technique and DNA sequencing applied on Enterococcus faecium isolates sampled from different sources and different parts of Romania A number of 32 isolates of E. faecium with different origins were analyzed with MLST and eBURST. The samples were taken from the wastewater (five isolates including one from Bucharest and four out of the 27 samples previously collected from Cluj-Napoca area), from birds (12 isolates collected from their faeces) and from humans (15 isolates collected from hospitalized persons: samples of secretions of infected wounds, blood or peritoneal fluid) (Table 1). Sequence types of E. faecium isolates were determined employing a MLST scheme, performing amplification of housekeeping genes by PCR and subsequent Sanger sequencing of the PCR products, as previously described by Homan et al. (2002). Table 1Table 1Table 1Table 1. List of Enterococcus faecium isolates used in MLST analysis E. faeciumE. faeciumE. faeciumE. faecium isolate isolate isolate isolate code code code code OriginOriginOriginOrigin SampledSampledSampledSampled SSSSampling ampling ampling ampling locationlocationlocationlocation Year of isolationYear of isolationYear of isolationYear of isolation G5 Birds Faeces Cluj 2007 G6 Birds Faeces Cluj 2007 G7 Birds Faeces Cluj 2007 G8 Birds Faeces Cluj 2007 G9 Birds Faeces Cluj 2007 G10 Birds Faeces Cluj 2007 G12 Birds Faeces Cluj 2007 G2-2 Birds Faeces Timişoara 2007 G3-2 Birds Faeces Timişoara 2007 G5-2 Birds Faeces Timişoara 2007 G9-1 Birds Faeces Timişoara 2007 G10-2 Birds Faeces Timişoara 2007 105/6 Human Blood Bucharest 2006 106/6 Human Blood Bucharest 2006 113/6 Human Blood Bucharest 2006 129/6 Human Blood Bucharest 2006 156/6 Human Blood Bucharest 2006 40/7 Human Blood Bucharest 2007 43/7 Human Blood Bucharest 2007 283/7 Human Blood Cluj 2007 284/7 Human Blood Cluj 2007 Cj117 Human Wound secretion Cluj 2008 Cj197 Human Peritoneal fluid Cluj 2008 Cj316 Human Peritoneal fluid Cluj 2008 18/7 Human Blood Oradea 2006 19/7 Human Wound secretion Oradea 2006 41/6 Human Blood Timişoara 2006 Cj14 Wastewater - Cluj 2008 Cj20 Wastewater - Cluj 2008 Cj22 Wastewater - Cluj 2008 Cj24 Wastewater - Cluj 2008 B176 Wastewater - Bucharest 2003 Băcilă I et al. (2023). Not Sci Biol 15(2):11546 5 The typing scheme used regions of seven structural genes (gdh - glucose-6-phosphate dehydrogenase; purK - 5-(carboxyamino) imidazole ribonucleotide synthase; pst - phosphate ATP-binding cassette transporter; atpA - ATP synthase, alpha subunit; gyd - glyceraldehyde-3-phosphate dehydrogenase; adk - adenylate kinase; ddl - D-alanine ligase). These genes are well-conserved structures in the bacterial genome, so the likelihood of mutations over time is rather low, allowing them to be successful tools in the accurate genetic characterization of certain pathogenic bacterial species, as precise identification of infectious agent strains is essential for both epidemiological studies and public health decisions. The sequences of each amplified gene fragment were compared to all previously identified (alleles) sequences for that locus, thus determining the number of alleles for all seven loci. Combining the number of alleles from all the loci defined the allelic profile of the strain. Each different allelic profile was considered a ST-sequence type. Such a ST represented a suitable and clear label for each isolate. Characterization of bacterial isolates was achieved by comparing the sequences obtained with the sequences existing in the databases. The PCR method is particularly precise, safe and time-efficient, allowing a large number of samples to be analyzed in a relatively short time. The difficulty encountered was the setting of optimal temperatures for the amplification of genetic material. Optimal boost temperatures varied in a fairly wide range, starting at 47 °C for the purK gene and reaching 59 °C for the ddl gene (Table 2). The primers required to amplify fragments belonging to the seven structural genes are also presented in Table 2. Table 2Table 2Table 2Table 2. Primers used for gene amplification in MLST technique and their annealing temperature NameNameNameName Sequence (5’Sequence (5’Sequence (5’Sequence (5’----3’)3’)3’)3’) Annealing temperature (Annealing temperature (Annealing temperature (Annealing temperature (°C) gdh1 5’-GGC GCA CTA AAA GAT ATG GT-3’ 52 gdh2 5’-CCA AGA TTG GGC AAC TTC GTC CCA-3’ purK1 5’-GCA GAT TGG CAC ATT GAA AGT-3’ 47 purK2 5’-TAC ATA AAT CCC GCC TGT TTC/T-3’ pstS1 5’-TTG AGC CAA GTC GAA GCT GGA G-3’ 50 pstS2 5’-CGT GAT CAC GTT CTA CTT CC-3’ atpA1 5’-CGG TTC ATA CGG AAT GGC ACA-3’ 48 atpA2 5’-AAG TTC ACG ATA AGC CAC GG-3’ gyd1 5’-CAA ACT GCT TAG CTC CAA TGG C-3’ 58 gyd2 5’-CAT TTC GTT GTC ATA CCA AGC-3’ adk1 5’-TAT GAA CCT CAT TTT AAT GGG-3’ 58 adk2 5’-GTT GAC TGC CAA ACG ATT TT-3’ ddl1 5’-GAG ACA TTG AAT ATG CCT TAT G-3’ 59 ddl2 5’-AAA AAG AAA TCG CAC CG-3’ The amplification reaction was performed in a 50 μL volume including: 1x PCR buffer; 1.5 mM MgCl2; 0.2 mM of each dNTP; 1 μM forward primer; 1 μM reverse primer; GoTaq Polymerase 1.25 U; and 6 μL genomic DNA. The following amplification program was used: initial denaturation for 3 min at 94 °C, followed by 35 cycles of 30 sec at 94 °C, 30 sec at 47-59 °C, and 30 sec at 72 °C; and a final extension for 5 min at 72 °C. Prior to sequencing, a purification step with Promega Wizard SV Gel and PCR Clean-UP System (Promega Corporation, Madison, USA) was carried out according to the manufacturer’s manual. Sequencing of both strands was performed using DTCS Quick Start Master Mix; 3.2 μL primer (1 pmol/µL); 0.5-10 μL of template DNA; and 0-9.5 μL H2O UV/UP. The following thermal cycle parameters have been used: 96 °C - 2 min; 40 cycles of 96 °C - 20 sec, 50 °C - 20 sec, 60 °C - 4 min. Excess primers and labelled ddNTPs were removed by purification with DNA Clean & Concentrator™ (Zymo Research, Orange, USA) according to the standard recommended by the manufacturer. The samples (total volume 20 μL) were prepared prior to sequencing by Băcilă I et al. (2023). Not Sci Biol 15(2):11546 6 adding 20 μL of HiDi formamide and then were loaded onto CEQ 8000 Genetic Analyzer (Beckman-Coulter, Fullerton, USA). Data analysis After performing the sequencing reactions, the sequences were analyzed using BioEdit (Hall, 1999) and eBURST (http://www.eBurst.mlst.net). The purpose of eBURST algorithm was to identify, based on sequence data, groups of isolates with related genotypes in population, and, then, to establish the founding genotype of each group. In the next step, the algorithm determined also the descendants of the founding genotype, which would be able to represent foundational genotypes for other subgroups of bacterial isolates. The graphical form of the analyses performed by the eBURST algorithm was a radial diagram centered on an ancestral founding genotype. Results Results Results Results Species identification of enterococcal isolates sampled from wastewater in Cluj-Napoca area The identification of enterococcal species isolated from wastewater was quite difficult to achieve and required both a constant repetition of the performed analyses and a combined use of several analysis methods. 27 isolates were identified in the collected samples of wastewater and the results obtained by a classical method have been confirmed by at least one standardized method (Table 3). Table 3Table 3Table 3Table 3. Identification of enterococcal isolated sampled from wastewater IsolateIsolateIsolateIsolate Classical identification Classical identification Classical identification Classical identification (biochemical series)(biochemical series)(biochemical series)(biochemical series) APIAPIAPIAPI®®®® 20 Strep 20 Strep 20 Strep 20 Strep (BioMérieux)(BioMérieux)(BioMérieux)(BioMérieux) VITEKVITEKVITEKVITEK Cj01 E. casseliflavus E. casseliflavus - Cj02 E. faecalis E. faecalis 99.9% - Cj03 E. faecalis E. faecalis 99.9% - Cj04 E. faecalis E. faecalis 99.9% - Cj05 E. faecium E. faecium 94.3% - Cj06 E. faecalis E. faecalis 99.9% - Cj07 E. faecalis E. faecalis 99.9% - Cj08 E. faecalis E. faecalis 99.9% - Cj09 E. faecalis E. faecalis 99.9% - Cj10 E. faecalis E. faecalis 99.9% - Cj11 E. faecalis E. faecalis 99.9% - Cj12 E. faecalis E. faecalis 99.9% - Cj13 E. casseliflavus - E. casseliflavus 99.9% Cj14 E. faecium E. faecium 97.5% E. faecium 94.26% Cj15 E. casseliflavus E. casseliflavus 99.9% E. casseliflavus 99.9% Cj16 E. faecalis E. faecalis 99.8% E. faecalis 99% Cj17 E. casseliflavus - E. casseliflavus 99.9% Cj18 E. casseliflavus - E. casseliflavus 99.9% Cj19 E. casseliflavus - E. casseliflavus 91.56% Cj20 E. faecium E. faecium 94.8% E. faecium 99% Cj21 E. casseliflavus - E. casseliflavus 97.41% Cj22 E. faecium E. faecium 97.8% E. faecium 98.5% Cj23 E. durans E. durans 93.8% E. durans 99% Băcilă I et al. (2023). Not Sci Biol 15(2):11546 7 Cj24 E. faecium E. faecium 97.9% E. faecium 99% Cj25 E. gallinarum - E. gallinarum 99% Cj26 E. gallinarum - E. gallinarum 85.66% Cj27 E. casseliflavus - E. casseliflavus 98.06% The percent indicate the degree of confidence in accuracy of the results The dominance of E. faecalis and E. faecium in the wastewater system was documented by previous studies that have examined species distributions of Enterococcus in wastewater (Blanch et al., 2003; Moore et al., 2008). Moreover, Ferguson et al. (2013) in a study performed on marine water and wastewater, found that E. faecalis was the dominant species in a ratio of over 90%, while E. faecium was present in a ratio of about 5%, the last 5% being represented by other enterococcal species. Nevertheless, in the present case, the species hierarchy differed from expectance (Figure 1). E. faecalis was placed on the first position with 11 identified isolates (40.74% from the total of analyzed isolates), on the second place was the species E. casseliflavus represented by eight isolates (29.62%), E. faecium was represented by only five isolates (18.50%), whereas the rest of the species were identified as E. gallinarum - two isolates (7.40%) and E. durans - one isolate (3.70%). Figure 1.Figure 1.Figure 1.Figure 1. Enterococcal ratio of species (%) isolated from wastewaters Identification of the Enterococcus faecium species by PCR This genetic identification method for species of enterococcal isolates was based on the specific amplification of a D-alanine-D-alanine ligase-encoding fragment (ddl). Following migration of the PCR products on 1% agarose gel, in case of E. faecium isolates amplification of a 550 bp fragment was obtained, whereas no amplification was visible for E. faecalis isolates (Figure 2). Băcilă I et al. (2023). Not Sci Biol 15(2):11546 8 Figure 2Figure 2Figure 2Figure 2. Agarose gel pattern showing amplification of a ddl gene fragment which signaled the presence of E. faecium 1: E. faecium Cj14; 2: E. faecium Cj20; 3: E. faecium Cj22; 4: E. faecium Cj24; 5: E. faecium ATTC 35667 (positive control); M: molecular marker 1 kbp (Axygen Biosciences, Union City, USA); 6: E. faecium G5; 7: E. faecium 106/6; 8: E. faecium 283/7; 9: E. faecalis ATTC 51299 (positive control); 10: E. faecalis Cj02; 11: E. faecalis Cj12; 12: negative control (no DNA) The isolates codes are the same as in Tables 1 and 3 The amplification reaction was performed post species identification by classical, biochemical, and automated phenotypic methods (Table 3), as a confirmation of the previous results. Inference of the phylogenetic relationships between Enterococcus faecium isolates sampled from Romania using the MLST approach and eBURST algorithm Allelic variation and genetic diversity of Enterococcus faecium isolates A number of 32 isolates of E. faecium from different sources were included in the MLST assay (Table 1). For these isolates a variable number of unique locus alleles was obtained, ranging from four unique alleles for the adk gene to 12 unique alleles for the atpA gene (Table 4). Table 4Table 4Table 4Table 4. Allelic variation in seven structural genes from Enterococcus faecium isolates LocusLocusLocusLocus Length of sequence (bp)Length of sequence (bp)Length of sequence (bp)Length of sequence (bp) Number of unique allelesNumber of unique allelesNumber of unique allelesNumber of unique alleles gdh 530 9 purK 492 6 pst 583 6 atpA 556 12 gyd 395 7 adk 437 4 ddl 465 6 Following the analysis of 32 isolates of E. faecium, the existence of 31 STs was established. Isolates 41/6 and 105/6 exhibited the same allelic profile. Both isolates originated from humans, being sampled from the same medium - blood. Surprisingly, isolates 105/6 was isolated in Bucharest, while isolates 41/6 was isolated in Timișoara. A number of five isolates showed 100% similarity for six of the analyzed loci, with differences appearing only in atpA locus. These isolates were: Cj22, 106/6, 197, Cj316, and 113/6. The source of these isolates differed: isolate Cj22 originated from wastewater and was isolated in Cluj-Napoca, isolates 106/6 and 113/6 were of human origin and have been collected from blood in Bucharest, the other two isolates, Cj197 and Cj316 were sampled from peritoneal fluids collected from patients hospitalized in Cluj-Napoca. Băcilă I et al. (2023). Not Sci Biol 15(2):11546 9 Another number of five isolates of E. faecium presented 100% similarity for five out of seven analyzed loci. Of these, isolates 43/7 and 129/6 were sampled in Bucharest from blood. The other three isolates derived from Cluj-Napoca, Cj24 isolate being collected from wastewater, while Cj117 and 284/7 isolates are of human origin, collected from plague and blood. There was also a batch of eight isolates which exhibited a 100% similarity for four out of seven loci, despite the fact that sources of origin are quite heterogeneous: isolates 18/7 and 19/7 are of human origin, being collected in Oradea from blood, respectively plaque discharge; isolates Cj20 and 283/7 were taken from Cluj-Napoca, but from different sources (wastewater and human blood); isolate B176 was sampled from wastewater in Bucharest; the remaining three isolates are of animal origin, being sampled from poultry faeces collected in Timișoara (G3-2 and G9 isolates) and Cluj (G10). The rest of the isolates showed 100% similarity for three or less of the analyzed loci. Assortment of Enterococcus faecium isolates sampled from Romania Based on the MLST data (sequence data, allele profiles) and subsequent analysis with eBURST algorithm, the division of E. faecium isolates sampled from Romania revealed three groups. The first group included Cj316, 106/6, Cj197, Cj22, 129/6, Cj117, Cj24, 284/7, and 43/7 isolates (Figure 3a); the second group included G9, G10-2, G7, G3-2, and G9-1 isolates (Figure 3b); and the third group comprised G8, G6, and 40/7 isolates (Figure 3c). The rest of the isolates (G5-2, G5, G2-2, 283/7, B176, 156/6, 41/6, Cj20, 19/7, 18/7, Cj14, G10, and G12) were not joined by eBURST in a particular group, assuming the lack of a phylogenetic bond between these isolates. The lack of framing of these isolates was due to the small number of common loci. Figure 3Figure 3Figure 3Figure 3. Grouping of Enterococcus faecium isolates using eBURST algorithm. a) - isolates Cj316, 106/6, Cj197, Cj22, 129/6, Cj117, Cj24, 284/7, and 43/7; b) – isolates G9, G10-2, G7, G3-2, and G9-1; c) - isolates G8, G6, and 40/7 Inference of phylogenetic relationships between Enterococcus faecium isolates sampled from Romania and from other parts of the world using the MLST database and the eBURST algorithm The next step after determining the phylogenetic relationships between the E. faecium isolates sampled from Romania, was to assign these isolates to the general group of E. faecium isolates sampled from different parts of the world. The isolates collected from Romania were introduced into the MLST database, and after applying the eBURST algorithm, the sequent ratio revealed the division of the isolates into 20 distinct groups (data not shown). Subsequently, a radial dendrogram of the group 1 MLST was generated (Figure 4). Băcilă I et al. (2023). Not Sci Biol 15(2):11546 10 The following Romanian isolates were included in group 1 MLST: Cj316, 106/6, Cj197, Cj22, 129/6, Cj117, Cj24, 284/7, and 43/7. Isolates Cj24, 106/6, Cj22, Cj197, and Cj316 represented SLVs (single locus variants), whereas the other four isolates were DLV (double locus variants). The other Romanian isolates were distributed into two separate MLST groups, respectively group 5 MLST and group 11 MLST (data not shown). This grouping confirmed the results obtained by applying the eBURST algorithm exclusively for isolates from Romania, when these isolates formed groups 2 and 3 (Figure 3 b, c). DiscussionDiscussionDiscussionDiscussion Enterococcus species distribution in wastewater The dominance of E. faecium and E. faecalis is well documented in previous studies which examined the distribution of Enterococcus species in wastewater systems (Sinton and Donnison, 1994; Blanch et al., 2003; Moore et al., 2008). Similarly, the prevalence of E. casseliflavus is also consistent with clinical human fecal samples (Ruoff et al., 1990; Stern et al., 1994). Nevertheless, the enterococcal rate is very variable, depending on various factors such as: degree of urbanization, number of hospitals, and so on. While performing a study on a hospital wastewater treatment plant, Karimi et al., (2016) found that out of 315 enterococci isolates, 162 (51.42%) were identiZed as E. faecium, 87 (27.61%) as E. hirae, 35 (11.11%) as E. faecalis, 11 (3.5%) as E. gallinarum, 7 (2.22%) as E. casseliEavus, and 4 (1.26%) as E. avium. Băcilă I et al. (2023). Not Sci Biol 15(2):11546 11 Figure 4Figure 4Figure 4Figure 4.... Dendrogram of group 1 of Enterococcus faecium isolates from MLST database. Green color is used to highlight the isolates sampled from Romania included in this group. Assortment of Enterococcus faecium isolates sampled from Romania A firstfirstfirstfirst groupgroupgroupgroup of enterococcal isolates was formed after eBURST analysis around Cj316 isolate (Figure 3a). This isolate was sampled from peritoneal fluid of human source collected in Cluj-Napoca. Interestingly, among the SLVs belonging to this isolate, the isolate Cj197 (sampled from the peritoneal fluid collected from patients hospitalized in Cluj-Napoca) was noted, indicating a significant relationship between these two isolates. Another interesting connection was based on the observation that other two SLVs of the Cj316 isolate are collected from the wastewaters from Cluj-Napoca, Cj22 and Cj24. This observation led to the assumption of a connection between the isolates collected in the hospital environment and those taken from the environment (here represented by wastewaters). Băcilă I et al. (2023). Not Sci Biol 15(2):11546 12 Other close connection, in terms of origin, within this first group was also observed on the DLV level, involving two isolates sampled in Cluj, namely the isolate Cj117, sampled from human wound secretion, and isolate 284/7, collected from human blood. The other three SLVs of the group were represented by the isolates 43/7, 106/6, and 129/6. These isolates were also connected by the place of origin, all three isolates being taken from samples collected from patients from Bucharest, two from the blood (106/6 and 129/6) and 43/7 from sputum. A second groupsecond groupsecond groupsecond group of E. faecium consisted of five isolates, all of animal origin, sampled from poultry faeces taken from Cluj and Timișoara. According to the dendrogram (Figure 3b), the central isolate of this group was G9, sampled from Cluj. The group contained three SLVs, two of them being isolated in Timișoara and one in Cluj. The only DLV was represented by the isolate G9-1, collected in Timișoara. The central isolate of the third groupthird groupthird groupthird group (Figure 3c) is the G8 isolate, of animal origin. The group had two SLVs represented by the isolates G6 (sampled in Cluj from bird faeces) and 40/7 (collected from human blood in Bucharest). Based on data obtained by querying the MLST database and the use of the eBURST algorithm, some general conclusions regarding the local epidemiology could be drawn for the Romanian isolates of E. faecium. Thus, for 17 isolates, straightforward connections in terms of their origin were established. However, for 14 isolates (45.65%) the phylogenetic pattern could not been so easily explained. Some possible answers could include (but are not limited to) the high genetic variability of bacterial isolates (Martínez-Carranza et al., 2018), the high degree of genetic recombination (Shikov et al., 2022), the horizontal genetic transfer between bacterial isolates (Heuer and Smalla, 2007), and genetic dispersion (Oh et al., 2010). Inference of phylogenetic relationships between Enterococcus faecium isolates sampled from Romania and from other parts of the world Subsequent to the eBURST analysis of 1387 isolates of E. faecium sampled from different parts of the globe, 20 groups of isolates were obtained. The main group was constituted of 1135 isolates. The founder genotype of this group was an isolate with a ST=17 (Figure 4) and an allelic profile of 1 1 1 1 1 1 1. According to the database, the first isolate with ST=17 was collected from a patient in the United Kingdom in 1992. This lineage (ST=17) was characterized by ampicillin resistance, pathogenicity island, and was associated with hospital outbreaks (Willems et al., 2005). The distribution of these isolates (ST=17) was higher in Europe, probably due to the existence of more extensive studies: 47 isolates in Germany, 23 isolates in the UK, 20 isolates in Netherlands, 15 isolates in Italy, 4 isolates in France, 3 isolates in Greece, 2 isolates in Denmark and Poland, whereas only one isolate was isolated from Hungary, Norway, and Serbia. By comparison, few isolates with ST 17 were identified on other continents: 11 isolates on the Australian continent, 4 isolates in the North America (USA), while in South America (Brazil) only one isolate (Willems et al., 2005). The source of E. faecium ST=17 isolates was in an overwhelming proportion the hospitalized patients. Still, there were also some exceptions, in Australia (two cases) and in the UK and France (one case each), where the source of the isolates was represented by people in outpatient treatment. Furthermore, in two cases found in 2001 in the USA, the isolates were taken from the environment. Concerning the Romanian isolates, the following isolates were included in group 1 MLST: Cj316, 106/6, C197, C222, 129/6, C117, Cj24, 284/7, and 43/7. The isolates Cj24, Cj22 (both from Cluj wastewater), 106/6 (from human blood in Bucharest), Cj197, and Cj316 (both from human peritoneal fluid collected from patients hospitalized in Cluj) represent SLVs emerging directly from the founding genotype ST=17. The other four isolates represented DLVs, but were not joined in the same subgroup. The isolates Cj117 (human wound secretion, Cluj) and 129/6 (human blood, Bucharest) were part of the same DLV. Isolates 284/7 (human blood, Cluj) and 43/7 (human sputum, Bucharest) belonged to two different DLVs. Băcilă I et al. (2023). Not Sci Biol 15(2):11546 13 The presence of five lines of E. faecium sampled from Romania as SLV, pertaining directly to the ST=17 complex, is particularly alarming because ST=17 is a predictor of subsequent bacteremia in hospitalized patients carrying VREF and has been frequently isolated from nosocomial settings (Kim et al., 2021). ConConConConclusionsclusionsclusionsclusions Out of 27 Enterococcus isolates sampled from wastewater, 11 isolates were identified as E. faecalis (40.74%); 8 isolates as E. casseliflavus (29.62%); 5 isolates as E. faecium (18.50%); 3 isolates as E. gallinarum (7.40%) and one isolate as E. durans. By applying the eBURST algorithm for E. faecium isolates sampled from Romania, three groups were delineated mainly on origin source. Thus, group 1 joined isolates from human sources and environment (wastewater), whereas groups 2 and 3 comprised isolates originated mainly from birds. This pattern suggested the existence of at least two phylogenetic lines of E. faecium isolates in Romania: a line that had mainly the human host prevalence, and a second line with mainly animal host prevalence. The first phylogenetic line highlighted the existence of a clear connection between the isolates of human origin and those of wastewater origin. This connection was most likely determined by the uncontrolled spills in flowing water networks, which underlined a direct link between the hospital environment and the surrounding environment (wastewater in the present case). The results of our analyses pointed out that the degree of isolate differentiation revealed by MLST and eBURST was high enough to be used in epidemiological investigations. Authors’ ContributionsAuthors’ ContributionsAuthors’ ContributionsAuthors’ Contributions The contributions of authors to the manuscript are as follows: conceptualization: OP, IB; field work: IB, EJ; data curation: IB, DȘ, EJ; formal analysis: IB; funding acquisition: OP; investigation: IB, EJ; methodology: IB; project administration: OP; writing - original draft: IB; writing - review and editing: IB, DȘ, EJ, OP. All authors read and approved the final manuscript. Ethical approvalEthical approvalEthical approvalEthical approval (for researches involving animals or humans) Not applicable. AcknowledgementsAcknowledgementsAcknowledgementsAcknowledgements The authors are thankful to the staff of Someșeni Wastewater Treatment Laboratory from Cluj-Napoca and the Cantacuzino Institute from Bucharest for help in collecting and identifying enterococcal samples. This work was partly supported by CEEX 28/2005 from the National Authority for Scientific Research, Romania, as well as by the Core Project BIORESGREEN, subproject BioClimpact no. 7N/03.01.2023, code 23020401. Conflict of InterestsConflict of InterestsConflict of InterestsConflict of Interests The authors declare that there are no conflicts of interest related to this article. Băcilă I et al. (2023). Not Sci Biol 15(2):11546 14 ReferencesReferencesReferencesReferences Adegoke AA, Madu CE, Reddy P, Stenström TA, Okoh AI (2022). Prevalence of vancomycin resistant Enterococcus in wastewater treatment plants and their recipients for reuse using PCR and MALDI-ToF MS. Frontiers in Environmental Science 9:797992. https://doi.org/10.3389/fenvs.2021.797992 Băcilă I, Jakab E, Popescu O (2007). Molecular characterization of Enterococcus faecium strains isolated from sewage using the multilocus sequence typing (MLST) technique. Annals NSBMB 16(3):69-77. Blanch AR, Caplin JL, Iversen A, Kühn I, Manero A, Taylor HD, Vilanova X (2003). Comparison of enterococcal populations related to urban and hospital wastewater in various climatic and geographic European regions. Journal of Applied Microbiology 94(6):994-1002. https://doi.org/10.1046/j.1365-2672.2003.01919.x Correa-Martínez CL, Jurke A, Schmitz J, Schaumburg F, Kampmeier S, Mellmann A (2022). Molecular epidemiology of vancomycin-resistant enterococci bloodstream infections in Germany: A population-based prospective longitudinal study. Microorganisms 10(1):130. https://doi.org/10.3390/microorganisms10010130 Dutka-Malen S, Even S, Courvalin P (1995). Detection of glycopeptide resistance genotypes and identification to the species level of clinically relevant enterococci by PCR. Journal of Clinical Microbiology 35:24-27. https://doi.org/10.1128/jcm.33.1.24-27.1995 Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG (2004). eBURST: Infering patterns of evolutionary descends among clusters of related genotypes from multilocus sequence typing data. Journal of Bacteriology 186:1516-1530. https://doi.org/10.1128/JB.186.5.1518-1530.2004 Ferguson DM, Griffith JF, McGee CD, Weisberg SB, Hagedorn C (2013). Comparison of Enterococcus species diversity in marine water and wastewater using Enterolert and EPA Method 1600. Journal of Environmental and Public Health 2013:848049. https://doi.org/10.1155/2013/848049 Foulquié Moreno MR, Sarantinopoulos P, Tsakalidou E, de Vuyst L (2006). The role and application of enterococci in food and health. International Journal of Food Microbiology 106:1-24. https://doi.org/10.1016/j.ijfoodmicro.2005.06.026 Francisco AP, Bugalho M, Ramirez M, Carriço JA (2009). Global optimal eBURST analysis of multilocus typing data using a graphic matroid approach. BMC Bioinformatics 10:152. https://doi.org/10.1186/1471-2105-10-152 Franz CMAP, Huch M, Abriouel H, Holzapfel W, Gálvez A (2011). Enterococci as probiotics and their implications in food safety. International Journal of Food Microbiology 151:125-140. https://doi.org/10.1016/j.ijfoodmicro.2011.08.014 Giraffa G (2014). Enterococcus. In: Batt CA, Tortorello ML (Eds). Encyclopedia of Food Microbiology (Second Edition). Academic Press, Amsterdam pp 674-679. https://doi.org/10.1016/B978-0-12-384730-0.00098-7. Hall TA (1999). BioEdit: a user friendly biological sequence alignment editor and analysis program for Windows 95/95/NT. Nucleic Acids Symposium Series 41:95-98. Heuer H, Smalla K (2007). Horizontal gene transfer between bacteria. Environmental Biosafety Research 6(1-2):3-13. https://doi.org/10.1051/ebr:2007034 Homan WL, Tribe D, Poznanski S, Li M, Hogg G, Spalburg E, ... Willems RJL (2002). Multilocus sequence typing for Enterococcus faecium. Journal of Clinical Microbiology 40:1963-1971. https://doi.org/10.1128/JCM.40.6.1963- 1971.2002 Karimi F, Samarghandi MR, Shokoohi R, Godini K, Arabestani MR (2016). Prevalence and removal efficiency of enterococcal species and vancomycin-resistant enterococci of a hospital wastewater treatment plant. Avicenna Journal of Environmental Health Engineering 3(2):e8623. https://doi.org/10.5812/ajehe.8623 Kim SH, Cho SY, Kim HM, Huh K, Kang C-I, Peck KR, Chung DR (2021). Sequence type 17 is a predictor of subsequent bacteremia in vancomycin-resistant Enterococcus faecium-colonized patients: a retrospective cohort study. Antimicrobial Resistance and Infection Control 10(1):108. https://doi.org/10.1186/s13756-021-00980-1 Klein G (2003). Taxonomy, ecology and antibiotic resistance of enterococci from food and the gastro-intestinal tract. International Journal of Food Microbiology 88:123-131. https://doi.org/10.1016/S0168-1605(03)00175-2 Kühn I, Iversen A, Burman LG, Olsson-Liljequist B, Franklin A, Finn M, ... Möllby R (2003). Comparison of enterococcal populations in animals, humans, and the environment – a European study. International Journal of Food Microbiology 88:133-145. https://doi.org/10.1016/S0168-1605(03)00176-4 Băcilă I et al. (2023). Not Sci Biol 15(2):11546 15 Lebreton F, Willems RJL, Gilmore MS (2014). Enterococcus diversity, origins in nature, and gut colonization. In: Gilmore MS, Clewell DB, Ike Y, Shankar N (Eds). Enterococci: from commensals to leading causes of drug resistant infection. Massachusetts Eye and Ear Infirmary, Boston pp 5-63. Martínez-Carranza E, Barajas H, Alcaraz L-D, Servín-González L, Ponce-Soto G-Y, Soberón-Chávez G (2018). Variability of bacterial essential genes among closely related bacteria: the case of Escherichia coli. Frontiers in Microbiology 9:1059. https://doi.org/10.3389/fmicb.2018.01059 Moore DF, Guzman JA, McGee C (2008). Species distribution and antimicrobial resistance of enterococci isolated from surface and ocean water. Journal of Applied Microbiology 105(4):1017-1025. https://doi.org/10.1111/j.1365- 2672.2008.03828.x O’Driscoll T, Crank CW (2015). Vancomycin-resistant enterococcal infections: epidemiology, clinical manifestations, and optimal management. Infection and Drug Resistance 8:217-230. https://doi.org/10.2147/IDR.S54125 Oh PL, Benson AK, Peterson DA, Patil PB, Moriyama EN, Roos S, Walter J (2010). Diversification of the gut symbiont Lactobacillus reuteri as a result of host-driven evolution. The ISME Journal 4:377-387. https://doi.org/10.1038/ismej.2009.123 Pérez-Losada M, Browne EB, Madsen A, Wirth T, Viscidi RP, Crandall KA (2005). Population genetics of microbial pathogens estimated from multilocus sequence typing (MLST) data. Infection, Genetics and Evolution 6:97-112. https://doi.org/10.1016/j.meegid.2005.02.003 Ruoff KL, de la Maza L, Murtagh MJ, Spargo JD, Ferraro MJ (1990). Species identities of enterococci isolated from clinical specimens. Journal of Clinical Microbiology 28(3):435-437. https://doi.org/10.1128/jcm.28.3.435-437.1990 Shikov AE, Malovichko YV, Nizhnikov AA, Antonets KS (2022). Current methods for recombination detection in bacteria. International Journal of Molecular Sciences 23:6257. https://doi.org/10.3390/ijms23116257 Sievert DM, Ricks P, Edwards JR, Schneider A, Patel J, Srinivasan A, Fridkin S, National Healthcare Safety Network (NHSN) Team and Participating NHSN Facilities (2013). Antimicrobial-resistant pathogens associated with healthcare-associated infections: summary of data reported to the National Healthcare Safety Network at the Centers for Disease Control and Prevention, 2009-2010. Infection Control and Hospital Epidemiology 34(1):1- 14. https://doi.org/10.1086/668770 Sinton LW, Donnison AM (1994). Characterisation of faecal streptococci from some New Zealand effluents and receiving waters. New Zealand Journal of Marine and Freshwater Research 28(2):145-158. https://doi.org/10.1080/00288330.1994.9516603 Spratt BG, Hanage WP, Li B, Aanensen DM, Feil EJ (2004). Displaying the relatedness among isolates of bacterial species – the eBurst approach. FEMS Microbiology Letters 241:129-134. https://doi.org/10.1016/j.femsle.2004.11.015 Stefani S, Agodi A (2000). Molecular epidemiology of antibiotic resistance. International Journal of Antimicrobial Agents 13:143-153. https://doi.org/10.1016/S0924-8579(99)00103-X Stern CS, Carvalho M da G, Teixeira LM (1994). Characterization of enterococci isolated from human and nonhuman sources in Brazil. Diagnostic Microbiology and Infectious Disease 20(2):61-67. https://doi.org/10.1016/0732- 8893(94)90093-0 Tavanti A, Davidson AD, Johnson EM, Maiden MCJ, Shaw DJ, Gow NAR, Odds FC (2005). Multilocus sequence typing for differentiation of strain of Candida tropicalis. Journal of Clinical Microbiology 43:5993-5600. https://doi.org/10.1128/JCM.43.11.5593-5600.2005 Toc DA, Butiuc-Keul AL, Iordache D, Botan A, Mihaila RM, Costache CA, ... Junie LM (2022). Descriptive analysis of circulating antimicrobial resistance genes in vancomycin-resistant Enterococcus (VRE) during the COVID-19 Pandemic. Biomedicines 10:1122. https://doi.org/10.3390/biomedicines10051122 Willems RJL, Top J, van Santen M, Robinson DA, Coque TM, Baquero F, .... Bonten MJM (2005). Global spread of vancomycin-resistant Enterococcus faecium from distinct nosocomial genetic complex. Emerging Infectious Diseases 11(6):821-8. https://doi.org/10.3201/eid1106.041204 Zhong Z, Zhang W, Song Y, Liu W, Xu H, Xi X, ... Sun Z (2017). Comparative genomic analysis of the genus Enterococcus. Microbiological Research 196:95-105. https://doi.org/10.1016/j.micres.2016.12.009 Băcilă I et al. (2023). Not Sci Biol 15(2):11546 16 The journal offers free, immediate, and unrestricted access to peer-reviewed research and scholarly work. Users are allowed to read, download, copy, distribute, print, search, or link to the full texts of the articles, or use them for any other lawful purpose, without asking prior permission from the publisher or the author. License License License License ---- Articles published in Notulae Scientia BiologicaeNotulae Scientia BiologicaeNotulae Scientia BiologicaeNotulae Scientia Biologicae are Open-Access, distributed under the terms and conditions of the Creative Commons Attribution (CC BY 4.0) License. © Articles by the authors; Licensee SMTCT, Cluj-Napoca, Romania. The journal allows the author(s) to hold the copyright/to retain publishing rights without restriction. Notes:Notes:Notes:Notes:  Material disclaimer: The authors are fully responsible for their work and they hold sole responsibility for the articles published in the journal.  Maps and affiliations: The publisher stay neutral with regard to jurisdictional claims in published maps and institutional affiliations.  Responsibilities: The editors, editorial board and publisher do not assume any responsibility for the article’s contents and for the authors’ views expressed in their contributions. The statements and opinions published represent the views of the authors or persons to whom they are credited. Publication of research information does not constitute a recommendation or endorsement of products involved.