Linking molecular and population stress responses in Daphnia magna exposed to cadmium

DNA microarrays can be used to measure environmental stress responses. If they are to be predictive of environmental impact, we need to determine if altered gene expression translates into negative impacts on individuals and populations. A large cDNA microarray (14000 spots) was created to measure molecular stress responses to cadmium in Daphnia magna,the most widely used aquatic indicator species, and relate responses to population growth rate (pgr). We used the array to detect differences in the transcription of genes in juvenile D. magna (24 h old) after 24 h exposure to a control and three cadmium concentrations (6, 20, and 37 mu g Cd2+ L-1). Stress responses at the population level were estimated following a further 8 days exposure. Pgr was approximately linear negative with increasing cadmium concentration over this range. The microarray profile of gene expression in response to acute cadmium exposure begins to provide an overview of the molecular responses of D. magna, especially in relation to growth and development. Of the responding genes, 29% were involved with metabolism including carbohydrate, fat and peptide metabolism, and energy production, 31% were involved with transcription/translation, while 40% of responding genes were associated with cellular processes like growth and moulting, ion transport, and general stress responses (which included oxidative stress). Our production and application of a large Daphnia magna microarray has shown that measured gene responses can be logically linked to the impact of a toxicant such as cadmium on somatic growth and development, and consequently pgr.

Systems biology meets stress ecology: linking molecular and organismal stress responses in Daphnia magna

Background: Ibuprofen and other nonsteroidal anti-inflammatory drugs have been designed to interrupt eicosanoid metabolism in mammals, but little is known of how they affect nontarget organisms. Here we report a systems biology study that simultaneously describes the transcriptomic and phenotypic stress responses of the model crustacean Daphnia magna after exposure to ibuprofen. Results: Our findings reveal intriguing similarities in the mode of action of ibuprofen between vertebrates and invertebrates, and they suggest that ibuprofen has a targeted impact on reproduction at the molecular, organismal, and population level in daphnids. Microarray expression and temporal real-time quantitative PCR profiles of key genes suggest early ibuprofen interruption of crustacean eicosanoid metabolism, which appears to disrupt signal transduction affecting juvenile hormone metabolism and oogenesis. Conclusion: Combining molecular and organismal stress responses provides a guide to possible chronic consequences of environmental stress for population health. This could improve current environmental risk assessment by providing an early indication of the need for higher tier testing. Our study demonstrates the advantages of a systems approach to stress ecology, in which Daphnia will probably play a major role.

The ecological niche of Daphnia magna characterized using population growth rate

The concept of an organism's niche is central to ecological theory, but an operational definition is needed that allows both its experimental delineation and interpretation of field distributions of the species. Here we use population growth rate (hereafter, pgr) to de. ne the niche as the set of points in niche space where pgr. 0. If there are just two axes to the niche space, their relationship to pgr can be pictured as a contour map in which pgr varies along the axes in the same way that the height of land above sea level varies with latitude and longitude. In laboratory experiments we measured the pgr of Daphnia magna over a grid of values of pH and Ca2+, and so defined its "laboratory niche'' in pH-Ca2+ space. The position of the laboratory niche boundary suggests that population persistence is only possible above 0.5 mg Ca2+/L and between pH 5.75 and pH 9, though more Ca2+ is needed at lower pH values. To see how well the measured niche predicts the field distribution of D. magna, we examined relevant field data from 422 sites in England and Wales. Of the 58 colonized water bodies, 56 lay within the laboratory niche. Very few of the sites near the niche boundary were colonized, probably because pgr there is so low that populations are vulnerable to extinction by other factors. Our study shows how the niche can be quantified and used to predict field distributions successfully.

Chronic toxicity of ibuprofen to Daphnia magna: Effects on life history traits and population dynamics

