Biomonitoring along the Tropical Southern Indian Coast with Multiple Biomarkers

We assessed the spatial and temporal variations of pollution indicators and geochemical and trace metal parameters (23 in total) from water and sediment (144 samples) of three different eco-niches (beach, fishing harbor, and estuary) in larger coastal cities of southern India (Cuddalore and Pondicherry) for one year. A total of 120 marine Pseudomonas isolates were challenged against different concentrations of copper solutions and 10 different antibiotics in heavy metal and antibiotic resistance approaches, respectively. The study shows that 4.16% of the isolates could survive in 250 mM of copper; 70% were resistant to minimum concentrations. Strains were resistant (98.4%) to at least one antibiotic in Cuddalore compared to the Pondicherry (78.4%) region. Pollution index (PI) (0–14.55) and antibiotic resistance index (ARI) (0.05–0.10) ratio indicated that high bacterial and antibiotic loads were released into the coastal environment. The degree of trace metal contamination in sediments were calculated by enrichment factor (EF), contamination factor (CF), pollution load index (PLI), and geo-accumulation index (Igeo). Statistical parameters like two-way analysis of variance (ANOVA), correlation, factor analysis and scatter matrix tools were employed between the 23 parameters in order to find sources, pathways, disparities and interactions of environmental pollutants. It indicates that geochemical and biological parameters were not strongly associated with each other (except a few) and were affected by different sources. Factor analysis elucidated, ‘microbe–metal’ interaction (Factor 1–48.86%), ‘anthropogenic’ factor (Factor 2–13.23%) and ‘Pseudomonas–Cadmium’ factor (Factor 3–11.74%), respectively.


Introduction
Coastal pollution can be caused by activities in large urban centers that have insufficient cleaning and being without sewage and wastewater treatment facilities. This may be particularly problematic during high tourist seasons and mass public visits on weekends because the beach (S1), Cuddalore fishing harbor (S2), Cuddalore Uppanar river estuary (S3), Pondicherry beach (S4), Pondicherry fishing harbor (S5) and Pondicherry Sangarabarani river estuary (S6). The sampling locations were demarcated by using a geographical positioning system (GPS) (Garmin GPS 60).

Sampling
A total of 144 (water and sediment) samples from three different eco-niches in both cities (Cuddalore and Pondicherry) were collected during one year from December 2013 to November 2014 (sampling dates: Postmonsoon = December 24th; January 26th & February 25th; Summer = March 24th; April 25th & May 24th; Premonsoon = June 24th; July 25th & August 25th; Monsoon = September 25th; October 24th & November 25th) for three different analytical approaches (bacteriological, geochemical, and trace metal analysis). Single random sampling was performed at three times during each season. A sample of 2000 ml of marine water was collected from 0-20 cm below the surface with a 2500 ml sterile container at each location. Physicochemical parameters i.e., pH, electrical conductivity (EC), total dissolved solids (TDS), and salinity were measured using a field kit (Thermo Orion 5-Star pH Multi-Meter) at the site and, dissolved oxygen (DO), biological oxygen demand (BOD), total alkalinity (TA) and total hardness (TH) were measured in the laboratory using a standard procedure [12,13]. Approximately 250 g of surface sediment samples were collected (0-5 cm depth) with a sterile spatulum and stored in aseptic polyethylene bags. All samples were kept in iceboxes and processed within 10 h of collection. For heavy metal analysis of water, one liter of sea water was acidified immediately with concentrated nitric acid.

Bacteriological analysis
Bacterial populations from water and sediment samples were analyzed by the pure culture technique (spread plating method) on specific selective media (S1 Table) plates with 100 μL of suitable serial dilutions [14]. The media plates were incubated at 37 ± 1˚C for 24 to 48 h except M-FC plates (M-FC agar plates were incubated at 44.5˚C ± 1˚C for 24-48 h) when the final counts of colonies were noted. Bacterial colonies from the selective media were considered as similar ("-like") to known microorganisms (LO) and Rapid Microbial Limit Test kits were used for bacterial identification [3,15]. The pollution index (PI) was calculated using the ratio of fecal coliforms (FC)/ fecal Streptococci (FS) [2,16]. All culture media were obtained from Hi-Media Pvt. Ltd., Bombay, India.

