Next Article in Journal
Geomats Used to Control Erosion on Reservoir Margins in Brazilian Hydroelectric Power Plants
Next Article in Special Issue
Aspergillus niger Culture Filtrate (ACF) Mediated Biocontrol of Enteric Pathogens in Wastewater
Previous Article in Journal
Study on the Horizontal Distribution Law of Flood Water and Sediment Factors under the Effect of Vegetation on a Curved Beach
Previous Article in Special Issue
Groundwater Nitrate Removal Performance of Selected Pseudomonas Strains Carrying nosZ Gene in Aerobic Granular Sequential Batch Reactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Persistence of Enterobacteriaceae Drawn into a Marine Saltern (Saline di Tarquinia, Italy) from the Adjacent Coastal Zone

1
Dipartimento di Ecologia e Biologia, Università degli Studi della Tuscia, Largo Università snc, 01100 Viterbo, Italy
2
Laboratoro di Ecologia dei Funghi Marini CONISMA, Università degli Studi della Tuscia, Largo Università snc, 01100 Viterbo, Italy
3
Dipartimento di Scienze dell’Ambiente e della Terra, Università di Milano-Bicocca, Piazza della Scienza 1, 20126 Milano, Italy
4
Institute of Water Research, University of Granada, C/Ramón y Cajal, 4, 18071 Granada, Spain
5
Faculty of Pharmacy, Campus de Cartuja, University of Granada, s/n, 18071 Granada, Spain
6
Laboratorio di Microbiologia Marina Applicata, CONISMA, Università degli Studi della Tuscia, Largo Università snc, 01100 Viterbo, Italy
*
Author to whom correspondence should be addressed.
Water 2021, 13(11), 1443; https://doi.org/10.3390/w13111443
Submission received: 30 April 2021 / Revised: 14 May 2021 / Accepted: 18 May 2021 / Published: 21 May 2021
(This article belongs to the Special Issue Biotechnology Applications in Water and Wastewater Treatment)

Abstract

:
Enterobacteriaceae is present in various niches worldwide (i.e., the gastrointestinal tracts of animals, clinical specimens, and diverse environments) and hosts some well-known pathogens (i.e., salmonellas, shigellas and pathogenic coliforms). No investigation has focused on its occurrence in marine salterns, and it is not clear if these hypersaline environments could be a reservoir for these bacteria including some potentially harmful members. In this study, a two-year metabarcoding survey was carried out on samples collected from different ponds of the “Saline di Tarquinia” salterns and the nearby coastal waters. Enterobacteriaceae was recorded almost constantly in the seawaters feeding the saltern. Its abundance was generally higher in the sea than in the ponds, probably due to the higher anthropic impact. The same trend was evidenced for the key genus (Escherichia/Shigella) and OTU (OTU 5) of the Enterobacteriaceae community. Various parameters affected taxon/OTU abundance: Enterobacteriaceae, Escherichia/Shigella and OTU5 decreased with increasing salinity and rains; moreover, Escherichia/Shigella and OTU 5 were higher in autumn than in spring. Although Enterobacteriaceae did not seem to find the most favourable conditions for a high-abundance persistence in the saltern environment, it did not disappear. These observations suggested this environment as a potential reservoir for bacteria with possible important health implications.

1. Introduction