The non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) is a widely used pharmaceutical that can be found in several freshwater ecosystems. Acute toxicity studies with Daphnia magna suggest that the 48 h EC50 (immobilisation) is 10-100 mg IB l(-1). However, there are currently no chronic IB toxicity dataon arthropod populations, and the aquatic life impacts of such analgesic drugs are still undefined. We performed a 14-day exposure of D. magna to IB as a model compound (concentration range: 0, 20, 40 and 80 mg IB l(-1)) measuring chronic effects on life history traits and population performance. Population growth rate was significantly reduced at all IB concentrations, although survival was only affected at 80 mg IB l(-1). Reproduction, however, was affected at lower concentrations of IB (14-day EC50 of 13.4 mg IB l(-1)), and was completely inhibited at the highest test concentration. The results from this study indicate that the long-term crustacean population consequences of a chronic IB exposure at environmentally realistic concentrations (ng l(-1) to mu g l(-1)) would most likely be of minor importance. We discuss our results in relation to recent genomic studies, which suggest that the potential mechanism of toxicity in Daphnia is similar to the mode of action in mammals, where IB inhibits eicosanoid biosynthesis. (C) 2007 Elsevier Ireland Ltd. All rights reserved.

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favorability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB)versus Thr(HB)-Asn(nHB). Significantly (p < or = 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi1 rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi1 conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi1 and chi2 conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favorability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http://www.rubic.rdg.ac.uk/betapairprefsparallel/).

Amino acid pairing preferences in parallel beta-sheets in proteins

Statistical approaches have been applied to examine amino acid pairing preferences within parallel beta-sheets. The main chain hydrogen bonding pattern in parallel beta-sheets means that, for each residue pair, only one of the residues is involved in main chain hydrogen bonding with the strand containing the partner residue. We call this the hydrogen bonded (HB) residue and the partner residue the non-hydrogen bonded (nHB) residue, and differentiate between the favourability of a pair and that of its reverse pair, e.g. Asn(HB)-Thr(nHB) versus Thr(HB)-Asn(nHB). Significantly (p <= 0.000001) favoured pairings were rationalised using stereochemical arguments. For instance, Asn(HB)-Thr(nHB) and Arg(HB)-Thr(nHB) were favoured pairs, where the residues adopted favoured chi(1) rotamer positions that allowed side-chain interactions to occur. In contrast, Thr(HB)-Asn(nHB) and Thr(HB)-Arg(nHB) were not significantly favoured, and could only form side-chain interactions if the residues involved adopted less favourable chi(1) conformations. The favourability of hydrophobic pairs e.g. Ile(HB)-Ile(nHB), Val(HB)-Val(nHB) and Leu(HB)-Ile(nHB) was explained by the residues adopting their most preferred chi(1) and chi(2) conformations, which enabled them to form nested arrangements. Cysteine-cysteine pairs are significantly favoured, although these do not form intrasheet disulphide bridges. Interactions between positively and negatively charged residues were asymmetrically preferred: those with the negatively charged residue at the HB position were more favoured. This trend was accounted for by the presence of general electrostatic interactions, which, based on analysis of distances between charged atoms, were likely to be stronger when the negatively charged residue is the HB partner. The Arg(HB)-Asp(nHB) interaction was an exception to this trend and its favourability was rationalised by the formation of specific side-chain interactions. This research provides rules that could be applied to protein structure prediction, comparative modelling and protein engineering and design. The methods used to analyse the pairing preferences are automated and detailed results are available (http:// www.rubic.rdg.ac.uk/betapairprefsparallel/). (c) 2005 Elsevier Ltd. All rights reserved.

Expression of target and reference genes in Daphnia magna exposed to ibuprofen

Background: Transcriptomic techniques are now being applied in ecotoxicology and toxicology to measure the impact of stressors and develop understanding of mechanisms of toxicity. Microarray technology in particular offers the potential to measure thousands of gene responses simultaneously. However, it is important that microarrays responses should be validated, at least initially, using real-time quantitative polymerase chain reaction (QPCR). The accurate measurement of target gene expression requires normalisation to an invariant internal control e. g., total RNA or reference genes. Reference genes are preferable, as they control for variation inherent in the cDNA synthesis and PCR. However, reference gene expression can vary between tissues and experimental conditions, which makes it crucial to validate them prior to application. Results: We evaluated 10 candidate reference genes for QPCR in Daphnia magna following a 24 h exposure to the non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) at 0, 20, 40 and 80 mg IB l(-1). Six of the 10 candidates appeared suitable for use as reference genes. As a robust approach, we used a combination normalisation factor (NF), calculated using the geNorm application, based on the geometric mean of three selected reference genes: glyceraldehyde-3-phosphate dehydrogenase, ubiquitin conjugating enzyme and actin. The effects of normalisation are illustrated using as target gene leukotriene B4 12-hydroxydehydrogenase (Ltb4dh), which was upregulated following 24 h exposure to 63-81 mg IB l(-1). Conclusions: As anticipated, use of the NF clarified the response of Ltb4dh in daphnids exposed to sublethal levels of ibuprofen. Our findings emphasise the importance in toxicogenomics of finding and applying invariant internal QPCR control(s) relevant to the study conditions.

The use of image analysis to estimate population growth rate in Daphnia magna

1. Population growth rate (PGR) is central to the theory of population ecology and is crucial for projecting population trends in conservation biology, pest management and wildlife harvesting. Furthermore, PGR is increasingly used to assess the effects of stressors. Image analysis that can automatically count and measure photographed individuals offers a potential methodology for estimating PGR. 2. This study evaluated two ways in which the PGR of Daphnia magna, exposed to different stressors, can be estimated using an image analysis system. The first method estimated PGR as the ratio of counts of individuals obtained at two different times, while the second method estimated PGR as the ratio of population sizes at two different times, where size is measured by the sum of the individuals' surface areas, i.e. total population surface area. This method is attractive if surface area is correlated with reproductive value (RV), as it is for D. magna, because of the theoretical result that PGR is the rate at which the population RV increases. 3. The image analysis system proved reliable and reproducible in counting populations of up to 440 individuals in 5 L of water. Image counts correlated well with manual counts but with a systematic underestimate of about 30%. This does not affect accuracy when estimating PGR as the ratio of two counts. Area estimates of PGR correlated well with count estimates, but were systematically higher, possibly reflecting their greater accuracy in the study situation. 4. Analysis of relevant scenarios suggested the correlation between RV and body size will generally be good for organisms in which fecundity correlates with body size. In these circumstances, area estimation of PGR is theoretically better than count estimation. 5. Synthesis and applications. There are both theoretical and practical advantages to area estimation of population growth rate when individuals' reproductive values are consistently well correlated with their surface areas. Because stressors may affect both the number and quality of individuals, area estimation of population growth rate should improve the accuracy of predicting stress impacts at the population level.

Population growth rate and carrying capacity for springtails Folsomia candida exposed to ivermectin

Forecasting the effects of stressors on the dynamics of natural populations requires assessment of the joint effects of a stressor and population density on the population response. The effects can be depicted as a contour map in which the population response, here assessed by Population growth rate, varies with stress and density in the same way that the height of land above sea level varies with latitude and longitude. We present the first complete map of this type using as our model Folsomia candida exposed to five different concentrations of the widespread anthelmintic veterinary medicine ivermectin in replicated microcosm experiments lasting 49 days. The concentrations of ivermectin in yeast were 0.0, 6.8 28.83 66.4 and 210.0 mg/L wet weight. Increasing density and chemical concentration both significantly reduced the population growth rate of Folsomia candida, in part through effects on food consumption and fecundity. The interaction between density and ivermectin concentration was "less-than-additive," implying that at high density populations were able to compensate for the effects of the chemical. This result demonstrates that regulatory protocols carried out at low density (as in most past experiments) may seriously overestimate effects in the field, where densities are locally high and populations are resource limited (e.g., in feces of livestock treated with ivermectin).

Towards a population ecology of stressed environments: the effects of zinc on the springtail Folsomia candida

1. To understand population dynamics in stressed environments it is necessary to join together two classical lines of research. Population responses to environmental stress have been studied at low density in life table response experiments. These show how the population's growth rate (pgr) at low density varies in relation to levels of stress. Population responses to density, on the other hand, are based on examination of the relationship between pgr and population density. 2. The joint effects of stress and density on pgr can be pictured as a contour map in which pgr varies with stress and density in the same way that the height of land above sea level varies with latitude and longitude. Here a microcosm experiment is reported that compared the joint effects of zinc and population density on the pgr of the springtail Folsomia candida (Collembola). 3. Our experiments allowed the plotting of a complete map of the effects of density and a stressor on pgr. Particularly important was the position of the pgr= 0 contour, which suggested that carrying capacity varied little with zinc concentration until toxic levels were reached. 4. This prediction accords well with observations of population abundance in the field. The method also allowed us to demonstrate, simultaneously, hormesis, toxicity, an Allee effect and density dependence. 5. The mechanisms responsible for these phenomena are discussed. As zinc is an essential trace element the initial increase in pgr is probably a consequence of dietary zinc deficiency. The Allee effect may be attributed to productivity of the environment increasing with density at low density. Density dependence is a result of food limitation. 6. Synthesis and applications. We illustrate a novel solution based on mapping a population's growth rate in relation to stress and population density. Our method allows us to demonstrate, simultaneously, hormesis, toxicity, an Allee effect and density dependence in an important ecological indicator species. We hope that the approach followed here will prove to have general applicability enabling predictions of field abundance to be made from estimates of the joint effects of the stressors and density on population growth rate.

Joint effects of density and a growth inhibitor on the life history and population growth rate of the midge Chironomus riparius

Results of previous laboratory studies suggest that high population density often buffers the effects of chemical stressors that predominately increase mortality. Mortality stressors act to release more resources for the survivors and, therefore, produce less-than-additive effects. By contrast, growth stressors are expected to have opposite results or more-than-additive effects. We investigated the effects of a growth inhibitor (lufenuron) on larval growth and survival of Chironomus riparius and examined its joint effects with density on population growth rate (PGR). Exposure to 60 mu g/kg sediment or greater inhibited larval growth, and exposure to 88 mu g/kg or greater often resulted in mortality before reaching emergence. The effects of lufenuron, however, differed with population density. At 88 mu g/kg, mortalities and, to a lesser extent, reduced fecundity resulted in a reduction in PGR at low density. Conversely, when populations were initiated at high density, PGR was similar to that of controls, because the few survivors reached maturity sooner and started producing offspring earlier. The effect of density as a growth stressor therefore was stronger than the effect of lufenuron, which had effects similar to those of a mortality stressor and produced less-than-additive effects. Longterm studies under field conditions, however, are needed before less-than-additive effects are considered to be the norm.

The genetic basis of resistance to anticoagulants in rodents

Anticoagulant compounds, i.e., derivatives of either 4-hydroxycoumarin (e.g., warfarin, bromadiolone) or indane-1,3-dione (e.g., diphacinone, chlorophacinone), have been in worldwide use as rodenticides for > 50 years. These compounds inhibit blood coagulation by repression of the vitamin K reductase reaction (VKOR). Anticoagulant-resistant rodent populations have been reported from many countries and pose a considerable problem for pest control. Resistance is transmitted as an autosomal dominant trait although, until recently, the basic genetic mutation was unknown. Here, we report on the identification of eight different mutations in the VKORC1 gene in resistant laboratory strains of brown rats and house mice and in wild-caught brown rats from various locations in Europe with five of these mutations affecting only two amino acids (Tyr139Cys, Tyr139Ser, Tyr139Phe and Leu128Gln, Leu128Ser). By recombinant expression of VKORC1 constructs in HEK293 cells we demonstrate that mutations at Tyr139 confer resistance to warlarin at variable degrees while the other mutations, in addition, dramatically reduce VKOR activity. Our data strongly argue for at least seven independent mutation events in brown rats and two in mice. They suggest that mutations in VKORC1 are the genetic basis of anticoagulant resistance in wild populations of rodents, although the mutations alone do not explain all aspects of resistance that have been reported. We hypothesize that these mutations, apart from generating structural changes in the VKORC1 protein, may induce compensatory mechanisms to maintain blood clotting. Our findings provide the basis for a DNA-based field monitoring of anticoagulant resistance in rodents.