1 Introduction

Water is an essential commodity to all life and there is not a substitute for water (Stikker, 1998). Anthropological activities, industrial discharges, sewage open outlets, hospital waste, and aquaculture effluents have been identified as major contaminants particularly in surface water in Sri Lanka (Mahagamage et al., 2016; Liyanage & Manage, 2018). The quality of the water is the major indicator to determine human health, mainly in populated areas where have a high risk to receive a multitude of contaminants such as antibiotics, heavy metals, organic wastes, and xenobiotics (Qadir et al., 2008; Mahagamage & Manage, 2015). In Sri Lanka, 29% populations only get pipe-borne water facilities where day by day water requirement including drinking is fulfilled from groundwater, river, springs, and streams (CEA, 2017; Mahagamage & Manage, 2019).

In Sri Lanka, there are 103 rivers starting from the central highlands and flowing through the country with increasing fertility (MOFE 2001; Mahagamage & Manage, 2018). The meandering part of the Kelani river basin acts as home to more than 20% of the Sri Lankan population and it fulfills more than 80% of drinking water requirement in Colombo district (Danish Hydraulic Institute, 1999). Compared to the other river basins, the highest urban population is recorded in the Kelani river basin (Mahagamage et al., 2020). Therefore, the pollution rate in the Kelani river is increasing as a result of poor local authority service delivery, industrial discharges, weak environmental management, and governance (CEA, 2017; Manage et al., 2020) coupled with inadequate awareness and education (Mahagamage & Manage, 2018). Finally, it leads to an increased risk of disease transmission to people who live in the Kelani river basin (Mahagamage et al., 2020).

Greater attention has been received from national and international agencies about the antibiotic contamination of surface and groundwater along with antibiotic-resistant bacteria resulting in the development of serious pathogenic diseases in public (WHO, 2018). The use of a wide variety of antibiotics in large amounts, including non-biodegradable antibiotics in medical, aquaculture, and livestock settings, ensures that antibiotic residues remain in the aquatic environment for a long time without degradation (Russell, 2004). According to the World Health Organization (WHO), antibiotics resistant genes (ARGs) are most critical in human health challenges in the next century as ARGs serve as pathways of carrying out genetic manipulation and generate pressure for the development of antibiotic resistance in a susceptible microbial strain through horizontal gene transfer (Manage, 2018; Manage & Liyanage, 2019; Liyanage & Manage, 2016a).

The high adaptive capacity of microorganisms habitually associated with the spread of antibiotic resistance in the environment. Bacteria are easily susceptible to genetic changes and acquisition of mobile genetic elements which are related to antibiotic resistance (Zhang & Li, 2011), dissolved DNA (Hernandes et al., 2013) through horizontal gene transfer (HGT) (Hernandes et al., 2013). Rapid acquisition of antibiotic resistance by pathogenic bacteria may occur via cell to cell contacts with other bacteria, which carry resistance genes on mobile genetic elements (MGEs) (Bennett, 2008; Huddleston, 2014; Manage & Liyanage, 2019). ARGs might exist and persist in the absence of selection and challenge the dynamic equilibrium of ecosystems and enter into humans through food chains (Huddleston, 2014). Thus, the development of ARGs in the environment may create a pandemic situation in global health shortly (Hernandes et al., 2013).

Tetracycline and penicillin are widely used in livestock farming and healthcare. The major mechanisms for tetracycline resistance are prevention of drug uptake (Efflux) into the cells, ribosomal protection, and enzyme inactivation (Roberts, 2005). Efflux genes of tetracycline have been reported as tet A, tet B, tet C, tet D, tet K, and tet L (Roberts, 2005) where the tet M and tet O are common genes for ribosomal protection proteins while only tet S, tet X, tet 34, and tet 37 are coding for a protein capable of enzymatic inactivation of tetracycline (Roberts, 2005). Penicillin is one of the common antibiotics and major penicillin resistance mechanisms are antibiotic hydrolysis mediated by the bacterial enzyme beta-lactamase (bla TEM, bla OXA), changes in penicillin binding proteins (PBP) (amp a), and decreasing of porin channel formation (OPR D).

At present, European Union legislation does not include any specific regulations regarding the presence of antibiotic-resistant bacteria and ARGs in water nor their concentration thresholds. The screening of ARGs and antibiotic-resistant bacteria in the clinical and aquatic environment would therefore be the first step in tackling the rapidly growing resistance to antibiotics (Berendonk et al., 2015). The complexity of the processes and the relative scarcity of studies done, knowledge on ARGs and their role in the environment is still poorly known, and therefore, environmental studies are given priority.