Enterobacteriaceae (phylum Proteobacteria) is a complex family under the taxonomic point of view that includes 33 genera and 134 species (List of Prokaryotic Names with Standing in Nomenclature, https://lpsn.dsmz.de/family/enterobacteriaceae-1 (accessed on 17 May 2021), [1]). It hosts members, distributed worldwide, usually isolated from various niches, such as the gastrointestinal tracts of animals (from insects to humans), clinical specimens, and diverse environments [2,3].
The gastrointestinal tract of animals is one of their most common habitats [4], representing the environment that is historically associated with this family.
Most genera of Enterobacteriaceae have an intimate association as commensals with animals (including humans), being part of the gut microbiota. However, various pathogenic variants, associated with infections, diseases, and syndromes, have been recognised [5]. In this context, the prototypical example is that of Escherichia coli. The vast majority of the E. coli population in humans has a commensal role and is adapted to colonise the host without causing disease. However, some pathogenic strains can colonise the host, often resulting in clinically significant pathologies [6]. The current classification of pathogenic E. coli includes enteroaggregative, enteropathogenic, enteroinvasive, enterotoxigenic, diffusely adherent and Shiga toxin-producing strains [7].
Species of Enterobacteriaceae are responsible for numerous animal infections, also with economic implications. For instance, Klebsiella pneumoniae, E. coli, and Citrobacter freundii, which cause bovine mastitis, reduce animal health and may result in considerable economic losses to the dairy industry [8,9]. E. coli, Citrobacter spp., Enterobacter spp., Klebsiella spp., Kluyvera spp. and Salmonella spp. are involved in farmed fish infections/outbreaks and are responsible for monetary losses in aquaculture [10,11,12].
Many Enterobacteriaceae are pathogens or opportunistic pathogens involved in several human diseases, from foodborne and waterborne diseases to nosocomial infections (i.e., septicaemia, meningitis, respiratory and urinary infections) [5,13,14,15]. Salmonella spp. and enterohaemorrhagic E. coli are among the most common foodborne pathogens, affecting millions of people annually [16]. In particular, salmonella represents the second most common causative agent of food-transmitted diseases in Europe and the United States [17,18], and the third most common cause of death among foodborne diseases worldwide [19]. Moreover, various species of Enterobacteriaceae (such as Escherichia spp., Salmonella spp., Shigella spp., Klebsiella spp., and Enterobacter spp.) are highly significant causes of intestinal and extraintestinal nosocomial infections [4]. Furthermore, some of them represent a serious public health threat due to their multidrug resistance [20,21]. Cephalosporin- and carbapenem-resistant Enterobacteriaceae was ranked in the category of critical priority pathogens by the WHO [22]. Infections caused by carbapenem-resistant Enterobacteriaceae are associated with high morbidity and mortality of hospitalised patients [23].
In addition to the clinical/health implications, another important issue is the detection of Enterobacteriaceae in an increasing number of environmental niches [5]. Members of this family have been found in soil [24,25,26,27,28], aquatic habitats [24,29,30,31], and associated with vegetation (fruits, vegetables, grains, flowering plants, and trees) [29,32,33,34,35,36,37,38]. For instance, various species (belonging to different genera, such as Enterobacter, Gibbsiella, Klebsiella, Kosakonia, Mangrovibacter, Phytobacter and Pluralibacter) have been reported in relationships with plants, as endophytes (involved in plant growth promotion, nutrient uptake, nodulation, defence against phytopathogens and increase in stress tolerance) [38,39,40,41,42,43,44] or phytopathogens [45,46,47].
Soil, water, and vegetation have proven to be environmental reservoirs for Enterobacteriaceae. In this context, some work has investigated their relationship with the environmental habitats and/or possible implications for public health [48].
Although Enterobacteriaceae members are found in an increasing number of environmental niches, their occurrence in some habitats is still poorly surveyed, and some environments are underestimated as potential reservoirs of these bacteria. In particular, their occurrence in marine salterns and possible variations according to the relative environmental stress conditions (such as intense salinity and temperature fluctuations) seem not to have been investigated. Moreover, it was not evaluated if these environments could represent a possible reservoir for this group, hosting bacteria with important health implications.
The “Saline di Tarquinia” (ST) are marine salterns, located ca. 80 km northwest of Rome (Italy), along a low sandy coast of the North Tyrrhenian Sea. ST, converted into a Nature Reserve in 1980, ceased salt production in 1997. This led to a dramatic reduction in anthropic activities, and also affected the management and maintenance of the pond water balance over time. However, ST is an extreme environment, featuring large and sudden variations of environmental parameters (in particular, salinity and water temperature) in the majority of the ponds. The ST salterns are located downstream from a large floodplain, bounded to the northeast by the hilly area of the city of Tarquinia and to the northwest and southeast by the mouths of the rivers Marta and Mignone, respectively. They cover an area of about 135 ha and consist of ca. 100 shallow interconnected ponds that allow the establishment of an increasing salinity gradient, starting from the water pumped from the nearby sea which is transferred through low-salinity ponds (concentrators), to ponds where salt precipitation occurs (crystallisers) [49]. Some studies have been carried out to investigate the microbial diversity of the site [50,51,52,53,54,55,56,57]. The site algal microflora was studied both from the taxonomic point of view [50,56] and for some aspects regarding their pigment production, which is particularly interesting at the biotechnological level [52,55]. Moreover, ST bacteria were studied by culture-dependent and culture-independent methods, showing a very diversified community with members adapted to broad salinity and temperature variations [51,53,54]. The ST bacterial diversity showed increasing simplification along the gradient according to salinity increases. Moreover, the communities were structured by various factors: salinity, sampling site, and sampling year and month [57] However, investigations related to some particular bacterial groups are quite scarce. The only study considering a specific bacterial group was focused on inspecting the presence of Vibrio within the ST system, which was demonstrated to be a possible reservoir for members of this genus even if at very low levels of abundance [49].
The aim of the current study was to investigate the presence and distribution of Enterobacteriaceae within the ST marine salterns, in order to understand if this hypersaline environment can represent a possible reservoir of potentially harmful bacteria. The work (carried out by a metabarcoding approach) surveyed the Enterobacteriaceae occurrence over a broad time span (2 years), analysing samples taken from the coastal waters (representing the saltpan seawater inputs) and various ponds along the salinity gradient. In addition, the influence of environmental parameters on Enterobacteriaceae distribution within the ST system was also studied.

2. Materials and Methods

2.1. Sample Collection and Characterisation

Samplings within the ST (North Tyrrhenian Sea, Italy; 42°12′07.8″ N 11°43′17.8″ E) were carried out monthly over two years (from May 2012 to April 2014). Water samples were collected from the sea (S) and three different ponds with increasing salinities (P5, low-salinity concentrator pond; P24, intermediate-salinity concentrator pond, and P37, crystallisation pond), selected according to their annual range of salinity in order to cover the whole ST gradient (Figure 1).
Sample collection, measurement, and the determination of environmental parameters were carried out as reported by Gorrasi et al. [57].
All samples, obtained by pooling three different sub-samples [58], were stored in sterile bottles and kept refrigerated (4 °C) during transport to the ST laboratory to be processed within one hour. Then, 1 L of water was vacuum-filtered on sterile membranes (0.22 µm, Millipore, Burlington, MA, USA), which were washed twice with sterile saline (NaCl) solutions (with concentrations roughly corresponding to the salinity recorded on the sites) to remove possible nutrients from the filters [59]. The filters were kept frozen until DNA extraction. The remaining sample aliquots were used to determine BOD5 and chlorophyll pigment concentration, which were measured by standard methods [60].
Environmental parameters (salinity, pH, conductivity, and water temperature) were recorded during samplings by common handheld probes. Daily rainfall data were obtained from a local government institution [61]. All parameters are reported in Table S1.

2.2. DNA Extraction, 16S rDNA Amplicon Libraries and Sequencing

Total DNA was extracted from filters using a GeneMATRIX Bacterial & Yeast Genomic DNA Purification Kit (EURx Sp. z o.o., Gdansk, Poland) as previously reported [58].
The multiplexed amplicon libraries were prepared using a dual PCR amplification protocol, amplifying the V5–V6 hypervariable regions of the 16S rDNA. The primers used, 783F (5′-CAGGATTAGATACCC) and 1046R (5′-CGACRRCCATGCANCACCT) [62,63], were modified by adding external barcodes to allow the parallel processing of multiple samples. The strategy and related protocols used to prepare the libraries have already been described in detail by Gorrasi et al. [57].
The second step of library preparation, with the addition of standard Nextera indexes (Illumina, San Diego, CA, USA), and sequencing were carried out at Nuova Genetica Italiana SRL (Monza-Brianza, Italy). Amplicon libraries were sequenced by Illumina MiSeq (Illumina, San Diego, CA, USA) using a 2 × 250 bp paired-end protocol.

2.3. Sequence Processing and Data Analysis

Sequence processing and data analyses were carried out as reported by Gorrasi et al. [57].
Reads from the sequencing were demultiplexed according to the indices and internal barcodes.
Reads were processed using the UPARSE pipeline [64]. Forward and reverse reads were merged with perfect overlapping and quality, filtered with default parameters. Suspected chimeras and singleton sequences (i.e., sequences appearing only once in the whole dataset) were removed. OTUs were defined on the whole dataset clustering the sequences at 97% similarity and defining a representative sequence for each cluster. The abundance of each OTU was estimated by mapping the sequences of each sample against the representative sequence of each OTU at 97% of similarity.
Taxonomic classification of the OTU representative sequences was obtained by RDP classifier [65], using a 50% confidence cut-off as suggested for short sequences [66].

2.4. Statistical Methods

Principal component analysis (PCA) was performed to summarise sample distribution according to Enterobacteriaceae genera.
Redundancy analysis (RDA) was performed to investigate the influence of environmental parameters on the Enterobacteriaceae community. OTU abundance data for the RDA analysis were square-root-transformed. The environmental parameters taken into account were: sampling site, salinity, water temperature, pH, chlorophyll, BOD5, rains (cumulative data relative to the seven days preceding the sampling) and sampling month (entered as Fourier series transformed data, sin(Month) and cos(Month), to account for seasonality) [67]. Sampling site, rains and sin(Month) were used as predictors for the RDA analysis. The other variables were removed, being collinear (Pearson |r| > 0.6) with the predictors or not significant according to the forward selection. The forward selection was performed for the stepwise selection of the explanatory variables, correcting the significance according to the false discovery rate (FDR) in order to avoid Type I error inflation [68]. The RDA significance was assessed by the Monte Carlo permutation test, using 9999 permutations.
The abundance variation of Enterobacteriaceae, the most abundant genera, and OTUs in relation to the key environmental parameters was investigated by generalised linear models (GLMs), assuming a Poisson distribution and with a log link function. Additionally, in the GLM analysis, the periodic regression based on Fourier series transformation of months was used to account for taxa seasonal variation. The Pearson matrix was calculated using the statistical software Systat 8.0 (Systat Software Inc., Point Richmond, CA, USA), whereas the PCA, RDA and GLM analyses were performed using the CANOCO v. 5.1 software package (Microcomputer Power, Ithaca, NY, USA).
Post hoc tests were used to assess pairwise differences in taxa abundance variation among sampling sites. Non-parametric tests (Kruskal–Wallis, one-way analysis of variance, and Kolmogorov–Smirnov two-sample tests) were performed, based on data assessments obtained by Shapiro–Wilk and Levene tests (carried out to test data normal distribution and homogeneity of variance assumptions, respectively). All these tests were run using Systat 8.0.
Spearman’s rank correlation coefficient (ρ) was calculated to determine associations among the most abundant OTUs (Systat 8.0).

3. Results and Discussion

3.1. Spatio-Temporal Variation of Enterobacteriaceae

Across the 96 samples, a total of 2,566,189 bacterial sequences were obtained after quality control. Among them, 45,778 (resulting in 26 different OTUs, by sequence clustering at 97% of similarity) were assigned to Enterobacteriaceae using a 50% confidence cut-off.
Over the two-year survey, Enterobacteriaceae was always detected in the sea, except in April and June 2013. Moreover, in 13 months out of 24 it was a major taxon, with relative abundance (Ra) in the range 2.5–18.1% (Figure 2; Table S2).
Abundance of this taxon was generally higher in the sea (S) than in the various ponds, and it decreased along the salinity gradient (Table S2). The exceptions were May 2012, September 2012, December 2012, November 2013, and March 2014, when it was higher in P5 than in S. In particular, in May 2012 in the P5 and in S, it was 19.7 and 7.2%, respectively.
A high presence of Enterobacteriaceae in coastal waters, in particular if close to anthropized areas, could be due to various causes, possibly related to environmental contamination, such as the discharge of sewage effluents, recreational activities, and aquaculture.
As for ST, they border a rather high-impact touristic settlement that, based on seasonal issues, may contribute to the increased presence of enteric bacteria. In addition, the mouths of the two above-mentioned rivers bordering the floodplain release their freshwaters close to the site, possibly contributing with further inputs of bacteria from the upstream urbanised areas, where there are also various agricultural and animal farming activities. Sea hydrodynamics (currents and tides) can facilitate the diffusion of these bacteria toward the saltern, and their loading into the pond system by seawater pumping. These are probably the reasons for the generally higher levels of this bacterial group in our sea samples. Moreover, the higher Enterobacteriaceae occurrence in the sea than in the ponds could also be due to more suitable (less stressing) conditions characterising the marine coastal environment.
In general, it was evident that Enterobacteriaceae driven into the ST system by the seawater influx did not seem to find the most favourable conditions for its persistence in the hypersaline environment. However, despite this, these bacteria generally did not disappear, as shown by monthly data. The presence of Enterobacteriaceae in the ponds could also be due to the contribution of the avifauna, which in some periods could be quite abundant.
To obtain an overview of Enterobacteriaceae occurrence in the ST system and in the neighbouring sea, the accumulated data were analysed (Figure 3A). Over the two years, Enterobacteriaceae showed mean abundances of 4.12, 2.88, 0.41 and 0.60% in S, P5, P24 and P37, respectively. However, in some periods, its abundances were notably above the average: in July 2012 in S (Ra = 18%) and in May 2012 in P5 (Ra = 19.7%).
Considering the global data, Enterobacteriaceae abundance was higher in S and P5 (low-salinity concentrator pond, showing salinities similar to those of the sea) than in P24 (intermediate-salinity concentrator pond) and P37 (crystallisation pond) (p < 0.01).
Among the 26 Enterobacteriaceae OTUs, 17 were assigned at the genus level. The genera detected were Escherichia/Shigella, Siccibacter, Buttiauxella, Atlantibacter, Yokenella, Pseudescherichia and Metakosakonia.
The genera abundances in all samples are reported in Table S3.
Metakosakonia were definitely underrepresented, being found only in two ST samples (May 2012—S and March 2014—P24) out of 96, and with very low abundances (Ra~0.01%).
Atlantibacter, Yokenella and Pseudescherichia, when detected, always represented rare taxa (with Ra values in the range 0.001–0.06%, 0.001–0.3% and 0.002–0.01%, respectively).
Over the two-year survey, Buttiauxella presence was revealed only in May 2012; it was found in all sites as a rare taxon (0.005–0.08%), except in P5 where it showed a Ra = 2.1%.
Siccibacter was found in 26 samples, but it showed Ra ≥ 1% (1.2 and 12% in May 2012—S and May 2012—P5, respectively) in only two samples.
Escherichia and Shigella species have essentially identical 16S rRNA sequences. Thus, the reference databases for taxonomic identification (including RDP) use the combined genus name Escherichia/Shigella [69]. Therefore, it is not possible to uniquely attribute the relative sequences to one of these two genera.
Escherichia/Shigella was the most represented taxon among the Enterobacteriaceae (82% of the Enterobacteriaceae sequences across all samples). It was found in 75 samples and was recorded with Ra ≥ 1% in 30 samples. Its mean abundances in S, P5, P24 and P37 were 3.45, 1.90, 0.34 and 0.46%, respectively. However, it showed above-average abundances in S in July 2013 (15.6 %) and September 2013 (10.7%) and in P24 in May 2012 (Ra = 1.2%) and August 2013 (Ra = 2%).
Considering the overall two-year data, as already observed for the whole Enterobacteriaceae, Escherichia/Shigella was more abundant in S and P5 (where it also showed much broader abundance variations) than in P24 and P37 (p < 0.01) (Figure 3B).
In addition, the PCA of sampling sites, according to the abundance of the Enterobacteriaceae genera, revealed Escherichia/Shigella as the most indicative taxon for sample variability. The PCA explained 99.58% of the total variance along the first two axes (82.56% explained by the first axis) (Figure 4).
From a general aspect, it is important to note that members of these two genera are mainly animal-adapted (including humans) enteric bacteria and, among them, some members are known human pathogens [7,70]. It is possible to speculate that their frequent detection within the ST, and in particular, with higher abundance in the neighbouring sea, could be mainly due to contamination of anthropogenic and/or animal origin in the coastal area [71].
Among the 26 Enterobacteriaceae OTUs, two were found as the most abundant (having a total number of mapping sequences on the whole ST dataset > 1000): OTU 5 (37,725 sequences) and OTU 206 (6714 sequences).
Most (82%) of the Enterobacteriaceae sequences found across all samples were ascribed to OTU 5, which was the most abundant and represented OTU (Table S4). It was assigned to Escherichia/Shigella and represented the main genus fraction (97.8–100% of this taxon across all samples). Moreover, in 73 samples it represented a consistent proportion (50–100%) of Enterobacteriaceae. OTU 5 showed mean abundances of 3.45, 1.90, 0.34 and 0.46% in S, P5, P24 and P37, respectively. It was recorded with above-average abundances in S in July 2012 (Ra = 15.6%) and September 2013 (10.7%), and in P24 in May 2012 (Ra = 1.2%) and August 2013 (Ra = 2%).
The cumulative data evidenced that OTU 5 was more abundant in S and P5 than in P24 and P37 (p < 0.01) (Figure 5A). Overall, OTU 5 abundance variation among sampling sites essentially coincided with those of Enterobacteriaceae and Escherichia/Shigella because it was the predominant fraction of both taxa.
OTU 206 (belonging to Siccibacter) was detected in 69 samples, but in 63 it represented a rare taxon. It was detected with Ra ≥ 1% (1.1–2.5%) in 6 samples: four collected in the sea (May 2012—S, July 2012—S, September 2013—S, and October 2013—S), and two in P5 (May 2012—P5 and November 2012—P5). Considering the overall data, its highest mean abundance (~0.6%) was recorded in the sea, where in July 2012 the OTU reached 2.5% (Table S4). Overall, OTU 206 was significantly (p < 0.01) more abundant in the sea than in P24 and P37 (Figure 5B).

3.2. Environmental Parameters Driving Enterobacteriaceae Occurrence (RDA) and Variation of Key Taxa and Most Abundant OTUs in Relation to Them

The RDA was performed to investigate the key parameters affecting Enterobacteriaceae occurrence and distributions and to explore the relationships between the most abundant OTUs and these parameters. Among the various environmental parameters taken into account, water temperature, pH, and conductivity were excluded from the RDA, due to these factors being collinear with other parameters (Pearson |r| > 0.6). The tested parameters were: sampling site, salinity, chlorophylls, BOD5, rains and sampling month (entered as sin(Month) and cos(Month) to account for seasonality). The RDA analysis showed that the Enterobacteriaceae community changed significantly with sampling site, rains and sin(Month) (Table 1).
Among the three factors, the sampling site had the main effect on the Enterobacteriaceae community, explaining the highest percentage of total variance (17.79%, sum of variance percentages explained by single sampling site; Table 1).
As shown by the RDA triplot (Figure 6), a clear separation between samples collected from S and P5 and those collected from P24 and P37 was observed. A partial overlap of S and P5 sample clusters was evidenced, while a large overlap was observed between those of P24 and P37. The sample ordering was based on the OTU composition, which means that the compositions of S and P5 were rather similar, and those of P24 and P37 were very similar. Therefore, the observed overlapping of clusters was mainly due to OTU 5 contribution (represented the preponderant fraction of the Enterobacteriaceae), and secondly to that of OTU 6, which showed similar abundances in S/P5 and P24/P37 (box and whisker plots, Figure 3).
The RDA showed that the P24 and P37 samples collected in spring and winter were those mostly affected by rains. Actually, rainfalls had a great impact on P24 (intermediate-salinity concentrator pond) and P37 (crystallisation pond): during the rainiest periods, salinity dramatically dropped in these two ponds. This was particularly evident in January–April 2013 and February 2014, when the recorded salinity values were close to those of the sea or even lower (Table S1). As for the abundance variation of Enterobacteriaceae, Escherichia/Shigella, and the most abundant OTUs in relation to the rains, they changed significantly (p < 0.02) according to this parameter, decreasing with rainfall increase (Figure 7A). The decrease in these bacteria in the ponds during rainy periods was probably due to the strong brine dilution caused by the precipitation. In fact, both concentrations and the crystallisation pond were very shallow (~15 cm), and their water volume was strongly affected by rains.
In addition, samples were ordered along the sin(Month) vector according to the season, with those collected in spring clearly projected towards the vector tip, whereas those collected in autumn were on the opposite side. Significant abundance variations in relation to this parameter were observed only for OTU 5, and more generally for the whole Escherichia/Shigella taxon, being higher in autumn than in spring (p = 0.048) (Figure 7B).
As for the most abundant OTUs, the RDA showed that they were projected towards S and P5 and on the opposite side of the rain tip vector. In fact, as discussed above, they showed the highest abundances in S and P5 and decreased with rainfall increase. In the RDA triplot, the two vectors representing OTU 5 and OTU 206, besides having the same direction, were very close to each other. Actually, these OTUs showed a high positive correlation (as indicated by the Spearman’s rank correlation coefficient, ρ = 0.962), suggesting their high association and thus their co-occurrence.
As mentioned, salinity was not significant for the RDA model, based on the contribution of each single OTU. However, the GLM analyses were performed on the whole family, its key genus, and OTUs to understand whether the abundance variation in relation to this parameter could still be significant. For Enterobacteriaceae, Escherichia/Shigella and OTU 5 significant abundance variations were evidenced (Figure 7C). They showed, as expected, the same trends with an abundance decrease at increasing salinity. This was congruent with the above results, which showed that their abundance decreased from the near the sea to the ST ponds and along the salinity gradient.

4. Conclusions

To the best of our knowledge, this was the first study investigating the presence and distribution of Enterobacteriaceae (based on monthly data recorded over two years) within marine salterns, in order to understand if these hypersaline environments can represent possible reservoirs of potentially harmful bacteria.
In the studied system, Enterobacteriaceae was almost always detected in the waters sampled from the nearby sea (the feeding point for the ST system). Its abundance and those of its key genus and OTU (Escherichia/Shigella and OTU5, respectively) were higher in the sea than in the ponds, probably due to a combination of factors: the high impact of human activities in the coastal area and the more suitable (less stressing) environmental conditions characterising the sea.
It was evident that Enterobacteriaceae loaded into the pond system by seawater pumping did not seem to find the most favourable conditions for its persistence at high levels of abundance in the hypersaline environment. However, despite this, it generally did not disappear, suggesting that this hypersaline environment could represent a possible reservoir for this group, hosting bacteria with important health implications.
This work also showed the importance of some parameters in determining taxon/OTU abundance. Enterobacteriaceae, Escherichia/Shigella and OTU5 decreased with increasing salinity and rain. Moreover, a seasonal pattern was evidenced for Escherichia/Shigella and OTU 5, which were higher in autumn than in spring.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/w13111443/s1, Table S1: Environmental parameters recorded for the ST samples from May 2012 to April 2014, Table S2: Relative abundances of Enterobacteriaceae across the ST samples, Table S3: Relative abundances of Enterobacteriaceae genera across the ST samples, Table S4: Relative abundances of Enterobacteriaceae OTUs across the ST samples.

Author Contributions

Conceptualisation, S.G., A.F. and M.F.; methodology, S.G., A.F., M.P. and M.F.; software, A.F.; validation, S.G., A.F., M.P., A.G.-M., J.G.-L., B.M.-P. and M.F.; formal analysis, S.G., A.F., M.P. and M.F.; investigation, S.G. and M.F.; resources, A.F., M.P. and M.F.; data curation, S.G. and A.F.; writing—original draft preparation, S.G., M.P., A.F., A.G.-M., J.G.-L., B.M.-P. and M.F.; writing—review and editing, S.G., M.P., A.F., A.G.-M., J.G.-L., B.M.-P. and M.F.; visualisation, S.G., A.F., M.P. and M.F.; supervision, A.F. and M.F.; project administration, M.F.; funding acquisition, A.F. and M.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequence data are available at the European Nucleotide Archive (ENA), study accession number PRJEB38856 (http://www.ebi.ac.uk/ena/data/view/PRJEB38856, accessed on 12 January 2021).

Acknowledgments

The authors wish to thank the staff of ‘Comando Unità per la Tutela Forestale Ambientale e Agroalimentare Carabinieri, Posto fisso di Tarquinia UTB Carabinieri’, Riserva Naturale Statale Saline di Tarquinia, Località Saline, 01016 Tarquinia-Viterbo, for their kind support during the sampling.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. LPSN-List of Prokaryotic names with Standing in Nomenclature, Family Enterobacteriacveae. Available online: https://lpsn.dsmz.de/family/enterobacteriaceae-1 (accessed on 3 March 2021).
  2. Brenner, D.J.; Farmer, J.J.  Enterobacteriaceae. In Bergey’s Manual of Systematics of Archaea and Bacteria; Whitman, W.E.B., Ed.; John Wiley & Sons: New York, NY, USA, 2015; pp. 1–24. [Google Scholar]
  3. Xu, Z.; Xia, M.; Huo, Y.X.; Yang, Y. Intestinirhabdus alba gen. nov., sp. nov., a novel genus of the family Enterobacteriaceae, isolated from the gut of plastic-eating larvae of the Coleoptera insect Zophobas atratus. Int. J. Syst. Evol. Microbiol. 2020, 70, 4951–4959. [Google Scholar] [CrossRef] [PubMed]
  4. Farmer, J.J.; Farmer, M.K.; Holmes, B. The Enterobacteriaceae: General characteristics. In Topley & Wilson’s Microbiology and Microbial Infections, 10th ed.; Borriello, S.P., Murray, P.R., Funke, G., Eds.; John Wiley & Sons: Chichester, UK, 2010; Volume 2, pp. 1317–1359. [Google Scholar]
  5. Janda, J.M.; Abbott, S.L. The Changing Face of the Family Enterobacteriaceae (Order: “Enterobacterales”): New Members, Taxonomic Issues, Geographic Expansion, and New Diseases and Disease Syndromes. Clin. Microbiol. Rev. 2021, 34, e00174-20. [Google Scholar] [CrossRef] [PubMed]
  6. Chattopadhyay, S.; Sokurenko, E.V. Evolution of pathogenic Escherichia coli. In Escherichia coli, 2nd ed.; Donnenberg, M.S., Ed.; Academic Press: Cambridge, MA, USA, 2013; pp. 45–71. [Google Scholar]
  7. Scheutz, F.; Strockbine, N.A.  Escherichia. In Bergey’s Manual of Systematics of Archaea and Bacteria; Whitman, W.E.B., Ed.; John Wiley & Sons: New York, NY, USA, 2015; pp. 1–49. [Google Scholar]
  8. Gao, J.; Barkema, H.W.; Zhang, L.; Liu, G.; Deng, Z.; Cai, L.; Shan, R.; Zhang, S.; Zou, J.; Kastelic, J.P.; et al. Incidence of clinical mastitis and distribution of pathogens on large Chinese dairy farms. J. Dairy Sci. 2017, 100, 4797–4806. [Google Scholar] [CrossRef] [PubMed]
  9. Nonnemann, B.; Lyhs, U.; Svennesen, L.; Kristensen, K.A.; Klaas, I.C.; Pedersen, K. Bovine mastitis bacteria resolved by MALDI-TOF mass spectrometry. J. Dairy Sci. 2019, 102, 2515–2524. [Google Scholar] [CrossRef] [PubMed]
  10. Elsherief, M.F.; Mousa, M.M.; El-Galil, H.A.; El-Bahy, E.F. Enterobacteriaceae associated with farm fish and retailed ones. Alex. J. Vet. Sci. 2014, 42, 99–104. [Google Scholar] [CrossRef] [Green Version]
  11. Das, A.; Behera, B.K.; Acharya, S.; Paria, P.; Chakraborty, H.J.; Parida, P.K.; Das, B.K. Genetic diversity and multiple antibiotic resistance index study of bacterial pathogen, Klebsiella pneumoniae strains isolated from diseased Indian major carps. Folia Microbiol. 2019, 64, 875–887. [Google Scholar] [CrossRef]
  12. Preena, P.G.; Dharmaratnam, A.; Raj, N.S.; Raja, S.A.; Nair, R.R.; Swaminathan, T.R. Antibiotic-resistant Enterobacteriaceae from diseased freshwater goldfish. Arch. Microbiol. 2021, 203, 219–231. [Google Scholar] [CrossRef]
  13. Fleckenstein, J.M.; Bartels, S.R.; Drevets, P.D.; Bronze, M.S.; Drevets, D.A. Infectious agents of food-and water-borne illnesses. Am. J. Med. Sci. 2010, 340, 238–246. [Google Scholar] [CrossRef]
  14. European Food Safety Authority and European Centre for Disease Prevention and Control (EFSA and ECDC). The European Union summary report on trends and sources of zoonoses, zoonotic agents and food-borne outbreaks in 2017. EFSA J. 2018, 16, e05500. [Google Scholar]
  15. Tack, D.M.; Ray, L.; Griffin, P.M.; Cieslak, P.R.; Dunn, J.; Rissman, T.; Jervis, R.; Lathrop, S.; Muse, A.; Duwell, M.; et al. Preliminary incidence and trends of infections with pathogens transmitted commonly through food-Foodborne Diseases Active Surveillance Network, 10 US Sites, 2016–2019. MMWR Morb. Mortal. Wkly. Rep. 2020, 69, 509–514. [Google Scholar] [CrossRef]
  16. World Health Organization. Available online: https://www.who.int/news-room/fact-sheets/detail/food-safety (accessed on 22 April 2021).
  17. European Food Safety Authority and European Centre for Disease Prevention and Control. The European Union One Health 2019 Zoonoses Report. EFSA J. 2021, 19, e06406. [Google Scholar]
  18. Lee, H.; Yoon, Y. Etiological Agents Implicated in Foodborne Illness World Wide. Food Sci. Anim. Resour. 2021, 41, 1–7. [Google Scholar] [CrossRef] [PubMed]
  19. World Health Organization. WHO Estimates of the Global Burden of Foodborne Diseases: Executive Summary (No. WHO/FOS/15.02.); WHO Press: Geneva, Switzerland, 2015; pp. 1–2. [Google Scholar]
  20. Bassetti, M.; Peghin, M.; Pecori, D. The management of multidrug-resistant Enterobacteriaceae. Curr. Opin. Infect. Dis. 2016, 29, 583–594. [Google Scholar] [CrossRef] [PubMed]
  21. van Duin, D.; Doi, Y. The global epidemiology of carbapenemase-producing Enterobacteriaceae. Virulence 2017, 8, 460–469. [Google Scholar] [CrossRef] [PubMed]
  22. World Health Organization. Global Priority List of Antibiotic-Resistant Bacteria to Guide Research, Discovery, and Development of New Antibiotics; WHO Press: Geneva, Switzerland, 2017; pp. 1–7. [Google Scholar]
  23. Zhu, C.; Liyanapathirana, V.; Li, C.; Pinto, V.; Hui, M.; Lo, N.; Wong, K.T.; Dissanayake, N.; Ip, M. Characterizing mobilized virulence factors and multidrug resistance genes in carbapenemase-producing Klebsiella pneumoniae in a Sri Lankan Hospital. Front. Microbiol. 2018, 9, 2044. [Google Scholar] [CrossRef] [Green Version]
  24. Izard, D.; Ferragut, C.; Gavini, F.; Kersters, K.; De Ley, J.; Leclerc, H. Klebsiella terrigena, a new species from soil and water. Int. J. Syst. Evol. Microbiol. 1981, 31, 116–127. [Google Scholar] [CrossRef] [Green Version]
  25. Shelobolina, E.S.; Sullivan, S.A.; O’Neill, K.R.; Nevin, K.P.; Lovley, D.R. Isolation, characterization, and U (VI)-reducing potential of a facultatively anaerobic, acid-resistant bacterium from low-pH, nitrate- and U (VI)-contaminated subsurface sediment and description of Salmonella subterranea sp. nov. Appl. Environ. Microbiol. 2004, 70, 2959–2965. [Google Scholar] [CrossRef] [Green Version]
  26. Ishii, S.; Ksoll, W.B.; Hicks, R.E.; Sadowsky, M.J. Presence and growth of naturalized Escherichia coli in temperate soils from Lake Superior watersheds. Appl. Environ. Microbiol. 2006, 72, 612–621. [Google Scholar] [CrossRef] [Green Version]
  27. Manter, D.K.; Hunter, W.J.; Vivanco, J.M. Enterobacter soli sp. nov.: A lignin-degrading γ-Proteobacteria isolated from soil. Curr. Microbiol. 2011, 62, 1044–1049. [Google Scholar] [CrossRef]
  28. Zhang, H.; Guo, S.H.; Sun, B.; Zhang, J.; Cheng, M.G.; Li, Q.; Hong, Q.; Huang, X. Mangrovibacter yixingensis sp. nov., isolated from farmland soil. Int. J. Syst. Evol. Microbiol. 2015, 65, 2447–2452. [Google Scholar] [CrossRef]
  29. Bagley, S.T.; Seidler, R.J.; Brenner, D.J. Klebsiella planticola sp. nov.: A new species of Enterobacteriaceae found primarily in nonclinical environments. Curr. Microbiol. 1981, 6, 105–109. [Google Scholar] [CrossRef]
  30. Baek, C.; Shin, S.K.; Yi, H. Limnobaculum parvum gen. nov., sp. nov., isolated from a freshwater lake. Int. J. Syst. Evol. Microbiol. 2019, 69, 1826–1830. [Google Scholar] [CrossRef] [PubMed]
  31. Niu, S.; Ma, W.; Jin, M.; Chen, J.; Li, S.; Zou, X. Complete genome sequence of Kosakonia sp. strain CCTCC M2018092, a fucose-rich exopolysaccharide producer. Microbiol. Resour. Announc. 2019, 8, e00567-19. [Google Scholar] [CrossRef] [Green Version]
  32. Martínez, L.; Caballero-Mellado, J.; Orozco, J.; Martínez-Romero, E. Diazotrophic bacteria associated with banana (Musa spp.). Plant Soil 2003, 257, 35–47. [Google Scholar] [CrossRef]
  33. Stephan, R.; Van Trappen, S.; Cleenwerck, I.; Vancanneyt, M.; De Vos, P.; Lehner, A. Enterobacter turicensis sp. nov. and Enterobacter helveticus sp. nov., isolated from fruit powder. Int. J. Syst. Evol. Microbiol. 2007, 57, 820–826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Peng, G.; Zhang, W.; Luo, H.; Xie, H.; Lai, W.; Tan, Z. Enterobacter oryzae sp. nov., a nitrogen-fixing bacterium isolated from the wild rice species Oryza latifolia. Int. J. Syst. Evol. Microbiol. 2009, 59, 1650–1655. [Google Scholar] [CrossRef]
  35. Brady, C.; Cleenwerck, I.; Venter, S.; Coutinho, T.; De Vos, P. Taxonomic evaluation of the genus Enterobacter based on multilocus sequence analysis (MLSA): Proposal to reclassify E. nimipressuralis and E. amnigenus into Lelliottia gen. nov. as Lelliottia nimipressuralis comb. nov. and Lelliottia amnigena comb. nov., respectively, E. gergoviae and E. pyrinus into Pluralibacter gen. nov. as Pluralibacter gergoviae comb. nov. and Pluralibacter pyrinus comb. nov., respectively, E. cowanii, E. radicincitans, E. oryzae and E. arachidis into Kosakonia gen. nov. as Kosakonia cowanii comb. nov., Kosakonia radicincitans comb. nov., Kosakonia oryzae comb. nov. and Kosakonia arachidis comb. nov., respectively, and E. turicensis, E. helveticus and E. pulveris into Cronobacter as Cronobacter zurichensis nom. nov., Cronobacter helveticus comb. nov. and Cronobacter pulveris comb. nov., respectively, and emended description of the genera Enterobacter and Cronobacter. Syst. Appl. Microbiol. 2013, 36, 309–319. [Google Scholar]
  36. Behera, P.; Ramana, V.V.; Maharana, B.; Joseph, N.; Vaishampayan, P.; Singh, N.K.; Shouche, Y.; Bhadury, P.; Mishra, S.R.; Raina, V.; et al. Mangrovibacter phragmitis sp. nov., an endophyte isolated from the roots of Phragmites karka. Int. J. Syst. Evol. Microbiol. 2017, 67, 1228–1234. [Google Scholar]
  37. Almasia, R.; Henríquez, M.; Levican, A.; Poblete-Morales, M. Genome sequence of a potentially new Buttiauxella species, Strain B2, isolated from rhizosphere of olivillo trees (Aextoxicon punctatum). Microbiol. Resour. Announc. 2020, 9, e01351-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Madhaiyan, M.; Saravanan, V.S.; Blom, J.; Smits, T.H.; Rezzonico, F.; Kim, S.J.; Weon, H.Y.; Kwon, S.W.; Whitman, W.B.; Ji, L. Phytobacter palmae sp. nov., a novel endophytic, N2 fixing, plant growth promoting Gammaproteobacterium isolated from oil palm (Elaeis guineensis Jacq.). Int. J. Syst. Evol. Microbiol. 2020, 70, 841–848. [Google Scholar] [CrossRef] [PubMed]
  39. Kämpfer, P.; Ruppel, S.; Remus, R. Enterobacter radicincitans sp. nov., a plant growth promoting species of the family Enterobacteriaceae. Syst. Appl. Microbiol. 2005, 28, 213–221. [Google Scholar] [CrossRef]
  40. Rameshkumar, N.; Lang, E.; Nair, S. Mangrovibacter plantisponsor gen. nov., sp. nov., a nitrogen-fixing bacterium isolated from a mangrove-associated wild rice (Porteresia coarctata Tateoka). Int. J. Syst. Evol. Microbiol. 2010, 60, 179–186. [Google Scholar] [CrossRef]
  41. Wei, C.Y.; Lin, L.; Luo, L.J.; Xing, Y.X.; Hu, C.J.; Yang, L.T.; Li, Y.R.; An, Q. Endophytic nitrogen-fixing Klebsiella variicola strain DX120E promotes sugarcane growth. Biol. Fertil. Soils 2014, 50, 657–666. [Google Scholar] [CrossRef]
  42. Kandel, S.L.; Joubert, P.M.; Doty, S.L. Bacterial endophyte colonization and distribution within plants. Microorganisms 2017, 5, 77. [Google Scholar] [CrossRef] [Green Version]
  43. Li, Q.; Guo, R.; Li, Y.; Hartman, W.H.; Li, S.; Zhang, Z.; Tringe, S.G.; Wang, H. Insight into the bacterial endophytic communities of peach cultivars related to crown gall disease resistance. Appl. Environ. Microbiol. 2019, 85, e02931-18. [Google Scholar] [CrossRef] [Green Version]
  44. Sharma, S.; Chen, C.; Navathe, S.; Chand, R.; Pandey, S.P. A halotolerant growth promoting rhizobacteria triggers induced systemic resistance in plants and defends against fungal infection. Sci. Rep. 2019, 9, 4054. [Google Scholar] [CrossRef]
  45. Dickey, R.S.; Zumoff, C.H. Emended description of Enterobacter cancerogenus comb. nov. (formerly Erwinia cancerogena). Int. J. Syst. Evol. Microbiol. 1998, 38, 371–374. [Google Scholar] [CrossRef]
  46. Chung, Y.R.; Brenner, D.J.; Steigerwalt, A.G.; Kim, B.S.; Kim, H.T.; Cho, K.Y. Enterobacter pyrinus sp. nov., an organism associated with brown leaf spot disease of pear trees. Int. J. Syst. Evol. Microbiol. 1993, 43, 157–161. [Google Scholar] [CrossRef] [Green Version]
  47. Brady, C.; Denman, S.; Kirk, S.; Venter, S.; Rodríguez-Palenzuela, P.; Coutinho, T. Description of Gibbsiella quercinecans gen. nov., sp. nov., associated with Acute Oak Decline. Syst. Appl. Microbiol. 2010, 33, 444–450. [Google Scholar] [CrossRef]
  48. Jang, J.; Hur, H.G.; Sadowsky, M.J.; Byappanahalli, M.N.; Yan, T.; Ishii, S. Environmental Escherichia coli: Ecology and public health implications-a review. J. Appl. Microbiol. 2017, 123, 570–581. [Google Scholar] [CrossRef] [Green Version]
  49. Gorrasi, S.; Pasqualetti, M.; Franzetti, A.; Pittino, F.; Fenice, M. Vibrio Communities Along a Salinity Gradient Within a Marine Salterns Hypersaline Environment (Saline di Tarquinia, Italy). Environ. Microbiol. 2020, 22, 4356–4366. [Google Scholar] [CrossRef] [PubMed]
  50. Alfinito, S.; Iberite, M.; Fumanti, B. The Algal Microflora of the Saltworks of Tarquinia (Italy). Hydrobiology 1990, 203, 137–146. [Google Scholar] [CrossRef]
  51. Barghini, P.; Silvi, S.; Aquilanti, A.; Marcelli, M.; Fenice, M. Bacteria from Marine Salterns as a Model of Microorganisms Adapted to High Environmental Variations. J. Environ. Prot. Ecol. 2014, 15, 897–906. [Google Scholar]
  52. Barghini, P.; Giovannini, V.; Fenice, M.; Gorrasi, S.; Pasqualetti, M. High Lutein Production by a Halo-tolerant Strain of Dunaliella sp. (Chlorophyceae) Isolated from Solar Salterns in Central Italy. J. Environ. Prot. Ecol. 2018, 19, 704–712. [Google Scholar]
  53. Barghini, P.; Pasqualetti, M.; Gorrasi, S.; Fenice, M. Bacteria from the “Saline di Tarquinia” Marine Salterns Reveal Very Atypical Growth Profiles with Regards to Salinity and Temperature. Mediterr. Mar. Sci. 2018, 19, 513–525. [Google Scholar] [CrossRef] [Green Version]
  54. Barghini, P.; Pasqualetti, M.; Gorrasi, S.; Fenice, M. Study of Bacterial Diversity of a Saltern Crystallisation Pond (‘Saline di Tarquinia’, Italy) and Its Correlation with Salinity Variations. J. Environ. Prot. Ecol. 2018, 19, 139–145. [Google Scholar]
  55. Pasqualetti, M.; Bernini, R.; Carletti, L.; Crisante, F.; Tempesta, S. Salinity and Nitrate Concentration on the Growth and Carotenoids Accumulation in a Strain of Dunaliella salina (Chlorophyta) Cultivated Under Laboratory Conditions. Transit. Water. Bull. 2010, 4, 94–104. [Google Scholar]
  56. Tempesta, S.; Paoletti, M.; Pasqualetti, M. Morphological and Molecular Identification of a Strain of the Unicellular Green Alga Dunaliella sp. Isolated from Tarquinia Salterns. Transit. Water. Bull. 2010, 4, 60–70. [Google Scholar]
  57. Gorrasi, S.; Franzetti, A.; Ambrosini, R.; Pittino, F.; Pasqualetti, M.; Fenice, M. Spatio-Temporal Variation of the Bacterial Communities along a Salinity Gradient within a Thalassohaline Environment (Saline di Tarquinia Salterns, Italy). Molecules 2021, 26, 1338. [Google Scholar] [CrossRef]
  58. Gorrasi, S.; Pesciaroli, C.; Barghini, P.; Pasqualetti, M.; Fenice, M. Structure and Diversity of the Bacterial Community of an Arctic Estuarine System (Kandalaksha Bay) Subject to Intense Tidal Currents. J. Mar. Syst. 2019, 196, 77–85. [Google Scholar] [CrossRef]
  59. Gorrasi, S.; Pesciaroli, C.; Barghini, P.; Pasqualetti, M.; Giovannini, V.; Fenice, M. Metagenetic Profiling of the Bacterial Communities of an Intertidal Pool in Kandalaksha Bay (White Sea, Russia). J. Environ. Prot. Ecol. 2019, 20, 1317–1324. [Google Scholar]
  60. Rice, E.W.; Baird, R.B.; Eaton, A.D.; Clesceri, L.S. Standard Methods for the Examination of Water and Wastewater, 22nd ed.; American Public Health Association: Washington, DC, USA, 2012. [Google Scholar]
  61. Centro Funzionale Regionale, Regione Lazio, Annali. Available online: http://www.idrografico.regione.lazio.it/annali/index.htm (accessed on 1 March 2021).
  62. Huber, J.A.; Welch, D.B.M.; Morrison, H.G.; Huse, S.M.; Neal, P.R.; Butterfield, D.A.; Sogin, M.L. Microbial Population Structures in the Deep Marine Biosphere. Science 2007, 318, 97–100. [Google Scholar] [CrossRef] [Green Version]
  63. Wang, Y.; Qian, P.Y. Conservative Fragments in Bacterial 16S rRNA Genes and Primer Design for 16S Ribosomal DNA Amplicons in Metagenomic Studies. PLoS ONE 2009, 4, e7401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Edgar, R.C. UPARSE: Highly Accurate OTU Sequences from Microbial Amplicon Reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef] [PubMed]
  65. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naive Bayesian Classifier for Rapid Assignment of rRNA Sequences into the New Bacterial Taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef] [Green Version]
  66. Claesson, M.J.; O’Sullivan, O.; Wang, O.; Nikkilä, J.; Marchesi, J.R.; Smidt, H.; de Vos, W.M.; Ross, R.P.; O’Toole, P.W. Comparative Analysis of Pyrosequencing and a Phylogenetic Microarray for Exploring Microbial Community Structures in the Human Distal Intestine. PLoS ONE 2009, 4, e6669. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Legendre, P.; Legendre, L. Numerical Ecology, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2012. [Google Scholar]
  68. Verhoeven, K.J.F.; Simonsen, K.L.; McIntyre, L.M. Implementing false discovery rate control: Increasing your power. Oikos 2005, 108, 643–647. [Google Scholar] [CrossRef]
  69. Edgar, R. Taxonomy annotation and guide tree errors in 16S rRNA databases. PeerJ 2018, 6, e5030. [Google Scholar] [CrossRef] [Green Version]
  70. Strockbine, N.A.; Maurelli, A.T.  Shigella. In Bergey’s Manual of Systematics of Archaea and Bacteria; Whitman, W.E.B., Ed.; John Wiley & Sons: New York, NY, USA, 2015; pp. 1–26. [Google Scholar]
  71. Regione Lazio, Amministrazione Trasparente, Piano Tutela Acque Regionale Aggiornamento. Available online: http://www.regione.lazio.it/binary/rl_amministrazione_trasparente/tbl_contenuti/Piano_Tutela_Acque_Regionale_Aggiornamento.pdf (accessed on 22 April 2021).
