Influence of convolution filtering on coronary plaque attenuation values: observations in an ex vivo model of multislice computed tomography coronary angiography

Attenuation variability (measured in Hounsfield Units, HU) of human coronary plaques using multislice computed tomography (MSCT) was evaluated in an ex vivo model with increasing convolution kernels. MSCT was performed in seven ex vivo left coronary arteries sunk into oil followingthe instillation of saline (1/∞) and a 1/50 solution of contrast material (400 mgI/ml iomeprol). Scan parameters were: slices/collimation, 16/0.75 mm; rotation time, 375 ms. Four convolution kernels were used: b30f-smooth, b36f-medium smooth, b46f-medium and b60f-sharp. An experienced radiologist scored for the presence of plaques and measured the attenuation in lumen, calcified and noncalcified plaques and the surrounding oil. The results were compared by the ANOVA test and correlated with Pearson’s test. The signal-to-noise ratio (SNR) and contrast-to-noise ratio (CNR) were calculated. The mean attenuation values were significantly different between the four filters (p < 0.0001) in each structure with both solutions. After clustering for the filter, all of the noncalcified plaque values (20.8 ± 39.1, 14.2 ± 35.8, 14.0 ± 32.0, 3.2 ± 32.4 HU with saline; 74.7 ± 66.6, 68.2 ± 63.3, 66.3 ± 66.5, 48.5 ± 60.0 HU in contrast solution) were significantly different, with the exception of the pair b36f–b46f, for which a moderate-high correlation was generally found. Improved SNRs and CNRs were achieved by b30f and b46f. The use of different convolution filters significantly modifief the attenuation values, while sharper filtering increased the calcified plaque attenuation and reduced the noncalcified plaque attenuation.
Introduction
Multislice computed tomography (MSCT) of the coronary arteries is becoming an established method to exclude the presence of significant stenosis because of the high negative predictive value [1–6]. Stenosis, however, occurs in a later stage of coronary artery disease (CAD); consequently, the atherosclerotic plaque composition and morphology may represent a better predictor of plaque stability [7–9]. Until recently, only invasive techniques, such as intracoronary ultrasound (ICUS), have been able to provide consistently reproducible quantitative and qualitative information on plaque composition [10]. ICUS, however, cannot be used for routine evaluation of plaque characteristics because of its invasiveness and related increased risk, additional time and cost.
Magnetic resonance imaging may be used to characterize in vivo the atherosclerotic plaque morphology and composition in carotid [11] and coronary arteries [12, 13].
Several recent in vivo studies have stressed that MSCT can be used to characterize coronary atherosclerotic plaque and has the potential to become the noninvasive modality of choice for plaque imaging [14–17]. This technique is able to provide information on noncalcified coronary plaques, which may be useful to identify a high-risk plaque [18–20].
MSCT plaque assessment is based upon the variable X-ray attenuation of the tissue components. However, several parameters, such as lumen attenuation, convolution filtering, body mass index of the patient and contrast-to-noise ratio (CNR) of the images, are able to modify the attenuation values that are being used to define the composition of coronary plaques.
The aim of the present study was to investigate the effect of increasing convolution filtering on plaque components with varying coronary attenuation in an ex vivo coronary model.
Material and methods
Specimens
Seven ex vivo left anterior descending (LAD) coronary arteries were obtained at autopsy (Fig. 1) from seven deceased patients: four male and three female, with a mean age of 68 ± 10 years: five (four males and one female with an age between 62 and 71 years) jad died of noncardiovascular diseases, and two (one male and one female with an age of 65 and 73 years, respectively) had died of ischemic heart disease. The Institutional Review Board approved the study protocol.
Fig. 1Example of proximal left coronary artery specimen. CX Left circumflex, LAD left anterior descending artery, LM left main coronary artery
The specimens were prepared and examined separately. Each coronary artery was fitted with two cannulas fixed with surgical thread in the proximal (in the left main) and distal ends (in the LAD). The circumflex artery had been closed earlier at its end with thread. Only the LAD was used because of the major length (segments considered with a mean length of 7.57 cm) and higher prevalence of atherosclerotic disease and for technical reasons (ability to fix a surgical thread in the left main, without hampering the evaluation of the proximal LAD).
Contrast material
Two solutions were used: a saline and a 1/50 dilution of contrast material (400 mgI/ml Iomeprol; Bracco, Italy). The attenuation values (Hounsfield Units, HU) of the two solutions, measured in a 10-ml syringe after dilution, were 15.9 ± 1.8 HU (1/∞; defined as saline; no contrast material was diluted) and 414.8 ± 5.6 HU (1/50; defined as contrast solution). We preferred to use a 1/50 dilution to achieve a mean coronary in-lumen density >300 HU – which is not typically achieved in in vivo cardiac computed tomography angiography (CTA) studies – in order to emphasize the difference with the saline solution. Our study benefits from two consecutive scans of the same vessel (e.g. with and without contrast media) for a comparative evaluation of plaques, which is not truly feasible during in vivo studies due to high radiation exposure.
Experimental settings
A box was filled with olive oil. Prior to positioning the specimen in the oil, saline was injected through the sheaths to wash out the air bubbles in the lumen as much as possible.
Once the specimen had sunk into the oil, which simulated epicardial fat, the solution was instilled through the sheaths using a 10-ml syringe from the proximal end of the specimen. The injection was finished when the solution was observed leaking out of the distal end of specimen. The leaking solution was removed from the specimen using a syringe. The same procedure was used to fill the LAD with both solutions. The specimen was kept for all investigations in the same longitudinal position (head-LM-to feet-distal end of LAD ) to obtain direct cross-section images.
Scan parameters
A MSCT scan (Somatom Sensation 16; Siemens, Germany) was performed following the intra-coronary injection of two solutions.
Scans were performed at the following parameters: slices/collimation, 16/0.75 mm; rotation time, 375 ms; feed per rotation 3.0 mm (pitch: 0.25 mm); 120 kv; 400 mAs; effective slice thickness, 1 mm; reconstruction increment, 0.5 mm; field of view (FOV), 100 mm. Four convolution filters were used: b30f (smooth), b36f (medium smooth), b46f (medium) and b60f (sharp). Based on our clinical experience in routine evaluation, b30f was used instead of b20f because of a better delineation of the coronary arteries and atherosclerotic lesions.
The scan geometry was based on a retrospective electrocardiogram (ECG)-gated protocol (the same used for the in vivo examination). This protocol is based on a low pitch that allows a retrospective reconstruction of multiple phases within the cardiac cycle. A demo-ECG was switched on and a heart rate of 71 occurred during the scan. The half-rotation reconstruction algorithm brought the effective temporal resolution down to 187 ms.
Data collection and analysis
A radiologist with a level 3 expertise in cardiac CT [21] loaded the data sets into a dedicated workstation (Leonardo; Siemens, Germany) and performed orthogonal views of all specimens for each convolution filter with two solutions.
One operator performed all of the measurements. Each specimen was evaluated for the presence of atherosclerotic plaque. The coronary wall thickening of plaque (minimum size observed: 0.3 mm2) was clearly visible from the surrounding oil with low attenuation and from the lumen after injection of contrast material. Plaques were targeted regardless of their size. The operator loaded the two data sets for each solution of the same specimen into a workstation screen, scrolling the data sets in parallel with standard soft-tissue window settings (window width:700 HU; window center:140 HU).
Once a plaque was detected in the two solutions by means of a comparative evaluation, the operator drew four regions of interest (ROIs) for the contrast material: in the lumen of the vessel (defined as lumen), the soft tissue of the coronary plaque (defined as noncalcified plaque), the calcification within the coronary artery wall (defined as calcified plaque) and the oil surrounding the plaque (defined as surrounding). The ROIs (minimum size: 0.1 mm2) were drawn as large as possible while avoiding the borders of each structure in order to limit the partial volume effect. Once the four ROIs were drawn, the operator could copy and paste these into other stacks of images. The two data sets for each solution were then evaluated using four convolution filters (b30f, b36f, b46f and b60f). For each ROI the mean and the standard deviation of the attenuation value were collected. Overall, 3200 attenuation values were obtained.
Statistical analysis
The attenuation values are presented as means and standard deviations. The signal-to-noise ratio (SNR) was calculated as the mean attenuation value of the ROI/standard deviation of the surrounding oil attenuation value. The standard deviation HU value of the surrounding air is often used for the estimation of noise. However, we did not include the air in our images. The standard deviation of the surrounding oil can estimate the noise because of its reliable signal. SNRs are presented as means for each kernel.
The contrast-to-noise ratio (CNR) was calculated as (mean attenuation A – mean attenuation B)/image noise (defined as the standard deviation of the surrounding oil). CNRs of (lumen – noncalcified plaque), (noncalcified plaque – surrounding), (calcified plaque – noncalcified plaque) and (calcified plaque – lumen) were calculated and presented as average for each kernel and solution.
Statistical evaluation was performed with dedicated software (SPSS ver. 10.1; SPSS, Chicago, Ill.). A one-way ANOVA test was performed to analyze the statistical difference between kernels in each structure for two solutions. A p < 0.0001 was considered for statistical significance. The Tukey HSD test evaluated the statistical absolute difference between any two groups of values. The Pearson’s correlation was also tested between all the different convolution filters.
After clustering for kernel and structure, the attenuation values obtained with the two solutions were compared with Student’s t-test and correlated with Pearson’s test.
Results
A total of 100 sections containing mixed plaques (n = 20) were evaluated in the seven coronary specimens. Two solutions and four convolution kernels were available (800 slices) at each level, and four ROIs (lumen, surrounding, calcified plaque, noncalcified plaque) were sampled in each slice. A total of 3200 measurements were performed. The results are summarized in Table 1. The attenuation values were all significantly different between the convolution filters in each structure with both solutions (p < 0.0001). The results of the Tukey HSD to test the difference between pairs of four filters showed that the noncalcified plaque attenuation values were all significantly different using both solutions (p < 0.01), with the exception of pair b36f–b46f. The attenuation values between b30f and b46f were significantly different at a high concentration in calcified plaques, whereas for the same pair of filters the attenuation values of the lumen and the surrounding were not significantly different. After clustering the paired attenuation values of the four structures for the convolution filter, a moderate-high correlation was generally found, with the exception of pairs b30f–b46f in surrounding (r = 0.150), b30f–b60f in lumen (r = 0.206) and surrounding (r = 0.311), b36f–b60f in lumen (r = 0.132) and b46f-b60f (r = 0.201) in lumen in the saline solution.
Table 1Summary of the attenuation values measured in each structure in the two solutions (e.g. saline and contrast solutions) with increasing convolution filter b30fb36fb46fb60fMeanSDSNRMeanSDSNRMeanSDSNRMeanSDSNRLumenSaline58.140.429.952.140.59.448.236.019.123.56.01.8Contrast329.193.7209.8312.685.544.0334.191.7169.7331.3138.317.8Noncalcified PlaqueSaline20.839.18.014.235.82.414.032.04.03.232.40.3Contrast74.766.639.368.263.39.266.366.525.948.560.02.4Calcified Plaque Saline740.5392.0327.0758.5360.3133.6785.7388.5294.91145.8517.487.8Contrast795.4333.8498.0838.7364.2116.0885.6382.6488.71194.8520.761.0SurroundingSaline127.12.7−54.7−128.03.1−22.6−128.72.7−47.6−134.811.3−10.5Contrast−128.22.6−84.5−126.84.3−17.6−128.52.4−7.2−123.710.0−6.4
Whereas the convolution filter did not significantly affect the attenuation values of the lumen and the surrounding, sharper filters decreased the attenuation of the plaque and, conversely, increased the attenuation of the calcification with both solutions (Fig. 2).
Fig. 2Mean attenuation of the four structures in the two solutions (e.g. saline and contrast) with increasing convolution filters. Whereas the convolution filter does not significantly affect the attenuation values of the lumen and the surrounding, sharper filters decrease the attenuation of the noncalcified plaque and, conversely, increase the attenuation of the calcified plaque with both solutions
The SNR values appear to have a common trend with increasing convolution filtering and both solutions. However, the best SNR was achieved with the b30f and b46f filters. Plaque values follow the lumen values, while the surrounding and calcification patterns are similar (Fig. 3).
Fig. 3The SNR values appear to have a common trend with increasing convolution filtering and both solutions (e. g. saline and contrast). However, the best SNR is achieved with b30 and b46 filters. Noncalcified plaque values follow the lumen values, while the surrounding and the calcified plaque patterns are similar
The CNR values were calculated between structures with increasing convolution filters and different solutions (Table 2). Improved CNRs were achieved with the b30f and b46f filters. The lumen-noncalcified plaque CNR was particularly improved with the b30f and b46f filters using a high solution of contrast media. The noncalcified plaque-surrounding CNR was improved with the b30f and b46f filters in contrast solution. Calcified plaque-noncalcified plaque CNR increased with the contrast solution with the b30f and b46f filters, while it decreased using the b36f and b60f filters. The b46f filter was found to be ideal for calcified plaque-lumen contrast since it improved the CNR at contrast solution.
Table 2Summary of CNR values calculated between structures with increasing convolution filters and different solutions b30fb36fb46fb60fSalineLumen – noncalcified plaque21.87.115.11.5Non – calcified plaque-surrounding62.725.051.610.8Calcified plaque – noncalcified plaque318.9131.3290.887.4Calcified plaque – lumen297.1124.2275.785.9Contrast solutionLumen – noncalcified plaque170.634.8143.815.4Noncalcified plaque – surrounding123.726.997.18.8Calcified plaque – noncalcified plaque458.7106.7462.858.5Calcified plaque – lumen288.271.9319.043.1
After clustering the paired attenuation values obtained for the four structures in order to test the statistical difference between saline and contrast solution for each convolution filter, values were all significantly different (p < 0.05) apart from surrounding with b46f and calcification with b60f. Pearson’s test has shown good correlation only between the attenuation values of calcification (r = 0.738, r = 0.826, r= 0.759, r = 0.723 with b30f, b36f, b46f, b60f, respectively).
Discussion
The evaluation of the coronary artery wall is a data set that has become available since the introduction of coronary MSCT. Several studies have reported the ability of MSCT to detect coronary plaques and subsequently characterize their tissue composition (i.e. calcified, mixed, noncalcified) [14, 20]. In an in vivo study, Kopp et al. showed an excellent correspondence between ICUS criteria (soft, intermediate and calcified plaques) and MSCT attenuation values [14]. Schroeder et al. subsequently found a strong correlation between tissue density measurements within the plaque and the qualitative ICUS criteria. The CT attenuation results were 14 ± 26 HU for plaques classified by ICUS as soft (i.e. predominantly lipid) plaques, 91 ± 21 HU for plaques classified by ICUS as intermediate (i.e. predominantly fibrous) plaques and 419 ± 194 HU for plaques classified by ICUS as calcified. These researchers reported no overlap in the mean attenuation among the three groups [15]. Leber et al. showed that lesion echogenity correlated well with MSCT attenuation measurements in coronary plaque, with MSCT correctly classifying 78% of sections containing hypoechoic plaque areas (soft plaques), 78% of sections containing hyperechoic plaque areas (fibrous plaques) and 95% of sections containing calcified plaque tissue. The MSCT density measurements within coronary lesions revealed significantly different values for hypoechoic, hyperechoic, and calcified plaques [17]. The close correlation between ICUS and MSCT in terms of the detection of calcified and noncalcified coronary plaques was also shown by Achenbach et al., although this group found a sensitivity of only 53% for exclusively noncalcified plaque [16].
Becker et al. compared the atherosclerotic lesions of 11 human cadaver heart specimens detected by MSCT with the histopathological macroscopic findings according to American Heart Association criteria. They concluded that MSCT is a promising tool for the characterization of atherosclerotic coronary plaques [22].
Such information is of great value since the risk of acute coronary syndromes caused by plaque disruption and thrombosis depends on plaque composition (i.e. noncalcified, predominantly lipid plaques with positive remodeling) rather than stenosis severity [18–20]. Reliable and reproducible non-invasive methods for the assessment of plaque constitution could be important in risk stratification of patients with CAD.
Nevertheless, several aspects (i.e. the impact of lumen attenuation, convolution filters, body mass index of the patient, CNR of the images and coronary calcification) of the methodology need to be better addressed in order to validate the accuracy of attenuation values measured within the coronary plaques using MSCT.
It is controversial if the absolute HU values are valid measurements for determining the plaque composition and identifying the corresponding tissue. A phantom study revealed that intracoronary attenuation can significantly affect the measured plaque attenuation so that accurate plaque tissue differentiation may be influenced [23]. Furthermore, Cademartiri et al. recently reported that the intravascular attenuation modifies significantly the attenuation of the coronary atherosclerotic plaques in an ex vivo coronary specimen [24]. Therefore, the characterization of the plaque on the basis of absolute attenuation values should be reported with caution, and the intraluminal attenuation should also be reported. The calibration of tissue attenuation values with contrast in-lumen values may improve accurate tissue component identification.
New possibilities for coronary artery plaque imaging may be offered by dual-source CT, which provides robust diagnostic image quality independent of the heart rate. Potential applications of dual-energy CT include tissue characterization and Ca quantification of atherosclerotic plaques [25, 26].
In our experimental study we addressed the issue of the variability of coronary plaque attenuation by varying the convolution filters (Fig. 4).
Fig. 4Example of plaque with increasing filtering using a saline and contrast solution
The direct back projection of CT attenuation profiles results in unsharp images. To prevent unsharp and approximate images, each projection has to be convolved with a predetermined mathematical function, the convolution kernel; this kernel represents a filtering procedure that is able to enhance or reduce object boundaries. Convolution influences image characteristics by the choice and design of the kernel (i.e. smooth, sharp). A smooth filter reduces spatial resolution as well as image noise; a sharp filter increases spatial resolution as well as image noise. The noise is caused by the fluctuations in the number of X-ray quanta registered by the detector and is affected by convolution and dose. The SNR evaluates the ratio between the intensity of the signal and its statistical fluctuation. Suzuki et al. have demonstrated that the accuracy of a vascular diameter measurement with an automated software is affected by the choice of convolution filtering [27]. Seifarth et al. recently compared the effect of different convolution filtering on the visualization of coronary artery stents. Using the edge-enhancing b46f, they observed a significant improvement in lumen assessment compared to the conventional b30f [28].
Our study shows how the use of different convolution filters significantly affects the measured coronary plaque attenuation (Fig. 3). The use of sharper convolution filters produces an increase in the mean attenuation value of the calcified component of plaque. Increasingly sharper filtering produces a reduction in the mean attenuation values of the noncalcified plaque component. On the basis of our results, it is still challenging to rely on absolute HU values since the protocol used in coronary MSCT has not been standardized. Furthermore, absolute plaque attenuation values determined in the various studies are not comparable. It is reasonable to use a proper convolution filtering with improved SNR and CNR to better focus on plaque components. The b36f intermediate filter obtained a slight decrease of noncalcified plaque attenuation and an increase of calcified plaque attenuation. However, improved SNRs and CNRs were also achieved by the b30f and b46f filters. The b36f underestimated the noncalcified plaque attenuation and overestimated the calcified plaque attenuation in the likely attempt to weigh both components.
There are several limitations to our study. The first one is inherent to the lack of an evaluation of intra- and inter-observer variability. A second limitation is related to the fact that the coronary specimens were collected from elderly patients with a high prevalence of calcification. The blooming artifacts of calcification may impair coronary noncalcified plaque visualization and characterization based on HU values [29]. Furthermore, we did not investigate the influence of mAs and voltage (kV) on the blooming effect. Nevertheless, it has been reported that the voltage setting employed (i.e.120 kV) might be reasonably chosen during in vivo studies [29]. In vivo studies should also better address the issue of body mass, which affects coronary plaque measurements [30]. Moreover, according to data mined from the literature, it is reasonable to expect that plaque size (especially in calcified plaques) and lumen area would be influenced by the use of different convolution filters due to variable partial volume effect. Edge-enhancing sharper reconstruction filters are reported to increase the visible in-stent lumen diameter [28, 31].
The pathological correlation has not been provided since the aim of our study was not to compare pathology with MSCT, but to assess the influence of convolution filtering on plaque attenuation measurement.
In conclusion, the use of different convolution filtering significantly modified plaque attenuation values. Therefore, convolution filtering should be reported when attenuation values are shown in order to provide a standardization of the methodology. Sharper convolution kernels increased the attenuation of the calcium within the coronary plaques and reduced the attenuation of soft plaque tissues. Improved SNR and CNR were achieved by using the b30f and b46f filters. A smooth filter (e.g. b30f) may be used for clinical routine evaluation, while a medium filter (e.g. b46f) may be considered the best choice for the assessment of highly calcified or stented vessels. The use of a proper filtering according to plaque type could give a more reliable assessment of plaque attenuation values in terms of HU.

