Introduction
In early 1993 media reports alleged clusters of anophthalmia and microphthalmia in England, and postulated that these might be linked to exposure to the pesticide (fungicide) Benomyl (DuPont, Wilmington, DE).[1] This was not the first time that alleged clusters of anophthalmia and microphthalmia had been reported in relation to environmental exposure: earlier reports dating from 1986 investigated clusters of microphthalmia in connection with two high temperature waste incinerators in Wales and Scotland.[2 3]
Immediate assessment of the English clusters was complicated by two factors. Firstly, there was a lack of diagnostic detail in the reports, which meant that an expected prevalence of the conditions could not be confidently asserted.[4] Secondly, there was the so called Texas sharpshooter problem (a Texan fires randomly at a barn door and subsequently draws a bull's eye in the densest cluster of bullet holes), a usual dilemma in the post-hoc assessment of clusters. Even on a random, and therefore non-uniform, spatial pattern it is possible to draw boundaries around apparent clusters of cases in such a way that the density of the cases in that area far exceeds expectation.[5] However, this has no statistical validity without reference to the background pattern of all cases from whom the dusters were identified.
We report on the results of a study set up in response to public concern to establish the presence or absence of any geographical variation in anophthalmia and microphthalmia, including large scale regional variation, excess prevalence in rural areas, or localised clustering. The analysis is based on a register of all cases of anophthalmia and microphthalmia born in England in 1988-94, which was established for the purposes of this study.[6] The time period was chosen to overlap as far as possible with the period of concern and to collect enough cases for geographical analysis, but not to go back too far and risk major underascertainment of cases in earlier years.
Subjects and methods
Data
Our study is based on a register of all cases of anophthalmia and microphthalmia born in England in 1988-94. The methodology by which the register was established is described elsewhere.[6] There were 444 cases registered, excluding cases of trisomy 13 and holoprosencephaly or cyclops.[6]
As there is no formal cut off point between mild microphthalmia and the normal eye it was expected that mild microphthalmia might show considerable geographical variation due to diagnostic and reporting differences alone. Severity of microphthalmia was based on the information supplied by clinicians in specially designed questionnaires.[6] A subgroup of 237 (53.4%) severe cases was defined by excluding 113 (25.5%) mild cases and 94 (21.2%) cases where the severity of microphthalmia was unknown.[6]
A further subgroup of 189 (79.7%) severe cases of unknown aetiology was defined by excluding 48 (20.3%) cases with known aetiology from the severe subgroup. The aetiology was established independently by a medical geneticist (Robin Winter) based on information obtained in questionnaires. The 48 cases included those with a syndrome of genetic origin, a family history strongly suggestive of genetic origin, bilateral conditions with parental consanguinity, and maternal infection.
The Office for National Statistics supplied postcoded birth registration data for the entire study period. There were 4 538 790 live births during 1988-94.
Postcodes of cases and all births allowed geographical localisation (grid reference with 100 m resolution) and assignment to census enumeration districts or wards where necessary. Cases were postcoded to residence at birth but in 33 (7.4%) cases (including 10 (4.2%) severe cases) the address at birth could not be obtained from the notification process.[6] Of all births from the Office for National Statistics' database, 12 616 (0.3%) were without valid postcodes.
Cases and births were grouped by regional health authority and district health authority of residence at birth according to the health authorities' boundaries before 1994 (table 1). In one (0.23%) case the region was unknown and in 14 (3.2%) cases the district was unknown. The urban or rural nature of the area of residence was determined using the population density of the 1991 census ward of residence. Wards were classed into five groups by dividing at quintiles of population density (table 2). Three per cent of births were in the group with the lowest population density, the proportion in this rural group in individual regions ranging up to 11% in East Anglia. As a validation of the main urban and rural analysis, an indicator of the urban or rural (dichotomous) nature of 1991 census enumeration districts was obtained from the Office for National Statistics. According to this indicator, 9% of births occurred in a rural setting. The state of socioeconomic deprivation of the enumeration district of residence was assessed according to the deprivation index developed by Carstairs.[7] This was standardised to Great Britain, and five socioeconomic groups derived corresponding to fifths of enumeration districts.
Table 1 Number of cases (all and severe) of anophthalmia and microphthalmia, prevalence per 10 000 births, and proportion of cases of unknown severity by region of residence in England, 1988-94
(*) Prevalence estimated assuming that proportion of severe cases in additional cases of unknown severity is same as proportion in cases of known severity.
Table 2 Number of cases (all, severe, severe of unknown aetiology) of anophthalmia and microphthalmia and prevalence per 10 000 births, by urban or rural residence (measured by population density)
(*) Prevalence estimated assuming that proportion of severe cases in additional cases of unknown severity is same as proportion in cases of known severity.
The estimated prevalence of severe cases was calculated by assuming that the proportion of such cases among cases of unknown severity was the same as the proportion of severe cases among those of known severity. For statistical analysis the denominator (number of births) rather than the numerator (number of cases) was manipulated in these estimations to avoid false inflation of statistical significance.
Statistical analysis
Comparison of prevalence of anophthalmia and microphthalmia across regions and groups defined by population density and deprivation was informed by the Pearson and Armitage [chi square] tests, and where necessary adjusted for confounding using Poisson regression. Variation in underlying prevalences across regions and across districts was estimated by the Martuzzi-Hills method,[8] which removes the random sampling variation expected in small numbers. This variation is expressed as the 5th to 95th centile range of prevalence ratios, relative to the overall prevalence. Expected numbers were calculated in each district of residence (for all cases and severe cases), stratified for region and population density fifth.
We used two tests of localised clustering: Cuzick-Edwards[9] and Diggle-Chetwynd."[10] For both these tests we used three controls per case randomly chosen from all births. For each case the Cuzick-Edwards test counts the number of other cases in its k nearest neighbours among all remaining cases and controls (we repeated with k= 1, 2, 3 ... 8). For each case the Diggle-Chetwynd test counts the number of cases within k kilometres (we repeated with k = 1, 2, 3 ... 50). Evidence for localised clustering exists if the sum of these counts is significantly greater than would be expected if cases were an independent sample of births.
The data included two pairs of siblings (all four cases were severe) who strongly influenced the results at short distances (1 nearest neighbour or 1 km distance). One of each pair was therefore removed from the clustering analyses on the basis that two members of a pair were probably not independent events. Both pairs had been excluded from the most restrictive subgroup of severe cases of unknown aetiology as most probably of genetic aetiology.
Results
The overall prevalence of anophthalmia and microphthalmia in England was 1.0 per 10 000 births. Table 1 shows the regional prevalence of all cases. Regional variation in prevalence did not reach statistical significance (P = 0.07 for all cases, P = 0.76 for severe cases). The 5th to 95th centile range of regional to national prevalence ratios for all cases estimated by the Martuzzi-Hills method was 0.83 to 1.22. There was no statistically significant heterogeneity in prevalence of all cases across districts (P [is greater than] 0.20), with an estimated 5th to 95th centile range of the ratios of observed to expected numbers of cases of 0.83 to 1.19.
The prevalence of anophthalmia and microphthalmia increased with decreasing population density--that is, it was higher in rural areas (table 2). The relative risk was 1.79 (95% confidence interval 1.15 to 2.81) in the most rural group compared with the most urban group for all cases, and 2.37 (1.38 to 4.08) for severe cases (table 2). Among the 22 most rural cases, 36% (eight cases) were bilateral compared with an average of 35% (157 cases) overall. Poisson regression showed that urban and rural variation was reduced after controlling for region, both for all cases (relative risk for most rural fifth compared with most urban fifth was 1.61 (95% confidence interval 1.01 to 2.58)) and for severe cases (relative risk for most rural fifth compared with most urban fifth was 2.25 (1.29 to 3.95)). The validation analysis using the Office for National Statistics' urban and rural indicator showed that 20 out of 22 cases within the areas of lowest population density lived in enumeration districts classified as rural. Forty (9%) cases were classified as rural. The odds ratio for living in a rural enumeration district was 1.21 (0.85 to 1.72).
Prevalence of anophthalmia and microphthalmia varied little by socioeconomic deprivation of enumeration districts (table 3). Relative risk for the most deprived group compared with the least deprived group was 0.92 (0.68 to 1.24). Tests for both heterogeneity and trend showed that the variation was easily explained by chance (P [is greater than] 0.20).
(*) Groups defined by dividing at quintiles of population of Great Britain: 1=most affluent, 5=most deprived.
Neither the Cuzick-Edwards (table 4) nor the Diggle-Chetwynd tests showed overall evidence of localised clustering at a national level (after removal of one of each of two pairs of siblings) whether or not controls were frequency matched to cases by population density fifth. Both tests were also performed for three regions separately: Trent, Northern, and Oxford. These regions had been selected a priori as having probably the most complete case ascertainment, or being areas of prior interest for clustering. The Cuzick-Edwards test showed no statistically significant clustering in these regions, while the Diggle-Chetwynd test showed statistically significantly more cases than expected within 2 km of an index case in Trent and within 50 km of an index case in the Oxford region. There was no evidence of clustering by the Cuzick-Edwards test within the most rural fifth or the two most rural fifths (table 4). The Diggle-Chetwynd test showed significant clustering in the most rural fifth only in the subgroup of severe cases of unknown aetiology, with more cases than expected within 7-9 km of other cases but based on only one case-case pair in Yorkshire. There was no significant clustering by this test when the two most rural fifths were considered.
Table 4 Observed number of cases of anophthalmia and microphthalmia within three nearest neighbours, number expected, and probability of exceeding the expected number by at least observed number, according to Cuzick-Edwards test
Table 5 gives an empirical idea of the number of cases that occurred close together. Results are given for the most rural fifth and second most urban fifth. Altogether, 48 cases had another case within 1 km and one case had two cases within 1 km. The maximum number of cases within 5 km of another case was eight, with 15 cases having 6-8 other cases within 5 km. The maximum number of cases within a 20 km radius urban area was 50. We may infer from the results of the tests for clustering that these numbers are close to those expected given the distribution of births and the prevalence of anophthalmia and microphthalmia.
Table 5 Counts of index cases of anophthalmia and microphthalmia according to number of other cases within specified distances(*) (by rural or urban residence of index case)
([dagger]) First fifth of population density (22 cases).
([double dagger]) Fourth fifth of population density (124 cases).
Discussion
Aetiological significance of the findings
We found little or no evidence of large scale geographical variation or localised clustering in anophthalmia and microphthalmia during 1988-94. However, we did find a gradient in prevalence from urban to rural, with a 80% excess prevalence (95% confidence interval 15% to 181%) in the least densely populated wards compared with the most densely populated. Further analysis, whether by controlling for region or using a dichotomous indicator of urban or rural only, tended to reduce the strength of this association, but given the a priori hypothesis of a rural excess that we set out to test we cannot dismiss this finding.
Whether a rural excess, if real, is associated with pesticide exposure, as postulated in the media,[1] or one or more other environmental or genetic risk factors would need further investigation. The evidence for a causal link between exposure to pesticide and anophthalmia and microphthalmia is currently weak. Case reports of anophthalmia and microphthalmia do exist where mothers were exposed to pesticides in early pregnancy.[11 12] Since exposure to pesticides is common the significance of the case reports to date cannot be assessed. Two epidemiological studies could not find a relation between markers of Benomyl exposure and anophthalmia and microphthalmia. The first was an ecological study which correlated regional use of the product in Italy to the prevalence of anophthalmia and microphthalmia.[13] The second study identified a cohort of farming families in Norway in 1960-89 and found only one report of anophthalmia and microphthalmia in the medical birth register for the families that were potentially exposed.[14] Other epidemiological studies have examined congenital malformations in relation to markers of maternal exposure to pesticides, but the diversity of pesticides, methods of exposure measurement, and malformations in these studies make it difficult to draw valid conclusions. Animal experiments have induced anophthalmia and microphthalmia (as well as other congenital anomalies) with exposure to Benomyl at high doses,[15] but this has not been considered to constitute a risk to humans at lower doses.[16 17]
Probably the strongest evidence for an environmental aetiology for anophthalmia and microphthalmia concerns maternal infections, and lack of spatial clustering is of interest in this regard. The relation between maternal infections and anophthalmia and microphthalmia is well established for rubella, toxoplasmosis, cytomegalovirus, and varicella.[18 19] Other viruses are implicated but with less evidence, including parvovirus B19,[20 21] herpes simplex type 2,[22] Epstein-Barr,[23] and coxsackie A9.[24] Furthermore, there is evidence that influenza, fever, or hyperthermia may cause both malformations of the central nervous system and microphthalmia.[25-28] A specific effect of hyperthermia as a cause of anophthalmia and microphthalmia is supported by animal experiments.[25 29] Hyperthermia may be a potentiating coteratogen in association with other exposure (including chemicals).
Our finding of lack of spatial clustering suggests that if a large proportion of cases can be attributed to maternal infection then either this infection must not show strong spatial clustering or it must cluster in both. space and time in such a way that clustering cannot be detected at a spatial level only. Further study of the anophthalmia and microphthalmia register will investigate variation in prevalence over time and space simultaneously in relation to known variations in maternal infections and pesticide use.
Other risk factors described in the literature for anophthalmia and microphthalmia include solvent misuse and exposure to x rays or drugs such as thalidomide, isotretinoin, warfarin, and alcohol.[19 30] Anophthalmia and microphthalmia are also associated with a wide range of genetic syndromes,[31] both chromosomal and monogenic, and it is likely that when a genetic syndrome is not evident, genetic susceptibility through one or more genes nevertheless plays a part in the response to environmental exposure.
We could find no association between prevalence of anophthalmia and microphthalmia and socioeconomic deprivation. Whatever the main risk factors are in determining the occurrence of anophthalmia and microphthalmia it seems they are not strongly related to socioeconomic factors.
Data quality and potential artefacts
The geographical analysis that we performed depends for its validity on the quality of the underlying data. As reported elsewhere,[6] we believe that about 15% (80) of cases may not have been reported to the register that we established, and that it is probable that these included mainly mild microphthalmia and children who did not survive their first year of life. It is difficult to establish a register retrospectively although our use of multiple sources of notification overcame this to a large extent. No source, even paediatricians who had seen the majority of cases, reported more than one quarter of cases to the register.[6] Underascertainment can have two main effects: to reduce the overall case number and thus statistical power (where 15% underascertainment is not a serious problem), and to introduce spurious geographical variation due to geo graphical variation in ascertainment. Since we found little or no regional or district variation or clustering we doubt if these results would essentially have changed with complete ascertainment. The theoretical possibility that underascertainment might obscure geographical variation and clustering seems unlikely.
We suspected when designing this study that mild microphthalmia might be variably reported and that this would lead to significant geographical variation. In fact our results for all cases, including mild cases, suggest that this did not occur to any great extent. Restricting the analysis to severe cases was difficult because of the unknown severity status of a proportion of cases. Nevertheless, we believe that it was helpful in determining that any geographical variation we did find, including the urban-rural gradient, was not an artefact of the variable reporting of mild cases.
Implications of the study
A geographical and clustering analysis such as that reported here is very much a first stage in looking at the environmental epidemiology of a congenital anomaly. Generalised clustering analyses may not be sensitive to some particular geographical trends, as we showed with the urban-rural gradient, and it is therefore necessary to move forward to hypothesis testing about specific exposures, whether within an ecological or case-control context. Our results are reassuring in terms of the absence of strongly localised environmental exposures causing clusters of children to be born with anophthalmia and microphthalmia, but further study is needed of the aetiology of anophthalmia and microphthalmia to prevent this condition in future generations.
We thank all the members of our advisory committee, Paul Elliott who helped launch the study, Richard Collin and Barry Jones (ophthalmologist advisers), Beverley Botting, Peter Diggle, Tony Gattrell, Ruth Gilbert, Catherine Peckham, and Robin Winter; Chris Grundy and John Charlton for supplying rural indicators, and Martine Vrijheid for help on some of the analyses.
Contributors: HD coordinated the study (particularly its design and interpretation) and drafted the paper; she will act as guarantor for the paper. AB participated in the study design (particularly data collection methods and design of the questionnaire), collected the data, and participated in the statistical analysis (particularly analyses of prevalence rates). BGA supervised the statistical analysis and informed its interpretation. PW performed various statistical analyses particularly those concerning clustering.
Funding: Department of Health.
Conflict of interest: None.
[1] Paduano M, McGhie J, Boulton A. Mystery of babies with no eyes. Observer 1993 Jan 17:3.
[2] Welsh Office. The incidence of cangenital malformations in Wales, with special reference to the District of Torfaen, Gwent. Cardiff: Welsh Office, 1985.
[3] Scottish Office. Report of a working party on microphthalmos in the Forth Valley health board area 1988. Edinburgh: Scottish Information Office, 1988.
[4] Dolk H, Elliott P. Evidence for "clusters of anophthalmia" is thin. BMJ 1993;307:203.
[5] Alexander FE, Cuzick J. Methods for the assessment of disease clusters. In: Elliott P, Cuzick J, English D, Stern R, eds. Geographical and environmental epidemiology: methods for small area studies. Oxford: Oxford University Press;1992:238-50.
[6] Busby A, Dolk H, Collin R, Jones B, Winter R. Compiling a national register of babies born with anophthalmia and microphthalmia in England 1988-94. Arch Dis Child 1998 (in press).
[7] Carstairs V, Morris R. Deprivation and health in Scotland. Aberdeen: Aberdeen University Press, 1991.
[8] Martuzzi M, Hills M. Estimating the degree of heterogeneity between event rates using likelihood. Am J Epidemiol 1995;141:369-74.
[9] Cuzick J, Edwards R. Spatial clustering for inhomogeneous populations. JR Statist Soc B 1990;52:73-104.
[10] Diggle P, Chetwynd AG. Second-order analysis of spatial clustering for inhomogeneous populations. Biometrics 1991;47:1155-63.
[11] Romero P, Barnett PG, Midtling JE. Congenital anomalies associated with maternal exposure to oxydemeton-methyl. Environ Res 1989;50:256-61.
[12] Sherman JD. Chlorpyrifos (dursban) associated birth defects: report of four cases. Arch Environ Health 1996;51:5-8.
[13] Spagnolo A, Bianchi F, Calzolari E, Clementi P, Mastroiacovo P, Meli P, et al. Anophthalmia and benomyl in Italy: a multicenter study based on 940,615 newborns. Reprod Toxicol 1994;397-403.
[14] Kristensen P, Irgens L. Clusters of anophthalmia--no link with benomyl in Norway. BMJ 1994;308:205-6 (letter).
[15] Hoogenboom ER, Rausdell JF, Ellis WG, Kavlock RJ, Zeman FJ. Effects on the fetal rat eye of maternal benomyl exposure and protein malnutrition. Curr Eye Res 1991;10:601-12.
[16] International programme on chemical safety. Environmental health criteria 148. Benomyl. Geneva: World Health Organisation, 1993.
[17] Ministry of Agriculture, Fisheries, and Food pesticides safety directorate. Evaluation of benomyl. London: MAFF, 1992.
[18] Lambert S, Hoyt C. Ocular manifestations of intrauterine infection. In: Taylor D, ed. Paediatric ophthalmology. Boston: Blackwell Scientific; 1990:91-102.
[19] Warburg M. Update of sporadic microphthalmos and coloboma. Ophthalmic Paediatr Genet 1992;13:111-22.
[20] Hartwig NG, Vermeij-Keers C, Van-Elsacker-Niele AM, Fleuren GJ. Embryonic malformations in a case of intrauterine parvovirus B19 infection. Teratology 1989;39:295-302.
[21] Wetland HT, Vermeij-Keers C, Salimans MMM, Fleuren GJ, Verwey RA. Parvovirus B19 associated with fetal abnormality. Lancet 1987;i:682-3.
[22] Von Herzen JL, Bernischke K. Unexpected herpes simplex infection in a newborn. Obstet Gynecol 1976;50:728-30.
[23] Buys ML. Eye. Birth defects encyclopaedia. Dover, MA: Center for Birth Defects Information Services, 1990.
[24] Knox EG, Lancashire RJ. Epidemiology of congenital malformations. London: HMSO, 1991.
[25] Edwards MJ. Hyperthermia as a teratogen. Teratogenesis, carcinog, mutagen 1986;6:563-82.
[26] Jones KL. Smith's recognizable patterns of human malformation. Philadelphia: WB Saunders, 1988.
[27] Spraggett K, Fraser FC. Teratogenicity of maternal fever in woman--a retrospective study. Teratology 1982;25:78A.
[28] Fraser FC, Skelton J. Possible teratogenicity of maternal fever. Lancet 1978;ii:634.
[29] Sulik KK, Cook CS, Webster WS. Teratogens and craniofacial malformations: relationships to cell death. Development 1988;103(suppl):213-32.
[30] Donald JM, Hooper K, Hopenhayn-Rich C. Reproductive and developmental toxicity of toluene: a review. Envir Health Perspect 1991;94:237-44.
[31] Warburg M. An update on microphthalmos and coloboma. A brief survey of genetic disorders with microphthalmos and coloboma. Ophthalmic Paediatr Genet 1991; 12:57-63.
(Accepted 11 June 1998)
RELATED ARTICLE: Key message
* Clusters of anophthalmia and microphthalmia in England have been alleged in the media, with hypothesised links to environmental exposure such as pesticides
* To answer concerns about clustering a register has been established of all cases of anophthalmia and microphthalmia born in England in 1988-94
* There is no large regional or district variation in prevalence
* Rural areas have a roughly twofold excess in prevalence, which requires further confirmation and investigation
* There is very little evidence for localised clustering in England in 1988-94
Editorial by Mariman
Environmental Epidemiology Unit, Department of Public Health and Policy, London School of Hygiene and Tropical Medicine, London WC1E 7HT H Dolk, senior lecturer A Busby, research fellow B G Armstrong, senior lecturer in medical statistics P H Walls, senior computing scientist
Correspondence to: Dr Dolk h.dolk@lshtm.ac.uk
BMJ 1998;317:905-10
COPYRIGHT 1998 British Medical Association
COPYRIGHT 2000 Gale Group