- Research article
- Open Access
- Open Peer Review
Identifying environmental risk factors for human neural tube defects before and after folic acid supplementation
- Yilan Liao†1, 3,
- Jinfeng Wang1Email author,
- Xinhu Li†2,
- Yaoqin Guo†1 and
- Xiaoying Zheng3Email author
© Liao et al; licensee BioMed Central Ltd. 2009
- Received: 17 February 2009
- Accepted: 16 October 2009
- Published: 16 October 2009
Birth defects are a major cause of infant mortality and disability in many parts of the world. Neural tube defects (NTDs) are one of the most common types of birth defects. In 2001, the Chinese population and family planning commission initiated a national intervention program for the prevention of birth defects. A key step in the program was the introduction of folic acid supplementation. Of interest in the present study was to determine whether folic acid supplementation has the same protective effect on NTDs under various geographical and socioeconomic conditions within the Chinese population and the nature in which the influence of environmental factors varied after folic acid supplementation.
In this study, Heshun was selected as the region of interest as a surrogate for helping to answer some of the questions raised in this study on the impact of the intervention program. Spatial filtering in combination with GIS software was used to detect annual potential clusters from 1998 to 2005 in Heshun, and Kruskal-wallis test and multivariate regression were applied to identify the environmental risk factors for NTDs among various regions.
In 1998, a significant (p < 0.100) NTDs cluster was detected in the west of Heshun. After folic acid supplementation, the significant clusters gradually moved from west to east. However, during the study period, most of the clusters appeared in the middle region of Heshun where more than 95 percent of the coal mines of Heshun are located. For the analysis, buffer regions of the coal mine zone were built in a GIS environment. It was found that the correlations between environmental risk factors and NTDs vary among the buffer regions.
This suggests that the government needs to adapt the intervention measures according to local conditions. More attention needs to be paid to the poor and to people living in areas near coal mines.
- Normalize Difference Vegetation Index
- Coal Mine
- Environmental Risk Factor
- Folic Acid Supplementation
- Buffer Region
Birth defects, as defined by the March of Dimes Birth Defects Foundation, refer to any anomaly, functional or structural, that presents in infancy or later in life and is induced by events preceding birth, whether inherited, or acquired. Varying from minor cosmetic irregularities to life-threatening disorders, birth defects are the major cause of infant mortality and a leading cause of disability . NTDs (neural-tube birth defects) are one of the most common forms of birth defects, often occurring between the third and fourth weeks of gestational age. They result in structural defects that occur anywhere along the neuroaxis from the developing brain to the sacrum and often result in the exposure of neural tissue .
Birth defects have a substantial public health impact on mortality, morbidity, disability, and to the cost of health care provision. Fortunately, they can often be prevented and early intervention is an important component in the minimization of their consequences. However, for the vast majority of birth defects, the etiology is still unknown.
An environmental cause is any non-genetic factor that increases the risk of a birth defect for the exposed individual . Such factors include fetal infection , maternal illness , nutritional deficiencies [6, 7], drug ingestion , chemical exposures , air pollution , radiation , etc. Different birth defects may be caused by different risk factors. In this study, we have limited our research to NTDs.
In recent years, many countries have made efforts to reduce the incidence of birth defects. China has the largest population in the world, and also has the highest incidence of birth defects, with the highest levels compared with the national average being observed in the Shanxi Province . Shanxi also has the highest rate of NTDs in the world. In 1999, the State Family Planning Commission launched the Birth Defect Intervention Project across the country. The four counties of Pingyao, Heshun, Zhongyang and Pingding in Shanxi were designated as the pilot regions for the project carried out in 2001. The intervention work mainly involves dynamic monitoring before and after pregnancy and again following childbirth, and as a result is referred to as a three-tiered intervention. The first-tier intervention service covers newlywed couples and couples who were entitled to have a second child. Before pregnancy, couples were offered counseling and advised to take supplements, including folic acid and iodine. The second-tier intervention covered all pregnant women. In middle and west China, women in some experimental counties took a vitamin supplement containing 0.4 mg folic acid everyday from 3 months before pregnancy to 3 months after pregnancy. However, a 4-mg folic acid supplement daily was advised for women with a previous pregnancy involving a neural-tube defect and women who have taken antiepileptic drugs before. The third-tier intervention covered the effective follow-up treatment of all children who were already affected by birth defects. In this national intervention project, the main measure taken to reduce the risk of NTDs was to counsel women of child-bearing age to take folic acid supplements .
For this analysis, we included all live and still births occurring in Heshun from January 1, 1998 to December 31, 2005, born to women at the hospital or at home, and who were residents of the county during that time period. We also included all therapeutic abortions in residents in that area where the estimated date of delivery fell in the time period of interest. All NTDs cases, regardless of pregnancy outcome, were verified by doctors in the hospitals. Records of NTDs cases were collected from the local family planning department. The NTDs in the study included anencephaly, spina bifida, encephalocele, and holoprosencephaly, among others. In addition, the Ministry of Science and Technology of the People's Republic of China provided ethics approval for our study.
The local planning department declined to provide identifiers to link substantiated NTDs cases to births, so we were unable to conduct the study at an individual level; instead, we conducted an ecological study. That is, we examined NTDs cases in relation to the environmental characteristics of the villages in which their mothers lived. There were 326 villages and one town in the study area. Since the main object of this study was to find the changes in the relationships between environmental risk factors and NTDs, the town was not included as the environmental factors there are somewhat complex. In addition, birth defect registers in the town were not included.
In the study, the environmental factors of various villages were classified into socioeconomic and geographical factors. The socioeconomic factors reported useful information on medical conditions (the number of doctors), the per-capita incomes (per-capita net incomes), the agricultural chemical exposures (the use of fertilizers and pesticides), and the crop yields (vegetable and fruit productions) of every village during the 1998-2005 period. All socioeconomic data were provided by the Heshun statistical bureau. The geographical factors included elevation, gradient, vegetations coverage (normalized difference vegetation index), access condition (distance to main roads), water shortage condition (distance to rivers) and the geological background (distance to faults) of the villages. The normalized difference vegetation index (NDVI) is an index that provides a standardized method of comparing vegetation greenness between satellite images and can be used as an indicator of relative biomass and greenness. The source of the NDVI dataset used in the study is the VITO (Flemish Inst. Technological Research, Belgium), http://www.vgt.vito.be. In addition, Li et al.  found that the occurrence ratio of NTDs in Heshun had a significant negative correlation with increasing-distance from faults. So we chose the distance between villages and faults as an important influence factor.
Geographic Information System (GIS) is a computer package used to store, manage, analyze, and map geographical data. It plays a significant role in data processing and has insuperable advantages over traditional methods. GIS allows the addition of relevant layers which can be used for analyzing the spatial relationships among selected factors. Also, it offers database capabilities that can handle attributes data effectively. Attribute calculations are simple and relatively accurate. We used ARCGIS 9.0i as the GIS platform to locate the 326 villages and to quantify the selected factors.
We assumed that the pilot work of the intervention project began to show effects by 2002. Because NTDs are rare, the NTDs cases reported during the 1998-2001 period and during the 2002-2005 period were combined for the purpose of calculating the occurrence ratios before and after intervention. As a result of the Birth Defect Intervention Project, the rate of NTDs in Heshun reduced from 261.7 per 10,000 births to 215.5 per 10,000 births during the 1998-2005 period.
Spatial filtering method
The spatial analysis of birth-related morbidity provides a way of identifying inadequate health care access as well as potential environmental and behavioral causal factors. Cluster analysis is frequently used to identify an unusually high occurrence of morbidity that is clustered in space and time . The spatial filtering methodology employs non-parametric statistical techniques as an exploratory spatial data analysis tool. It can provide estimates for the unknown risk or relative risk without making distributional assumptions. The idea behind this technique is to create a series of circular moving windows or disks centered on point locations on a regular grid that covers the study area. For each of the disks, the rate is computed as the crude rate of events to population, or, alternatively, as the ratio of cases to controls. The resulting values at the grid point are then represented in smooth form by means of contour plots, a three-dimensional surface, or similar techniques . This method has been used previously to study clusters of congenital anomalies, infant mortality, and other forms of birth morbidity [15–17]. In the study, we assumed that the probability of a case resulting in a NTDs was equal to the proportion of all births in the region that resulted in NTDs. A 'smoothed' probability map could be drawn where the significance levels of high rates of NTDs by percentage for each individual circle was calculated and mapped in isarithmic form.
The data were statistically analyzed using SPSS Win 13.0 using a 95% level of significance. Kruskal-Wallis test is a nonparametric procedure that can be used to compare more than two populations in completely randomized design. So Kruskal-wallis test was used to compare the average rates of NTDs in the various buffer regions of coal mines. A linear regression analysis was used to investigate the relationship between environmental factors and NTDs. A bivariate correlation analysis was then applied to explore the relationship among different factors. Villages in which the number of births over a five-year period was less than 5 were excluded in the analysis to ensure the reliability of the conclusions. In addition, the crude ORs and 95% CI was applied to explore the association between residences near coal mines and risk of NTDs.
The detecting results of different spatial filter sizes
the number of clusters
the number of grids included in clusters
The approximately 4.8 km distance interval between grid intersections resulted in 204 grid points in Heshun. Meaningful NTDs rates were estimated for the grid points that had at least 5 births within a 4.8 km vicinity. The spatial filter circle surrounding each of these points is the area from which an estimate of the NTDs rates is made. We counted the normal births and NTDs within the circle and assigned the observed rate to that location. When we repeated this for a grid of such estimates, we were able to interpolate the NTDs rates as a continuous spatial distribution. Neighboring grid points share circular patterns that overlap, thereby sharing observations. Isarithmic maps with a constant range of values were constructed in GIS after the estimated rates were assigned to grid points.
As seen in Figure 3, clusters of NTDs were often seen in the buffer regions of 0-6 and 7-12 km during the study period. Although the government had initiated a birth defects intervention project two years previously, the hot spots for the NTDs continued to appear in these buffer regions into the years 2004 and 2005. Correspondingly, the buffer region of 31-42 km had only one cluster of NTDs in 1998.
Kruskal-Wallis test results of the average NTDs rates of the buffer regions before and after folic acid supplementation
Degrees of Freedom
Multivariate regression analysis
Multivariate regression analyses of the association between environmental factors and NTDs as the dependent variable
Model (R2) and Independent Variables
Rates for the 1998-2001 period in the buffer region of 7-12 km (0.415)
Rates for the 2002-2005 period in the buffer region of 7-12 km (0.242)
Distance to faults
Rates for the 1998-2001 period in the buffer region of 13-18 km (0.191)
Rates for the 2002-2005 period in the buffer region of 13-18 km (0.711)
distance to roads
distance to rivers
Rates for the 1998-2001 period in the buffer region of 25-30 km (0.486)
Rates for the 2002-2005 period in the buffer region of 25-30 km (0.362)
The study suggests that the effect of folic acid supplementation varied among the different subzones in Heshun. The NTDs rates in the buffer regions of 7-12, 13-18, and 25-30 km decreased significantly due to folic acid supplementation. However, folic acid supplementation did not play a role in reducing the rate in the buffer regions of 0-6, 19-24, and 31-42 km. In fact, the rate for the sub-region near the coal mines zone was actually found to increase. The environmental risk factors in these sub-regions are not the same. It is significant that we were able to identify environmental correlates of NTDs among the sub-regions.
The buffer region of 0-6 km is the area close to the coal mines. Spatial filtering detected clusters of NTDs in this area almost every year. Infants there had a significantly higher risk of NTDs compared with those living in other areas throughout the whole of the study period (OR = 1.347; 95% CI = 1.004-1.808). So, in addition to the three-tiered intervention program, the government should strengthen its local environmental governance.
Many of the villages in the buffer region of 7-12 km are located at high altitudes and located at significant distances from rivers. This area is an economically deprived area of Heshun and lacks agricultural production materials such as fertilizer and pesticide. Before intervention, there was a poor state of health amongst the population and many villages had no doctor. People generally saw a doctor at the rural clinics of the neighboring villages. Some of the records of NTDs cases were collected from clinics and hospitals, so the number of doctors was in direct proportion to the NTDs incidence during the 1998-2001 period. With the implementation of interventions, the number of rural doctors increased substantially. The production of fruit, distance to factories and to faults were correlated with NTDs in these areas. It was unexpected that the number of NTDs was found to increase with increasing distance from factories. Taking the spatial interaction into account, we qualified the environmental factors of neighboring villages within a specified distance and then explored the relationships between them and the distance to factories. The number of doctors (r = -0.364, p = 0.029), NDVI (r = -0.420, p = 0.011), vegetable production (r = -0.502, p = 0.002) and per capita net income (r = -0.420, p = 0.011) were all negatively correlated with this variable. So the factor variable 'distance to factories', could be taken as a composite indicator of the peripheral environment. Faults indicate intensive geological structural movement. And the concentration of radon in soil, water, and air near fault zones were found to be far higher than away from faults . However, it needs to approve in the future whether a fetus in utero is at particular risk of damage from radiation of radon. In previous work, the higher rates of certain types of NTDs, such as anencephaly and spina bifida, were found to be associated with living at lower altitude . However, in our study, the elevation was found to be in direct proportion to the NTDs for the 1998-2001 period in the 13-18 km buffer region. Similarly, we tried to understand the relationship between the neighboring environmental characters of villages where NTDs cases occurred and their elevations. In the neighborhood of villages with higher altitudes, the number of rural doctors (r = -0.473, p = 0.015), the use of chemical fertilizers (r = -0.821, p = 0.000), the fruit yield (r = -0.645, p = 0.000), and the per capita net income (r = -0.535, p = 0.005) of neighbors were relatively small. The difficult living environments made it more likely for the infants in these villages to contract NTDs. This is different from the 7-12 km buffer region. The 13-18 km buffer region was covered by dense vegetation, and the vegetable yield there continued to increase during the 2002-2005 period.
Nevertheless, limited medical access, and conditions of water shortage were the main handicaps in reducing NTDs prevalence in these areas.
Although no selected environmental risk factors related directly to NTDs in the 19-24 km buffer, we found that the use of fertilizers and pesticides was the highest in the county. Moreover, the production of fruit in this buffer was far less than the one in the 13-18 km buffer. This may be used to explain why the rate curve of NTDs there suddenly became higher.
The 25-30 km buffer region used more fertilizers than the other areas before 2002. The multivariate regression analysis found that the use of fertilizers was directly related to NTDs rates there at that time. Fertilizers are organic or inorganic substances, either natural or synthetic, used to supply elements (such as nitrogen, phosphate and potash) essential for plant growth. However, over application of nitrogen fertilizers and loss of phosphate fertilizer due to soil erosion may result in some undesirable environmental problems such as increase of nitrate in water and eutrophication of water bodies . Croen et al.  found that exposure to nitrate in drinking water at high concentrations was associated with increased risk for anencephaly. Moreover, the surrounding environment of the villages with higher elevations was also badly affected. So elevation became the other factor which associated with NTDs in this region. After 2001, this sub-region underwent a period of strong economic growth and the production of vegetables increased year by year. Infants in the villages with a higher slope had a higher risk of NTDs. From our investigation, it was found that most of these villages are in mountainous regions or regions lacking in water. There is significant scope to further reduce NTDs in this sub-zone by improving the quality of life of these residents.
However, we have not found any selected influence factors of NTDs in the buffer region of 31-42 km yet. So we cannot give any explanation for no change in NTDs there and attempt to deal with this problem in further work.
Although we found that folic acid supplementation had different effects on the various buffer regions, there are some limitations associated with our choice of methodology. A concern is that there may have been a significant amount of underreporting in this area, as some pregnant women would have chosen to have home births rather than hospital births. Therefore, we were limited to the use of data from hospital records and investigations undertaken in the villages. Also, some of the women may have relocated during their pregnancies, so there could be a migration bias in risk factor identification. In addition, our study is limited to a ecological analysis since we were unable to obtain any individual information otherwise required. The measures of environmental exposure are limited. Causes of individual NTDs cases cannot be identified, and the maternal characteristics of each NTDs is unknown.
Given that the causes of NTDs vary among sub-regions, it is sensible to adopt different intervention methodologies for each region. More attention should be paid to areas near coal mines. Needless-to-say, there is also an urgent need to try to improve the quality of life of people living in remote mountainous areas. Although in this study we only had the opportunity to investigate the effects of the birth defects intervention project in Heshun, a rural region in the north of China, the analysis methodology may provide some good ideas to Chinese government to improve treatment policies and intervention measures for birth defect in the future.
This work was supported by the Project of the National Natural Science Foundation of China (70571076&40471111), the Hi-Tech Research and Development Program of China (2006AA12Z15), the National Basic Research Priorities Program (2001CB5103) of the Ministry of Science and Technology of the People's Republic of China, and Knowledge Innovation Program of the CAS (KZCX2-YW-3-8). And the authors thank Ms Li Zhao for collecting the socio-economic data of the study area.
- Carmona RH: The global challenges of birth defects and disabilities. Lancet. 2005, 366: 1142-1144. 10.1016/S0140-6736(05)67459-4.View ArticlePubMedGoogle Scholar
- Frey L, Hauser WA: Epidemiology of neural tube defects. Epilepsia. 2003, 44: 4-13. 10.1046/j.1528-1157.44.s3.2.x.View ArticlePubMedGoogle Scholar
- Dolk H: Epidemiologic approaches to identifying environmental causes of birth defects. Semin Med Genet. 2004, 125C: 4-11. 10.1002/ajmg.c.30000.Google Scholar
- Hansen M, Kurinczuk JJ, Bower C, Webb S: The risk of major birth defects after intracytoplasmic sperm injection and in vitro fertilization. N Engl J Med. 2002, 346: 25-730. 10.1056/NEJMoa010035.View ArticleGoogle Scholar
- Gregg NM: Congenital cataract following German measles in the mother. Trans Ophthalmol Soc. 1942, Austr3: 35-46.Google Scholar
- Hook EB, Czeizel AE: Can terathanasia explain the protective effect of folic-acid supplementation on birth defects?. Lancet. 1997, 350: 513-515. 10.1016/S0140-6736(97)01342-1.View ArticlePubMedGoogle Scholar
- Brender JD, Zhan FB, Suarez L, Langlois P, Gilani Z, Delima I, Moody K: Linking environmental hazards and birth defects data. Int J Occup Environ Health. 2006, 12: 126-133.View ArticlePubMedGoogle Scholar
- Schardein JL: Chemically induced birth defects. 2000, New York, Marcel Dekker, 3View ArticleGoogle Scholar
- Brender JD, Suarez L, Felkner M, Gilani Z, Stinchcomb D, Moody K, Henry J, Hendricks K: Maternal exposure to arsenic, cadmium, lead, and mercury and neural tube defects in offspring. Environ Res. 2006, 101: 132-139. 10.1016/j.envres.2005.08.003.View ArticlePubMedGoogle Scholar
- Gilboa SM, Mendola P, Olshan AF, Harness C, Loomis D, Langlois PH, Savitz DA, Herring AH: Comparison of residential geocoding methods in population-based study of air quality and birth defects. Environ Res. 2006, 101: 256-262. 10.1016/j.envres.2006.01.004.View ArticlePubMedGoogle Scholar
- Araneta MRG, Destiche DA, Schlangen KM, Merz RD, Forrester MB, Gray GC: Birth defects prevalence among infants of Persian gulf war veterans born in Hawaii, 1989-1993. Teratology. 2000, 62: 195-204. 10.1002/1096-9926(200010)62:4<195::AID-TERA5>3.0.CO;2-5.View ArticlePubMedGoogle Scholar
- Hu HT: Methods to prevent Shanxi birth defect. 2003, China.org.cn, [http://www.china.org.cn/english/2003/sep/74927.htm]Google Scholar
- Wu JL, Wang JF, Meng B, Chen G, Pang LH, Song XM, Zhang KL, Zhang T, Zheng XY: Exploratory spatial data analysis for the identification of risk factors to birth defects. BMC Public Health. 2004, 4: 23-33. 10.1186/1471-2458-4-23.View ArticlePubMedPubMed CentralGoogle Scholar
- Li XH, Wang JF, Liao YL, Meng B, Zheng XY: A geo-analysis for the environmental cause of human birth defects. ET&C. 2006, 88 (3): 551-559.Google Scholar
- Esra O, Bryan LW, Su YK, Melina SM: Comparison of spatial scan statistic and spatial filtering in estimating low birth weight clusters. Int J Health Geogr. 2005, 4: 19-28. 10.1186/1476-072X-4-19.View ArticleGoogle Scholar
- Rushton G, Lolonis P: Exploratory spatial analysis of birth defect rates in an urban population. Stat Med. 1996, 15: 717-726. 10.1002/(SICI)1097-0258(19960415)15:7/9<717::AID-SIM243>3.0.CO;2-0.View ArticlePubMedGoogle Scholar
- Talbot TO, Kulldorff M, Forand SP, Haley VB: Evaluation of spatial filters to create smoothed maps of health data. Stat Med. 2000, 19 (Part 17/18): 2399-2408. 10.1002/1097-0258(20000915/30)19:17/18<2399::AID-SIM577>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Castilla EE, Lopez-Camelo JS, Campana H: Altitude as a risk factor for congenital anomalies. Am J Med Genet. 1999, 86: 9-14. 10.1002/(SICI)1096-8628(19990903)86:1<9::AID-AJMG3>3.0.CO;2-X.View ArticlePubMedGoogle Scholar
- Cao ZH: Environmental issues related to chemical fertilizer use in China. Pedosphere. 1996, 6 (4): 289-303.Google Scholar
- Croen LS, Todoroff K, Shaw GM: Maternal exposure to nitrate from drinking water and diet and risk for neural tube defects. AM J Epi. 2001, 153 (4): 325-331. 10.1093/aje/153.4.325.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2458/9/391/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.