coronary plaque
multislice computed tomography
coronary angiography
convolutions kernels


Development, standardization and refinement of procedures for evaluating effects of endocrine active compounds on development and sexual differentiation of Xenopus laevis

Xenopus laevis has been introduced as a model to study effects of endocrine-active compounds (EAC) on development and sexual differentiation. However, variable and inconsistent data have raised questions about the reliability of the test methods applied. The current study was conducted in two laboratories to develop, refine, and standardize procedures and protocols. Larvae were exposed in flow-through systems to 17β-estradiol (E2), at concentrations from 0.2 to 6.0 μg E2 L−1 in Experiment 1A, and 0.015 to 2.0 μg E2 L−1 in Experiment 1B. In both studies survival (92%, 99%) and percentage of animals that completed metamorphosis (97%, 99%) indicated reproducible biological performance. Furthermore, minor variations in husbandry led to significant differences in snout-to-vent length, weight, and gonad size. In Experiment 1A, almost complete feminization occurred in all E2 treatment groups whereas a concentration response was observed in Experiment 1B resulting in an EC50 of 0.12 μg E2 L−1. The final verified protocol is suitable for determining effects of EAC on development and sexual differentiation in X. laevis.
Introduction
During the last ten years, the number of chemicals evaluated for endocrine-active effects has increased. Potential exposure of wildlife populations, including fish, amphibians, birds, and mammals, to environmental estrogens has attracted much attention due to the feasible impact of these substances on survival and population sustainability [1–6]. Laboratory studies using a range of wildlife species have demonstrated that physiological processes including development and reproduction, are sensitive to environmental estrogens [7, 8]. Due to the diversity of wildlife, certain species have been used as representatives of the various taxa that are considered in eco-toxicological hazard assessment. One example of a well established standardized animal model is the frog embryo teratogenesis assay (FETAX) [9] using Xenopus laevis embryos to determine impacts of chemicals on early development. This species has been investigated focusing on larval stages for studying the hormonal regulation of metamorphosis, development, sexual differentiation, and the impact of endocrine active compounds (EAC) [10–14]. However, the results in the literature for basic biological parameters such as mortality, growth, development, and sexual differentiation vary considerably among laboratories [15–17] making interpretation difficult. For example, it has been reported that the herbicide atrazine feminizes and demasculinizes amphibians [18–21] but other investigators have not been able to replicate these findings [22, 23].
The published data on potential effects of atrazine on gonadal development in laboratory studies and reproductive effects in wildlife populations were evaluated by a United States Environmental Protection Agency (EPA) Scientific Advisory Panel (SAP) in 2003. The SAP identified a number of factors that adversely affected the quality and comparability of data in the open literature on atrazine. These factors included variability in study design, differences in animal husbandry conditions (e.g., feeding rates), excess biomass in treatment chambers, and the effects of these factors on water quality. Other factors that contributed to difficulty in interpreting the results were variability in the response to the positive control substance and the negative control, as well as differences in the use of gross and microscopic pathology terminology, limited statistical power and lack of compliance with Good Laboratory Practice standards. In consequence, the SAP published a White Paper [24] that highlighted inconsistencies and variability among studies and by so doing provided guidance for the development of a standardized study design to evaluate potential effects of EAC on development and sexual differentiation of X. laevis.
The purpose of this study was to develop, standardize, and refine a protocol meeting specifications of the 2003 SAP that could be used to assess estrogenic effects on development and sexual differentiation in X. laevis. A team of scientists representing endocrinology, toxicology, chemistry, pathology, and statistics was involved in the creation of the study design, selection of endpoints, and choice of data analysis tools. The team concluded that it was necessary to conduct two independent studies consecutively in two different locations to develop, refine and validate the study design by assessing effects of E2 on gonadal development of X. laevis exposed from NF (Nieuwkoop and Faber) stages 46–48 through stage 66 in a flow-through exposure system. Experiment 1A was conducted at Wildlife International Ltd, MD, USA (WLI) in order to develop and standardize conditions for rearing and handling X. laevis larvae, operation of a flow-through system and maintaining water quality and environmental parameters. A concentration range for E2 of 0.2, 1.5 and 6.0 μg L−1 was chosen to calculate the effective concentration needed to achieve a 50% effect (EC50) for gonadal feminization and to evaluate abnormal development of gonads. In addition, this preliminary study was used to establish the terminology for classifying gross and histological gonadal abnormalities, to cross-train laboratory personnel, and to develop the appropriate statistical analysis methodology. The procedures developed in Experiment 1A were refined and, implemented in Experiment 1B conducted at the Leibniz-Institute of Freshwater Ecology and Inland Fisheries, Germany (IGB).
Material and methods
Test organism
Tadpoles of X. laevis were obtained from the commercial supplier, Xenopus I, (Dexter, MI, USA). These larvae (approx. 1200 for each study) originated from ten pairs of adult breeders for each study. Xenopus I reared the eggs until day post fertilization (dpf) 3 and then shipped them to each laboratory.
Rearing
Upon arrival on dpf 4, tadpoles were acclimated to water and temperature conditions (21 ± 2 °C) over a six-hour time period. During the acclimation period larvae were maintained in two 40-L glass aquaria containing 30 L well water (ca. 600 larvae/tank, Experiment 1A) and three 40-L glass aquaria containing 30 L dilution water (ca. 400 larvae/tank, Experiment 1B) derived from the flow-through system. At the initiation of treatment, twenty-five larvae were randomly assigned to each glass aquarium (n = 200 larvae per treatment group) so that the biological loading rate was always less than 1g L−1 day−1. In Experiment 1A exposures were initiated at dpf 6 NF stage 45–46 and in Experiment 1B at dpf 8 when tadpoles had reached NF stage 47–48 as identified by visible hind limb buds [25]. Exposures continued until completion of metamorphosis or dpf 82, whichever came first. Throughout the acclimation period, water temperature was recorded twice daily whereas dissolved oxygen (DO) concentration and pH were measured daily. Experiment 1B was conducted in compliance with the local animal protection committee (G0048/05).
Feeding
Animals were fed Sera Micron (Lot No: 9698; Sera, Heinsberg, Germany), a commercially available food which was reported by the manufacturer to consist of 50.2% protein, 8.1% fat, 4.2% fibre, and 11.0% ash. A sample of Sera Micron was analysed by Lancaster Laboratories (Lancaster, PA, USA) for potential contaminants (metals and pesticides), and all tested contaminants were below their limits of quantification for analysis at that facility. A mixture of diet and dilution water was tested for estrogenicity and androgenicity using YES (yeast estrogen screen) [26] and YAS (yeast androgen screen) assays [27]. Based on these analyses no estrogenic or androgenic activity was detected.
During acclimation of larvae in Experiment 1A, 1 g Sera Micron was added to each tank daily. This amount of food was found to be excessive. In experiment 1B, after temperature acclimation, larvae were transferred to the acclimation tanks where they promptly attached to the glass aquaria walls. Therefore, no food was added on dpf 4. The following day some of the tadpoles were free swimming and 150 mg Sera Micron was added that evening by rinsing through a sieve. The amount of feed was chosen based on the intake of food by the larvae, and care was taken to ensure that feeding did not adversely affect water quality. From dpf 6 to 8 the free swimming tadpoles were fed 200 mg Sera Micron in each tank twice a day. Throughout the remaining course of both studies, tadpoles were fed a suspension of Sera Micron (77 mg L−1) three times daily. The total quantity of food added per tank each day increased from 300 mg/tank on dpf 8 to 900 mg/tank on dpf 33. The quantity of food added to each tank was adjusted periodically in Experiment 1A and on a daily basis in Experiment 1B to account for mortality and removal of tadpoles as they completed metamorphosis.
Study design
Because an almost completely feminized phenotype was observed in Experiment 1A in all treatment groups, the nominal E2 concentrations in Experiment 1B were adjusted to 0.015, 0.2, and 1.5 μg E2 L−1, respectively. Mid and high concentrations were selected to correspond to experiment 1A and allow comparison of results from the two laboratories.
Replicates
Each study consisted of at least one negative control group plus three E2-treated groups, each comprising eight replicated glass aquaria/treatment groups. Eight glass aquaria were arranged in two clusters of four tanks each. The clusters of negative control tanks and E2 treated tanks were randomly distributed within the environmental chamber in Experiment 1A and within water baths in Experiment 1B. Each cluster was supplied water via the same water flow meter and mixing chamber. Experiment 1A included 16 control tanks that were initially divided into two 8-tank groups designated as negative control and reference control, respectively. The animals from all treatment groups were handled in a blinded fashion so that the biologists and histologists involved in the experiments did not know the treatment applied to frogs being observed. An additional non-blind reference control group was included to provide reference information on rates of spontaneous incidence of gonadal abnormalities for blind histological assessment of gonads from the E2 treatment groups. Histological analysis was conducted on all frogs of Experiment 1A. At the end of Experiment 1A, it was discovered that slightly elevated temperature (1 °C) in one cluster of four negative control tanks had resulted in accelerated growth of the frogs in the tanks. Therefore, measurements of weight, length, gonad image area, and age at metamorphosis from those tanks were dropped from analysis, and the four remaining negative control tanks were combined with the reference control group to provide these measurements.
Flow-through system
Both studies were conducted using continuous flow-through systems operating at a flow rate of at least 50 L dilution water per tank per day. This was equivalent to a water exchange rate of approximately seven tank volumes per day. Test tanks consisted of 9-L glass aquaria (30 × 20 × 14.5 cm) containing 7 L water. Dilution water for Experiment 1A was obtained from a well of approximately 40 m depth located on the WLI site. In Experiment 1B the dilution water used was non-chlorinated municipal water. Dilution water samples from Experiment 1A were analysed for possible metal or pesticides by Lancaster Laboratories and water samples of Experiment 1B by Environmental Chemistry and Pharmanalytics (RCC Ltd., Itingen, Switzerland).
The flow-through systems in the two laboratories were essentially identical. Prior to use, the water was passed through several filters (5, 1, 0.5, and 0.45 μm) to remove particulate material, a UV-sterilizer, and a final particle filter (0.45 μm). After being filtered and temperature conditioned, the dilution water was divided into 16 streams. Each stream passed through a flow-control assembly, consisting of a rotameter flow gauge and a needle valve, to a mixing chamber. The flow from each mixing chamber was then split to supply a cluster of four tanks. Each tank received 140 mL min−1 dilution water (±5%).
Primary E2 solutions
The primary E2 stock solutions were continuously delivered to each mixing chamber (0.140 mL min−1), where it was vigorously mixed with the dilution water. Primary E2 stock solutions were metered into the mixing chambers by an eight-channel peristaltic pump (MasterflexL/S, Model 752455, Cole–Parmer, IL, USA). The flow-through system (temperature, water, filter pressure, flow-meter), and peristaltic pumps were calibrated four weeks prior to study initiation, and inspected visually twice daily throughout the study.
E2 solutions
Once every seven days a primary E2 stock solution was prepared in amber bottles by dissolving 6 mg E2 (Lot No: 11121AB, CAS: 50–28–2, purity 99.9%, Sigma Aldrich; WLI, Allentown, PA, USA; IGB: Taufkirchen, Germany) in alkaline solution (0.01 mol L−1 NaOH) applying moderate heating and constant mixing of the solution overnight similar to procedures employed by the Duluth Laboratory of the USEPA, Office of Research and Development (personal communication by J. Tietge). In Experiment 1A, the working stock solutions were prepared at nominal concentrations of 200, 1500, and 6000 μg E2 L−1, and in Experiment 1B at 15, 200, and 1500 μg E2 L−1. The quantity of NaOH in E2 test solutions delivered to each tank was negligible, and therefore no NaOH was added to the negative control.
Analysis of E2 concentration
Sampling
In Experiment 1A, water samples were collected from each tank in all treatment groups to confirm the operation of the diluter two days prior to study initiation. In Experiment 1B, E2 concentrations were monitored twice weekly for a two-week period prior to study initiation. After study initiation, water sampling was carried out twice weekly from alternating tanks in two tanks per cluster until study termination. Throughout the studies water samples were collected routinely two days after the freshly prepared E2 solution was connected to the flow-through system and a second time in the same week. At mid-depth from each test tank 190 mL water was collected using a glass pipette and 10 mL methanol (MeOH) was added. Of this solution 20 mL was decanted into a 20 mL amber scintillation vial for E2-enzyme-linked immunosorbent assay (ELISA) analysis. The remaining volume of 180 mL was stored as backup at 4 °C until E2 analysis was completed. Water samples from Experiment 1B were shipped to WLI for E2 analysis. The stability of E2 was confirmed by shipping samples of E2 (0.1 μg E2 L−1) that were prepared and analysed in parallel with the water samples.
Sample analysis
Rapid analysis was performed using ELISA kits (Ecologiena, Abraxis LLC, Warminster, PA, USA) to ensure the absence of the test substance in the negative control tanks and to verify the proper E2 concentration in each E2 tank. The ELISA method was based upon the E2 kit manual provided by the manufacturer. The method detection range of this assay was 0.05–0.5 μg E2 L−1. Parallel analyses of water samples for E2 concentration were performed by direct aqueous injection using an Agilent Series 1100 High-Performance Liquid Chromatograph coupled with a MDS Sciex API 3000 tandem mass spectrometer (LC–MS–MS) (Applied Biosystems, Foster City, CA, USA) and MDS Sciex API heated nebulizer ion source (Applied Biosystems) operated in the multi-reaction mode (MRM). Chromatographic separations were achieved using a Keystone Betasil C-18 column (50 mm × 2 mm, 3 μm particle size) (Thermo Hypersil, Madison, WI, USA). The method limit of quantitation (LOQ) for HPLC–MS–MS estradiol analyses was defined as 0.005 μg L−1.
Environmental conditions
Temperature
In Experiment 1A, all glass aquaria were situated in a walk-in environmental chamber (21 ± 2 °C). In Experiment 1B, glass aquaria were placed in temperature controlled water baths (22 ± 1 °C). In both experiments water temperatures were continuously recorded, and temperatures in all glass aquaria were measured weekly using a liquid-in-glass thermometer.
Light
The target light intensity at tank level was 100–500 lux and was measured weekly using a Sper Scientific (AZ, USA) light meter. A photoperiod of 16 h of light and 8 h of dark was used in Experiment 1A, and was adjusted to a 12 h light, 12 h dark cycle in Experiment 1B. A 30-min low light (around 30% intensity of full light) transition period was installed in both laboratories.
Water quality
Water quality parameters were assessed according to ASTM (2003) [28] and US EPA guidelines. Water hardness, pH, conductivity, and the concentrations of nitrate, and ammonia were measured weekly in one tank from each cluster throughout the experiments. Hardness and alkalinity were measured by titration based on procedures described in “Standard Methods for the Examination of Water and Wastewater” [29]. Nitrate concentrations were analysed with a Hach DR/700 colorimeter (Hach Company, Colorado, USA) and specific conductance was detected using a YSI conductivity meter (WLI: Yellow Springs Instruments Model 33, IGB: 3200, Yellow Springs, OH, USA). pH and ammonia concentrations were measured with an Orion pH/ISE Meter 720Aplus (Thermo Electron, USA). DO concentrations were measured in alternate tanks per treatment three times a week with an Orion Model 850plus dissolved oxygen meter (Thermo Electron). Commencing at dpf 43 in Experiment 1A, and at dpf 36 in Experiment 1B, gentle external aeration was applied by a glass pipette in each test tank to ensure maintenance of an adequate DO concentration.
Maintaining water quality
The tanks required daily cleaning to minimize microbial growth and to maintain adequate water quality. Bio-film that accumulated on the inner walls and bottom of the tanks was scraped off daily, and detritus was siphoned from the tanks. Daily cleaning proved inadequate to maintain appropriate water quality throughout the entire exposure period. Hence, in Experiment 1A all tanks were replaced with clean tanks weekly and in Experiment 1B twice throughout the experiment when siphoning did not lead to satisfactory results.
In-life assessment of tadpoles
Tadpoles were monitored for changes in general health, swimming behaviour and morphological appearance daily. The numbers of dead or moribund tadpoles and tadpoles completing metamorphosis were recorded daily for each tank.
Snout-to-vent length and weight
Each tadpole that had either completed metamorphosis by NF stage 66, or had failed to achieve metamorphosis by dpf 82, was removed from the test tank and euthanized by immersion in MS 222 (tricaine methanesulfonate) (WIL: Sigma Aldrich, Allentown, PA, USA; IGB: Sigma–Aldrich, Taufkirchen, Germany) (2 g L−1 buffered to pH 7.5 with NaHCO3). Immediately thereafter, the snout-to-vent length was measured to the nearest millimetre. Each frog was blotted dry and then weighed to the nearest milligram. A unique identifier was assigned to each frog that included: study number, in-life laboratory ID, colour code, replicate, and animal number.
Gross pathology
The pleuroperitoneal cavity was opened to expose the visceral organs. The liver, stomach, and intestine were examined using a dissecting stereo microscope (WIL: Olympus SZ61-TR, Olympus America, Melville, NY, USA; IGB: Olympus, SZX7, Olympus, Hamburg, Germany), and unusual variations in size, colour, and/or abnormal structure were recorded. The gastrointestinal tract was excised to expose the ventral surfaces of the kidney and gonads, which were similarly assessed for gross abnormalities. Tadpoles that did not complete metamorphosis at study termination (Experiment 1A: 25, Experiment 1B: 7) were staged according to NF and, the gonad morphology was assessed (if possible), and the carcasses were archived in fixative. Data derived from such animals were excluded from further analysis.
Assessment of gonad morphology
Using a dissecting stereo microscope, gross evaluations of gonad morphology were performed on all animals that reached NF stage 66. To enhance visualization of the gonads, which presented as thin, pale-tan strips of tissue on the ventromedial margin of the kidney, several drops of Bouin’s solution (WIL: Sigma Aldrich, Allentown, PA, USA; IGB: Sigma–Aldrich, Taufkirchen, Germany) were applied to partially fix these organs. Based on gross observation, the gonads of each animal were identified as testes, ovaries, or malformed. Examples of observed malformations included mixed sex, intersex, pearling and segmental aplasia. Mixed sex was defined as the co-occurrence of both ovarian and testicular tissue in a single gonad. Intersex was assessed as ovarian and testicular tissue in the same individual as separate gonads (left/right). The term pearling was characterized by the presence of multiple, prominent segmental enlargements and/or attenuations along the length of one or both gonads, whereas gonads that presented as nodular islands of testicular or ovarian tissue, with either intervening membranous connections or no connections at all, were classified as segmental aplasia. Gross findings for the gonads of each frog were recorded separately for the left and right gonads. The gonadal findings for each frog were verified independently by a second biologist. If the two biologists differed in their interpretations, the disparate opinions were discussed and a consensus finding was recorded.
Prior to the initiation of the first E2 experiment a list of standardized terms for changes of potential morphological gonadal observations was developed. Subsequent to the completion of both experiments the applied terminology was refined in order to better characterize and group different types of gonad changes observed and which in turn could be analysed. For example, the subcategories narrow, slightly narrow, truncated, slightly truncated, and margin entire were integrated into the main category hypoplasia. This term describes all gonads which appear smaller than those characterised as “typical/normal”. Since the term “typical/normal” is rather general, further sub-categories allow the biological variability of gonadal morphology to be described. However, during the investigations it was revealed that too few observations in each sub-category ruled out reasonable statistical analyses. Therefore, a refined list of qualitative descriptive terminology for gonadal features was created as presented in Table 1.
Table 1Glossary of terms used to describe features observed during gross morphological examination of gonads in X. laevis NF stage 66 FeatureDefinition1.Adhesion (other tissue)Gonads joined to other abdominal tissue2.Aplasia (agenesis)Complete lack of gonad development3.Segmental AplasiaGonad is longitudinally discontinuous 3-A. Tissue separationOne or more areas along the length of the gonad is poorly developed, attenuated with an essentially complete lack of gonadal tissue 3-B. Extraneous gonadal tissueA small disjunct cluster of gonad-like tissue4.BifurcationDivision of the gonad oriented longitudinally, along the cranial–caudal axis 4-A. Protuberance(s)Small projection(s) of the gonad along the margins 4-B. SymmetricA “Y” -shaped single gonad with relatively equal division5.Angular deformityOne or more bends in the gonad to an excessive degree 5-A. Gonad foldedThe gonad is bent over or doubled up so that one part lies on another6.DisplacedThe gonad, or a section of it, is not located in the typical location 6-A. Displaced laterallyThe gonad, or a section of it, is not located in the typical location. Instead it is located further away from the frog’s midline 6-B. Displaced mediallyThe gonad, or a section of it, is not located in the typical location. Instead it is located closer to or crosses over the frog’s midline 6-C. Displaced craniallyThe gonad is not located in the typical location. Instead it is located more towards the anterior 6-D. Displaced caudallyThe gonad is not located in the typical location. Instead it is located more towards the posterior7.FusedLeft and right gonads joined to varying degrees at one or more locations along cranial–caudal axis8.HypertrophySize of the gonad is larger than typical 8-A. WideGonad appears broad. Lateral width is large relative to the length 8-B. ElongatedGonad appears long (cranial–caudal). Length is large relative to the length of the kidneys and the torso of the frog 8-C. ThickGonad relatively increased in the distance between the dorsal and ventral surfaces 8-D. EnlargedA generalized increase in the size of the gonad along all axes9.Segmental hypertrophyA gonad where one or more areas are excessively large 9-A. Enlargement (mass)A section of the gonad substantially bigger in width and/or thickness than the remaining tissue 9-B. PearlingPronounced multiple segmental enlargement of the gonad 9-C. Partly thickA section or sections of the gonad is relatively increased in the distance between the dorsal and ventral surfaces10.HypoplasiaSize of the gonad is decreased 10-A. NarrowGonad appears strap-like. Lateral width is markedly reduced relative to the length 10-B. Slightly narrowLateral width is somewhat reduced relative to the length 10-C. TruncatedGonad appears short (cranial–caudal). Length is small relative to the length of the kidneys and the torso of the frog 10-D. Slightly truncatedGonad appears somewhat short (rostal–caudal). Length is somewhat reduced relative to the length of the kidneys and the torso of the frog 10-E. ThinGonad relatively reduced in distance between the dorsal and ventral surfaces 10-F. Margin entireAll of the ovary or indeterminate gonad has a smooth edge; no scalloping or lobes11.Segmental hypoplasiaA gonad where one or more areas are excessively reduced, attenuated, or poorly developed but not separated 11-A. Partly narrowA section or sections of the gonad has a small lateral width relative to the length and the width of the other sections 11-B. Partly thinA section or sections of the gonad is relatively reduced in distance between the dorsal and ventral surfaces 11-C. Margin slightly sinuateGonad has only shallow waved or scalloped edge (typically associated with a narrow or partly narrow ovary) 11-D. Margin partially entireA section or sections of the gonad has a smooth edge (typically associated with a narrow or partly narrow ovary) 11-E. PearlingPronounced multiple segmental attenuation or narrowing of the gonad12.IntersexOvarian and testicular tissue present as separate structures (i.e. not contained in the same gonad) (left/right)13.Mixed sexOvarian and testicular tissue present in the same gonad14.TranslucentGonad appears not so dense, light able to pass through diffusely 14-A. Tissue slightly translucentTo a small extent the gonad appears not so dense, a small amount of light able to pass through diffusely15.Segmental translucenceA section or sections of the gonad appears not so dense, light able to pass through diffusely16.Melanophores decreasedApparently fewer pigment-containing cells than typically seen. Applies only to females, mixed sex frogs, and frogs of undetermined sex 16-A. ApigmentationNo melanophores in the internal tissue. Applies only to females, mixed sex frogs, and frogs of undetermined sex17.Internal melanophoresPigment-containing cells within the gonadal tissue. Applies only to males, mixed sex frogs, and frogs of undetermined sex18.OtherAny finding not listed. Describe as needed. For any “other” finding that study personnel consider not incidental, notify the Study Director and other laboratories as soon as possible via email including a description and photograph
Photography
After gross gonad inspection was complete, each gonad was photographed in situ using a digital camera (WIL: Olympus DP12–2, Olympus America, Melville, NY, USA; IGB: Olympus DT5, Hamburg, Germany) attached to the stereo microscope. Each digitalized image included a millimetre measurement scale placed adjacent to the gonads to permit gonad size to be determined from the photograph. Following photography, each frog was placed in a labelled individual container of 30 mL of Bouin’s solution for approximately 48 h. At the end of the fixation period, each carcass was rinsed several times in 70% ethanol and was placed in 30 mL 10% neutral buffered formalin (WIL: VWR, Westchester, PA, USA; IGB: Histofix, Roth, Karlsruhe, Germany). A complete gonadal histological evaluation was performed on each frog of experiment 1A and results are presented in a separate paper (Wolf et al., in preparation). In some cases the growth and histological findings differ and a better understanding of results is gained by considering both types of observation.
Gonad measurement
The biological variability of morphological features of the gonads and the subjective nature of the assessment prompted us to develop a more quantitative metric of gonad size. Thus, individual measurements of gonad image area were derived from photographs of the gonads. Gonad image areas were obtained using image-processing software (Image Pro Plus, Version 5.1, Media Cybernetics, Silver Spring, MD, USA) that calculated the combined area of the left and right gonad whose outline had been manually traced around the digital image. The program was calibrated by the millimetre measurement scale prior to measurement for the gonad of each tadpole.
Statistical analysis
For both experiments, observations were collected for each animal. However, consistent with the experimental design, the tank is considered the primary experimental unit. When tank differences are present, statistical analysis should accommodate this experimental structure directly using nested “random effects” analyses or indirectly by analysing tank means or tank percentages. Neither approach is completely satisfactory for all endpoints. Nested models have specific distributional requirements and, in the case of incidence data, become unstable when frequencies are at or near zero. Analyses of tank means can also be problematic for these experiments. Estradiol-induced feminization causes the number of males per tank to be smaller for treated groups than for the control group. This, in turn, results in greater variation between estradiol-treated tanks than between control tanks. Most parametric and non-parametric statistical procedures, however, require equal variability between tanks within every group. When tank differences are small or absent, however, the individual data from all tanks within the same treatment group can be pooled and analysed more simply using the frog as the basic experimental unit.
For Experiments 1A and 1B, endpoints were tested for differences between tanks using a test of homogeneity. A chi-square homogeneity test was used for frequency endpoints (males, females, mixed sex) and a Kruskal–Wallis test was used for measurement endpoints (snout-to-vent length, body weight, age at completion of metamorphosis, gonad image area). This homogeneity test was calculated for each experimental group and the p-values combined using Tippett’s minimum-p method. Tank differences were found for several frequency endpoints in Experiment 1A. In addition, there was evidence that the presence of a humidifier caused a temperature-induced impact on measurement endpoints for two clusters of tanks. In contrast, Experiment 1B did not reveal any statistically significant (p < 0.05) tank effects. Consequently, all subsequent statistical analyses for Experiment 1B can be conducted on individual animals with all tanks pooled. Thus, analyses of Experiment 1A are quite complex and varied while Experiment 1B results can be analysed rather simply. In order to summarize the results consistently to allow comparisons of both experiments, only animal-level means and standard deviations are presented. The means and standard deviations given in this report for both experiments were calculated from individual animals using Excel Version 2000. For Experiment 1B these simple statistics are identical with those obtained from the statistical analyses described below. For Experiment 1A these means (and standard deviations) should be considered only approximate. They may differ slightly from means obtained using more sophisticated analyses that incorporate a tank structure.
In Experiment 1B, the estradiol-treated groups were compared with controls in a step-wise manner that preserved power and protected against excessive false positives. An overall test of group differences and a test for trend with dose were both conducted. If either of these two overall tests showed statistical significance at the 5% level (i.e., p < 0.05) then comparisons of each dose group to the control were made. For frequency endpoints a chi-square homogeneity test and a Cochran–Armitage test for trend were used for the overall homogeneity and trend tests, respectively. Comparisons of dose groups to control were conducted using Fisher’s exact test. For measurement data essentially all endpoints showed statistically significant deviations from normality (p < 0.01, Shapiro–Wilk test). Consequently, a Kruskal–Wallis nonparametric test was used for group homogeneity and a Jonckheere–Terpstra test was used for trend. Comparisons of dose groups to control were conducted using the nonparametric Wilcoxon–Mann–Whitney U-test. All statistical analyses of Experiment 1B were calculated using Release 9.1 of the SAS statistical package (2004).
Results
It is emphasized that Experiment 1A was performed first in order to standardize conditions for rearing X. laevis and for operation of the flow-through system. Improvements of these procedures were implemented in Experiment 1B resulting in optimised growth and development conditions for X. laevis tadpoles. Those parameters that were successfully achieved during Experiment 1B are summarised in Table 2 and may be used as guidance for similar experiments. The following sections focus on results obtained in Experiment 1B, presenting results from Experiment 1A, as needed, to allow comparison of the two experiments.
Table 2Final study design specifications. Brief summary of the most important parameters for exposing X. laevis tadpoles in a flow-through exposure systemParameterCharacteristicWater-quality parameters  Supply water: filtered, contaminant freeFiltered and UV sterilized; tested for chemical contaminants  Dissolved oxygen≥60% of saturation  pH7.9 to 8.3  Ammonia0.0 to 0.35 mg L−1  Nitrate0.04 to 2.29 mg L−1  Hardness90 to 130 mg CaCO3 L−1  Alkalinity133 to 172 mg CaCO3 L−1  Specific conductance726 to 817 μg S cm−1  Supplemental aerationInitiated not later than dpf 36Animal husbandry  Animal supplyXenopus, Dexter, MI, USA  Shipped from SupplierDay 3 post fertilization (dpf)  Acclimationdpf 4–7  Treatment perioddpf 8–82  Feed: Sera MicronSera Micron (contaminants and estrogenic potential evaluated) Contents: 50% protein, 8% fat; 4% fibre; 11% ash.  Feeding rate—acclimationdpf 4–5; per 300 larvae: 200 mg; if free swimming, dpf 6–8: 200 mg twice daily (based on the intake of food by the larvae)  Feeding rate—NF stages 46–6677 mg L−1: fed three times daily, 300 mg/tank from dpf 8, increasing to 900 mg/tank at dpf 33, then adjusted during experiment as animals were removed (dead, froglets)  Exposure chamber cleaningDaily  Exposure chamber replacementAs needed  Water-quality parametersHardness, ammonia, nitrate, pH, specific conductivity twice a week; DO three times a week  Temperature22 ± 1°C (verified twice daily, measured once a week)  Light intensity100–500 lux (measured once weekly, at level of surface of water)  Light:dark cycle12:12; 30-min transition (30% of full light)Experimental design  Performed in compliance with Good Laboratory Practice StandardYes  Colour coded (blind) treatmentYes  Randomization of tadpolesTadpoles distributed in rotation one or two at a time to transport vessels to contain 25, then vessels were randomly assigned to test tanks  Randomization of tank clustersCompletely randomized placement; one cluster consists of four tanks  Negative controlWater (eight replicates, two clusters)  Treatment 17β-estradiol concentrations0.015, 0.2, 1.5, 6.0 μg E2 L−1; eight replicates; two clusters for each group  Tank volume9 L (water volume 7 L), tanks covered  Tank flow rate50 L day−1  Number of animals200 larvae per treatment group  Animal load densityLess than 1 g L−1 water; 25 larvae/tank  HPLC–MS–MS verificationdpf 8 all tanks; dpf 16–81 weekly (alternating tanks)
Environmental conditions
In both experiments experimental conditions generally met the criteria listed in Table 2. In Experiment 1B the mean and SD for hardness and alkalinity were 103 ± 5.8 mg and 162 ± 6.8 mg CaCO3 L−1, respectively. Specific conductance was 794 ± 18.0 μS cm−1, pH 8.1 ± 0.09, nitrate 1.6 ± 0.5 mg L−1 and ammonia 0.08 ± 0.09 mg L−1. The measurements from negative control tanks were not significantly different from those from the E2 tanks. Technical problems in Experiment 1B prevented accurate determination of DO concentration, however, based on subsequent studies using the same procedures, DO was consistently above 60% saturation.
E2 exposure verification
In each experiment, E2 was not present in any of the negative control tanks. Measured concentrations of E2 stock solutions ranged from 96.6 to 106% in Experiment 1A and between 90.5 and 108% in Experiment 1B. In Experiment 1A, the mean concentration of E2 for three treatment groups began to decline after experiment initiation and fell below 80% of nominal E2 concentration by dpf 29. Between exposure dpf 43 and 50 mean concentration of E2 decreased to minima of 0.4, 43, and 38% of nominal concentrations of E2 in the 0.2, 1.5, and 6.0 μg E2 L−1 treatment groups, respectively. After dpf 50, concentrations of E2 began to increase in all treatment groups, and exceeded 60% of nominal concentration on dpf 64 (Fig. 1a). In Experiment 1B, concentrations of E2 declined after study initiation and dropped below 60% of nominal E2 concentration by dpf 23. On dpf 37 concentrations of E2 fell to 47, 46, and 50% in the 0.015, 0.2, and 1.5 μg E2 L−1 treatment groups, respectively. The concentration of E2 increased to 60% of nominal concentration on dpf 65 (except the 1.5 μg E2 L−1 with 44.4%) (Fig. 1b).
Fig. 1Comparison of measured mean E2 concentrations in water samples of E2 treatments in the flow-through system. Mean values represent the concentrations of E2 of weekly sampling intervals from dpf 6 or 8 through the last sampling interval on termination of the experiments on dpf 82. a, E2 concentrations of Experiment 1A; b, E2 concentrations of Experiment 1B. The limit of quantification was 0.00500 μg E2 L−1, calculated as the product of the concentration of the lowest calibration standard (0.100 μg E2 L−1) and the dilution factor of the matrix blank samples (0.0500 μg E2 L−1)
In-life assessment of tadpoles
Appearance and behaviour
No effects of E2 treatment on tadpoles condition, behaviour, or external appearance was observed in any of the treatment groups at any time throughout the entire course of Experiment 1A or Experiment 1B.
Survival
Excellent survival was achieved in both experiments in controls as well as in E2 treated groups. The total survival was 92% in Experiment 1A and 99% in Experiment B. In Experiment 1A, survival in the combined negative control, 0.2, 1.5, and 6.0 μg E2 L−1 treatment groups was 90.8, 94.5, 95.0, and 88.0%, respectively. In Experiment 1B survival rates were 100% in the negative control and 97.5, 97.5, and 99.5% in the 0.015, 0.2, and 1.5 μg E2 L−1 exposure groups, respectively.
Snout-to-vent length and weight
In Experiment 1A, for females the average length values were, 22.2 ± 0.1, 22.1 ± 0.1, and 22.8 ± 0.1 mm in the 0.2, 1.5, and 6.0 μg E2 L−1 treatment group compared with 22.0 ± 0.2 mm in the combined control group. Males treated with E2 had an average snout-to-vent length of 21.8 ± 0.1, 22.1 ± 0.6, and 23.0 ± 0.1 mm in the combined control, 1.5, and 6.0 μg E2 L−1 treatment group, respectively. The increase in the mean weight of males was more pronounced than the increase in the mean weight of females (Fig. 2).
Fig. 2Body weight (means ± SD) of a, females, and b, males, after exposure to E2 from dpf 6 or 8 in a flow-through system until completion of metamorphosis or day 82 post-fertilization. 0.015 μg E2 L−1 was not used in Experiment 1A and 6.0 μg E2 L−1 was not used in Experiment 1B. Neg Ctrl, negative control group
In Experiment 1B no statistically significant differences between the snout-to-vent length for females and males were observed. The snout-to-vent length for females and males in the negative control group was 17.2 ± 1.5 and 17.0 ± 1.6 mm, respectively. For females and males treated with 1.5 μg E2 L−1 snout-to-vent length was 17.5 ± 1.5 and 16.3 ± 1.3 mm, respectively. No statistically significant treatment-related effect was observed for body weight either in females or in males. At completion of metamorphosis, frogs from Experiment 1A weighed more than twice as much as frogs from Experiment 1B, regardless of the gender (Fig. 2).
Time course of metamorphosis
On average, the process of metamorphosis for tadpoles took longer in Experiment 1A than in Experiment 1B. In Experiment 1A, the age at which the first frog in the experiment completed metamorphosis was delayed by 6 days (dpf 44) in each E2 treatment group compared with animals in Experiment 1B (dpf 38) (Fig. 3). In Experiment 1A exposure of females to 0.2, 1.5, and 6.0 μg E2 L−1 resulted in a slight delay of metamorphosis. In Experiment 1B, the time course of metamorphosis for females in the 0.015 μg E2 L−1 treatment group was similar to that for females in the negative control group but exposure to 0.2 and 1.5 μg E2 L−1 caused a delay of completion of metamorphosis (Fig. 3a,b). A comparable picture for the cumulative time course of metamorphosis of males is shown in Fig. 3c,d. In Experiment 1A, males treated with 0.2 and 1.5 μg E2 L−1 had a delay in completing metamorphosis of approximately ten days compared with males in the negative control group. In Experiment 1B metamorphosis of male frogs in the 0.2 μg E2 L−1 treatment group was slightly postponed whereas for males treated with 1.5 μg E2 L−1 completion of metamorphosis was delayed by three days. The lowest E2 concentration used obviously did not affect the mean time to complete metamorphosis.
Fig. 3Cumulative portion of X. laevis frogs completing metamorphosis after exposure to different concentrations of E2 from day 6 or 8 post-fertilization in a flow-through system until completion of metamorphosis or day 82 post-fertilization. The number of individuals per treatment group (n) is given in the legend. a, portion of females in Experiment 1A; b, portion of females in Experiment 1B; c, portion of males in Experiment 1A (just one male was observed in the 6.0 μg E2 L−1 group); d, portion of males in Experiment 1B. Neg Ctrl, negative control group
Exposure to E2 had no effect upon the percentage of surviving tadpoles that completed metamorphosis on dpf 82 (97.4% in Experiment 1A, 99.1% in Experiment 1B). By dpf 82, 25 (out of 920) and 7 (out of 789) individuals failed to complete metamorphosis in Experiment 1A and Experiment 1B, respectively, and no treatment-related effect was observed. In Experiment 1B, 1, 2, and 4 individuals failed to complete metamorphosis by dpf 82 in the negative control, 0.2, and, 1.5 μg E2 L−1 treatment groups, respectively.
Mean age at completion of metamorphosis
The differences in the time course of metamorphosis described in the previous section are reflected in the mean ages of frogs at completion of metamorphosis in each of the treatment groups. In Experiment 1A, the mean ages at completion of metamorphosis of females treated with 0.2, 1.5, and 6.0 μg E2 L−1 were dpf 64.2 ± 7.6, 62.9 ± 7.9, and 65.1 ± 9.2, respectively, and substantially greater than the mean age in the negative control group (dpf 61.4 ± 8.2) (Fig. 4). In contrast, in Experiment 1B females treated with the different E2 concentrations showed a dose response with a significant increase noted in the 0.2 μg E2 L−1 group (dpf 50.7 ± 7.5) and a greater effect (dpf 51.5 ± 6.6) in the 1.5 μg E2 L−1 treatment group compared with dpf 47.8 ± 5.6 for females in the negative control group. For males in Experiment 1A, the mean age at completion of metamorphosis was increased from dpf 60.6 ± 9.4 in the negative control to dpf 70.9 ± 9.3 and dpf 72.3 ± 9.0 in the 0.2 and 1.5 μg E2 L−1 treatment groups, respectively. In Experiment 1B, the mean age at completion of metamorphosis of males was statistically significant greater in the 0.2 μg E2 L−1 (dpf 51.2 ± 6.6) and 1.5 μg E2 L−1 (dpf 58.0 ± 8.3) treatment groups than in the negative control group (dpf 48.5 ± 6.8) (Fig. 4).
Fig. 4Mean age at completion of metamorphosis (mean ± SD) by treatment group for female, a, and male, b, frogs at completion of metamorphosis or on day 82 post-fertilization. 0.015 μg E2 L−1 was not used in Experiment 1A and 6.0 μg E2 L−1 was not used in Experiment 1B. Only one male frog was observed in the 6.0 μg E2 L−1 group in Experiment 1A. Neg Ctrl, negative control group. For Experiment 1B significant differences are marked by asterisks (*p < 0.05, **p < 0.01, ***p < 0.001; two-sided Wilcoxon–Mann–Whitney test of equality with negative control group). A Kruskal–Wallis test of overall group differences was significant (females p < 0.0001, males p = 0.0009) and a Jonckheere–Terpstra test of trend was significant (females p < 0.0001, males p = 0.0062)
Appearance and examination of non-gonadal organs
No externally visible morphologic abnormalities were observed in each treatment group in either experiment. Likewise, inspection of the liver, stomach, intestine, kidney for tumours, lesions or any other remarkable features had no effect.
Determination of sex, percent female, male, and mixed sex by gross examination
The percentages of male, female, and mixed-sexed animals were calculated based on frogs that completed metamorphosis (NF stage 66) among all animals. In Experiment 1A, the percentage of female phenotype was increased in all E2 treated groups. In the 0.2, 1.5, and 6.0 μg E2 L−1 treatment group 88.4, 95.8, and 98.3% of frogs showed female phenotype. In Experiment 1B, the incidence of female phenotype in the E2 treated groups displayed a statistically significant distinct dose-responsive increase with a high percentage of females observed for the 0.2 (70.5%) and 1.5 μg E2 L−1 (92.3%) treatment groups, respectively (Fig. 5). The E2 concentration that results in an increase of female phenotype to 75% is defined as the effective concentration (EC50). Based on the feminization effect observed in Experiment 1B the EC50 was calculated to be 0.12 μg E2.
Fig. 5Individual bars represent frequencies of sex of X. laevis frogs among all animals according to treatment group based on gross morphological evaluation; the total number of frogs (n) that completed metamorphosis, is given on top of the bars. Tadpoles were exposed to E2, in a flow-through system, from day 6 or 8 post-fertilization through completion of metamorphosis or day 82 post-fertilization. The EC50 was calculated from the results of Experiment 1B and was 0.12 μg E2. The numbers of frogs in the negative control in Experiment 1A was higher for evaluation of gender because the animals which were excluded from calculation of weight and length were included in calculation of the portion of males, females, and mixed sex. Neg Ctrl, negative control group. For Experiment 1B asterisks denote significant differences (**p < 0.01, ***p < 0.001; two-sided Fisher’s exact test for comparison to negative control). Both an overall exact chi-square test of equality of all experimental groups and a Cochran–Armitage trend test were statistically significant (p < 0.05) in every case
In Experiment 1A, the percentage of mixed sex individuals were 6.9% (13/189), 2.1% (4/191) and 1.1% (2/176) in the 0.2, 1.5, and 6.0 μg E2 L−1 treatment groups, respectively. In Experiment 1B, the incidences of mixed-sex animals were statistically significant with 6.2% (12/193) and 4.1% (8/195) in the 0.2 and 1.5 μg E2 L−1 treatment groups, respectively (Fig. 5). In both experiments the percentage of frogs displaying mixed sex gonads was higher in E2-treated groups than in the negative control groups. When results from the two experiments are compared, the percentage of mixed sex animals was seen to be similar in equal E2 concentrations.
Segmental aplasia was detected in Experiment 1A in one of nine male individuals out of the 0.2 μg E2 L−1 treatment group. In Experiment 1B, segmental aplasia was observed in four of 101 (4.0%), nine of 92 (9.8%), eighteen of 45 (40.0%) and four of seven (57.1%) male frogs in the 0.015, 0.2, and 1.5 μg L−1 treatment groups, respectively. Thus, the frequency of segmental aplasia in the 0.2 and 1.5 μg E2 L−1 groups was significantly elevated. Pearling was observed in two of nine and one of five male frogs the 0.2, and 1.5 μg E2 L−1 groups, respectively, in Experiment 1A. In Experiment 1B, the occurrence of pearling was determined in two frogs out of 45 within the 0.2 μg E2 L−1 treatment group. Intersex was not found in any of the 1593 Xenopus evaluated for gonad morphology in both studies. Figure 6a,b demonstrate ordinary developed ovarian and testes tissue that was assumed to be “normal gonad morphology” for gross morphological evaluation. Figure 6c to e illustrate gonads showing mixed sex, pearling, and segmental aplasia and clarify the terminology which was applied in these experiments to describe morphological changes in gonadal appearance of X. laevis due to E2 treatment. To allow the variability of the morphological appearance of gonads to be described, an extensive terminology to include the biological variability was developed and used in the two experiments (Table 1). In Experiment 1A there were apparent differences in the occurrence of ovaries that were truncated, narrow entire margin, or had enlarged masses. However, no differences in morphological characteristics were observed between E2-treated and control frogs in Experiment 1B.
Fig. 6a. Ovary of X. laevis revealing segmented, lobular structure and clearly visible internal melanocytes in the tissue; the length of the ovary is slightly shorter than the length of the kidney. b. In contrast, a testis is shorter relative to the length of the kidney (half of the kidney), straight shaped, and does not contain internal melanocytes. In general testis tissue appears much denser than ovary tissue. c. Mixed sex gonads were defined as the co-occurrence of both ovarian and testicular tissue in a single gonad. The ovary tissue in the presented gonad shows a reduced lobular structure, just a few melanophores, whereas the testis tissue appears bulbous. Arrows indicate the different gonad structures, f, ovary-like tissue, and m, male like tissue. d. Pearling (p) was characterized by the presence of multiple, prominent segmental enlargements and/or attenuations along the length of one or both gonads. e. Segmental aplasia (segA) was classified as nodular islands of testicular or ovarian tissue, with either intervening membranous connections or no connections. The kidney and the attached gonad tissue were fixed with Bouin’s solution
Gonadal measurement
Because of the subjective nature of gross gonad morphological assessment, quantitative measurements of gonad area were performed on digital images.
In Experiment 1A, mean gonad areas of female and male frogs were approximately twice as large as in animals of the corresponding gender in Experiment 1B. The mean gonad image areas for females in Experiment 1A were 2.20 ± 0.82 mm2, 2.21 ± 0.89 mm2, 2.09 ± 0.92 mm2, and 2.23 ± 0.92 mm2 in the combined negative control, 0.2, 1.5, and 6.0 μg E2 L−1 treatment groups, respectively. In Experiment 1B, for females treated with 0.015, 0.2, and 1.5 μg E2 L−1 mean gonad image areas were 1.21 ± 0.41 mm2, 1.27 ± 0.44 mm2, and 1.32 ± 0.41 mm2, respectively, compared with 1.23 ± 0.42 mm2 in the negative control group. Mean gonad areas in females treated with 1.5 μg E2 L−1 were significantly increased (p < 0.05). In Experiment 1A, the measured mean gonad areas of males in the combined control, 0.2 and 1.5 E2 L−1 were 1.28 ± 0.2 mm2, 1.39 ± 0.2 mm2, and 1.70 ± 0.3 mm2, respectively. In Experiment 1B, males treated with 0.015, 0.2, and 1.5 μg E2 L−1 had a slight increase of mean gonad areas 0.82 ± 0.20 mm2, 0.85 ± 0.19 mm2, and 0.89 ± 0.14 mm2, respectively, compared with the negative control (0.77 ± 0.22 mm2). A Kruskal–Wallis test of overall group differences was marginally non-significant (p = 0.079) but a Jonckheere–Terpstra test of trend was significant (p = 0.0092).
Discussion
Gonadal development in anuran amphibians has been proposed as a valuable model to evaluate the physiological consequences of exposure to EAC with (anti)estrogenic and/or (anti)androgenic activity [30]. It has long been known that treatment of tadpoles with gonadal steroids can disrupt normal gonadal development in several anuran species and might even lead to complete sex reversal [31, 32]. However, studies on the potential effects of EAC on anuran gonadal development have produced variable and, in some cases, inconsistent results. From analysis of published studies it appears that some of the variability between EAC studies may be because of differences between the experimental conditions employed. One example to illustrate this situation is the outcome of three different exposure studies assessing the effects of bisphenol A (BPA) on gonadal development in X. laevis [15–17]. By using a semi-static exposure system, feminization of tadpoles due to BPA treatment was reported by Kloas [15] and Levy [17]. Attempts to replicate these findings in a flow-through exposure system by Pickford [16] were not successful. A more recent example includes inconsistent results observed in different studies on the putative effects of atrazine on survival, growth, metamorphosis, and sexual development of X. laevis [20, 21, 23]. The failure to reproduce the original findings leaves regulatory authorities with many uncertainties and in consequence, may prevent sound assessment of potential risks associated with the use of the chemicals of concern.
Contradictory results from several studies of atrazine, and the criticism associated with each of these studies, made it impossible for the US EPA to perform an unambiguous evaluation of the actual impact of atrazine on development and sexual differentiation in X. laevis. The US EPA therefore requested further independent studies to be conducted using a standardized test protocol including validated positive control treatments and an optimized study design [33]. To address uncertainties resulting from inconsistencies of previously reported data, the current multi-site study was undertaken in order to develop and refine an optimized testing protocol for the use in developmental studies with X. laevis tadpoles. The current paper presents the results of the two independent exposure experiments that were performed by two different laboratories. A first experiment (Experiment 1A) served as a pilot study to identify appropriate testing conditions for the development of a standardized testing protocol. Accordingly, the test protocol used in the follow-up experiment (Experiment 1B) included some modifications reflecting the experience gained in Experiment 1A. In combination, the current multi-site study aimed at a characterization of the general performance of the tadpoles using the optimized testing protocol and a detailed evaluation of concentration-dependent effects of E2 on tadpole growth, larval development and sexual differentiation.
The testing protocol under investigation specified the use of a flow-through exposure system to achieve E2 exposure of tadpoles. The main advantages of flow-through systems versus static renewal systems are the maintenance of constant testing conditions (e.g., water quality) and the continuous delivery of the test substance to the exposure vessels to achieve stable exposure of the test species to the test substance. In the two exposure experiments with E2, however, difficulties were encountered with regard to the maintenance of the target test concentration of E2 even under flow-through conditions. Although E2 concentrations in the stock solutions were within acceptable ranges in both exposure studies, and the quantity of E2 metered into the tanks was correct, actual E2 concentrations in the exposure tanks decreased to levels well below nominal test concentrations during the entire course of the experiments (Fig. 1). These observations clearly indicate that in future studies care must be taken to maintain constant nominal concentrations of the test compound throughout exposure. It is recommended the flow rate of the stock solutions be adjusted on a daily basis according to the actual test substance concentration present in the tanks.
Comparison of the time-course of E2 concentration changes with estimates of the biomass present in E2 treatment tanks at different times during Experiments 1A and 1B suggests that E2 concentrations were inversely related to changes in total biomass. Aqueous E2 concentrations decreased concurrent with the increase in body size and body weight of tadpoles during premetamorphosis (NF stages up to 52), prometamorphosis (NF stages 53 to 57), and climax (58 to 66). Concentrations of E2 began to increase as tadpoles that completed metamorphosis were removed from the test tanks, and as remaining tadpoles began to lose weight during tail resorption. Figure 7 illustrates that E2 concentrations were inversely related to changes in total biomass of the tadpoles. Furthermore, it should be noted that the biomass was greater in Experiment 1A than in Experiment 1B and the decline in aqueous E2 concentration was more marked in Experiment 1A. These observations, combined with the relatively high octanol–water partition coefficient of E2 (log Kow 3.94) suggest that additional measures (e.g., gradual increase of amounts of E2 delivered to the test tanks) are required to maintain a constant exposure to E2 in terms of actual aqueous E2 concentrations. In semi-static systems E2 disappears quickly, which suggests the animal may absorb the E2. Only few studies are available that have addressed this aspect. Results from Hayes and Licht [34] suggest that exogenous E2 is readily metabolized by X. laevis and other anuran species, and there was no evidence of bioaccumulation of E2 in anuran tadpoles. However, metabolic enzymes and the sites of E2 metabolism in anuran tadpoles are largely unknown. Possible candidate metabolites might include various E2 conjugates, but neither the test solutions nor the test animals have been examined analytically for the presence of such metabolites. To the best of our knowledge, there are also no data available regarding possible developmental profiles of E2-metabolising capacities that might explain the transient reduction of aqueous E2 concentrations in the present studies. It is, further, not known whether any of the putative metabolites still possess estrogenic or other hormonal activity. Finally, it is possible that factors other than changes in tadpole biomass may have contributed to the observed changes of E2 concentrations. For example, a small fraction of E2 may have been absorbed either in the biofilm or directly on the surface of the glass tanks. Another possibility could be that E2 partition from water to fat in the administered Sera Micron feed (~8% fat) occurred. Notably, daily food amounts administered to each tank were adjusted to total larvae biomass. Therefore, it is an interesting possibility that at least a portion of the E2 dose may reach biological targets in the larvae as a result of oral ingestion of the feed.
Fig. 7Correlation of changes in E2 concentration and the biomass of growing tadpoles during the exposure period in Experiments 1A and 1B
An important objective underlying the development of a standardized testing protocol for exposure of X. laevis tadpoles was to maintain appropriate water quality throughout the exposure experiments. Thus special care was taken to standardize and refine water quality parameters to ensure optimal growth and development of X. laevis tadpoles. Except for a transient decline of DO concentration below desirable levels, all water-quality parameters were within acceptable limits. The decrease of DO concentrations through dpf 43 in Experiment 1A and inability to maintain the DO concentration above 60%, may have been due to insufficient cleaning. Based on this observation, the cleaning procedure was enhanced, and additional aeration was provided by gently bubbling air through glass pipettes in each tank.
An interesting finding in these two experiments was the effects on growth and development of X. laevis tadpoles that occurred due to slight differences in rearing procedures between Experiment 1A and Experiment 1B. Tadpoles in Experiment 1A required more time to complete metamorphosis and the metamorphic froglets (NF stage 66) were larger and heavier than those in Experiment 1B. Whereas the relationship between age and size is in accordance with established amphibian life-history models [25], the factors that are responsible for the inter-laboratory differences are not entirely clear in this case. In Experiment 1A, the photoperiod was 16 hours light:8 hours dark compared with 12 hours light:12 hours dark in Experiment 1B. The mean temperature was approximately 1 °C lower and more variable (±1 °C) in Experiment 1A. The differences in photoperiod and water temperature are certainly prime candidates to explain the slightly slower development in Experiment 1A. On average, tadpoles in Experiment 1A consumed food 13 days longer than those in Experiment 1B. Subsequently, the longer food consumption of tadpoles in Experiment 1A led to greater weight compared with that of the froglets in Experiment 1B. Based on the results achieved in experiments 1A and 1B, and to consider the reported differences on growth and development of X. laevis in the literature, it is obviously important to standardize the test systems.
Results from Experiment 1B indicated that most of the tadpoles (92%) from the negative control group completed metamorphosis within a relatively short time (18 days), whereas metamorphosis was delayed (26 days) in a few tadpoles only (7.5%; 15/199). In contrast, time to metamorphosis in Experiment 1A, was much more variable, because 78.5% (216/275) of the tadpoles completed metamorphosis within twenty-two days and 21.5% (59/275) in thirteen days. Furthermore, tadpoles in Experiment 1A completed metamorphosis later than in Experiment 1B (Fig. 8). Thus, the use of the improved testing protocol in Experiment 1B was associated with a reduction in the time needed for completion of metamorphosis, which resulted in much more homogenous development of the test organisms. The latter finding is particularly important given the fact that still only few data are available upon which one could build a set of performance criteria to assess the validity of exposure experiments with X. laevis tadpoles. Homogenous growth and development in the negative control group is certainly one important aspect that deserves special attention. Generally, it can be proposed that homogenous rates of tadpole growth and development will reflect appropriate culture and testing conditions. Furthermore, any reduction in the inter-individual variability of growth and development will enable more sensitive assessment of the possible effects of a test compound on the development of tadpoles. To this end, the robust demonstration of altered rates of development due to test compound exposure will also provide additional information for interpretation of the effects of the test compound on gonadal development (e.g., identification of systemic toxicity).
Fig. 8Proportion of male and females frogs (n) of the negative control groups completing metamorphosis during the course of the studies, based on the total number of animals surviving in Experiments 1A and 1B
In contrast, it is still difficult to assess the validity of an exposure experiment based on the time required for completion of metamorphosis. According to Nieuwkoop and Faber [25] the average time to reach stage 66 is dpf 58. In Experiment 1B, the average time to stage 66 was still shorter, being dpf 47.8 ± 5.6 and 48.5 ± 6.8 for female and male frogs, respectively. This reduction of the time needed for tadpoles to complete metamorphosis is likely a result of the various efforts to optimize an array of biotic and abiotic factors known to affect tadpole growth and development.
An important element of experimental approaches used to assess the estrogenic activity of chemicals is the choice of an appropriate positive control treatment to confirm the validity of the experimental system. Therefore, a major focus of the current multi-site study was to establish a concentration–response relationship for E2 treatment on sexual differentiation of X. laevis. The results from Experiment 1A were not appropriate to demonstrate this relationship because all E2 concentrations (0.2, 1.5, and 6.0 μg E2 L−1) employed in Experiment 1A caused a very high percentage of phenotypic females. Based on these observations, the E2 test concentrations were modified in Experiment 1B, and the use of E2 concentrations of 0.015, 0.2 and 1.5 μg E2 L−1 revealed a clear concentration-dependent increase in the number of phenotypic females. The data from Experiment 1B allowed the calculation of an EC50 of 0.12 μg E2 L−1 resulting in 75% gonadal female phenotype.
Literature information about E2 concentration-dependent effects on growth, development, and sexual differentiation in X. laevis is rare. In most studies single E2 concentrations ranging from 0.3 to 100 μg L−1 were used to study effects on gonadal development [15–17, 22, 31, 35, 36] and resulted in percentage of female phenotype between 60, and 100%. Hu et al. [14] applied 1, 10, and 100 μg E2 L−1 and reported an increase of female phenotype to 66, 99, and 100%, respectively. Rearing X. laevis larvae under optimized conditions in Experiment 1B demonstrated a distinctive sensitivity for gonadal differentiation. This is demonstrated by a clear concentration-responsive effect, and a distinct effective concentration of 0.2 μg E2 L−1, as shown by the frequency of female/male individuals and the occurrence of substantial mixed sex and other abnormalities. Comparing the results of Experiments 1A and 1B regarding gonadal effects of 0.2 and 1.5 μg E2 L−1, a similar effect was observed, indicating high reproducibility of the testing protocol. Based on the endogenous E2 concentrations measured in adult X. laevis being approximately 28 ng mL−1 for females and 3 ng mL−1 for males, [37] and in froglets approximately 0.002 μg E2/fresh body weight [38] the calculated EC50 appears to represent a physiological E2 concentration in X. laevis.
Another question of high interest is how estrogenic EAC might impact gonad gross morphology. While determination of phenotypic sex was a core endpoint in the assessment of the effects of E2 on sexual differentiation, a number of other gross abnormalities were also evaluated. To that end, a glossary of gonadal gross morphological terminology and criteria was included in the study protocol (Table 1). The importance of this glossary cannot be understated, as a substantial degree of diagnostic variability, within and among studies, can be attributed to differences in terminology [39]. For example the word intersex is characterized differently by several authors and is often used interchangeably with the term mixed sex, whereas in the glossary, these terms have different meanings. Therefore, it is difficult to describe morphological changes triggered by EAC in a way that is clear and generally accepted. Furthermore, assessments of gonadal morphology involve subjective judgement, which make studies difficult to interpret and compare. In the current multi-site study the uncertainty arising from the subjective judgement that was necessarily involved in evaluation of the gross morphological assessment each gonad was reduced by having each finding verified by a second biologist. Additionally, digital photographs of each gonad were taken to allow later verification of assessments, and to allow measurement of gonadal image area using image-processing software. Measurement of gonad image area enables objective verification of some of the subjective size-related observations made during the gross assessments. Trend analysis revealed for both sexes in Experiment 1B an increase of mean gonadal area that coincides with an increased age at completion of metamorphosis.
In both experiments, gross morphological evaluation of the gonads revealed the presence of mixed sex, segmental aplasia, and pearling at 0.2 μg E2 L−1. It is noteworthy that intersex, characterized as ovarian and testicular tissue in the same individual as separate gonads (left/right), was not found in any of the 1593 frogs evaluated in these two experiments. Interestingly, in both experiments the highest incidence of these testicular abnormalities was observed in tadpoles exposed to 0.2 μg E2 L−1. This was most likely due to the greater degree of complete sex reversal at the two higher E2 concentrations, as suggested by the significantly increased proportion of phenotypic females compared with phenotypic males in the 1.5 and 6.0 μg E2 L−1 treatment groups. Based on gross morphological evaluation no differences for the terms of main-categories such as hyper- or hypoplasia were observed compared with individuals in the negative control in Experiment 1B.
In conclusion, this multi-site study has developed and established an optimized study design and a standardized test protocol to determine optimal growth and development of X. laevis larvae. High survival rates, homogenous rates of development, relatively short times to completion of metamorphosis of control animals, and concentration-dependent effects of E2 treatment on gonadal development were achieved using appropriate testing conditions. Based upon these results, it is concluded that the protocol provides a basis for evaluation on potential effects of chronic EAC exposure on gonadal development in X. laevis using a flow-through system. Furthermore, it is our opinion that the protocol could also be successfully applied to investigate the effects of other types of EAC, such as (anti)thyroidal chemicals, on the development of X. laevis.

development
sexual differentiation
xenopus laevis
17β-estradiol
flow-through exposure system
standardized protocol