Determination of antibiotic resistance
A total of 120 marine Pseudomonas were isolated from water and sediment samples by using selective media (cetrimide agar) plates. Twenty colonies of Pseudomonas sp. were isolated from each region and challenged against ten standard antibiotics on Mueller Hinton agar (MHA) by the disc diffusion method for multiple antibiotic resistance analysis [2]. The results were interpreted based on the recommendations of the National Committee for Clinical Laboratory Standards for antimicrobial susceptibility testing [17]. Ten standard antibiotic discs represented seven different chemical structural classes of antibiotics: aminoglycosides (gentamycin 10 mcg), ß-lactams (amoxicillin 10 mcg, ampicillin 10 mcg, methicillin 10 mcg and penicillin-G 10 mcg), glycopeptides (vancomycin 10 mcg), macrolides (erythromycin 10 mcg), quinolones (ciprofloxacin 10 mcg), tetracyclines (tetracycline 10 mcg) and others (chloramphenicol 10 mcg). The antibacterial resistance index (ARI) of each location was calculated according to Hinton et al. [18], using the formula ARI = y / nx, where 'y' represents the actual number of resistance determinants recorded from a population of size 'n' and 'x' as the total number of antibacterials tested in the sensitivity test.

Assessment of metal toxicity
A total of 120 Pseudomonas isolates were challenged against five different concentrations (5, 10, 50, 100 and 250 mM) of copper metal salt solution for metal resistance studies by plate diffusion and tube dilution methods. A stock solution of copper (CuSO 4 .5H 2 0) was prepared in triple distilled water and was sterilized at 121˚C for 15 min. For a plate diffusion assay, 500 μL of copper solution was added to a central well (1 cm in diameter and 4 mm in depth) of each nutrient agar plate to allow metal diffusion for one day. Eight isolates were inoculated in each plate by the radial streaking method and were incubated at 37±1˚C for 48 h. All the trials were performed in triplicate. After incubation, the percentage of bacterial tolerance was calculated in terms of the ratio: length of growth in mm vs length of the total inoculated streak. For the tube dilution method, the appropriate amount of metal solution and 200 μL of standard culture (test culture suspension prepared in sterile 0.85% saline matching an optical density of 0.5 McFarland standards corresponding to 108 CFU/mL) were added into nutrient broth containing test tubes, and the final volume of 10 ml was made-up with sterile nutrient broth. Tubes were read after incubation at 37±1˚C for 5 days [19].

Assessment of heavy metals
One liter of sea water was filtered through a 0.45 μm nitrocellulose membrane filter paper and adjusted to pH 2 with HNO 3 taken in a separatory funnel. From a freshly prepared solution of amino-pyrolidine dithiocarbamate (APDC) 10 ml (3% w/v) were added into the funnel, and the mixture was shaken by a mechanical shaker for 10 min. Furthermore, 25 ml of methyl-isobutyl-ketone (MIBK) was added to this mixture and shaken for 15 minutes. The phases were allowed to separate and the top organic phase was collected. The bottom aqueous phase was again shaken with 25 ml of MIBK, and the organic phase was obtained and pooled with the previous phase. The pooled organic phase was mixed with 2 ml of 50% HNO 3 , and shaken vigorously for 10 min to separate the bottom acid layer [20]. The marine surface sediment samples were air dried at room temperature and homogenized using an agate mortar and pestle, and sieved using a 63 μm sieve. The particles > 63 μm in size were retained in pre-cleaned and sterilized plastic bottles for trace metal analysis. The trace metals in marine sediment were analyzed as described by Thuy et al. [21] in a slightly changed procedure. One gram of dried sediment was added to 10 ml of 'aqua regia' (i.e. HCl:HNO 3 = 3:1) in Teflon PTFE vessels and treated at 100˚C in a hot air oven for 24 h. After incubation, the residue was dissolved with 20 ml of HNO 3 and filtered in a Millipore filtration unit with 0.45 μm nitrocellulose filter paper. The final volume was made up with HNO 3 to a volume of 25 ml in a volumetric tube.
The concentration of metals in each water and sediment samples were determined by atomic absorption spectrophotometry (GBC SensAA-AAS, Australia) by a flame method. The blanks and internal standards were run at regular intervals to verify the accuracy of the method. Percent recoveries of metals ranged from 70 to 88% and the relative standard deviation of triplicate measurements was < 5%. Four important parameters [22] were used to evaluate the degree of trace metal contaminations in marine sediments such as enrichment factor (EF), contamination factor (CF), pollution load index (PLI) and geo-accumulation index (I geo ).