ARGs and the levels, and resistance to antibiotic-resistant bacteria in the surface and groundwater of the meandering zone of the Kelani river basin in Sri Lanka are limited (Liyanage et al., 2015; Mahagamage et al., 2020). Even in Sri Lanka, tetracycline and penicillin are widely used for aquaculture and human medicine. WHO (2018) reported that antibiotic treatment for most of the waterborne diseases has created development of antibiotic resistance in pathogenic bacteria in the environment. Thus, studies are needed to find out the antibiotic resistance in the running water sources which are being used for human consumption that will help develop a strategic plan to protect surface and groundwater for the safeguard of water consumers. This is the first report on the screening and quantification of ARGs in the meandering zone of the Kelani river basin in Sri Lanka. Accordingly, in the present study, two categories of antibiotics (penicillin and tetracycline) and their resistance genes (OPR D, bla TEM, bla OXA, amp a, tet A, and tet M) were selected. The presence of these antibiotics and ARGs was identified and quantified using high-performance liquid chromatography (HPLC) and quantitative polymerase chain reaction (qPCR), respectively. Thus, the present study was aimed to evaluate the prevalence of penicillin and tetracycline resistance genes in surface and groundwater sources within the meandering zone of the Kelani river basin.

2 Methodology

2.1 Sampling Sites and Sample Collection

The Kelani river basin was physically divided into three zones: head, transition, and meandering based on the land-use practices, geology, geography, and anthropogenic activities of the river basin (Mahagamage & Manage, 2018). In the present study, the meandering zone of Kelani river basin which starts from Awissawala to Mattakkuliya which cover three districts (Kaluthara, Gampaha, and Colombo) was selected to collect surface and groundwater samples considering the population density, industrial setting (soap, detergents, food, beverages, and pharmaceuticals), and land-use practices (home gardens, coconut, paddy, and rubber cultivation) (CEA, 2014; Mahagamage et al., 2014). Forty (40) ground and surface water samples were collected (Fig. 1, Table 1) in triplicate using a Rutner sampler (2L) to sterile amber color glass bottles (1L) for antibiotics and ARGs analysis following the WHO sampling protocol (WHO, 2009). Then, the collected water samples were stored in the icebox at 4 °C during transportation and kept in refrigerated condition for subsequent analysis.

Fig. 1
figure 1

Ground and surface water sampling locations with land use in the meandering zone of the Kelani river basin

Table 1 Ground and surface water sampling locations in the meandering zone of the Kelani river basin

2.2 Quantification of Tetracycline and Penicillin in Water Samples

One liter of water sample was adjusted to pH 3 and filtered through 0.22 μm Millipore filters. Filtered water samples were spiked with antibiotics at a final concentration of 100 mg L−1. Sep-Pak Plus C18 cartridge was used to concentrate the antibiotics after conditioning with 100% methanol and following 5 ml deionized water (Liyanage et al., 2015). Prepared samples were passed through the C18 cartridges set up which was held with the SPE unit with maintaining approximately 1–2 ml/ml. The 100% methanol (mobile phase) is used for eluting analytes from cartridges (Liyanage and Manage, 2014; Liyanage and Manage, 2016b). The target antibiotics were quantified by using Agilent 1200 series HPLC equipped with a diode array fluorescence detector (Fernandez-Torres et al., 2010; Liyanage and Manage, 2016a).

2.3 Extraction of Environmental DNA from Surface and Groundwater Samples

For DNA extraction, 250 ml of water sample was collected and filtered through a 47 mm polycarbonate filter (0.22 μm pore size, Millipore). Each filter paper was soaked in 10 ml of 70% methanol and stored at − 20 °C until use. Extraction of DNA from the filter paper was carried out according to the modified method described by Kim et al. (2012) and purification of DNA was subjected for PCR analysis.

2.4 Detection of ARGs by PCR

PCRs were performed to detect antibiotic resistance genes and the PCR mixtures contained 0.5 µl of target primer (10 µm), 5 µl Go taq reaction buffer, 0.5 µl dNTPs, 2.0 µl from 25 mM MgCl2, and 0.1 µl of Gotaq DNA polymerase, adjusted to a total volume of 25 µl. Purified DNA (5 µl) was used as the PCR template. Optimized conditions used for the primers were shown in Table 2. PCR amplification (initial denaturation at 95 °C for 2 min, denaturation at 95 °C for 30 s, annealing for 30 s, and extension at 72 °C for 1 min/35 cycles) was performed using the BIOLAB PCR system thermal cycler (BYQ6078E-757, China). The final extension was performed at 72 °C for 5 min. The amplified product was analyzed using electrophoresis.

