The impact of physico-chemical water quality parameters on bacterial diversity in the Vaal River, South Africa

This study aimed to identify bacterial community structures in the Vaal River using PCR-DGGE (polymerase chain reaction denaturing gradient gel electrophoresis) and high-throughput sequencing. The impact of physico-chemical characteristics on bacterial structures was investigated through multivariate analysis. Samples were collected from 4 sampling stations along the Upper Vaal River during winter (June 2009) and summer (December 2010). Physico-chemical analysis was conducted on-site. Additional physico-chemical data were obtained from statutory bodies. DNA was directly isolated from water samples and PCR amplified using universal bacterial primer pairs. PCR products were subjected to DGGE fingerprinting and high-throughput sequencing, followed by Shannon-Weaver diversity calculations, cluster analysis and multivariate analysis. Physico-chemical parameters did not exceed the prescribed South African water quality standards for domestic use, aquatic ecosystems, livestock watering and irrigation. DGGE banding patterns revealed similar bacterial community structures for 3 of the 4 sampling stations. PCA and RDA indicated that pH, water temperature and inorganic nutrient concentrations could be used to explain changes in bacterial community structures. High-throughput sequencing data showed that bacterial assemblages were dominated by common freshwater groups: Cyanobacteria, Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Bacteroidetes and Actinobacteria. Other freshwater phyla such as Deltaproteobacteria, Epsilonbacteria, Acidobacteria, Verrucomicrobia, Firmicutes, Fusobacteria, Flavobacteria and Fibrobacteres were found in low proportions. This study provides an overview of the dominant bacterial groups in the Upper Vaal River and the impact of environmental changes on bacterial diversity.


INTRODUCTION
Socio-economic growth and development of the Vaal River require continuous augmentation of this water resource to meet the growing water requirements of communities in Gauteng, the Free State, North West and Northern Cape provinces (DWAF, 2009b).Water quality has drastically deteriorated due to constant disposal of industrial and domestic waste into the river.Salinisation, eutrophication and microbiological pollution are currently the main problems affecting the water quality (DWAF, 2009a).The Department of Water Affairs and Forestry (DWAF) of South Africa, in line with the South African National Water Act (NWA), Act No. 36 of 1998, stipulated regulatory guidelines and criteria a water system must meet to ensure that the country's water resources are fit for use.A structured biomonitoring programme was implemented by the DWAF in 2009 to determine the exact sensitivity and health status of the Vaal River (DWAF, 2009a).Criteria routinely monitored to ensure sustainability, optimal water use and protection of the water resource include: physico-chemical characteristics, stream flow, discharge loads and microbiological pollutants, in particular, Escherichia coli (DWAF, 2009a;2009b).The detection of Escherichia coli only indicates the presence of faecal contamination and not necessarily the degree of industrial pollution.Therefore, in-depth studies on the microbial communities in the Vaal River are essential to understand the microbial processes underlying secondary pollution and changes in the physico-chemical quality of water.
DGGE has been applied in numerous research studies involving the assessment of microbial diversity of rivers, streams, lakes and sediment, to determine the water quality of the resource (De Figueiredo et al., 2010;Essahale et al., 2010;De Figueiredo et al., 2011;Haller et al., 2011).This method opened up new avenues of research on the diversity of microorganisms present in complex aquatic environments.Currently, metagenomic analysis of microbial ecology, such as highthroughput sequencing (HTS), has been the focus of several environmental studies such as soil, (Lemos et al., 2011), freshwater lakes (Marshall et al., 2008) and deep sea microbiota (Sogin et al., 2006).Metagenomic analysis provides extensive information on community structure and composition (Kakirde et al., 2010).In addition, phylogenetic and functional analyses of microorganisms can be determined at community level (Cowan et al., 2005).
The objectives of this study were (i) to identify the bacterial community structures in the planktonic phase of the Vaal River using 16S rDNA PCR-DGGE and high-throughput sequencing, and (ii) determine the impact of physico-chemical characteristics on bacterial community structures using principle component analysis (PCA) and redundancy analysis (RDA).Samples were collected from the planktonic phase in sterile glass bottles and preserved on ice for not longer than 6 h prior to nucleic acid isolation.Physico-chemical analysis was conducted in situ.Additional physico-chemical data were obtained from the Department of Water Affairs (2012) and the South African Weather Service (2012).A summary of the physicochemical variables of all studied sampling sites is shown in Table 1.

Nucleic acid isolation
A hundred millilitres of water samples were filtered through a 0.45 μm nitrate cellulose membrane filter (Whatman, Missouri, USA) and subsequently lysed in a 1 mg/mℓ lysozyme solution  387 that contained 0.25-0.50mm glass beads (Sigma-Aldrich Co., Missouri, USA) for bacterial cell disruption.The lysis solution was incubated at 37°C for 10 min while agitated in a vortex.Proteinase K (1 mg/mℓ) was then added and the lysis solution was incubated at 56°C for an additional 30 min.DNA was isolated from the crude lysate using the PeqGold Bacterial DNA Kit (PEQLAB Biotechnologie GmbH, Erlangen, Germany).The quality and quantity of the isolated nucleic acids were determined using the Nanodrop ND1000 (NanoDrop Technologies, Delaware, USA) and agarose electrophoresis.

PCR amplification and DGGE analysis of bacterial community structures
The highly variable V3 region of the 16S rDNA gene fragments were PCR amplified using the universal primer pair 341F-GC and 907R (~ 500 bp) (Muyzer et al., 1993) PCR products were analysed by DGGE using a DCode Universal Detection System (Bio-Rad Laboratories, Hercules, California, USA).Four reference species, namely Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus and Streptococcus faecalis, were included in all DGGE studies.DGGE analysis was conducted at a denaturing gradient of 30-50% in 1 mm vertical polyacrylamide gels (8% (wt/vol) acrylamide in 1 × TAE).20 µℓ of amplification product were mixed with 5 µℓ of loading buffer (6× Orange Loading Dye, Fermentas Life Sciences, Maryland, USA) and loaded into the gel.Electrophoresis was performed at a constant temperature of 60°C for 16 h at 100 V in 1 × TAE buffer (40 mM Tris-acetate, 1 mM EDTA, pH 8.0).Polyacrylamide gels were stained with ethidium bromide (10 mg/ℓ) for 45 min and visualised with a Gene Genius Bio Imaging System (Syngene, Cambridge, UK) and GeneSnap software (version 6.00.22).None of the DGGE gels were digitally enhanced or modified.Bands of interest were only highlighted for better visualisation and not analytical purposes.Selected DNA bands of interest were excised from gels with a sterile scalpel and eluted in 20 μℓ of sterile nucleasefree water for 12 h at 4°C. 2 µℓ of the elute were used as DNA template in PCR amplification reactions with primer pair 341F and 907R (Muyzer et al., 1993) and conditions described above.PCR products were subsequently purified and sequenced using a BigDye ® Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, California, USA) and Genetic Analyzer 3130 (Applied Biosystems, California, USA).Sequences were aligned to 16S rRNA sequences in the National Center of Biotechnology Information Database (NCBI) using BLASTN searches to determine their identity.A total of 23 bacterial nucleotide sequences were submitted to the GenBank database under accession numbers JQ085826 -JQ085849.

High-throughput sequencing
HTS analysis was performed by Inqaba Biotech, South Africa, using the Roche 454 GS-FLXTM System.The V1-V3 region of the 16S rRNA gene was amplified using primer pair 27F (GAGTTTGATCCTGGCTCAG) and 518R (ATTACCGCGGCTGCTGG) to produce ~ 500 bp fragments.Subsequently, sequences were trimmed to remove GS tags and further analysed with the CLC Bio Genomics Workbench version 4.7.2 software (CLC Bio, Aarhus, Denmark).Sequences shorter than 200 bp in length were excluded from data sets.All remaining sequences were subjected to the National Center for Biotechnology Information (NCBI) database for BLAST analysis.Sequences were then submitted to Pintail version 1.0 to detect the presence of PCR artefacts.PCR products with chimeric properties were eliminated from data sets prior to phylogenetic analysis.The remaining 922 sequences were submitted to GenBank with accession numbers JN865256-JN866178.

Statistical analysis
Bacterial community diversity was calculated with the Shannon-Weaver diversity index (H'), based on DGGE profiles.The Shannon-Weaver indices (H') were calculated according to Zhang et al. (2011).Similarities between the banding patterns generated by PCR-DGGE of the various sampling sites were compared by cluster analysis as indicated by Gafan et al. (2005).Cluster analyses were displayed graphically as UPGMA dendrograms.
The distribution of samples according to environmental factors was analysed by PCA.The statistical significance of the relationships between bacterial community structures, DGGE banding profiles, high-throughput sequencing data and water quality was further assessed by RDA.Environmental variables selected are summarised in Table 2. Multivariate analysis was performed by a Monte Carlo permutations test using unlimited permutations.Analysis was carried out using CANOCO software version 4.5.

Physico-chemical characteristics
Selected physico-chemical parameters measured or obtained are listed in Table 1.These parameters showed all physicochemical values to fall within the prescribed South African water quality guidelines for domestic use (DWAF, 1996a), aquatic ecosystems (DWAF, 1996b), livestock watering (DWAF, 1996c), irrigation (DWAF, 1996d) and aquaculture (DWAF, 1996e) (Table 2).Water temperatures were between 10 and 13°C in June and December temperatures exceeded 20°C (24.4-28.7°C).The temperatures of inland aquatic ecosystems in South Africa generally range between 5 and 30°C but can fluctuate depending on the geographical features of the region and catchment area, seasonal changes and the impact of anthropogenic activities (DWAF, 1996b).In December, the flow velocity increased sequentially from Deneysville to downstream sampling stations (Scandinawieë Drift).This trend was not observed in June, when rainfall was low.

DGGE analysis
In this study, PCR-DGGE was able to give spatial information about the dominant bacterial communities in the Vaal River system (Fig. 2).Previous studies suggest that band intensity is related to the relative abundance of the corresponding phylotypes in the sample mixture (Murray et al., 1996;Riemann et al., 1999).Thus, bands with relatively high intensities were assumed to be dominant taxa.
DGGE profiles demonstrated high resolution and intensity at a denaturing gradient of 30-50%.Four bacterial species, Escherichia coli, Pseudomonas aeruginosa, Streptococcus faecalis and Staphylococcus aureus, were included in all DGGE studies, to determine the potential of using such an approach to establish the presence of these species in water samples.Corresponding bands for Staphylococcus aureus and Pseudomonas aeruginosa were detected for Vaal Barrage, Parys and Scandinawieë Drift.In addition, Parys illustrated a band with similar migration patterns to Escherichia coli.All corresponding bands were excised and sequenced but produced poor quality sequences with indefinite identification.Since sequence data could not confirm accurate identification of excised bands, results remain inconclusive.
Vaal Barrage, Parys and Scandinawieë Drift displayed similar DGGE patterns for the dominant bands in June and December (Fig. 2).However, DGGE profiles for Deneysville varied to some extent from the three other sites.Although some dominant bands showed similar migration patterns to Vaal   On the other hand, bacterial diversity for Deneysville was higher in December than in June.The Shannon-Weaver indices (Fig. 4), however, contradicted the DGGE diversity data.These showed a higher bacterial diversity for Vaal Barrage and Scandinawieë Drift during December compared to June.The Shannon-Weaver index calculation includes the presence and absence of bands, but also band intensity that could be used to explain the contradiction (Zhang et al., 2011).A total of 24 bacterial bands were excised, sequenced and compared to sequences in the NCBI database (Table 3).Approximately 75% of the bacterial sequences recovered displayed high sequence homologies (> 97%) with the known database sequences.However, 50% of these sequences showed the highest sequence similarity to uncultured bacteria obtained directly from freshwater samples.These results support the presence of many uncultured and potentially undescribed bacterial taxa in freshwater ecosystems.Taxonomic classifications of the partial 16S rDNA sequences obtained affiliated to Cyanobacteria (B4, B13-B15, B17, B23), Bacteroidetes (B6, B11, B22), Betaproteobacteria (B2, B12, B24) and uncultured bacteria (B1, B3, B5, B7-B10, B16, B18-B21).Bacterial communities for June displayed relative abundances of 8%, 17%, 17% and 58% for Cyanobacteria, Bacteroidetes, Betaproteobacteria and uncultured bacteria, respectively.In contrast, the relative abundance for Cyanobacteria increased to 42% in December, whereas Bacteroidetes, Betaproteobacteria, and uncultured bacteria respectively accounted for 8%, 8% and 42% of the four main phylogenetic groups.
While identification of the four indicator organisms employed in DGGE profiling remained inconclusive by Sanger sequencing, HTS analysis verified that two of the bands did in fact belong to the Pseudomonadaceae family and Esche richia spp.Additional opportunistic pathogens detected in low quantities at Vaal Barrage, Parys and Scandinawieë Drift included Roseomonas sp., Ralstonia sp., Serratia sp. and Stenotrophomonas sp.

Distribution of bacterial diversity in the Vaal River
The Shannon-Weaver diversity indices (H') were calculated from DGGE banding patterns as the number and relative intensity of bands (Fig. 4).Indices were used to compare the overall structure of bacterial communities among the four sampling sites.H' for June and December samples ranged from 0.27-0.46and 0.70-0.86,respectively.Bacterial diversity gradually increased from upstream to downstream sites, except for Parys in December which consisted of a lower diversity.Similar trends were also observed for HTS data.Cluster analysis was performed to gain an overview of the association of bacterial communities at the four sampling stations during June and December (Fig. 5).UPMGA dendrograms showed grouping of samples according to season.June samples showed high similarity (> 94%) among bacterial communities for Vaal Barrage, Parys and Scandinawieë Drift.A similar trend was observed for the December samples where Vaal Barrage and Scandinawieë Drift were defined by a 100% similarity.Noticeable was the grouping of the December Parys and Deneysville samples (100% similarity).Grouping of these two sampling sites may be attributed to similar banding patterns of a few dominant DGGE bands (Fig. 2).Diversity indices (H') and cluster analyses could be associated with DGGE profiles which reflected variations in the distribution, abundance and composition of bacterial taxa.The effect of different sampling periods is illustrated by the PCA analysis results (Fig. 6A).The June samples, with negative  and positive scores along the first axis, are separated from the December samples, which showed a positive score along the second axis.The first axis was mainly defined by ammonium, nitrate, phosphate, chloride, sulphate, TDS, conductivity and rainfall.The second axis was related to temperature, day length and flow rate.
RDA plots for high-throughput sequencing data (Fig. 6D) showed (i) positive correlations between the flow rate and abundances of Gammaproteobacteria, Deltaproteobacteria and Fibrobacteres along the first axis, (ii) positive correlations between rainfall, TDS, nitrate, ammonium, chloride and sulphate concentrations, and abundances of Acidobacteria and Actinobacteria along the second axis, and (iii) positive correlations between ammonium, chloride and phosphate concentrations, and abundances of Fusobacteria, Verrucomicrobia and Euglenoida along the second axis.Betaproteobacteria negatively correlated with Gammaproteobacteria.A high abundance of Betaproteobacteria was detected in June but decreased considerably in December.An opposite inclination was observed for Gammaproteobacteria.

Microbial community dynamics
Knowledge and insight into the diversity and function of freshwater microorganisms is an essential requirement for the sustainable management of freshwater resources.In addition, changes in bacterial community structures might be used as potential bio-indicators of environmental disturbances.The aim of this study was to examine bacterial community structures in a segment of the Vaal River, in response to environmental parameters, using a PCR-DGGE and high-throughput sequencing approach.High-throughput sequencing provided an overview of the dominant bacterial communities in the planktonic phase and marked shifts in composition, as attested to by PCA and RDA.
The composition of bacterial communities in a given environment depends on the interaction between various factors such as the geographic environment (Zhang et al., 2011), temperature (Hall et al., 2008), pH (Yannarell and Triplett, 2005), flow rate (Crump and Hobbie, 2005), light intensity (Sigee, 2005) and nutrient concentrations (Pomeroy and Wiebe, 2001).In this study of a segment of the Vaal River, the physicochemical parameters varied with sampling station and season of sampling.PCA and RDA analysis indicated that bacterial community structures were mainly influenced by pH, temperature and inorganic components.
The bacterial community structures were similar for the three sampling sites during each sampling period.However, the June bacterial community structures were different to the December assemblages.DGGE results suggested that bacterial diversity was higher during June compared to December.These results were, however, contradicted by the Shannon-Weaver indices.The latter analysis included presence-absence, as well as (abundance) band intensity data.This could be used to explain the contradiction (Zhang et al., 2011).Diversity index analysis of the high-throughput sequencing data showed similar trends to the Shannon-Weaver analysis of DGGE profiles.
Bacterial community structures could be correlated to inorganic nutrients as shown by PCA and RDA.The Vaal Barrage creates a buffering action that encapsulates organic and inorganic particles in the water-column for several weeks.This creates a relatively stable environment in which organisms can develop into a community.The planktonic bacteria then flow from here downstream to Parys and Scandinawieë Drift.Therefore, bacterial communities along this section of the Vaal River will be relatively similar.In addition, the dominant bacterial groups detected at these three points may be native species with broader niche capabilities, which allow them to grow and survive under a variety of environmental conditions (Anderson-Glenna et al., 2008).Recurrent native bacterial communities in aquatic ecosystems have been reported previously (Sekiguchi et al., 2002;Crump et al., 2003).It should be noted that the DNA amplification method used in this study did not discriminate between DNA derived from living cells versus DNA from dead cells and/or even naked or free DNA available in the water column.This aspect should be considered in future aquatic studies.
A feature highlighted in the present study was the relatively low bacterial diversity detected at Deneysville in June and December.Bacterial community structures at this sampling station largely consisted of Cyanobacteria, particularly Cyanophyta (Anabaena sp.), where pH and temperature were the main factors that affected the community structures.An alkaline pH was measured in June and December when temperatures in December were above 25°C.Optimum growth of Cyanophyta and the formation of surface algal blooms are the direct result of high nutrient concentrations (particularly phosphate) and physico-chemical characteristics (high pH, temperature and light intensity) (Sigee, 2005).In addition to these conditions, buoyancy also plays an important role in the development of Cyanophyta populations.Buoyancy allows algal populations to adopt an optimum position within the water column in relation to light and CO 2 availability (Sigee, 2005).This mechanism leads to changes in the water chemistry and light regime in the epilimnion that depress the growth of other phyto-and bacterioplankton groups (Sigee, 2005).
Although flow rate in this study was not shown to affect bacterial communities, previous studies have suggested that flow rate and hydraulic retention time have a substantial effect on community structures (Lindström and Bergström, 2004;Crump and Hobbie, 2005).Temporal variation in bacterial diversity was observed between the June and December samples.The Gauteng and North West Provinces received heavy rainfall in December 2010 that caused a drastic increase in flow rate, particularly at Parys and Scandinawieë Drift.The high flow rate resulted in flooding at these two sampling stations that likely changed the bacterial community structures.Bacterial communities in rivers with short hydraulic retention times would potentially remain undetected by DGGE due to high loss rates (wash-out effect) which in turn result in a lower bacterial density and diversity (Sommaruga and Casamayor, 2009).In contrast, rivers with an extended hydraulic retention time display an accumulation of nutrients which promotes a higher genetic diversity of bacteria.Although flow rate differences provide a reasonable explanation for the seasonal variation in bacterial, further investigations are needed to confirm this for the Vaal River.

Phylogenetic diversity of bacterial communities
Phylogenetic affiliation of the dominant groups retrieved from the freshwater samples by PCR-DGGE and high-throughput sequencing corresponded to Cyanobacteria, Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Bacteroidetes and Actinobacteria.Other freshwater phyla such as Deltaproteobacteria, Epsilonbacteria, Acidobacteria, Verrucomicrobia, Firmicutes, Fusobacteria, Flavobacteria and Fibrobacteres were found in low proportions.
Cyanobacteria accounted for a large proportion of bacterial diversity during December, which agrees well with the physico-chemical characteristics of the water samples.Several studies have indicated that Cyanobacteria tend to dominate phytoplankton communities in pristine freshwater systems (Anderson-Glenna et al., 2008;Foong et al., 2010) whereas other authors have reported an increase in the prevalence of Cyanobacteria in response to fluvial, organic and urban wastewater pollution (Douterelo et al., 2004;Ibekwe et al., 2012).Due to the trophic status of the Vaal River, cyanobacterial blooms usually occur during late spring and summer and often consist of Microcystis aeruginosa, Oscillatoria sp. and Anabaena floss-aqua (Cloot and Le Roux, 1997;DWAF, 2009a).In this study, Anabaena sp., Cymbella helvetica and Synechocystis sp. were in high abundance at Deneysville during December 2010.Anabaena spp.are among the most distributed toxin producers in eutrophic freshwater bodies (Berg et al., 1986).Their potential effects on aquatic ecosystems may be subtle or can cause major changes in the survival of sensitive species (DWAF, 2009a).In addition, these toxins may pose a serious health hazard for human and animal consumption.
Alphaproteobacteria, Betaproteobacteria, Gamma proteobacteria and Actinobacteria are ubiquitous groups in freshwater habitats (Gich et al., 2005;Anderson-Glenna et al., 2008) and are numerically important in river systems (Beier et al., 2008;Lemke et al., 2009).Members of Betaproteobacteria respond rapidly to organic and inorganic nutrient enrichment (Hahn, 2003;Simek et al., 2005) and have been isolated from various polluted and unpolluted freshwater bodies (De Figueiredo et al., 2011;Haller et al., 2011).Two important genera of this subphylum include Dechlorosomonas and Variovorax.Members of Dechlorosomonas are capable of oxidising aromatic compounds such as benzoate, chlorobenzoate and toluene (Coates et al., 2001), whereas Variovorax spp.are involved in plant growth and remediation of xenobiotics (Jamieson et al., 2009).Several opportunistic human pathogens of the Gammaproteobacteria group were detected at low abundance.Human diseases and infections are often associated with these pathogens (Berg et al., 2005;Mahlen, 2011) and have caused mortalities in immunocompromised individuals (Fergie et al., 1994;Paez and Costa 2008).Thus, although the opportunistic pathogens were present at low levels, their impact should not be underestimated.
RDA analysis revealed that nitrate, ammonium, chloride and sulphate were the four most influential inorganic factors responsible for shaping Actinobacterial and Acidobacterial communities.A few studies suggested that these two phyla participate in the nitrogen cycle in soils and sediments by reducing nitrate, nitrite and possibly nitric oxide (Gtari et al., 2007;Ward et al., 2009).Norris et al. (2011) also implicated some novel Actinobacteria from geothermal environments in growing autotrophically with sulphur as an energy source.Correlation between Verrucomicrobia and phosphate was also detected suggesting that this inorganic nutrient influenced the Verrucomicrobia community within the total bacterial population.The association between Verrucomicrobia and phosphate levels has seldom been discussed in previous studies of microbial ecology of freshwater resources (Lindström et al., 2005;Liu et al., 2009).Very little is known about the physiology and ecological roles of Actinobacteria, Acidobacteria and Verrucomicrobia in these habitats and the impact of physico-chemical characteristics on their community composition.
Members of Bacteroidetes usually inhabit mesotrophic and eutrophic water bodies that have high nutrient levels (Xi et al., 2007;de Figueiredo et al., 2011).This group is known to degrade polymeric organic matter, and to play an important role in the turnover of organic matter (Cottrell and Kirchman, 2000), and is often isolated from humic waters (Anderson-Glenna et al., 2008;Stabili and Cavallo, 2011).The Bacteroidetes-Flavobacterium-like lineages are often present in high abundance following the growth and decline of cyanobacterial blooms (Eiler and Bertilsson, 2007;Newton et al., 2011).Their presence and distribution is mainly determined by resource availability and is favoured during periods of high heterotrophic activity and enhanced growth (Eiler and Bertilsson, 2007).This phenomenon was evident in the high abundance of Bacteroidetes in June following the December 2008 to February 2009 cyanobacterial blooms.

CONCLUSIONS
This study investigated the impact of physico-chemical water quality parameters on bacterial community structures in a segment of the Vaal River.The PCR-DGGE approach and high-throughput sequencing analysis presented useful data in the identification of dominant bacterial groups at the four sampling stations.Molecular analysis showed that (i) bacterial community structures for June were different to the December assemblages, (ii) bacterial community structures for Vaal Barrage, Parys and Scandinawieë Drift were similar, (iii) bacterial communities at Deneysville differed from the three other sites and were lower in diversity, and (iv) Cyanobacteria, Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, Bacteroidetes and Actinobacteria were the dominant bacterial groups detected and were shown to be impacted by physico-chemical water quality parameters.This study contributed to the identification of bacterial phylotypes, their spatial succession and the effect of physico-chemical characteristics on these freshwater bacterial communities.A detailed study on the relationships between the dominant bacterial taxa and specific physico-chemical water characteristics is required to improve our knowledge on how bacterial community structures in the Vaal River are affected.

Figure 1
Figure 1Geographical illustration of the Vaal River system.The four sampling stations are indicated on the map.

a
Figure 2 DGGE bacterial community analyses for 16S rDNA gene fragments from surface water during June 2009 and December 2010.Sampling sites selected along the Vaal River include Deneysville (D), Parys (P), Scandinawieë Drift (SD) and Barrage (B).Four indicator species were used as references: E.coli (E.c), Pseudomonas aeruginosa (P.a), Streptococcus faecalis (S.f) and Staphylococcus aureus (S.a).The DNA present in numbered bands was sequenced; identities are summarised in Table3.None of the DGGE gels were digitally enhanced or modified.Bands of interest were only highlighted for better visualisation and not for analytical purposes.

Figure 4
Figure 4 Shannon-Weaver diversity indices (H') for the Vaal River in June 2009 and December 2010 at Deneysville, Barrage, Parys, and Scandinawieë Drift.

Figure 5
Figure 5 Cluster analysis of DGGE band patterns obtained in June 2009 and December 2010 using Pearson correlation coefficient.DGGE profiles are graphically demonstrated as UPGMA dendrograms.
Figure 6 (A) PCA analysis of physico-chemical and microbial variables in the first-and second-axis ordination plots, (B) RDA tri-plot of DGGE bands (samples indicated using band [BN] numbers) and environmental variables (represented by arrows) in June 2009, (C) RDA tri-plot of DGGE bands (samples indicated using band [BN] numbers) and environmental variables (represented by arrows) in December 2010 and (D) RDA tri-plot of bacterial phyla and environmental variables (represented by arrows).

Table 3
. None of the DGGE gels were digitally enhanced or modified.Bands of interest were only highlighted for better visualisation and not for analytical purposes.http://dx.doi.org/10.4314/wsa.v39i3.7 Available on website http://www.wrc.org.zaISSN 0378-4738 (Print) = Water SA Vol.39 No 3 WISA 2012 Special Edition 2013 ISSN 1816-7950 (On-line) = Water SA Vol.39 No 3 WISA 2012 Special Edition 2013 389 Barrage, Parys and Scandinawieë Drift, a few distinct bands exhibited unique migration positions.A higher bacterial diversity, based on number of bands, was detected for Vaal Barrage and Scandinawieë Drift during June compared to December.