EF ¼ ðMe=FeÞsample=ðMe=FeÞbackground ð1Þ
where, (Me/Fe) sample is the metal to Fe ratio in the sample of interest; (Me/Fe) background is the natural background value of metal to Fe ratio.
where, C m Sample is the concentration of a given metal in coastal sediments, and C m Background is the value of the metal equals to the world surface rock average.
where, n is the number of metals (seven in the present study).
where, C m Sample is the measured concentration of elements n in the sediment sample and the C m Background is the geochemical background value. The factor 1.5 is introduced to include possible variation of the background values due to lithogenic effects. The enrichment factor (EF) categories, contamination level of CF and geo-accumulation index classes (seven) are given in S2 Table. Statistical analysis Pearson correlation coefficient and scatter matrix was employed for a better understanding of the relationship between the concentration of multiple variables (trace metals, geochemical, and bacteriological parameters) simultaneously by using the statistical package ORIGIN8.0. The two way analysis of variance (ANOVA) and scatter matrix was employed (ORIGIN8.0) to understand the variation in the variables between different stations, different locations and their interactions. A scatter matrix reveals a correlation pattern between two parameters and also indicates whether those parameters are supporting each other or not. A 45-degree line diagonal plot indicates a strong relation of two variables while the distribution of the plot explains a lower relation between the variables [13]. Principal component analysis (PCA)/ factor analysis (FA) was also performed (SPSS16.0) to trace the source/ pathway of the parameters with respect to the environment.

Geochemical analysis
The average mean value of pH was 7.4, 7.9, 7.6 and 7.1 in premonsoon, monsoon, postmonsoon and summer season, respectively (S1 Fig). The EC in water was due to ionization of dissolved inorganic solids and became a measure of calculated TDS. The EC, TDS and salinity ranged from 26000-67800 μS/cm, 12950-48525 mg/L and 20-38 ppt, respectively. In Cuddalore fishing harbor DO and BOD ranged from 3.3-8.3 mg/L and 3.2-7.9 mg/L, and the DO and BOD levels were highly fluctuating during all seasons. Mean level of TA and TH were 162 mg/L and 101.6 mg/L in monsoon; 106 mg/L and 79.1 mg/L in summer; 114 mg/L and 56 mg/ L in premonsoon; 92 mg/L and 61 mg/L in postmonsoon.

Indicators of bacterial populations
High levels of pollution indicating bacteria in coastal ecosystems were also indicating a common problem in all investigated areas that often lead to an impairment of beneficial uses and human health. Thus, we represent here the average mean values of all parameters in 144 samples (water and sediment) from the coastlines of both cities for all seasons showing a great variety of distribution patterns. The TVC for water and sediment samples were highest during the monsoon season and least during the postmonsoon. In water samples, the mean TVC ranged from 1. The TS did not show any location-wise trend and the mean densities were low during the postmonsoon. Similar observations were made by Vignesh and co-workers [1] at Tamil Nadu beaches and Kumarasamy et al [14] in the Cauvery estuary, Tamil Nadu. In water, FC counts were high at sampling site S2 during all seasons and the values were 1080 mL−1, 2200 mL−1, 440 mL−1 and 320 mL−1 during premonsoon, monsoon, postmonsoon and summer, respectively. Similar patterns were also noted in sediment samples. The VLO counts were higher at the Cuddalore coast than at the Pondicherry coast. In water, the mean values of VLO densities were 1268.33 mL−1, 1685.25 mL−1, 448.33 mL−1 and 1021 mL −1 during premonsoon, monsoon, postmonsoon and summer, respectively. In sediments, the VLO ranges were 2-5 folds higher than in water samples. During monsoon, the counts of VC and PA were generally higher at all sampling sites. Nagvenkar and Ramaiah [4] made similar observations in Mandovi and Zuari estuary at Goa. In water, Cuddalore fishing harbor (S2) showed high VC counts at all seasons i.e., 1700 mL−1, 2600 mL−1, 800 mL−1, 640 mL−1 during premonsoon, monsoon, postmonsoon, and summer, respectively. In sediments, the mean value of VC during monsoon were 1700 mL−1, 7700 mL−1, 3300 mL−1, 2000 mL−1, 3100 mL −1 and 1800 mL−1 at S1, S2, S3, S4, S5 and S6, respectively. In sediments, the counts of Pseudomonas aeruginosa (PA) i.e., 30 mL−1, 0 mL−1, 80 mL−1 and 30 mL−1 during postmonsoon, summer, premonsoon and monsoon, respectively, was found at Cuddalore Uppanar estuary (S3). The PI value was above 1 (> 1) indicating that the sampling sites were highly contaminated by human fecal matter rather than animal or other sources (Fig 3A and 3B).