Table 2 Primers, primer sequences, and annealing temperatures used to amplify selected resistant genes of the two antibiotics employed in the present study

2.5 Quantification of Antibiotic Resistance Genes

Quantitative PCR (qPCR) assays were employed to quantify five ARGs: bla TEM, bla OXA, amp a, tet A, and tet M, which give resistance to the penicillin and tetracycline. qPCR assays were performed using the SYBR Green (Applied Biosystems) using an Applied Biosystem 7500 real-time machine (Marti et al., 2013). Each gene was amplified using a specific primer set (Table 2) and the PCR conditions were optimized as initial denaturation at 95 °C for 3 min, followed by 40 cycles at 95 °C for 15 s and the annealing temperatures are given in Table 2 for the 20 s. A dissociation curve was constructed by increasing temperature from 65 to 95 °C to confirm the specificity of the amplified products.

2.6 Development of Standard Curves

Standard curves were generated by cloning the amplicon from positive controls into the PBR322 vector (Invitrogen, USA). The corresponding copy number was calculated using the following equation.

Copy number µl−1 = (A × 6.022 × 1023) (660 × B)−1.

where A is the plasmid DNA concentration (g µl−1), B is the plasmid length (bp) containing the cloned sequence, 6.022 × 1023 is the Avogadro’s number, and 660 is the average molecular weight of one base pair (Perini et al., 2011).

2.7 Water Quality Analysis

Water temperature, pH, and dissolved oxygen (DO) were measured using a portable multimeter (HACH–HQ 40D), and electrical conductivity (EC), total dissolved solids (TDS), and salinity were recorded using a portable conductivity meter (HACH–Sension EC5) at the site itself. Nitrates (as N-NO3), nitrite (as N-NO2), and total phosphate (TP) concentrations were measured by the Spectrophotometric (Spectro UV–VIS Double UVD 2960) methods (APHA 1999). Membrane Lactose Glucuronide Agar (MLGA) membrane filtration method was performed to determine the total coliform (TC) and E. coli count per 0.1 dm3 of the water samples (SLSI 2013).

2.8 Statistical Analysis

Comparison between average antibiotic concentration and ARGs concentrations among different sampling points were carried out using one-way ANOVA test. Differences were considered significant at p < 0.05. All statistical analysis was performed using Minitab 17 software. Pearson correlation tests were carried out for 10 physico-chemical and microbiological parameters for 40 sampling locations.

3 Results

3.1 Antibiotic Concentrations

In the present study, none of the penicillin and tetracycline group which targeted was detected at any sampling location except Kelani river mouth. However, the results of the study showed that the river mouth was contaminated with amoxicillin (AMX) in the penicillin group (0.003 ± 0.001 µg/ml).

3.2 Screening of Antibiotic Resistance Genes

Percentage of sampling sites from 5 to 15% were positive for penicillin resistance genes (bla TEM, bla OXA, OPR D, amp a, amp b) compare to the positive sampling sites percentage (~ 10%) against tetracycline resistance genes (tet A, tet M, tet S, tet B) in groundwater samples (Fig. 2). In groundwater samples, the highest number of positive sampling sites was recorded against bla TEM (15%), followed by descending order for bla OXA (10%), tet M (10%), amp a (5%), and amp b (5%), respectively (Fig. 2a). However, none of the groundwater sampling locations were positive for penicillin, OPR D, and tetracycline, tet A, tet B, and tet S.

Fig. 2
figure 2

Percentages of the a penicillin (bla TEM, bla OXA, OPR D, amp a, amp b) and b tetracycline (tet M, tet S, tet A, tet B) resistance gene in surface and groundwater sampling locations in the Kelani river basin

Similarly, resistance genes of tet A and tet M, belonging to the group tetracycline (OTC and TET), were detected as 5% and 15% sampling locations where penicillin resistance genes of bla TEM (30%), bla OXA (20%), amp a 10%), and amp b (10%) in surface water sampling locations, respectively (Fig. 2a and b). However, according to results of the present study tet B and tet S were not detected in any sampling location in the present study (Fig. 2b).

3.3 Quantification of Antibiotic Resistance Genes

Out of 5 genes that were initially studied during the pre-screening process of ground and surface water samples, 3 genes were selected from each resistance mechanism (changes in penicillin-binding protein, secretion of ß-lactamases, and decreasing of porin channel formation) based on the gene abundance in water samples. Two genes from 4 tetracycline resistance genes were selected for further analysis based on detection frequency and response mechanism. Subsequently, these selected genes were quantified using qPCR assays. All qPCR assays showed high R2 values (> 0.99) and high efficiencies from 96.84 to 107.71%.

