Geochemistry of ground water and the incidence of acute myocardial infarction in Finland

A Kousa1, E Moltchanova2, M Viik-Kajander2, M Rytkönen2, J Tuomilehto2, T Tarvainen3, M Karvonen2 for the Spat Study Group

1 Geological Survey of Finland, Kuopio, Finland
2 Department of Epidemiology and Health Promotion, National Public Health Institute, Helsinki, Finland
3 Geological Survey of Finland, Espoo, Finland

Correspondence to:
Correspondence to:
Anne Kousa
Geological Survey of Finland, POBox 1237, FIN-70211 Kuopio, Finland; anne.kousa@gsf.fi

Accepted for publication 13 August 2003



Study objective: To examine the association of spatial variation in acute myocardial infarction (AMI) incidence and its putative environmental determinants in ground water such as total water hardness, the concentration of calcium, magnesium, fluoride, iron, copper, zinc, nitrate, and aluminium.

Design: Small area study using Bayesian modelling and the geo-referenced data aggregated into 10 kmx10 km cells.

Setting: The population data were obtained from Statistics Finland, AMI case data from the National Death Register and the Hospital Discharge Register, and the geochemical data from hydrogeochemical database of Geological Survey of Finland.

Participants: A total of 18 946 men aged 35–74 years with the first AMI attack in the years 1983, 1988, and 1993.

Main results: One unit (in German degree °dH) increment in water hardness decreased the risk of AMI by 1%. Geochemical elements in ground water included in this study did not show a statistically significant effect on the incidence and spatial variation of AMI, even though suggestive findings were detected for fluoride (protective), iron and copper (increasing).

Conclusions: The results of this study with more specific Bayesian statistical analysis confirm findings from earlier observations of the inverse relation between water hardness and coronary heart disease. The role of environmental geochemistry in the geographical variation of the AMI incidence should be studied further in more detail incorporating the individual intake of both food borne and water borne nutrients. Geochemical-spatial analysis provides a basis for the selection of areas suitable for such research.

Keywords: Bayesian approach; acute myocardial infarction; geochemistry; ground water; small area analysis

Abbreviations: CVD, cardiovascular disease; AMI, acute myocardial infarction; CHD, coronary heart disease

Cardiovascular disease (CVD) is the major cause of death in most developed countries including Finland.1–3 The occurrence of coronary heart disease (CHD) varies between populations14 but also within populations inside a country.5 Already in 1947 Kannisto found that CHD mortality was much higher in the eastern part than in the western part of Finland. In the 1980s the CHD risk was still 40% higher in eastern Finland than that in western and southern parts of the country.6 The major CHD risk factors do not fully explain the geographical variation of CHD risk in Finland.7–9 Although the geographical differences have long been known, the reasons are still partly ambiguous.10 Besides a genetic predisposition11–13 several lifestyle and environmental factors have been implicated in the pathogenesis of CHD.14–16 Availability of trace elements in soil and ground water may be a cause of certain chronic ailments.17 Soils and rocks in the countries of northern Europe are poor sources of many essential trace elements.1718

Our recent study of the spatial distribution of the first acute myocardial infarction (AMI) event showed that despite the decreasing trend in AMI incidence, the geographical difference in incidence and high risk areas has remained within Finland.19 The presence of high risk areas for AMI suggests that genetic or environmental risk factors have accumulated in certain geographical locations in Finland. Our aim was to examine the possible association of spatial variation of AMI incidence with geochemical compounds in ground water.


*    METHODS
TOP
 ABSTRACT
 METHODS
RESULTS
 DISCUSSION
 REFERENCES
 
Finnish ground water is slightly acidic and very soft (1–4°dH) or soft (4–8°dH).20 Besides the geological factors affecting trace element composition, atmospheric, anthropogenic, and marine factors also contribute to the chemical composition of the ground water.21

The data on men aged 35–74 years with the first attack of AMI (18946 cases) were obtained from the nationwide Death Register and the Hospital Discharge Register. The national personal identification number was used to perform a computerised records linkage of the data for deaths and hospitalisation attributable to AMI (ICD-8 and ICD-9 codes 410–414). Both fatal and non-fatal events from the years 1983, 1988, and 1993 were included in the study. Cases with a previous hospitalisation for AMI were excluded. Data for these three years have been pooled. The data on population at risk, provided by coordinates of the place of residence, were obtained from Statistics Finland. The data were aggregated into 10 kmx10 km grid cells to ensure the protection of privacy of the individuals.