Antibiotic resistance
A comparison of antibiotic resistant Pseudomonas levels of two different coastal cities is provided in Table 1. In this study, antimicrobial resistance patterns showed that the resistance behavior was diverse between the sampling locations and types of samples owing to human fecal matter. In Cuddalore, more than 30% of isolates were resistant to single antibiotics whereas double and 3 antibiotic resistances were recorded in Pseudomonas sp. such as 18.3% and 16.6%, respectively. Among a total of 60 Pseudomonas isolates, 1.6% were susceptible to all antibiotics in Cuddalore while 21.6% were susceptible to all antibiotics in Pondicherry. The Cuddalore region showed a higher percentage (98.4%) of resistant strains than the Pondicherry region (78.4%) indicating that most of the strains were resistant to one or more of the tested antibiotics. With respect to antibiotic classes, the Pseudomonas isolates from the Cuddalore and Pondicherry coast were resistant to 10 types of antibiotics belonging to 7 different groups: β-lactam (23.9% and 17.3%), glycopeptide (56.6% and 46.6%), macrolides (18.3% and 11.6%), aminoglycosides (11.6% and 8.3%), tetracycline (8.3% and 15%), quinolone (0% and 0%) and others (8.3% and 5%). In sea water, antibiotic the resistance index (ARI) was higher at Cuddalore beach (0.10) and Cuddalore fishing harbor (0.10), followed by Cuddalore Uppanar river estuary (0.09), Pondicherry fishing harbor (0.08), Pondicherry Sangarabarani river estuary (0.07), and Pondicherry beach (0.05) ( Fig 3C). For both water and sediments, the ARI rate was higher in the Cuddalore region than in the Pondicherry coastal zone, indicating that the accumulation of sewage discharges and fecal matters are higher in the Cuddalore coastal region.

Metal toxicity
About 120 marine Pseudomonas isolates were investigated for their copper resistance from two coastal cities of southern India. The patterns of metal resistance were varying depending on their five different (5, 10, 50, 100 and 250 mM) concentrations. Average concentrations of heavy metals in seawater and sediments were about 50 mM [23]. In the present study, the tube dilution (or) minimal inhibitory concentration (MIC) method confirmed that only few of the 120 isolates could survive at higher (250 mM) concentrations of copper among the remaining isolates. More than 70% of the strains were resistant to minimum levels (5 mM) of metal concentrations whereas a maximum of 4.16% of the strains grew at 250 mM copper concentration. These strains could be used for metal extraction through biosorption. Table 2  percentages of microbial growth (which were producing very distinct effects) and their minimal inhibitory concentration levels against different metal concentrations. The Pseudomonas isolates were easily tolerating lower concentrations while higher concentrations strongly affected bacterial growth. A growth rate between 90-100% was observed for 5.8% of the bacterial populations at 5 mM of Cu, whereas no population was growing at a growth rate of 0-40% with 5 mM of Cu. At 10 mM of Cu, 40% of the populations showed a growth rate of 61-70%. In 50 mM, 56.6% of the populations showed growth rates between 51-60%, whereas 33.3% of the populations were observed with a 51-60% growth rate at 100 mM concentrations. At 250 mM Cu, 53.3% of the populations showed 41-50% growth rate. A higher metal resistant Pseudomonas sp. was observed at Cuddalore than at Pondicherry.

