|Home | About | Journals | Submit | Contact Us | Français|
A cohort of 1022 consecutive singleton births was generated during 1987–1988 in the Faroe Islands, where increased methylmercury exposure occurs from traditional seafood diets that include pilot whale meat. The prenatal exposure level was determined from mercury analyses of cord blood, cord tissue, and maternal hair. At age 14 years, 878 of 1010 living cohort members underwent detailed neurobehavioral examination. Eighteen participants with neurological disorders were excluded. Blood and hair samples obtained from the participants were analyzed for mercury. The neuropsychological test battery was designed based on the same criteria as applied at the examination at age 7 years. Multiple regression analysis was carried out and included adjustment for confounders. Indicators of prenatal methylmercury exposure were significantly associated with deficits in finger tapping speed, reaction time on a continued performance task, and cued naming. Postnatal methylmercury exposure had no discernible effect. These findings are similar to those obtained at age 7 years, and the relative contribution of mercury exposure to the predictive power of the multiple regression models was also similar. An analysis of the test score difference between results at 7 and 14 years suggested that mercury-associated deficits had not changed between the two examinations. In structural equation model analyses, the neuropsychological tests were separated into five groups; methylmercury exposure was significantly associated with deficits in motor, attention, and verbal tests. These findings are supported by independent assessment of neurophysiological outcomes. The effects on brain function associated with prenatal methylmercury exposure therefore appear to be multi-focal and permanent.
Methylmercury (MeHg) is a widespread contaminant of seafood and freshwater fish. The developing brain is considered the main target for MeHg toxicity, and the risk to consumers from MeHg exposure has been assessed by national and international agencies [11, 23, 28, 34]. Results from prospective epidemiological studies of birth cohorts from the Faroe Islands [18, 29] have contributed significantly to the data used for derivation of recommendations concerning human exposure limits, especially for pregnant women.
For the purposes of risk assessment, valid dose-effect relationships are essential, and the research in the Faroes has endeavored to represent both ends of the exposure-outcome relationship with the highest possible precision and sensitivity. Because of the need to obtain precise measures of the causative exposure [7, 14], the mercury concentration has been measured in cord blood and other biological samples from exposed participants. Detailed comparison of various prenatal exposure indices within the cohort has shown that the cord-blood concentration is consistently the most precise predictor of nervous system deficits determined during postnatal follow-up [7, 14].
At the same time, valid outcome variables must be sensitive to MeHg neurotoxicity and relatively robust to impacts of confounders. In choosing feasible effect parameters, an important consideration includes suitability for the age group and culture. As confounder-independent outcomes, neurophysiological measurements have provided evidence that developmental MeHg neurotoxicity is detectable through age 14 years . Clinical neurological tests have shown mercury-related adverse effects neonatally , but such tests may not be sensitive enough at school age . Standardized neuropsychological tests are useful to gain insight into functional domains and overall cognitive functioning and have been widely applied as sensitive indicators of central nervous system (CNS) dysfunction associated with exposure to neurotoxicants in both environmental and occupational exposure settings [9,39]. In choosing such measures for the study of children, it is important to consider the developmental stage of the child at both exposure and at the time of testing as well as cultural and psychometric parameters of the tests being administered [2,39].
Risk assessments have so far been based on functional outcomes up to early school age obtained in three prospective epidemiological studies [11, 23, 28, 34]. Data from older children and adolescents have been available only from less informative cross-sectional studies. Developmental exposure to inorganic lead is known to cause permanent CNS damage, thus suggesting that deficits attributable to early developmental neurotoxicant exposure may be irreversible [2, 40]. However, long-term effects of prenatal and early childhood exposures to neurotoxicants may become more difficult to demonstrate with time, because an increasing number of other factors play a role in test performance and it may be difficult, if not impossible, to define and measure these covariates. Thus, the increasing complexity of confounder adjustment may cause problems in the statistical analysis of cognitive test data. Despite these potential problems, prospective assessment of cohorts with well defined exposures occurring early in development is essential to understanding how neuropsychological outcome measures can be utilized as long-term indicators of early neurotoxicant insults.
We have prospectively followed a birth cohort from the Faroe Islands for 14 years. The Faroes are located in the North Atlantic between Norway, Shetland, and Iceland. In this Nordic fishing community, excess exposure to MeHg is mainly due to the traditional habit of eating meat from the pilot whale, while baseline exposures are due to frequent intake of other types of seafood that contains lower MeHg concentrations . In a study conducted when the cohort children were 7 years of age, the main finding was that decrements in specific functional domains were associated with prenatal MeHg exposure . Tests assessing the domains of attention, language, and verbal memory showed the most robust effects, while measures of motor speed and visuospatial function showed less consistent effects. Among several measures of prenatal and post-natal exposure, the strongest associations were obtained with the cord-blood mercury concentration as the exposure indicator [17–19]. These findings were robust in analyses controlled for age, sex and confounders, and they persisted after exclusion of high-exposure participants (i.e., maternal hair-mercury concentrations above 10 μg/g) . Ingestion of whale blubber causes exposure to lipophilic contaminants, notably polychlorinated biphenyls (PCBs), but the possible neurotoxic influence of this exposure did not explain the MeHg-associated neurobehavioral deficits [6, 16, 29]. The examination at age 14 years that is described in this paper included a clinical test battery similar to the one applied at age 7 years.
The cohort was assembled in the Faroe Islands during a 21-month period in 1986–1987 . The primary indicator of intrauterine exposure to MeHg was the mercury concentration in cord blood, and concentrations in maternal hair at parturition were also determined. Subsequently, mercury concentrations were measured in stored cord tissue (dry weight) from about half of the cohort members examined . MeHg exposure was found to vary considerably: 15% of the mothers had hair-mercury concentrations above 10 μg/g, while 4% were below 1 μg/g, a level that corresponds to the exposure limit recommended by the U.S. Environmental Protection Agency [28, 34]. Exposure to PCBs was assessed by analyzing cord tissue from 435 of the cohort subjects for major PCB congeners . The frequency of maternal intake of fish and of pilot whale meat was available from a questionnaire administered shortly after parturition . Cohort members were first invited for detailed examination at school age (7 years), when a total of 917 of eligible children (90.3%) participated .
At age 14 years, a total of 878 of 1,010 living cohort members (86.9%) were examined. The 438 boys and 440 girls examined had an average age of 13.83 (S.D., 0.32) years. Most examinations took place at the National Hospital in Tórshavn, the capital of the Faroe Islands, during April–June of 2000 and 2001. For families who had moved to Denmark, examinations were also offered in Odense in November, 2000. Transportation costs and other expenses were reimbursed in accordance with rates approved by the Faroese ethical review committee. Each day, four participants were examined during the morning and four during the afternoon. Each group of participants rotated between four stations staffed by health professionals, who administered the clinical tests, and who had no access to information on individual exposure levels. These stations each lasted approximately one hour each and covered physical examination, neurophysiological tests, and two series of neuropsychological tests, each of which was administered by the same examiner in the same sequence throughout.
The thorough pediatric examination included otoscopy and assessment of neurological function. None of the adolescents had current middle ear infection. A total of 18 participants examined had neurological disorders thought to be independent of MeHg exposure and were excluded from the data analysis. The MeHg exposure of these participants did not differ from that of other cohort members. Of 860 included participants that completed the 14-year examination schedule, 815 were examined also at age 7 years.
Blood specimens and hair samples were again obtained, and the proximal 2-cm hair segment was used for mercury analysis . Detailed quality assurance data showed acceptable analyses for both exposure biomarkers. Results in μg may be converted to nmol by multiplying by 5.0.
The study protocol was approved by the ethical review committee for the Faroe Islands and the U.S. institutional review boards, and written informed consent was obtained from all parents.
The test battery was designed based on the criteria previously described  as applied at the 7-year examinations . Because some of the tests used at age 7 were less feasible than desired and because other tests were thought to be appropriate for 14-year-old adolescents, some changes in the battery were made. Table 1 shows the tests used at the two examinations. Two Faroese-speaking board-certified psychologists administered the tests, each administering the same tests in the same sequence. Because test descriptions have been previously published, only the changes from the previous battery  will be outlined here.
Neuropsychological Examination System (NES2) Continuous Performance Test (CPT) : The total duration of the sustained attention test was increased to 10 minutes, and the inter-stimulus interval was reduced to 1000 ms as age-appropriate adjustments.
The Catsys® equipment (Danish Product Development, 3070 Snekkersten, Denmark) allows electronic assessment of simple reaction time using joint auditory and visual stimuli, as well as tapping speed, and speed of pronation-supination of the forearm .
Wechsler Intelligence Scale for Children – Revised (WISC-R) Block Designs: This visuospatial test has been standardized in the US up to age 17 years  and a translated Norwegian version up to 16 years . However, inspection of our data from the first year of testing showed a noticeable ceiling effect. To overcome this problem, we added the last three and most difficult test items from the Wechsler Adult Intelligence Scale – Revised (WAIS-R) version of the test . The usual discontinuation criteria and scoring rules for bonus points were applied.
Stanford-Binet Copying : Because this test had been successfully used in the field as an outcome parameter in MeHg-exposed children  and the scoring system is simpler, it replaced the Bender Visual Motor Gestalt Test, which was used as a measure of visual constructional abilities in the 7-year examinations. Administration of all 16 designs allowed equal numbers of test trials for all participants. Based on the manual , the last five (most difficult) items were identified as the most applicable to this group of participants given their age. Using these five items eliminated the noise from random errors on easier items. The outcome measures therefore included the total score for the full set items (16 designs) and the score for the last five items. As in our previous research involving other cohorts , an immediate memory condition was added as a measure of short-term memory for visuospatial material. This condition was administered immediately after the subject drew the 16 designs, with the instruction to draw all of the designs the subject could recall. The score for the number of correct designs drawn from memory was based on more lenient criteria, in which each drawing was scored as correct if the design was clearly recognizable, with no omissions of figure elements or gross distortions of proportions that would alter the perception of the gestalt .
Children’s Category Test : Level 2 of this test consists of 6 subtests and 83 items. The nonverbal stimulus material is designed to assess components of working memory, including concept formation, capacity to adopt strategies and flexibly invoke new strategies in problem solving, and the ability to benefit from experience. The subject receives immediate feedback as to the correctness of each response. The task does not require verbal or complex motor responses and is untimed. Score is total number of errors.
Wechsler Memory Scale - III (WMS-III) Spatial Span : For this task, the examiner taps with a finger a sequence on the blue blocks on a white form board. The subject has to repeat the demonstrated sequence. Series of increasing length are presented until the subject fails on both trials in a series of the same length. The series are administered first in forward and then in backward order. Score is the number of correct series repeated forward and backward. This test provides a visual analogue for WISC-R Digit Span. Both tests assess attention (forward spans) and working memory/executive function (backward spans).
Following descriptive analyses of the data, two approaches were taken to data analysis. First, multiple regression analysis was carried out as in the initial analyses of the 7-year data , and structural equation models were then developed, where test results were analyzed in larger, multi-test groupings .
The regression analysis was used to determine the associations between individual test outcomes, the cord-blood mercury concentrations and potential covariates. The covariates were selected to approximate the regression model used for analysis of the 7-year data set [8, 18]. Major obligatory confounders were therefore age, sex, maternal Raven score, and residence in town/village. Likewise, both maternal and paternal employment was again included, as was the time of the day of the examination (morning/afternoon), the language used (Faroese/Danish), and computer game experience (little/some/much; for computer-assisted tests only). Instead of the 7-year variable representing childcare before school entry, we used the participant’s grade in school.
An extended set of covariates was also examined, although some variables were available only for cohort members participating in the 7-year examinations (i.e., medical risk factors, maternal and paternal professional training, all of them included at age 7 years). These covariates also included family structure (whether the participant was living with both parents) and history of ever having tried tobacco smoking or alcoholic beverages. The participant’s school was rated as primary (up to 7th grade) or secondary (grades beyond 7th). The average Tanner score was obtained from the pediatric examination. All of these parameters showed only weak associations with MeHg exposures. A disadvantage of these supplementary analyses was that information was missing from several participants, thus affecting the statistical power of regression analyses that rely on complete cases. Additional analyses also incorporated the PCB exposure  and postnatal MeHg exposure parameters  as explanatory variables.
Contaminant concentrations entered the model after logarithmic transformation and the mercury regression coefficients were calculated to correspond to the change in the dependent variable associated with a doubling in MeHg exposure. Thus, these models assume a linear effect of log transformed mercury concentration on the neurobehavioral outcomes. We used generalized additive models to assess whether that assumption was appropriate, because these models do not require linearity assumptions while providing a smooth nonparametric dose-response curve . The regression coefficients were standardized by expressing them as percent of the standard deviation of the unadjusted outcome parameter. Some outcomes were transformed to obtain normally distributed residuals, but these transformations did not influence the significance of the mercury effect, and results of the regression analysis are therefore presented for the non-transformed data. However, transformed data were used for structural equation models (see below). The p values reported are two-sided, and for one-sided hypothesis testing these results should be divided by two.
The results from participants who were examined both at age 7 years and 14 years were further examined using summary statistics to determine whether the mercury effect had changed between the two examinations. In one approach, the difference between the two scores was used as the dependent variable and related to the prenatal MeHg exposure after confounder correction. As the test score difference may depend on the initial test performance level, a supplementary analysis was conducted with the log-transformed relative increase as the dependent variable. In these analyses, the set of covariates were supplemented so that it also included all potential confounders previously considered in the analysis of the 7 year data .
Because multiple regression analyses may not adequately address problems with exposure imprecision, multiple testing and missing data, structural equation models were applied in order to address these issues . These models incorporate two parts – measurement and structure. In the measurement part, the observed variables are considered to be manifestations of a smaller number of latent variables. The structural model then describes the relationship between the latent variables, possibly after adjustment for the effects of covariates. By including three independent measures of exposure (mercury in cord blood, mercury in maternal hair, and frequency of whale meat dinners during pregnancy), a latent exposure variable is generated that is likely to provide a closer approximation of true exposure than individual exposure measures considered alone .
Likewise, the neurobehavioral outcome variables were collected into a small number of latent variables determined from factor structure (Frodi Debes, unpublished results) and consistency with response modalities. According to this approach, test scores were assumed to reflect five broad functional groupings. These groupings are therefore empirically based and do not necessarily reflect functional domains typically used in neuropsychology. In order to include a sufficient number of tests in each group, outcomes were allowed to occur in more than one grouping. Motor outcomes included the three NES2 finger tapping scores, CATSYS supination/pronation score, CATSYS maximum finger tapping score. Attention outcomes were Digit Spans, Spatial Span, CATSYS mean reaction time, NES2-CPT mean reaction time, number of false responses, and number of missed responses. Spatial outcomes were the Children's category test, WISC-R Block Design, WISC+WAIS Block Designs, Copying total score, the Copying score for the last five designs, and number correctly recalled. Verbal outcomes comprised the two naming scores, the four California Verbal Learning Test (CVLT) scores, and the Similarities score. The memory outcomes were the four CVLT scores and the Copying recall score.
A separate model was developed for each of the neurobehavioral groupings considered as latent functions, and the mercury effect was expressed as percent of the standard deviation of the latent function result. The potential confounders previously considered in the primary regression model entered the structural equation model as covariates. Thus, we estimated the effect of a one-unit increase of the exposure after controlling for the potential confounders.
Because seafood also contains important nutrients, we conducted supplementary analyses, where maternal fish intake during pregnancy (number of fish dinners per week during pregnancy ) was included as a covariate that was allowed to affect both the MeHg exposure and the latent outcomes. Fig. 1 shows a path diagram that illustrates the initial model for verbal outcomes; similar model structures were developed for the four other sets of functional groupings.
The fit of these models was assessed in likelihood ratio testing against a more flexible model, which did not put any restrictions on the covariance matrix of the observed variables (known as the unrestricted model). If this test is rejected, then the proposed structural equation model does not fit the data as well as the unrestricted model. Models with a poor fit were modified until an acceptable fit was obtained by allowing for local dependence and differential item functioning [4,6]. Two outcomes are said to be locally dependent if the correlation between them cannot fully be explained by the underlying neurobehavioral function. Local dependence may arise because tasks resemble each other, and sub-scores on the same test may, therefore, show local dependence. In the present analysis, local dependence was anticipated for the two naming scores, the four verbal learning scores, the continuous performance scores, and the two copying scores, because tasks from the same test are interrelated, so that the scores are likely to be correlated even for participants with the same latent ability. Thus, by adjustment for local dependence, double-counting is avoided when highly correlated outcomes are included in the same model . Furthermore, an initial model did not allow for so-called differential item function and therefore assumes that the covariates affect the observed outcomes only through the latent outcome. This means that, except for scale differences, a one-unit increase in a given covariate is assumed to lead to the same increase in all outcomes representing the same latent variable . The appropriateness of this assumption was tested by also allowing covariates to affect the test scores directly.
The prenatal MeHg exposures of the participants of the 14-year examinations (Table 2) were similar to those of the cohort as a whole . Exposure levels at age 14 years averaged about one-fifth of those experienced prenatally, although exposures at age 7 years were slightly higher .
Table 3 shows associations with potential confounders. Continuous variables were trichotomized for the purpose of this table only. The expected pattern of associations reflects dietary habits and local differences in availability of whale meat during 1986–1987. Accordingly, participants from fishing villages had higher prenatal exposures than those living in towns or in Denmark. Participants tested in Denmark or in the Danish language had lower exposure. Lower exposure was also associated with a higher maternal Raven score. Adolescents from smaller primary schools had higher exposure, as had participants examined in the afternoon, at an older age, and those living with both parents.
The neuropsychological test results were generally within expected limits, based on Scandinavian and similar norms for this age range (Table 4). However, the Faroese adolescents performed above expectation given available norms for Block Designs, with a marked ceiling effect during the first year of examinations. Adding the three last designs from the adult version (WAIS-R) corrected for this imperfection in resolution at the high end, but was available for less than half of the participants.
All multiple regression analyses were based on complete cases (Table 5). Thus, the number of participants included depended upon the choice of exposure biomarker and confounders. Higher prenatal MeHg exposure was associated with lower finger tapping scores, increased reaction time, and lower cued naming scores; the cord-blood mercury concentration was a significant predictor of outcomes of all three types. The maternal hair mercury concentration showed significant or near-significant associations with deficits only on the three conditions of finger tapping, and the two measures of reaction time. The cord tissue mercury concentration showed no clear relation to finger tapping and reaction time results, but was associated with deficits on the naming and verbal learning results. One surprising find was that higher MeHg exposure was associated with better scores on WMS-III Spatial Span.
The estimated mercury effects remained virtually unchanged after adjustment for the extended set of confounders, except that, in this reduced data set, the p-value for the effect on the cued naming score was slightly above a 5% significance level. Some outcomes were transformed to obtain normally distributed residuals, but this change barely affected the mercury effects. Table 5 therefore shows the effect on the untransformed variables. The logarithmic transformation of the exposure variables assume a linear effect of each doubling of the exposure. The appropriateness of this assumption was confirmed in non-parametric models (Figs 2 and and33).
In contrast to the prenatal exposure variables, markers of postnatal exposures were generally only weakly related to cognitive test scores at 14 years, with many coefficients suggesting effects in the direction opposite to expectation. Only one significant association was seen, i.e., for the hair mercury concentration at age 14 years on the NES2 finger tapping score on the preferred hand (effect, −5.18; p = 0.019). This effect appeared less clear when corrected for the full set of confounders (effect, −3.18; p = 0.17) and decreased further when adjusted for the prenatal exposure level (effect, −2.74; p = 0.26). On the other hand, prenatal effects changed only slightly after adjustments for postnatal exposure markers. However, due to limited power in a regression model with both prenatal and postnatal exposure variables, the p-values generally increased slightly. Thus, the p values for the three outcomes significantly associated with the cord blood mercury concentration without adjustment for postnatal exposure increased to levels ranging from 0.061 (CPT reaction time) to 0.11 (finger tapping), after adjustment for postnatal exposure. On the other hand, the Copying test (5 most difficult designs) now showed a significant p-value of 0.03 for the effect of prenatal exposure effects after adjustment for postnatal exposure.
Likewise, PCB exposure, although only available for half of the subjects, showed only weak associations with the outcomes, none of them reaching statistical significance. The p values for the three outcomes showing significant deficits associated with the cord blood mercury concentration (Table 5) varied from 0.15 (CPT reaction time) to 0.98 (finger tapping with both hands).
In general, a high percentage of the outcome variation could not be explained in the regression models and most covariate associations were weaker than at age 7 years . Age was not a general predictor of test performance at this developmental stage and within the limited age range covered, and was only weakly related to scores on verbal abstraction, attention, and manual motor speed tasks. Sex was a stronger predictor of test performance, where boys performed faster on tasks assessing manual motor speed, while girls earned higher scores on most other tasks. No sex difference was seen in vocabulary, nonverbal abstraction, and visuospatial memory. The maternal Raven score was again an important predictor for most outcomes. Residence in towns was positively associated with naming ability and, to some degree, with other abilities. Computer experience was associated with better performance on the continuous performance tasks and finger tapping. However, the full complement of predictors explained only about 10% of the total variation (Table 6). Although prenatal MeHg exposure seemed a less important predictor of neuropsychological performance at age 14 than at age 7 years, this risk factor appeared to represent about the same proportion of the total variance explained by the regression model at the two occasions.
Assuming that prenatal MeHg neurotoxicity causes permanent changes, the mercury-associated effects observed at age 14 years could represent the deficits previously observed at age 7 years, as modified by learning and maturation. An analysis was therefore conducted based on outcomes available from both examinations and with a significant mercury effect at age 7 or 14 years (Table 7). For most outcomes, prenatal MeHg exposure showed only weak effects on the test score difference. The sign of the coefficients generally suggested that highly exposed participants may have been able to compensate in part for deficits seen at age 7 years. However, this effect was weak and only statistically significant for one of the CVLT scores. For two finger tapping scores (non-preferred hand and both hands), this analysis indicated that, between ages 7 and 14 years, highly exposed participants had fallen even further behind those less exposed, but again this tendency was not statistically significant. A supplementary analysis of the relative difference in test scores yielded similar results, but these models provided a poorer fit to the data.
To take into account multiple statistical testing, exposure imprecision and missing data, we then applied structural equation models. The crude, initial models fitted the data poorly, but after adjustment for local dependence and differential item functioning, models with an acceptable fit were obtained for all five test groups (Table 8). The strongest mercury effect was found in regard to the motor and attention groupings. The effect for the verbal grouping is close to being significant at the 5% level, while the effect on spatial and memory groups is far from significant.
Maternal fish intake during pregnancy was then added to the structural equation models (Table 8). Although fish intake appeared to exert a weak beneficial effect, it was not statistically significant. Inclusion of this adjustment changed the mercury coefficients in the direction anticipated.
Significant mercury effects corresponded to a decrease corresponding to 6.5–8.5% of the standard deviation for each doubling in the exposure (Table 8). These results are in accordance with the strongest regression coefficients in Table 5. For motor function, a doubling of the true MeHg exposure corresponded to 1.3 fewer taps on NES2 finger tapping with the preferred hand. As a result, at age 14 years, this effect is slightly stronger than at 7 years. Other effects were about the same or less. In the regression analysis, a strong positive mercury effect had been obtained for Spatial Span, and this tendency led to poorer fit of the structural equation for attention. A better fitting model with a stronger effect would be obtained by excluding this score, but this possibility was not further pursued.
As expected, local dependence were observed for the two naming scores, the four CVLT scores, the NES2 continuous performance scores, and the two Copying scores. Likewise, differential item functioning occurred for some outcomes, mostly because test scores depended on sex. Allowance for these factors improved the fit of the models, but affected the mercury coefficients only minimally. For the spatial outcomes the final model did not fit data as well as the unrestricted model and the results should be cautiously interpreted.
This study presents results on neuropsychological performance of adolescents with widely differing degrees of prenatal exposure to MeHg from maternal seafood diets during pregnancy and lower postnatal exposures to this neurotoxicant. Adverse effects were identified in regard to motor speed, attention, and language. These effects were apparent both in multiple regression analyses and in structural equation models that take into account multiple testing, exposure imprecision, and incomplete data. Potential confounders were included in the analyses. These findings are in accordance with the results obtained from examination of the same cohort members at age 7 years and suggest that prenatal MeHg exposure is associated with enduring neurotoxic effects on CNS function . Of particular interest, the Boston Naming Test results continued to show effects of prenatal MeHg exposure. At age 7, prenatal MeHg exposure was associated with the capacities to name pictured objects both under cued and uncued conditions, thus suggesting that lexical development was more limited at high exposures. The findings for age 14 further support an exposure-related limitation of lexical development, because participants with higher exposure had a poorer ability to recall names of objects, even when cues were provided, than did participants with lower exposures.
While the effect at age 14 years appeared to be somewhat weaker than seven years before, the overall explanatory power of the statistical models was also more limited than it was previously. A multitude of events and incidents, beneficial and adverse, may have occurred during the time between the two examinations and may have attenuated the correlation between two sets of tests and decreased the association with early-life covariates. For example, differences in time of onset of puberty, together with very high growth rates in this phase of development, may contribute unadjusted variability that overshadows the association with defined predictors, including MeHg exposure. In addition, the covariates and methylmercury exposure might particularly affect the rate at which skills are acquired at the younger age, while innate ability might be a more important predictor of performance at 14 years. Thus, although the neuropsychological performance results at age 14 were more affected by variability that could not be captured by adjustments using known covariates, the relative importance of MeHg exposure was about the same at the two examinations. In addition, analyses that incorporated test results at both examinations showed that the mercury-associated deficits had changed only little between age 7 and 14 years.
When considering the estimated effect size, comparison of regression coefficients with the standard deviations of the test scores is appropriate . Because of the logarithmic transformation of the exposure scale, effects must be expressed in terms of log-scale increases. We found that a doubling of the prenatal MeHg exposure at age 7 years corresponded to a decrease in performance in the range of about 5–10% of the standard deviation [18,19]. The effect at age 14 years was similar in size (Table 8). The regression coefficients obtained for mercury may be compared to those obtained for age. At age 7 years, a doubling in MeHg exposure corresponded to a loss in development by about 1–2 months . At age 14, age relationships for the neuropsychological tests are less steep, and a doubling in the MeHg exposure therefore corresponded to greater delays in development.
While most tests showed a mercury association in the direction predicted, Spatial Span showed an unanticipated positive effect of increased prenatal MeHg exposures. This tendency caused a poor fit of one of the structural equation models. Although a chance finding cannot be excluded, one must also consider the possibility that MeHg exposure may be associated with intake of essential nutrients from seafood, and that, for example, long-chain n-3 fatty acids may have beneficial effects on brain development [1,22]. Partial adjustment for this factor was provided by adding maternal fish intake during pregnancy as a covariate. This addition changed the mercury regression coefficients in the direction anticipated to better reflect the true adverse effect of methylmercury per se. This confounding had previously been considered as an explanation of the weak tendency of decreased (improved) visual evoked potential latencies at higher prenatal MeHg exposure levels . It therefore seems likely that the outcomes of the present study are affected by both contaminant toxicity and beneficial nutrient effects in the opposite direction, although the different seafood components may not affect different brain functions to the same degree. Furthermore, it is apparent that any beneficial nutrients would have compensated only partially for the neurotoxic effects of MeHg in this cohort. However, adjustment for fish intake was based only on maternal recall of the frequency of fish dinners during pregnancy, thus not allowing any precise assessment of essential nutrient intakes. If more precise indicators of nutrient supply had been available, statistical control for their beneficial effects might have revealed an even stronger MeHg effect.
The results of the multiple regression analyses must be interpreted in light of various uncertainties, in particular exposure misclassification [7, 14]. Thus, while our previous studies have shown that the cord-blood mercury concentration is the exposure parameter that best predicts test performance, closely followed by the mercury concentration in cord tissue , both parameters are participant to substantial variability that exceeds the normal laboratory error. The mercury concentration in maternal hair is substantially more imprecise and may be severely affected by, e.g., permanent hair treatment . Such exposure imprecision is likely to be random and would therefore bias the regression coefficients toward null. However, the imprecision of the exposure markers was taken into account in the structural equation models by viewing the observed parameters as reflections of a true underlying exposure variable. These more advanced statistical analyses revealed that the grouped test score outcomes, labeled as motor, attention, and verbal functions, had the clearest dose-effect relationships for developmental MeHg exposure.
In both the multiple regression and structural equation analyses, tests or groups of tests that showed the best dose-effect relationships with regard to prenatal MeHg exposure tended to be highly sensitive tests with wide score ranges or extremely accurate outcome measurements. For example, reaction time and manual motor speed as measured by computer are highly accurate. The Boston Naming Test has 60 items, allowing much better discrimination of the subject’s abilities than tests with fewer items. Therefore, the mercury-associations of the tests with superior psychometric properties strengthen the argument that MeHg has a long-term adverse effect on brain development and function. The tests assessing visuospatial skills appeared to have a low ceiling in this population at 14 years, while responses on tests of short-term memory generally had a narrow range, thereby limiting their ability to reveal exposure-effect relationships. Even given these provisos, the confluent finding of effects on tasks assessing attention, manual motor speed, and language/semantic memory suggests that the effects of MeHg at these levels of brain development are multi-focal or diffuse, affecting widely divergent brain structures. The existence of sustained neurotoxic effects in adolescents is also supported by the neurophysiological examinations that were carried out in parallel . Thus, auditory evoked potential latencies at age 14 years increased at higher prenatal MeHg exposures, and the heart rate variability decreased. The notion of permanent effects is also consistent with the known neuropathological effects of prenatal MeHg exposure at much higher levels producing clinical Minamata disease [15,28], but clearly as a much more subtle picture.
Postnatal exposures were generally weakly related to cognitive test scores at 14 years, although postnatal exposures were poorly characterized and averaged only about one-fifth of the prenatal exposures. Significant effects were seen on the finger tapping scores, but these effects could be due to chance and were less certain, when corrected for the full set of confounders. Because the neuropsychological outcomes may be affected by both prenatal and postnatal exposures, and the exposures at different times were interrelated , these effects would be difficult to separate. In addition, using the postnatal exposure to adjust the regression coefficients for prenatal exposure is problematic, as this approach will involve two error-prone and associated predictors, thereby causing complex biases they may be difficult to predict . Given these caveats, postnatal exposure in this population did not appear to affect the neurobehavioural outcomes to any appreciable degree. However, our results on peak V of the evoked potentials showed a possibly specific effect of postnatal (recent) MeHg exposures , thus supporting the finding that recent exposure could conceivably cause a decrease in motor speed. This effect is perhaps surprising, because the average MeHg exposure of the 14-year-olds was reflected by a hair-mercury concentration of about 1 μg/g, a level that corresponds to the Reference Dose used by the U.S. Environmental Protection Agency . However, these findings are in accordance with several cross-sectional studies carried out in 7-to-12-year-old children, who also showed decreased motor speed associated with increased hair-mercury concentrations . While the hair-mercury concentration in the latter studies reflected the child’s own recent exposure, an association with prenatal exposures could not be ascertained.
PCB is another neurotoxicant present in seafood, and increased exposures to methylmercury and PCBs are correlated . At age 7 years, prenatal PCB exposure was associated with deficits on CPT reaction time, naming, and long-term verbal memory, although these tendencies all but disappeared after adjustment for mercury exposure [16,18]. This result was confirmed in a structural equation analysis that allowed for measurement errors in biomarkers of both MeHg and PCB . No association with PCB exposure was seen at age 14 years.
The findings of the present study add to the substantial evidence on developmental methylmercury neurotoxicity [11,15,23,28]. One previous study carried out in the Seychelles  found no significant association between prenatal methylmercury and performance on neurobehavioral tests at age 8 years. However, the two studies are not contradictory, because the Seychelles study cannot exclude an effect at the magnitude observed at the Faroes . In addition, the mercury effect was assessed only from analyses of maternal hair collected six months after parturition and after adjustment for the child’s postnatal mercury exposure . Joint analyses based on three major prospective studies are in excellent agreement with the findings in the Faroes [23,28].
The long-term effects of developmental neurotoxicity are of substantial medical and societal importance. Significant exposure-associated deficits were observed at adolescence, although subtle neurotoxicity, e.g., due to lead exposure , may become more difficult to document in adolescents, who may have had a longer and partially unknown exposure record and who may also differ in regard to the timing and speed of puberty development and a variety of other cofactors. The present study differs from most previous studies by emphasizing domain-related brain functions, rather than formal IQ testing. However, some WISC-R subtests were included in this study, and MeHg neurotoxicity can be expressed in terms of approximate IQs . Such IQ deficits are known to affect the likelihood that an adolescent completes high school and the chance of completing an academic or vocational education. The associated costs to society are therefore substantial .
This study was supported by grants from the US National Institute of Environmental Health Sciences (ES09797) and the Danish Medical Research Council. The contents of this paper are solely the responsibility of the authors and do not represent the official views of the NIEHS, NIH or any other funding agency. Advice on neurobehavioral test selection was contributed by the following members of the steering committee convened by NIEHS: David Bellinger, Kim Dietrich, Annette Kirshner, and David Otto. Arne Ludvig administered one of the test panels to all participants; Nicolina Sørensen and Flemming Juul Hansen conducted the pediatric examinations; Anne Sofie Veyhe administered the CATSYS tests; and Mari Ann Ellendersen was in charge of cohort maintenance and clinical scheduling.