Geochemical data were obtained from the hydrogeochemical database of the Geological Survey of Finland.21 The data on total water hardness (°dH), Ca, Mg, Fe, F-, NO3-(mg/l) and Cu, Zn, and Al (µg/l) were available. Element concentrations were determined with different methods, for example, ICP-MS, ICP-AES, iconography, and AAS. The original data contained from 3621 up to 12 407 ground water samples.

The geochemical data were interpolated into a regular grid by using the ALKEMIA software developed at Geological Survey of Finland.22 In the ALKEMIA Smooth interpolation method, the nearest samples to the grid cell receive greater weight. The value of the cell is a weighted median of sample values.23–26

Bayesian spatial conditional autoregressive model (CAR) with covariates, which is currently in wide use in the field of the disease mapping, was applied in this study.27–30 Because Finland is sparsely inhabited, we propose one modification, which is pertinent to the sparsely populated areas. In the case of the 10 kmx10 km grid over Finland (excluding Lapland), some grid cells are empty and have to be omitted from the analysis; thus 5% of cells would be omitted. However, once we take environmental factors into account, assuming that the disease risk is influenced by both demographic factors (that is, people who actually live within the grid cell) and environmental factors in each cell whether or not it is inhabited, the omission of unpopulated cells results in a loss of information. The covariates included in the model were the age of onset of AMI and the levels of geochemical compounds in the ground water. The following modification is thus proposed.

Let Yik denote the number of cases in the cell i and age group k. Furthermore, let Nik denote the respective population at risk. The proposed probability distribution is then as follows:



that is, the Poisson distribution is assigned to the inhabited cells and the uninhabited cells naturally have no cases of the disease with the unit probability. Also we assign common regression structure to the µik:



where

{alpha}, is the baseline risk

{lambda}i, is the local unexplained spatial random effect

ßk, is the effect of age group k on the risk level

K, is the age group, k = 0,...,K

{xi}, is a vector of environmental covariate effects

Zi, is a vector of environmental covariates for area i

In this analysis, the age axis was divided into eight, five year age groups: 35–39, 40–44, 45–49, 50–54, 55–59, 60–64, 65–69, and 70–74. A non-proportional hazard model described their effect, which for AMI is more appropriate than the proportional hazards.

As outlined in the preceding section seven geochemical covariates were included in the analysis.

The regression coefficients ß and {xi} were given non-informative Normal priors N(0,0.00001), the background level {alpha} was given an improper flat prior



and the {lambda} were given a CAR structure:



where

{lambda}-i, are spatial variation parameters in the neighbourhood of i

mi, is the number of neighbors for cell i

{tau}, is the overall level of spatial precision (inverse spatial variation)

In the CAR models a neighbourhood structure needs to be defined. The neighbours were defined to be all those cells adjacent to the cell i through side or corner. Thus each cell could have at most eight (8) neighbours.

The model was fitted using WinBUGS. A total of 10000 iterations with 5000 burn-in were run. "Burn-in" denotes iterations, which were discarded because of non-convergence of the model at the early stages of the algorithm. The evaluation of the test results showed that a satisfactory convergence was reached.

The posterior joint and marginal distributions of the parameters of interest were estimated and summarised. The p% highest density regions (HDR), defined as most compact set of parameter values the posterior density mass over which is p/100, is used in Bayesian statistics to describe the variability of the estimate. It is thus by its nature somewhat similar to the frequentist confidence interval.


*    RESULTS
TOP
 ABSTRACT
 METHODS
 RESULTS
DISCUSSION
 REFERENCES
 
Age group and the total water hardness, Ca, Mg, Fe, F-, Cu, Al, Zn, and NO3- concentrations in the ground water were included in the analyses as covariates. The overall age adjusted incidence of AMI among men aged 35–74 year was 480/100 000/year (posterior 95% HDR 473, 487). Table 1Go gives information on the chemical contents of ground water. Table 2Go illustrates the number of AMI cases, population at risk, and AMI incidence by age and water hardness. One unit (°dH) increment in water hardness decreased the risk of AMI by 1% (table 3Go). The levels of other geochemical elements included in this study did not have any additional effect on the spatial variation of the incidence of AMI.


View this table:
[in this window]
[in a new window]
  		Table 1 Geochemical concentrations in ground water in Finland
 
View this table:
[in this window]
[in a new window]
  		Table 2 Number of AMI cases, population at risk, and the AMI incidence per year by age and water hardness among 35 to 74 year old men in Finland in 1983, 1988, and 1993 (pooled data)
 
View this table:
[in this window]
[in a new window]
  		Table 3 Estimated effects of the geochemical covariates on the incidence of the first AMI among 35–74 year old Finnish men in 1983, 1988, and 1993 (pooled data)
 