Assessment of trace metals
Concentrations of various trace metals in different water and sediment samples are shown in    were less than 1 at all stations except Cd. This indicates a low level of sediment contamination in the proposed study area. The CF values were suggesting that these stations were highly polluted with Cd, which was due to its location nearest to the industrial effluents and the aquaculture sewage outfall point of Cuddalore coastal area. The heavy metal enrichment factor is a measure of untreated effluents and sewage discharges in the aquatic environment. The enrichment factor has been used to examine the anthropogenic contribution of metals in sediments and to differentiate the metal source. This holds for anthropogenic or natural based elements against background levels. At some of the locations, the EF value of Cd was higher than 40, indicating a higher Cd pollution in these regions and, Pb and Zn were crossing the EF value of 2, suggesting that these areas were only moderately polluted (Table 3). However, the enrichment factor values obtained for most heavy metals (Cu, Ni, Fe, Cr & Pb) being less than 1 reveals that these elements were originating from non-anthropogenic heavy metals in the sediments at most of the sampling sites. The I geo values calculated for Cd (1.69) showed higher values throughout the study period in the industrial zone of Cuddalore and Pondicherry. Interestingly, I geo values greater than 1 revealed that a significant portion of metals originated from geological sources.

Two-way analysis of variance (ANOVA)
The two-way analysis of variance (ANOVA) test examines the influence of different categorial independent variables on one dependent variable. It can not only determine the main effect of contributions of each independent variable but also identifies if there is a significant interaction effect between the independent variables. Differences in the concentration of parameters between the sampling locations were statistically significant ( Table 4) which indicated that strong signals of seasonal variations and local variations could be discerned. It was revealed   that significant (ANOVA, P < 0.001 and P < 0.05) variations were observed between different seasons and locations, and the interactions of water and sediments. Descriptive statistics of microbial, physiochemical and trace metal parameters are represented in S3 Table. Results from two-way ANOVA demonstrate that copper in sediment samples had no significant effects in all divisions such as between seasons, between locations and their interactions. In sediment samples, insignificant effects were also observed in PA and Ni at divisions of their interactions (seasons vs locations). Similar to this, an insignificant effect was observed with pH (between locations and their interactions) and salinity (between seasons and their interactions) of water samples. Most of the parameters had a significant effect between seasons and between locations and their interactions. No significant (ANOVA, P > 0.05) effect was observed in PA, salinity and TH between seasons.

Correlation and scatter matrix
Correlations between the parameters can provide information on the relationship of the sources and pathways of the factors. Correlations between the variables are listed in Table 5.
High positive correlations were observed between bacterial parameters except PA with other bacterial parameters. This indicates that bacterial groups were strongly associated with each other ('p'value-<0.0001) due to their same pollution origin/ sources. Interestingly, the geochemical parameters were not associated with bacterial parameters except TA-bacterial parameters. At the same time, heavy metal parameters were positively associated with each other and their parameters were highly significant ('p'value-<0.05), except Cd-Cr, Ni & Zn-Ni. A poor correlation between TH, TA and heavy metals may have resulted from a tremendous deviation in the content of those parameters which proved that they were affected by different sources. A great relocation difference was also observed with physiochemical parameters. Commonly, less positive and negative correlations were observed between salinity and all other parameters which confirmed that they were not mutually associated with each other. A similar observation was also observed between pH and all other parameters. However, bacterial parameters were mutually associated with heavy metal parameters confirming that the minimum level of trace metals were acting as micronutrients and were supporting microbial growth with high significances ('p'value-<0.05).
The scatter matrix (Fig 5) explained that the distribution pattern of parameter levels and the linear line distribution indicated that parameters were mutually associated with each other. This pattern was also used to estimate the co-variance matrix for example of the normal distribution of multi-variance. The pair-wise plots are explaining the univariate distribution of the variables and also indicate the deviated values which are not plotted on the 45-degree liner line. The scatter matrix showed a strong association between the bacterial parameters except PA and a similar pattern is also followed between the Cu and Fe; EC and TDS. While dispersed values were observed between the PA, pH, DO, BOD, TA and TH = all parameters (except TVC, TC, TS, FC, FS, VLO and VC). But, a moderate level in the 45˚diagonal plot was observed between the remaining parameters especially between trace metals.