The highest ARGs ranged from 1.572 × 102 to 700.576 × 102 copy/ml for penicillin was recorded in surface water samples, where from 0.372 × 102 to 312.700 × 102 copy/ml were detected in groundwater samples (Fig. 3a). Among the penicillin resistance genes, the bla TEM (700.576 × 102 copy/ml) was recorded as the highest concentration where the amp a was ranged from 0.372 × 102 to 371.67 × 102 copy/ml. OPR D gene was recorded as 1.572 × 102 copy/ml (Fig. 3a).

Fig. 3
figure 3

a Variations of total copy numbers of antibiotic resistance genes against ß-lactam resistance genes. b Variations of total copy numbers of antibiotic resistance genes against tetracycline resistance genes; when error bars are not shown, the standard deviation was less than the width of the symbol.

However, among the tetracycline-resistant genes, both tet M and tet A genes were detected in surface water samples ranging from 1.346 × 102 to 439.875 × 102 copy/ml, and only tet M was recorded as 215.987 × 102 copy/ml in groundwater (Fig. 3b).

3.4 Surface and Groundwater Quality

The surface water pH varied between 5.72 and 7.92 and the lowest pH was recorded from the Kaduwela sampling location where the EC was ranged between 52 and 2330 μS/cm and the highest EC value was recorded in the Hemilton canal. N-nitrate, N-nitrite, and TP concentration were ranged between 1.21 to 11.08 mg/L, 1.03 to 27.66 μg/L, and 55.0 to 191.4 μg/L, respectively. The groundwater pH varied between 4.21 and 7.30 and the lowest was recorded in the Dekatana sampling location. The EC of groundwater was ranged between 79.4 and 610.7 μS/cm and the highest EC was recorded in the Habarakada sampling point. N-nitrate, N-nitrite, and TP concentrations were ranged between 0.25 to 11.67 mg/L, < 1.00 to 29.33 μg/L, and 44.0 to 237.6 μg/L, respectively. Further, it was found that the entire lower part of the Kelani river basin (surface and groundwater) contaminated with TC and FC (98%). Pearson correlation showed a positive low correlation in temperature, TP, and FC; pH and FC; DO and nitrate where a negative low correlation was found in EC and DO.

4 Discussion

At present, greater attention has been paid to the heavy use of antibiotics in hospitals and the aquaculture industry by national and international health agencies due to the development of antibiotic resistance (Aminov and Mackie, 2007; WHO, 2018). The U.S. Centers for Disease Control and Prevention (CDC) has estimated due to antibiotic resistance approximately 23,000 deaths occur annually in the USA and estimated that the direct cost for treatment is about $20 billion (CDC, 2013). However, the value for economic losses in developing countries has not been calculated so far due to a lack of published information (Laxminarayan et al., 2013).

In Sri Lanka, ß-lactams and tetracycline are the most intensively used antibiotic classes for human therapies according to 2015, annual report of the Ministry of Health, Sri Lanka (Ministry of Health, 2015). However, the same antibiotics which are used for human infections are being prescribed for different diseases of farm animals and as growth hormone in the aquaculture industry (NAQDA, 2015). In Sri Lanka, information regarding antibiotic production and usage patterns makes nervous due to limited published data.

Liyanage and Manage (2016a) have reported contaminations of amoxicillin (water: 0.001–0.023 ppm, sediment: 0.001–0.005 ppm), tetracycline (water: 0–0.001 ppm, sediments: not detected), sulfardiaoxine (SDI) (water: 0.001–0.003 ppm, sediments: 0.001–0.003 ppm), sulfarmethaxasol (SMX) (water: 0.001–0.018 ppm, sediments: 0.001–0.002 ppm), and erythromycin (ERM) (water: 0.001–0.008 ppm, sediments: 0.001–0.003 ppm) in hospital effluent water in Sri Lanka. Further, Liyanage et al. (2015) documented that the concentration of TET in aquaculture farms was ranged between 0.001 ± 0.011 and 0.112 ± 0.017 ppm whereas OTC ranged from 0.008 ± 0.012 to 0.234 ± 0.014 ppm.