*    DISCUSSION
TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
REFERENCES
 
The large geographical variation and changes in the incidence of AMI in Finland cannot be explained by individual lifestyle or genetic factors alone; environmental exposures must also contribute to the development of the disease. The classic risk factors and socioeconomic status provide only a partial explanation for the excess CHD risk in eastern Finland.731 The age distribution of the population did not have an effect on the geographical variation of the incidence of AMI. The results support the early observations of the inverse relation between the AMI incidence and total water hardness. An inverse relation between water hardness and CVD mortality has been detected in several studies.32–37 They have suggested that CHD mortality can be related to the amount of magnesium and calcium in drinking water.3638–44 In some studies an association between CVD and water hardness was not found.45–48 Much of the disagreement in earlier studies may be related to the complexity of the ecological analysis and the difficulty to apply results from ecological studies at the individual level.

In the general population, the magnesium intake has decreased over the years especially in the western world.49 Some previous studies have shown that a large number of subjects had a lower intake of magnesium than the recommended dietary amount (350 mg/day).42 It has been suggested that magnesium in water, appearing as hydrated ions, has a higher bioavailability than magnesium in food, which is bound in different compounds that are less easily absorbed.50

Fluoride concentrations of around one mg/l in household water may be beneficial.4041 Recent studies have also provided evidence that high serum iron and copper concentrations are associated with the CHD.5152 In this study one mg/l increment in the fluoride concentration in the drinking water was associated with a 3% decrease in the risk of AMI. In our study one µg/l increment in copper and one mg/l increment in iron on average increased the risk of AMI by 4% and 10%, respectively. The differences were not, however, statistically significant. The non-significant results in our study may be attributable to excessive smoothing technique. Thus, our study provides further supportive evidence for the importance of the ground water fluoride, iron and, copper concentrations for the risk of AMI.

CHD has a multifactorial aetiology. The method of spatial analysis used in this study is especially useful for testing the impact of several factors simultaneously. The validity of the Bayesian method used in this study has been also demonstrated earlier studies.192753 Additional simulations have been run to check the validity of the proposed changes to it regarding the inclusion of the uninhabited cells in the analysis.

Ground water reflects the contents of trace elements in soil and bedrock2154 but only a small proportion of the population use locally produced food supplies, cereals, and vegetables. Individual studies on the role of intake of both food and water-borne nutrients should incorporate environmental exposure or control for it.


*    FOOTNOTES
 
Funding: this work was partly supported by Academy of Finland (no 78422), by the Yrjö Jansson Foundation and Juho Vainio Foundation.

Conflicts of interest: none declared.


*    REFERENCES
TOP
 ABSTRACT
 METHODS
 RESULTS
 DISCUSSION
 REFERENCES
 
  1. Thom TJ, Epstein FH, Feldman JJ, et al. Total mortality and mortality from heart disease, cancer, and stroke from 1950 to 1987 in 27 countries: highlights of trends and their interrelationships among causes of death. Washington, DC: US DHHS PHS, National Institutes of Health, NIH Publication, 1992:92–3088.
  2. Murray CJL, Lopez AD. Evidence-based health policy-lessons from the global burden of disease study. Science 1996;274:740–3.[Free Full Text]
  3. Uemura K, Pisa Z. Trends in cardiovascular disease mortality in industrialized countries since 1950. World Health Stat Q 1988;41:155–78.[Medline]
  4. Tunstall-Pedoe H, Kuulasmaa K, Amouyel P, et al. Myocardial infarction and coronary deaths in the World Health Organization MONICA Project. Registration procedures, event rates and case fatality in 38 populations from 21 countries in 4 continents. Circulation 1994;90:583–612.[Medline]
  5. Kannisto V. The causes of death as demographical factors in Finland. [In Finnish, English summary]. Helsinki: Kansantaloudellisia tutkimuksia-Economic studies XV, 1947.
  6. Tuomilehto J, Arstila M, Kaarsalo E, et al. Acute myocardial infarction (AMI) in Finland—baseline data from the FINMONICA AMI register in 1983–1985. Eur Heart J 1992;13:1153.[Abstract/Free Full Text]
  7. Jousilahti P, Vartiainen E, Tuomilehto J, et al. Role of known risk factors in explaining the difference in the risk of coronary heart disease between eastern and southwestern Finland. Ann Med 1998;50:481–7.
  8. Vartiainen E, Jousilahti P, Alfthan G, et al. Cardiovascular risk factor changes in Finland, 1972–1997. Int J Epidemiol 2000;29:49–56.[Abstract/Free Full Text]
  9. Forsen T, Eriksson J, Qiao Q, et al. Short stature and coronary heart disease: a 35-follow-up of the Finnish cohorts of the Seven Countries Study. J Intern Med 2000;248:326–32.[CrossRef][Medline]
  10. Koskinen S. Origins of regional differences in mortality from ischaemic heart disease in Finland. Jyväskylä: NAWH, Research reports 41, 1994:204.
  11. Kontula K, Koivisto U- M, Miettinen H, et al. Suomalaiset sepelvaltimotautigeenit. [The Finnish IHD-genes, in Finnish]. Duodecim 1994;110:721–9.[Medline]
  12. Koskenvuo M, Kaprio J, Romanov K. Twin studies in metabolic diseases. Ann Med 1992;24:379–81.[Medline]
  13. MacCluer J, Kammerer C. Dissecting the genetic contribution to coronary heart disease. Am J Hum Genet 1991;49:1139–44.[Medline]
  14. Tuomilehto J, Puska P, Korhonen H, et al. Trends and determinants of ischaemic heart disease mortality in Finland: with special reference to a possible levelling off in the early 1980s. Int J Epidemiol 1989;18 (3 suppl 1) :109–17.
  15. Vartiainen E, Puska P, Pekkanen J, et al. Changes in risk factors explain changes in mortality from ischaemic heart disease in Finland. BMJ 1994;309:23–7.[Abstract/Free Full Text]
  16. Jousilahti P, Tuomilehto J, Vartiainen E, et al. Importance of risk factor clustering in coronary heart disease mortality and incidence in eastern Finland. J Cardiovasc Risk 1995;2:63–70.[CrossRef][Medline]
  17. Masironi R. Geochemistry, soils and cardiovascular diseases. Experientia 1987;43:68–74.[CrossRef][Medline]
  18. Koljonen T. Bedrock and associated geophysical and other features. In: Koljonen T, ed. The geochemical atlas of Finland. Vol part 2. Till: Geological Survey of Finland, 1992.
  19. Karvonen M, Moltchanova E, Viik-Kajander M, et al. Regional inequality in the risk of acute myocardial infarction in Finland: a case study of 35- to 74-year-old men. Heart Drug 2002;2:51–60.
  20. Tarvainen T, Lahermo P, Hatakka T, et al. Chemical composition of well water in Finland—main results of the "One thousand wells" project. In: Autio S, ed. Special paper 31. Espoo: Geological Survey of Finland, 2001:57–76.
  21. Lahermo P, Ilmasti M, Juntunen R, et al. The geochemical atlas of Finland, Part 1. The hydrogeochemical mapping of Finnish ground water. Espoo: Geological Survey of Finland 1990:66.
  22. Ahlsved C, Lampio E, Tarvainen T. ALKEMIA—a VAX minicomputer database and program package for geochemical exploration. J Geochem Explor 1991;41:23–8.[CrossRef]
  23. Björklund A, Lummaa M. Representation of regional, local and residual variability of geochemical data by means of filtering techniques. Proceedings 2nd international symposium on methods of prospecting geochemistry, Irkutsk, USSR, 1983:25–34.
  24. Björklund A, Gustavsson N. Visualization of geochemical data on maps: new options. J Geochem Explor 1987;29:89–103.[CrossRef]
  25. Tarvainen T. Environmental applications of geochemical databases in Finland. [PhD Thesis]. Espoo: Geological Survey of Finland, 1996.
  26. Gustavsson N, Lampio E, Tarvainen T. Visualization of geochemical data on maps at the Geological Survey of Finland. J Geochem Explor 1997;59:197–207.[CrossRef]
  27. Ranta J, Penttinen A. Probabilistic small area risk assessment using GIS-based data: a case study on Finnish childhood diabetes. Stat Med 2000;19:2345–59.[CrossRef][Medline]
  28. Osnes K, Aalen OO. Spatial smoothing of cancer survival: a Bayesian approach. Stat Med 1999;18:2087–99.[CrossRef][Medline]
  29. Clayton DG, Bernardinelli L, Montomoli C. Spatial correlation in ecological analysis. Int J Epidemiol 1993;22:1193–202.[Abstract/Free Full Text]
  30. Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Annals of the Institute of Statistical Mathematics 1991;43:1–59.
  31. Salomaa V, Niemelä M, Miettinen H, et al. Relationship of socioeconomic status to the incidence and prehospital, 28-day, and 1-year mortality rates of acute coronary events in the FINMONICA Myocardial Infarction Register Study. Circulation 2000;101:1913–18.[Medline]
  32. Crawford MD, Gardner MJ, Morris JN. Mortality and hardness of local water-supplies. Lancet 1968;i:827–31.
  33. Masironi R, Pisa Z, Clayton D. Myocardial infarction and water hardness in European towns. Environ Pathol Toxicol 1980;4:77–87.
  34. Karppanen H, Pennanen R, Passinen L. Minerals, coronary heart disease and sudden coronary death. Adv Cardiol 1978;25:9–24.[Medline]
  35. Nerbrand CH, Svärdsudd K, Ekland J, et al. Cardiovascular mortality and morbidity in seven counties in Sweden in relation to water hardness and geological settings. Eur Heart J 1992;13:721–7.[Abstract/Free Full Text]
  36. Rylander R, Bonevik H, Rubenowitz E. Magnesium and calcium in drinking water and cardiovascular mortality. Scand J Work Environ Health 1991;17:91–4.[Medline]
  37. Piispanen R. Water hardness and cardiovascular mortality in Finland. Environmental Geochemistry and Health 1993;15:201–8.[CrossRef]
  38. Punsar S, Karvonen MJ. Drinking water quality and sudden death: observations from West and East Finland. Cardiology 1979;64:24–34.[Medline]
  39. Anderson TW, Neri LC, Schreiber GB, et al. Ischemic heart disease, water hardness and myocardial magnesium. [Letter]. Can Med Assoc J 1975;113:199–203.[Abstract]
  40. Luoma H, Helminen SK, Ranta H, et al. Relationships between the fluoride and magnesium concentrations in drinking water and some components in serum related to cardiovascular diseases in men from four rural districts in Finland. Scand J Clin Lab Invest 1973;32:217–24.[Medline]
  41. Luoma H, Aromaa A, Helminen S, et al. Risk of myocardial infarction in Finnish men in relation to fluoride, magnesium and calcium concentration in drinking water. Acta Med Scand 1983;213:171–6.[Medline]
  42. Rubenowitz E, Molin I, Axelsson G, et al. Magnesium in drinking water in relation to morbidity and mortality from acute myocardial infarction. Epidemiology 2000;11:416–21.[CrossRef][Medline]
  43. Rubenowitz E, Axelsson G, Rylander R. Magnesium and calcium in drinking water and death from acute myocardial infarction in women. Epidemiology 1999;10:31–6.[Medline]
  44. Punsar S, Erämetsä O, Karvonen MJ, et al. Coronary heart disease and drinking water. J Chron Dis 1975;28:259–87.[CrossRef][Medline]
  45. Huel G, Thomazeau R, Derriennic F, et al. Water hardness and cardiovascular mortality. A study of 947 Alsatian communities. Rev Epidemiol Sante Publique 1978;26:381–90.[Medline]
  46. Smith WC, Crombie IK. Coronary heart disease and water hardness in Scotland—is there a relationship? J Epidemiol Community Health 1987;41:227–8.[Abstract]
  47. Maheswaran R, Morris S, Falconer S, et al. Magnesium in drinking water supplies and mortality from acute myocardial infarction in north west England. Heart 1999;82:455–60.[Abstract/Free Full Text]
  48. Reunanen A, Knekt P, Marniemi J, et al. Serum calcium, magnesium, copper and zinc and risk of cardiovascular death. Eur J Clin Nutr 1996;50:431–7.[Medline]
  49. Saris NE, Mervaala E, Karppanen H, et al. Magnesium. An update on physiological, clinical and analytical aspects. Clin Chim Acta 2000;294:1–26.[CrossRef][Medline]
  50. Durlach J. The importance of magnesium in water. In: Magnesium in clinical practice. London: John Libbey, 1988:221–2.
  51. Salonen JT, Salonen R, Seppänen K, et al. Interactions of serum copper, selenium, and low density lipoprotein cholesterol in atherogenesis. BMJ 1991;302:756–60.
  52. Salonen JT, Nyyssönen K, Korpela H, et al. High stored iron levels are associated with excess risk of myocardial infarction in eastern Finnish men. Circulation 1992;86:803–11.[Medline]
  53. Rytkonen M, Ranta J, Tuomilehto J, et al. The SPAT Study Group. The Finnish Childhood Diabetes Registry Group. Bayesian analysis of geographical variation in the incidence of type I diabetes in Finland. Diabetologia 2001;44 (suppl 3) :B37–44.
  54. Kousa A, Nikkarinen M. Geochemical environment in areas of low and high coronary heart disease mortality. In: Autio S, ed. Special paper 23. Espoo: Geological Survey of Finland, 1997:137–48.