Figure 1. Detailed map of the “Saline di Tarquinia” (Viterbo, Italy) marine salterns with indication of the sampling sites. 1 = S, sea; 2 = P5, low-salinity concentrator pond; 3 = P24, intermediate-salinity concentrator pond; and 4 = P37, crystallisation pond. The map was generated using Google Earth Pro version 7.3.1 and graphically edited.
Figure 1. Detailed map of the “Saline di Tarquinia” (Viterbo, Italy) marine salterns with indication of the sampling sites. 1 = S, sea; 2 = P5, low-salinity concentrator pond; 3 = P24, intermediate-salinity concentrator pond; and 4 = P37, crystallisation pond. The map was generated using Google Earth Pro version 7.3.1 and graphically edited.
Water 13 01443 g001
Figure 2. Occurrence of Enterobacteriaceae among the ST sampling sites during the two-year period, May 2012–April 2014.
Figure 2. Occurrence of Enterobacteriaceae among the ST sampling sites during the two-year period, May 2012–April 2014.
Water 13 01443 g002
Figure 3. Box and whisker plots showing the abundance variation of Enterobacteriaceae (A) and Escherichia/Shigella (B) across the sampling sites. The thick lines represent the median; the box upper and lower limits indicate the 25th and the 75th percentiles, respectively; the whiskers indicate the data that lie beyond the 5th percentile (lower whisker) and the 75th percentile (upper whisker); and the dots represent the outliers. The different letters indicate differences between the mean values of different groups.
Figure 3. Box and whisker plots showing the abundance variation of Enterobacteriaceae (A) and Escherichia/Shigella (B) across the sampling sites. The thick lines represent the median; the box upper and lower limits indicate the 25th and the 75th percentiles, respectively; the whiskers indicate the data that lie beyond the 5th percentile (lower whisker) and the 75th percentile (upper whisker); and the dots represent the outliers. The different letters indicate differences between the mean values of different groups.
Water 13 01443 g003
Figure 4. Principal component analysis (PCA) of sampling sites according to the abundance of the Enterobacteriaceae genera. The panels show genera according to their occurrence in ST samples (A) and sample ordering according to genera contribution (B). In panel B, symbol size is indicative of the Escherichia/Shigella abundance in each sample. The percentages of explained variation by the first and the second axes are 82.56 and 17.02%, respectively.
Figure 4. Principal component analysis (PCA) of sampling sites according to the abundance of the Enterobacteriaceae genera. The panels show genera according to their occurrence in ST samples (A) and sample ordering according to genera contribution (B). In panel B, symbol size is indicative of the Escherichia/Shigella abundance in each sample. The percentages of explained variation by the first and the second axes are 82.56 and 17.02%, respectively.
Water 13 01443 g004
Figure 5. Box and whisker plots showing the abundance variation of OTU5 (A) and OTU 206 (B) across the sampling sites. The thick lines represent the median; the box upper and lower limits indicate the 25th and the 75th percentiles, respectively; the whiskers indicate the data that lie beyond the 5th percentile (lower whisker) and the 75th percentile (upper whisker); and the dots represent the outliers. The different letters indicate differences between the mean values of different groups.
Figure 5. Box and whisker plots showing the abundance variation of OTU5 (A) and OTU 206 (B) across the sampling sites. The thick lines represent the median; the box upper and lower limits indicate the 25th and the 75th percentiles, respectively; the whiskers indicate the data that lie beyond the 5th percentile (lower whisker) and the 75th percentile (upper whisker); and the dots represent the outliers. The different letters indicate differences between the mean values of different groups.
Water 13 01443 g005
Figure 6. Redundancy analysis (RDA) triplot showing the relationships among the ST samples, the Enterobacteriaceae OTUs, and the environmental variables affecting the Enterobacteriaceae community. The first and the second axes explained 30.79% and 0.92% (F = 8.4, p = 0.001) of the variance in the Enterobacteriaceae community, respectively. The red arrows represent the variables sin(Month) and rains. The polygons enclose samples collected from the same site (light blue polygon: S samples; orange polygon: P5 samples; red polygon: P24 samples; brown polygon: P37 samples). The blue arrows represent OTUs; only the names of the two most abundant OTUs are reported.
Figure 6. Redundancy analysis (RDA) triplot showing the relationships among the ST samples, the Enterobacteriaceae OTUs, and the environmental variables affecting the Enterobacteriaceae community. The first and the second axes explained 30.79% and 0.92% (F = 8.4, p = 0.001) of the variance in the Enterobacteriaceae community, respectively. The red arrows represent the variables sin(Month) and rains. The polygons enclose samples collected from the same site (light blue polygon: S samples; orange polygon: P5 samples; red polygon: P24 samples; brown polygon: P37 samples). The blue arrows represent OTUs; only the names of the two most abundant OTUs are reported.
Water 13 01443 g006
Figure 7. GLM plots showing the significant trends of the Enterobacteriaceae family, its key genus, and OTUs in relation to the environmental parameters. The abundance variation of Enterobacteriaceae, Escherichia/Shigella, OTU 5 and OTU 206 in relation to the rains (A); abundance variation of Escherichia/Shigella and OTU 5 in relation to the sin(Month) (B); and abundance variation of Enterobacteriaceae, Escherichia/Shigella and OTU 5 in relation to the salinity (C).
Figure 7. GLM plots showing the significant trends of the Enterobacteriaceae family, its key genus, and OTUs in relation to the environmental parameters. The abundance variation of Enterobacteriaceae, Escherichia/Shigella, OTU 5 and OTU 206 in relation to the rains (A); abundance variation of Escherichia/Shigella and OTU 5 in relation to the sin(Month) (B); and abundance variation of Enterobacteriaceae, Escherichia/Shigella and OTU 5 in relation to the salinity (C).
Water 13 01443 g007
Table 1. Descriptive statistics of RDA.
Table 1. Descriptive statistics of RDA.
VariableExplains %Contribution %Fpp(adj) 1
Sampling site: S11.2632.912.00.000940.00848
Rains8.6625.310.20.000940.00424
Sampling site: P55.9917.57.50.004710.01414
sin(Month)5.3415.67.20.004710.01060
Sampling site: P370.270.80.40.678550.69180
Sampling site: P240.270.80.40.691800.69180
1 Corrected according to the FDR.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gorrasi, S.; Pasqualetti, M.; Franzetti, A.; Gonzalez-Martinez, A.; Gonzalez-Lopez, J.; Muñoz-Palazon, B.; Fenice, M. Persistence of Enterobacteriaceae Drawn into a Marine Saltern (Saline di Tarquinia, Italy) from the Adjacent Coastal Zone. Water 2021, 13, 1443. https://doi.org/10.3390/w13111443

AMA Style

Gorrasi S, Pasqualetti M, Franzetti A, Gonzalez-Martinez A, Gonzalez-Lopez J, Muñoz-Palazon B, Fenice M. Persistence of Enterobacteriaceae Drawn into a Marine Saltern (Saline di Tarquinia, Italy) from the Adjacent Coastal Zone. Water. 2021; 13(11):1443. https://doi.org/10.3390/w13111443

Chicago/Turabian Style

Gorrasi, Susanna, Marcella Pasqualetti, Andrea Franzetti, Alejandro Gonzalez-Martinez, Jesus Gonzalez-Lopez, Barbara Muñoz-Palazon, and Massimiliano Fenice. 2021. "Persistence of Enterobacteriaceae Drawn into a Marine Saltern (Saline di Tarquinia, Italy) from the Adjacent Coastal Zone" Water 13, no. 11: 1443. https://doi.org/10.3390/w13111443

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop