Identification and application of bile metabolites to assess the exposure of fish to pharmaceuticals in the environment

Thousands of tons of pharmaceuticals are consumed yearly worldwide. Due to the continuous and increasing consumption and their incomplete elimination in wastewater treatment plants (WWTP), pharmaceuticals and their metabolites can be detected in receiving waters, although at low concentrations (ng to low μg/L). As bioactive molecules the presence of pharmaceuticals in the aquatic environment must be considered potentially hazardous for the aquatic organisms. In this thesis, the biotransformation and excretion of pharmaceuticals in fish was studied. The main biotransformation pathways of three anti‐inflammatory drugs, diclofenac, naproxen and ibuprofen, in rainbow trout were glucuronidation and taurine conjugation of the parent compounds and their phase I metabolites. The same metabolites were present in fish bile in aquatic exposures as in fish dosed with intraperitoneal injection. Higher bioconcentration factor in bile (BCFbile) was found for ibuprofen when compared to diclofenac and naproxen. Laboratory exposure studies were followed by a study of uptake of pharmaceuticals in a wild fish population living in lake contaminated with WWTP effluents. Of the analyzed 17 pharmaceuticals and six phase I metabolites, only diclofenac, naproxen and ibuprofen was present in bream and roach bile. It was shown, that diclofenac, naproxen and ibuprofen excreted by the liver can be found in rainbow trout and in two native fish species living in the receiving waters. In the bream and roach bile, the concentrations of diclofenac, naproxen and ibuprofen were roughly 1000 times higher than those found in the lake water, while in the laboratory exposures, the bioconcentration of the compounds and their metabolites in rainbow trout bile were at the same level as in wild fish or an order of magnitude higher. Thus, the parent compounds and their metabolites in fish bile can be used as a reliable biomarker to monitor the exposure of fish to environmental pharmaceuticals present in water receiving discharges from WWTPs.

Persistent organic pollutants in sediments and fish from Lake Victoria, Uganda

This thesis describes the occurrence and sources of selected persistent organic pollutants (POPs) such as polychlorinated dibenzo-p-dioxins (PCDDs), polychlorinated dibenzofurans (PCDFs), polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs) and hexachlorocyclohexanes (HCHs) in the northern watershed of Lake Victoria. Sediments and fish were collected from three highly polluted embayments (i.e. Murchison Bay, Napoleon Gulf and Thurston Bay) of the lake. The analysis for PCDD/Fs, PCBs and PBDEs was done using a high resolution mass spectrometer coupled to a gas chromatograph (GC), and a GC equipped with an electron capture detector was used for HCHs. Total (Σ) PCDD/Fs, PCBs and PBDEs in sediments ranged from 3.19 to 478, 313 to 4325 and 60.8 to 179 pg g-1 dry weight (dw), respectively. The highest concentrations of pollutants were found at sites close to industrial areas and wastewater discharge points. The maximum concentrations of PCDD/Fs, PCBs, PBDEs and HCHs in fish muscle homogenates were 49, 779, 495 and 45,900 pg g-1 wet weight (ww), respectively. The concentrations of the pollutants in Nile perch (Lates niloticus) were significantly greater than those in Nile tilapia (Oreochromis niloticus), possibly due to differences in trophic level and dietary feeding habits among fish species. World Health Organization-toxic equivalency quotient (WHO2005-TEQ) values in the sediments were up to 4.24 pg g-1 dw for PCDD/Fs and 0.55 pg TEQ g-1 dw for the 12 dioxin-like PCBs (dl-PCBs). 23.1% of the samples from the Napoleon Gulf were above the interim sediment quality guideline value of 0.85 pg WHO-TEQ g-1 dw set by the Canadian Council for Ministers of the Environment. The WHO2005-TEQs in fish were 0.001-0.16 pg g-1 for PCDD/Fs and 0.001-0.31 pg g-1 ww for dl- PCBs. The TEQ values were within a permissible level of 3.5 pg g−1 ww recommended by the European Commission. Based on the Commission set TEQs and minimum risk level criteria formulated by the Agency for Toxic Substances and Disease Registry, the consumption of fish from Lake Victoria gives no indication of health risks associated to PCDD/Fs and PCBs. Principal component analysis (PCA) indicated that anthropogenic activities such as agricultural straw open burning, medical waste incinerators and municipal solid waste combustors were the major sources of PCDD/Fs in the watershed of Lake Victoria. The ratios of α-/γ-HCH varied from 0.89 to 1.68 suggesting that the highest HCH residues mainly came from earlier usage and fresh γ-HCH (lindane). In the present study, the concentration of POPs in fish were not significantly related to those in sediments, and the biota sediment accumulation factor (BSAF) concept was found to be a poor predictor of the bioavailability and bioaccumulation of environmental pollutants.

Diel patterns and tissue-specificity of environmental responses in fish

Humans are profoundly changing aquatic environments through climate change and the release of nutrients and chemicals. To understand the effects of these changes on natural populations, knowledge on individuals’ environmental responses is needed. At the molecular level, the environmental responses are partly mediated by chances in messenger RNA and protein levels. In this thesis I study messenger RNA and protein responses to an assortment of environmental stressors in fish. As daily (diel) rhythms are known to be ubiquitous in different tissues, I particularly focus on diel patterns in the responses. The studied species are the three-spined stickleback (Gasterosteus aculeatus L.) and the Arctic char (Salvelinus alpinus L.), both of which have circumpolar distribution in the Northern hemisphere. In the first two studies, three-spined sticklebacks were exposed to both the non-steroidal anti-inflammatory drug diclofenac and low-oxygen conditions (hypoxia), and their responses measured at separate time points in the liver and gills. The results show how the seemingly unrelated environmental stressors, hypoxia and anti-inflammatory drugs, can have harmful combined effects that differ from the effects of each stressor alone. Moreover, both stressors disturbed natural diel patterns in gene expression. In the third study, I studied the responses of three-spined sticklebacks to two test chemicals: one used in hormonal medicine (17α-ethinyl-oestradiol) and one used as a plasticizer and solvent chemical (di-n-butyl phthalate). The results suggest that the phthalate can affect genes related to spermatogenesis in fish testes, while estrogen-mimicking compounds can lead to numerous disturbances in the endocrine system. In the final study, the temperature-dependence of diel rhythms in messenger RNA levels were evaluated in the liver tissue of the Arctic char, a cold-adapted salmonid. The results show that cold acclimation repressed diel rhythms in gene expression compared to warm-acclimated fish, in which the expression of hundreds of genes was rhythmic, suggesting the circadian clock of the Arctic fish species can be sensitive to temperature. Overall, the results of the thesis indicate that fishes’ responses to abiotic factors interact with their diel rhythms, and more studies on the consequences of these interactions are needed to comprehensively understand human impacts on ecosystems.

Conservation Genetics of Exploited Finnish Salmonid Fishes

Genetic diversity is one of the levels of biodiversity that the World Conservation Union (IUCN) has recognized as being important to preserve. This is because genetic diversity is fundamental to the future evolution and to the adaptive flexibility of a species to respond to the inherently dynamic nature of the natural world. Therefore, the key to maintaining biodiversity and healthy ecosystems is to identify, monitor and maintain locally-adapted populations, along with their unique gene pools, upon which future adaptation depends. Thus, conservation genetics deals with the genetic factors that affect extinction risk and the genetic management regimes required to minimize the risk. The conservation of exploited species, such as salmonid fishes, is particularly challenging due to the conflicts between different interest groups. In this thesis, I conduct a series of conservation genetic studies on primarily Finnish populations of two salmonid fish species (European grayling, Thymallus thymallus, and lake-run brown trout, Salmo trutta) which are popular recreational game fishes in Finland. The general aim of these studies was to apply and develop population genetic approaches to assist conservation and sustainable harvest of these populations. The approaches applied included: i) the characterization of population genetic structure at national and local scales; ii) the identification of management units and the prioritization of populations for conservation based on evolutionary forces shaping indigenous gene pools; iii) the detection of population declines and the testing of the assumptions underlying these tests; and iv) the evaluation of the contribution of natural populations to a mixed stock fishery. Based on microsatellite analyses, clear genetic structuring of exploited Finnish grayling and brown trout populations was detected at both national and local scales. Finnish grayling were clustered into three genetically distinct groups, corresponding to northern, Baltic and south-eastern geographic areas of Finland. The genetic differentiation among and within population groups of grayling ranged from moderate to high levels. Such strong genetic structuring combined with low genetic diversity strongly indicates that genetic drift plays a major role in the evolution of grayling populations. Further analyses of European grayling covering the majority of the species’ distribution range indicated a strong global footprint of population decline. Using a coalescent approach the beginning of population reduction was dated back to 1 000-10 000 years ago (ca. 200-2 000 generations). Forward simulations demonstrated that the bottleneck footprints measured using the M ratio can persist within small populations much longer than previously anticipated in the face of low levels of gene flow. In contrast to the M ratio, two alternative methods for genetic bottleneck detection identified recent bottlenecks in six grayling populations that warrant future monitoring. Consistent with the predominant role of random genetic drift, the effective population size (N_e) estimates of all grayling populations were very low with the majority of N_e estimates below 50. Taken together, highly structured local populations, limited gene flow and the small N_e of grayling populations indicates that grayling populations are vulnerable to overexploitation and, hence, monitoring and careful management using the precautionary principles is required not only in Finland but throughout Europe. Population genetic analyses of lake-run brown trout populations in the Inari basin (northernmost Finland) revealed hierarchical population structure where individual populations were clustered into three population groups largely corresponding to different geographic regions of the basin. Similar to my earlier work with European grayling, the genetic differentiation among and within population groups of lake-run brown trout was relatively high. Such strong differentiation indicated that the power to determine the relative contribution of populations in mixed fisheries should be relatively high. Consistent with these expectations, high accuracy and precision in mixed stock analysis (MSA) simulations were observed. Application of MSA to indigenous fish caught in the Inari basin identified altogether twelve populations that contributed significantly to mixed stock fisheries with the Ivalojoki river system being the major contributor (70%) to the total catch. When the contribution of wild trout populations to the fisheries was evaluated regionally, geographically nearby populations were the main contributors to the local catches. MSA also revealed a clear separation between the lower and upper reaches of Ivalojoki river system – in contrast to lower reaches of the Ivalojoki river that contributed considerably to the catch, populations from the upper reaches of the Ivalojoki river system (>140 km from the river mouth) did not contribute significantly to the fishery. This could be related to the available habitat size but also associated with a resident type life history and increased cost of migration. The studies in my thesis highlight the importance of dense sampling and wide population coverage at the scale being studied and also demonstrate the importance of critical evaluation of the underlying assumptions of the population genetic models and methods used. These results have important implications for conservation and sustainable fisheries management of Finnish populations of European grayling and brown trout in the Inari basin.

Molecular Evolution of Metazoan Hypoxia-inducing Factors

Most metazoans rely on aerobic energy production, which is dependent on adequate oxygen supply. In the case of reduced oxygen supply (hypoxia), the most profound changes in gene expression are mediated by transcription factors named hypoxia-inducible factors (HIF alpha). These proteins are post-translationally regulated by prolyl-4-hydroxylase (PHD) enzymes that are direct “sensors” of cellular oxygen levels. This thesis examines the molecular evolution of metazoan HIF systems. In early metazoans the HIF system emerged from pre-existing PHD oxygen sensors and early bHLH-PAS transcription factors. In invertebrates our analysis revealed an unexpected diversity of PHD genes and HIF alpha sequence characteristics. An early branching vertebrate, the epaulette shark (Hemiscyllium ocellatum) was chosen for sequencing and hypoxia preconditioning studies of HIF alpha and PHD genes. As no quantitative PCR reference genes were available, this thesis includes the first study of reference genes in cartilaginous fish species. Applying multiple statistical analysis we also discoveredthat commonly used reference gene software may perform poorly with some data sets. Novel reference genes allowed accurate measurements of the mRNAlevels of the studied target genes. Cartilaginous fishes have three genomic duplicates of both HIF alpha and PHD genes like mammals and teleost fishes. Combining functional divergence and selection analyses it was possible to describe how sequence changes in both HIF alpha and PHD duplicates may have contributed to the differential oxygen sensitivityof HIF alphas. Additionally, novel teleost HIF-1 alpha sequences were produced and used to reveal the molecular evolution of HIF-1 alpha in this lineage rich with hypoxia tolerant species.

Biogeochemistry in acid sulphate soil landscapes and small urban centres in Western Finland

The overall purpose of this thesis was to increase the knowledge on the biogeochemistry of rural acid sulphate (AS) soil environments and urban forest ecosystems near small towns in Western Finland. In addition, the potential causal relationship between the distribution of AS soils and geographical occurence of multiple sclerosis (MS) disease was assessed based on a review of existing literature and data. Acid sulphate soils, which occupy an area of approximately 17–24 million hectare worldwide, are regarded as the nastiest soils in the world. Independent of the geographical locality of these soils, they pose a great threat to their surrounding environment if disturbed. The abundant metal-rich acid drainage from Finnish AS soils, which is a result of sulphide oxidation due to artificial farmland drainage, has significant but spatially and temporally variable ecotoxicological impacts on biodiversity and community structure of fish, benthic invertebrates and macrophytes. This has resulted in mass fish kills and even eradication of sensitive fish species in affected waters. Moreover, previous investigations demonstrated significantly enriched concentrations of Co, Ni, Mn and Al, metals which are abundantly mobilised in AS soils, in agricultural crops (timothy grass and oats) and approximately 50 times higher concentrations of Al in cow milk originating from AS soils in Western Finland. Nevertheless, the results presented here demonstrate, in general, relatively moderate metal concentrations in oats and cabbage grown on AS soils in Western Finland, although some of the studied fields showed anomalous values of metals (e.g. Co and Ni) in both the soil and target plants (especially oats), similar to that of the previous investigations. The results indicated that the concentrations of Co, Ni, Mn and Zn in oats and Co and Zn in cabbage were governed by soil geochemistry as these metals were correlated with corresponding concentrations extracted from the soil by NH4Ac-EDTA and NH4Ac, respectively. The concentrations of Cu and Fe in oats and cabbage were uncorrelated to that of the easily soluble concentrations in the soils, suggesting that biological processes (e.g. plant-root processes) overshadow geochemical variation. The concentrations of K and Mg in cabbage, which showed a low spread and were strongly correlated to the NH4Ac extractable contents in the soil, were governed by both the bioavailable fractions in the topsoil and plant-uptake mechanisms. The plant´s ability to regulate its uptake of Ca and P (e.g. through root exudates) seemed to be more important than the influence of soil geochemistry. The distribution of P, K, Ca, Mg, Mn and S within humus, moss and needles in and around small towns was to a high degree controlled by biological cycling, which was indicated by the low correlation coefficients for P, K, Ca, Mg and S between humus and moss, and the low spread of these nutrients in moss and needles. The concentration variations of elements in till are mainly due to natural processes (e.g. intrusions, weathering, mineralogical variations in the bedrock). There was a strong spatial pattern for B in humus, moss and needles, which was suggested to be associated with anthropogenic emissions from nearby town centres. Geogenic dust affected the spatial distribution of Fe and Cr in moss, while natural processes governed the Fe anomaly found in the needles. The spatial accumulation patterns of Zn, Cd, Cu, Ni and Pb in humus and moss were strong and diverse, and related to current industry, the former steel industry, coal combustion, and natural geochemical processes. An intriguing Cu anomaly was found in moss. Since it was located close to a main railway line and because the railway line´s electric cables are made of Cu, it was suggested that the reason for the Cu anomaly is corrosion of these cables. In Western Finland, where AS soils are particularly abundant and enrich the metal concentrations of stream waters, cow milk and to some extent crops, an environmental risk assessment would be motivated to elucidate if the metal dispersion affect human health. Within this context, a topic of concern is the distribution of multiple sclerosis as high MS prevalence rates are found in the main area of AS soils. Regionally, the AS soil type in the Seinäjoki area has been demonstrated to be very severe in terms of metal leaching, this area also shows one of the highest MS rates reported worldwide. On a local scale, these severe AS soil types coincide well with the corresponding MS clustering along the Kyrönjoki River in Seinäjoki. There are reasons to suspect that these spatial correlations are causal, as multiple sclerosis has been suggested to result from a combination of genetic and environmental factors.

Tehokalastuksen vaikutukset Lappajärvessä : Kalastus elinkeinoksi Lappajärvellä -hankkeen vuosien 2001-2006 seurannan loppuraportti

Lappajärven, Etelä-Pohjanmaan suurimman järven, ekologinen tila on heikentynyt ja se luokitellaan nykyisin tyydyttäväksi. Rehevöityminen, 1990-luvun muikkukato ja yhteiskunnan rakennemuutos ovat vähentäneet kalastajien ja kalansaaliin määrää. Ammattikalastuksen edellytysten parantamiseksi toteutettiin vuosina 2001–2004 ”Kalastus elinkeinoksi Lappajärvellä” – hanke. Hankkeen vaikutuksia kalastoon ja järven ekosysteemiin seurattiin vuosina 2001–2006. Kokonaissaalis vuosina 2001–2004 oli 1655 tonnia (n. 29 kg/ha/vuosi). Vaikka saalistavoitetta (2100 tn) ei täysin saavutettu, oli hankkeella runsaasti myönteisiä, joskin osin vähäisiä tai epävarmoja vaikutuksia järven tilaan ja kalastoon. Kalastajien yksikkösaaliit etenkin kuhan ja isokokoisen ahvenen osalta nousivat hankkeen alussa. Loppuvuosina yksikkösaaliiden kehitys oli osittain ristiriitaista. Ulapan taloudellisesti vähempiarvoinen kalasto, varsinkin kuorekanta harveni selvästi. Muikku puolestaan lisääntyi ja pysyi runsaana. Kuhan ja isokokoisen ahvenen lisääntymistä selittivät pääosin lämpimät kesät ja kuhan onnistuneet istutukset. Petokalat, etenkin kuha, voivat osaltaan säädellä Lappajärven ulapan planktonsyöjäkalojen kantoja ja toisaalta estää särjen siirtymistä ulapalle. Tehokalastuksen loputtua alkoi näkyä merkkejä kalakantojen palautumisesta takaisin kohti lähtötilannetta. Rantavyöhykkeellä hankkeen vaikutukset jäivät selvästi vähäisemmiksi kuin ulapalla. Eläinplankton lisääntyi hieman ylemmissä vesikerroksissa, mutta tilanne alkoi palautua ennalleen hankkeen jälkeen. Sinileväkukinnot vaikuttivat vähentyneen hankkeen aikana, tosin kukintojen voimakas vaihtelu vaikeutti arviointia. Tehokalastus ja vesiensuojelun tehostuminen vähensivät fosforipitoisuuksia, tosin myös vähäsateiset vuodet vaikuttivat vähenemiseen. Pitoisuuksien lasku jatkui hankkeen jälkeen. Hanketta edeltänyt klorofyllipitoisuuksien kasvusuuntaus pysähtyi ja pitoisuustaso vakiintui. Tehokalastus pienensi klorofylli/fosfori-suhdetta ja vähensi sisäistä kuormitusta, mutta vaikutus alkoi heiketä hankkeen loputtua. Typpipitoisuuksiin, veden väriin tai näkösyvyyteen kalastus ei vaikuttanut. Tehokalastuksella voidaan parantaa Lappajärven tilaa, mutta kalastuksen on oltava jatkuvaa. Saalistavoitteeksi voisi riittää n. 20-30 kg/ha vuodessa, mutta samalla on jatkettava ulkoisen kuormituksen vähentämistä. Kalastuksen tulisi kohdistua arvokalojen lisäksi myös vähempiarvoisiin kalalajeihin ja samalla on tuettava petokalakantojen vahvistumista.

Automatic recognition of fish from video sequences

The problem of automatic recognition of the fish from the video sequences is discussed in this Master’s Thesis. This is a very urgent issue for many organizations engaged in fish farming in Finland and Russia because the process of automation control and counting of individual species is turning point in the industry. The difficulties and the specific features of the problem have been identified in order to find a solution and propose some recommendations for the components of the automated fish recognition system. Methods such as background subtraction, Kalman filtering and Viola-Jones method were implemented during this work for detection, tracking and estimation of fish parameters. Both the results of the experiments and the choice of the appropriate methods strongly depend on the quality and the type of a video which is used as an input data. Practical experiments have demonstrated that not all methods can produce good results for real data, whereas on synthetic data they operate satisfactorily.

Spatial and temporal variation in fish populations and assemblages in coastal waters of the northern Baltic Proper

The distribution and traits of fish are of interest both ecologically and socio-economically. In this thesis, phenotypic and structural variation in fish populations and assemblages was studied on multiple spatial and temporal scales in shallow coastal areas in the archipelago of the northern Baltic Proper. In Lumparn basin in Åland Islands, the fish assemblage displayed significant seasonal variation in depth zone distribution. The results indicate that investigating both spatial and temporal variation in small scale is crucial for understanding patterns in fish distribution and community structure in large scale. The local population of Eurasian perch Perca fluviatilis L displayed habitat-specific morphological and dietary variation. Perch in the pelagic zone were on average deeper in their body shape than the littoral ones and fed on fish and benthic invertebrates. The results differ from previous studies conducted in freshwater habitats, where the pelagic perch typically are streamlined in body shape and zooplanktivorous. Stable isotopes of carbon and nitrogen differed between perch with different stomach contents, suggesting differentiation of individual diet preferences. In the study areas Lumparn and Ivarskärsfjärden in Åland Islands and Galtfjärden in Swedish east coast, the development in fish assemblages during the 2000’s indicated a general shift towards higher abundances of small-bodied lower-order consumers, especially cyprinids. For European pikeperch Sander lucioperca L., recent declines in adult fish abundances and high mortalities (Z = 1.06–1.16) were observed, which suggests unsustainably high fishing pressure on pikeperch. Based on the results it can be hypothesized that fishing has reduced the abundances of large predatory fish, which together with bottom-up forcing by eutrophication has allowed the lower-order consumer species to increase in abundances. This thesis contributes to the scientific understanding of aquatic ecosystems with new descriptions on morphological and dietary adaptations in perch in brackish water, and on the seasonal variation in small-scale spatial fish distribution. The results also demonstrate anthropogenic effects on coastal fish communities and underline the urgency of further reducing nutrient inputs and regulating fisheries in the Baltic Sea region.