Factor analysis
To evaluate the possible relationships among variables, we performed statistical analysis such as principle component analysis (PCA)/ factor analysis (FA) with varimax normalization (PCA-V) to explore the sources and pathways of the variables. In this study, a factor with an eigen value > 1 was considered for subsequent discussion ( Table 6). The factor I explained 48.86% of the total variance and was best represented by bacterial parameters (except PA), heavy metals (expect Cd and Zn) and TA. This factor represents a pollution source from sewage and fecal contaminations throughout the coastal regions and it could be considered as a  'microbe-metal' interaction factor. An overwhelming 13.23% of the total variance was contributed by factor II, showing higher loading for geochemical parameters such as EC (0.929), TDS (0.929), and salinity (0.625) (except pH, DO, BOD, TA and TH), which is considered as an 'anthropogenic' factor. Factor 3 explained 11.74% of the total variation and has a positive loading on PA and Cd, a moderate positive loading on BOD, whereas a high negative loading on pH and TH is considered as a 'Pseudomonas-Cadmium' factor.

Discussion
A rapid increase of anthropogenic activities in coastal areas has caused a significant increase of pollutants and direct impacts on coastal ecosystems, particularly in vulnerable regions. Accumulation of pollutants in coastal ecosystems degrade environmental quality and deplete environmental filtering and self-purification capacity [26]. Microbial indicators have been used worldwide to indicate whether human/ animal wastes have contaminated a water body. The present study was conducted to evaluate the microbial (pollution indicator loads and antibiotic resistant and trace metal resistant strain loads) and metal pollution in populated coastal cities (Cuddalore and Pondicherry). This study was also focusing on the relationship between pollutants, both season-and location-wise. Moreover, the two-way ANOVA, correlation and PCA were helpful to determine whether the variables were significantly different between season and location, as well as their interactions. These statistical tools were shown to be useful for an effective understanding and interpretation of the different variables belonging to three different parameters from both samples. Total coliforms comprise a large group of organisms of which a major proportion detected in water samples might be of fecal origin, indicating that total coliforms might be a reliable indicator for fecal contamination [27]. The World Health Organization [28] recommended that fecal coliform counts must not exceed 1000/ 100 ml for intermediate contact and 130/ 100 ml for full recreational contact. In the present study from Indian waters, the coliform and Streptococci were lowest during post monsoon and highest during monsoon. A similar pattern was also reported in Indian coastal zones such as Tamil Nadu beaches [1], Cauvery river estuary [14], Tamirabarani river estuary [13], Mandovi and Zuari estuary [4], Gujarat waters [29], the Visakhapatnam coast [3], and Mumbai waters [30].
The fishing harbors (S2 and S5) of both cities received more pollution due to dense human populations and the confluence of a variety of discharges and fecal contaminations. Counts of coliforms observed from Cuddalore and Pondicherry areas were at least 2-5 orders of magnitude lower than in Mumbai harbor waters [30], and 1-2 fold lower at Tamil Nadu beaches [1]. The counts of pathogens from Uppanar (S3) and Sangarabarani (S6) estuaries were at least 1-3 orders lower than in the Seine estuary [31]. However, the counts of pathogens were higher in the Cauvery estuary [14,15] than in the S3 and S6 regions. When compared to the Mandovi and Zuari estuary [4], bacterial counts were 1-3 orders lower in the present report. Compared to the Uppanar and Sangarabarani estuary, total coliform and pathogen counts were 5-7 orders higher in the Coovum estuary in Chennai [2].
However, sediments get resuspended due to various biogenic and physical disturbances and sediments may contain 100 to 1000 fold higher bacterial counts than the overlying water layers [5]. The higher bacterial counts in sediments are due to autochthonous action [4], bacterial aggregations by dynamic flocculation [1], rich organic content [13], waves, wind, temperature, and tidal activities [2]. Co-sedimentation, i.e., the joint deposition of sediment particles and bacteria is a common way of bacterial settlement in moving waters. In addition to microorganisms, urine and feces containing organic matter as well as eutrophying substances in the form of phosphorus and nitrogen compounds [32]. Similar results were observed in rivers and coastal waters [1,13]. Thus, higher accumulation of physiochemical and trace metal parameters could not enhance microbial growth. The fecal matters from humans and animals can easily settle in sediment and this protects microorganisms and also enhances their survival [33]. Other nonpoint sources of contamination (malfunctioning sewage, septic systems, storm water drainage, and urban runoff) could also settle in the sediment, leading to bioaugmentation. By limits set by the United States Environmental Protection Agency (USEPA), recreational waters with concentrations exceeding the maximum contaminant limits of 33 or 200 cfu/100 ml for FC, respectively, presents a health risk [34]. Fecal coliform and Streptococci ratios of the study sites showed a pollution index (PI) ratio of more than one [35] which indicated that the sampling sites were possibly contaminated by human excretions. The average PI ratio of the entire study area further indicated specific variations among the locations with all-time occurrences of indicator organisms and it proved that the sampling locations were contaminated by human feces during all seasons.
Marine isolates were highly resistant to antibiotics when compared to clinical isolates [36]. The improper and overuse of antibiotics has promoted the emergence of resistant strains that provided a selective pressure on bacterial strains, and caused an imbalance between susceptible and resistant bacteria [2]. In the recent past, unnecessary usage of antibiotics by humans and animals for various purposes led to changes in microbial genetic ecology through mutation, selection, and the flux of genetic information (horizontal gene transfer) among microbes [37]. The presence of multi-drug resistant strains is alarming because infection with such strains leads to a higher fatality rate than with antibiotic-sensitive strains [38]. Highest antibiotic resistance was observed in marine bacterial isolates from different water and sediment samples of both regions investigated here. More than 85% were resistant to one or more of the tested antibiotics, whereas nearly 10% of the strains were susceptible to the antibiotics tested.
The Cuddalore region got higher resistant strains (98%) than the Pondicherry region (78%) due to municipal discharges, fishing harbor waste and high population density. In both regions, a low frequency (9%) of antibiotic resistance to 6-10 antibiotics was observed in marine Pseudomonas isolates. The highest frequency (33%) of antibiotic resistance to 1-3 antibiotics was observed in Cuddalore strains, whereas 27% of antibiotic resistance to 1-3 antibiotics was observed in Pondicherry strains. A minimal percentage of these organisms showed resistance to 4-5 antibiotics. The antibiotic resistance index (ARI) indicated that the Cuddalore region was higher exposed than the Pondicherry region due to an accumulation of municipal sewage, as well as open defecation from human settlements and animal husbandries. In the present study, a higher ARI ratio corresponded to high-risk antibiotic-exposed sites. This indicates that a system understanding is needed in order to prevent the spread of antibiotic resistant bacteria in the marine environment.
Aquatic bacteria can develop resistance against organic matter, inorganic matter, metals and antibiotics by adapting themselves to survive in extreme environments. Microorganisms exposed to selection pressures by toxic compounds develop resistance rapidly [39] and the resistance may vary depending on metal concentrations. In the present study, plate diffusion and tube dilution (MIC) methods confirmed higher sensitivity of Pseudomonas sp. growth against high copper concentrations (250 mM). Heavy metals generally inhibit microbial growth by blocking essential functional groups, displacing essential metal ions, or modifying the active conformation of biological molecules [23]. The Pseudomonas isolates from fishing harbor areas were much more resistant than in other areas which received different waste materials from shipping activities, domestic waste from slums, paintings, cargo/ ballast water discharges. The mechanism of metal resistance may be affected in two different ways: proteinmetal associations-biomagnification or blockage of cell walls and the systems of membrane transportation [40]. However, in the Cuddalore and Pondicherry regions, the percentage of Cu sensitive strains at 5, 10, 50, 100 and 250 mM concentrations were 11.9%, 23.4%, 46.7%, 78.4% and 95.84% and, 20.0%, 32.5%, 65.0%, 89.2% and 0%, respectively. Metal resistant bacteria detoxify metals with the aid of chaperones such as metalloproteins which reduce metal toxicity in water [23]. Above 60% of the strains were resistant to minimum levels of metal concentrations whereas a maximum of 4.16% of the strains grew at 250mM of copper. Heavy metal resistant strains were high in fishing harbors followed by estuaries and beaches in both cities due to the accumulation of several trace metal wastes in coastal sediments.
Increased accumulation of anthropogenic toxic metals in marine environments (fishing harbour, estuary and beach) of both cities from industrial and municipal sources stimulated their extreme persistence and high toxicity [41,42]. In the present study, variations of metal concentrations were compared between both water and sediment samples and were also compared season-and location-wise. The metal concentrations were low in most of the places when compared to the permissible limits of international standards [43,44,45] with some exceptions. The metal concentrations at all the sites from both cities did not show any particular trend and it was indicated that the sources might be diffuse. The abundance of metals in the study regions were in the order Fe > Zn > Cu > Cd > Pb > Ni > Cr. One can observe that the areas influenced by fishing have higher concentrations of toxic metals than other areas throughout the year. The metal concentrations in fishing harbor areas were probably related to nonpoint source influxes from various discharges such as dredging, cargo handling, dumping of ship wastes, spilling of cargo, chemicals, metal ores, ballast water and paint waste directly into the sea [2]. Sewage and urban waste water from town and several settlements, constitute other discharges that are negatively affecting water quality [46]. Similar to the fishing harbors, the Cuddalore estuary (S3) recorded higher values of all metals. The higher concentrations at this site might be due to contributions from agricultural farms and fishing communities dotted along the banks of this river that contained agrochemicals, industrial wastes, and municipal sewage.
Cadmium is extremely toxic and high levels of Cd in water could cause adverse health effects to consumers such as renal diseases and cancer [47]. The WHO has established a human health based guideline of 0.003 mg/L for direct usage [48]. These guidelines were exceeded at some of the sampling sites. The WHO recommended value for Zn in water for domestic supplies is 3 mg/L [48] and the Zn level was not crossing the limits in most places of our study. Copper concentrations were exceeding permissible limits (0.02 mg/L) at certain sites due to surface runoff, fishing boat paints and contributions of river and pipeline discharges to the coastal system. It was shown that concentrations of Ni and Pb were high in some of the places at all seasons and exceeded maximum permissible limits (0.01 and 0.1 mg/ L, respectively) at a few sampling sites.
In the factor analysis, factor 1 (48.86%) was influenced by bacterial parameters (except PA), heavy metals (expect Cd and Zn) and TA. Similar observations were recorded along the Chennai coast (Santhiya et al., 2012) and our earlier studies also proved that microbial and physiochemical parameters were not influenced by the same factor [1,13]. Factor 2 (13.23%) was contributed by geochemical parameters. The strong loading of these parameters could be explained by anthropogenic activities through sewage/ industrial/ agricultural/ shipping/ fishcatching activities, clearing of lands, runoff, and erosive processes taking place near the study area [49]. Factor 3 (11.74%) explains the PA and Cd loading with negative pH indicating that PA and Cd levels were fluctuating by precipitation of inorganic substances in this transient coastal aquatic system [50]. These parameters are indicators of organic pollution and the BOD implies that microorganisms utilized DO as biodegradable organic matter [51]. This is also reflecting the control of different sources like anthropogenic discharges with simple hardness association [4]. The DO and BOD did not show either high positive nor high negative values in all the factors. Elevated concentrations of organic substances strengthened anoxia or hypoxia by raising the BOD, the oxygen consumed during microbial decomposition and other organic reduced components [52].
High concentrations of trace metals caused severe health hazards to marine micro-and macrofauna. It may affect their morphology, physiology, biochemistry, enzyme activity, behavior, reproductive system, and causing skin diseases among others. The variation of the enrichment factor (EF) for Cd was significant and the EF values ranged from 9.34-60.84; most of the remaining metals have EF values below 1. This indicates that the sampling stations were more contaminated by cadmium than by other metals. The same pattern was also followed in the geo-accumulation index (I geo ), the cadmium level (range = 0.00-1.74) indicated that the sites were moderately polluted. Most of the contamination factor (CF) values for Cd were above 1 (except S4 during the monsoon season) (range = 0.87-4.83), demarcating that the sampling sites were considerably contaminated. The pollution load index (PLI) (range = 0.03-0.18) clearly explained that the sampling sites were below the baseline pollution category (Tomlinson et al., 1980), providing subsequent bioaccumulation in marine biota.

Conclusion
The studied levels of different pollutants (fecal indicators, antibiotic/ heavy metal resistant strains, physiochemical factors, and trace metals) from water and sediment samples of southern India endowed with major sources of pollution and the levels of bacteria and trace metals were crossing permissible limits of WHO and TNPSC. Different statistical approaches confirmed multiple heterogeneous but also independent sources and pathways of geochemical, microbial, and trace metal parameters. This antibiotic resistance and trace metal study revealed that most of the marine Pseudomonas isolates were resistant to at least one antibiotic. The present study in particular calls for the monitoring and subsequent management of several anthropogenic problems that will be necessary to solve in the frame of sustainable coastal management.