Many researchers have reported the presences of the potential of ARGs in river water samples all over the world (Pruden et al., 2006; Berendonk et al., 2015) including pristine environment (Pruden et al., 2013; Hsieh et al., 2011; Liyanage & Manage, 2016a). According to the authors’ knowledge, no records are available regarding ARGs (tet M, tet S, tet A, amp a, amp b, amp c, bla TEM, bla OXA) in the Kelani river basin which is the second-largest watershed in Sri Lanka. The river basin covers seven districts, namely, Colombo (446 km2, 19%), Gampaha (343 km2, 14.7%), Kaluthara (12 km2, 0.5%), Kegalle (1,028 km2, 44%), Rathnapura (70 km2, 3%), Nuwara Eliya (430 km2, 18.4%), and Kandy (4 km2, 0.2%). However, since the river flows through the urban centers, farmlands, and industrial zones, it is more vulnerable to water pollution (Mahagamage & Manage, 2014; Mahagamage et al., 2016) and it is the most polluted river in Sri Lanka due to the rapid growth of industries located in close vicinity of the river and passes through the most populated capital city of the country (Mahagamage et al., 2014; CEA, 2015). Thus, the present study is the first report on the screening of the potential of ARGs with special emphasis on penicillin and tetracycline in surface and groundwater samples collected from the meandering zone of the Kelani river basin, Sri Lanka.

The penicillin resistance genes (bla TEM, bla OXA, OPR D, amp a, amp b) were detected in a high percentage (5–30%) comparable to the resistance genes recorded for tetracycline (tet A, tet M, tet S, tet B) in surface water. Similarly, a high number of locations were positive for penicillin resistance genes (5–15%) compared to the resistance genes of tet A, tet M, tet B, and tet S belonging to the group tetracycline (OTC and TET) (~ 10%) in groundwater. Thus, results were suggesting usage of specific antibiotics in particular aspects causes ARGs development in the environment.

The results of the study clearly showed that ARGs are present in various environmental compartments. Moreover, the presence of antibiotics in a sub-inhibitory concentration may induce horizontal gene transfer within bacterial communities (Shah et al., 2008), which further increases the prevalence of ARGs.

In this study, a potent rapid molecular technique, real-time PCR, was applied to quantify the concentrations of ARGs (tet M, tet A, amp a, OPR D, bla TEM) in ground and surface water samples. The method proved to be very efficient, reproducible, and sensitive for the quantification of selected resistance genes in water samples with high efficiency values (96.84–107.71%).

Likewise, Mao et al., (2015) suggested that water in river water and river mouth is a significant point source of several ARGs into the receiving environment. In the case of penicillin and tetracycline resistance genes, concentration in water samples of the Kelani river may become severe due to detection of the respective gene without detection of any selected penicillin (AMX, ampicillin, cloxacillin) or tetracycline (TET, OTC) in water samples. Thus, it will be confirmed that the spreading of ARGs and ARB was not totally dependent on antibiotic concentration. Altogether, these observations undoubtedly demonstrate the contribution of effluent water discharges to the spread of antibiotic resistance in the natural environment.

Further, data obtained from the present study was used to analyze the correlation between the absolute concentrations of selected genes vs the antibiotics. However, no significant correlation (p > 0.05) was found between antibiotic concentration in the environment and antibiotic resistance genes. Further, it was found that the resistance genes were detected with the absence of tetracycline and penicillin in water. Thus, the results of the study suggest that the increase in the prevalence of resistance genes is caused by the persistence of these genes even the absence of selection pressure (Kim et al., 2012). However, the results presented in this study are consistent with previous studies suggesting that exposure to antibiotics could be the main factor that leads to selective pressure for ARGs (Franje et al., 2010).

The highest EC values were recorded in Mattakkuliya river mouth (1040 μS/cm) and the Hemilton canal (2330 μS/cm) which exceeded the standard given by SLSI for drinking water quality (750 μS/cm). All other water quality parameters remained the SLSI values given for drinking water except microbiological parameters (TC and FC). It was found that all the surface and groundwater samples were contaminated with TC and FC bacteria and exceeded the values given by SLSI for drinking water standards.

Thus, based on the results obtained from the present study and the available literature, the authors suggest that further studies are needed to update which concentrations of antibiotic resistance are developed, disseminated, and variations of that concentration with a different type of environment in the future.

5 Conclusion

The present study is the first study on antibiotics and ARGs in the meandering zone of the Kelani river basin in Sri Lanka. The results of the study revealed that surface and groundwater are possible reservoirs of ARGs. Further, pathogen-associated taxonomic groups in the Kelani river had implications for human health is suggested. The study highlighted that the great impact of heavy and unregulated use of antibiotics along with the discharge of partially or untreated wastewater into the aquatic environment might cause significant contamination of both antibiotics and ARGs. Furthermore, development of antibiotic resistance create huge problem in treating bacterial infections in near future. Thus, great concern should be given to the amount of antibiotic resistance in the aquatic environment that can serve as a harbor of resistance for the pathogens as ARGs and ARB may be creating a potential risk to human and animal health.