Association between CFL1gene polymorphisms and spina bifida risk in a California population
© Zhu et al; licensee BioMed Central Ltd. 2007
Received: 26 October 2006
Accepted: 12 March 2007
Published: 12 March 2007
CFL1 encodes human non-muscle cofilin (n-cofilin), which is an actin-depolymerizing factor and is essential in cytokinesis, endocytosis, and in the development of all embryonic tissues. Cfl1 knockout mice exhibit failure of neural tube closure at E10.5 and die in utero. We hypothesized that genetic variation within the human CFL1 gene may alter the protein's function and result in defective actin depolymerizing and cellular activity during neural tube closure. Such alterations may be associated with an increased risk for neural tube defects (NTDs).
Having re-sequenced the human CFL1 gene and identified five common single nucleotide polymorphisms (SNPs) in our target population, we investigated whether there existed a possible association between the genetic variations of the CFL1 gene and risk of spina bifida. Samples were obtained from a large population-based case-control study in California. Allele association, genotype association and haplotype association were evaluated in two different ethnicity groups, non-Hispanic white and Hispanic white.
Homozygosity for the minor alleles of the SNPs studied (rs652021, rs665306, rs667555, rs4621 and rs11227332) appeared to produce an increased risk for spina bifida. Subjects with the haplotype composed of all minor alleles (CCGGT) appeared to have increased spina bifida risk (OR = 1.6, 95% CI: 0.9~2.9), however, this finding is not statistically significant likely due to limited sample size.
The sequence variation of human CFL1 gene is a genetic modifier for spina bifida risk in this California population.
Neural tube defects (NTDs) are a group of severe congenital malformations characterized by a failure of neural tube closure during early embryonic development. NTDs are complex birth defects with a multi-factorial pattern of inheritance, requiring both genetic and environmental factors to contribute to their etiology . The development and closure of the neural tube is usually completed within 28 days post-conception, in a process that is tightly regulated yet prone to environmental perturbation . Periconceptional folic acid supplementation has been repeatedly reported to prevent 50~70% of NTDs [3–5]. Mutations and polymorphisms in folate pathway genes such as methylenetetrahydrofolate reductase (MTHFR), methionine synthase (MTR), methionine synthase reductase (MTRR), and betaine-homocysteine methyltransferase (BHMT), have been intensively investigated and in some studies have shown an association with NTD risk [6–9]. However, none of these factors individually contributes substantially to the population burden of NTDs.
In addition to folate, other developmental mechanisms have been postulated as contributors to abnormal neural tube development. Animal models have provided crucial mechanistic information and possible candidate genes to explain susceptibility to NTDs . More than 80 genetic mouse models exhibit NTDs, with new ones emerging from gene targeting studies and large-scale mutagenesis screens on a regular basis . A survey of the genes whose disruption causes NTD indicates multiple key signaling pathways and cellular functions that are essential for neural tube closure. One such gene candidate involves non-muscle cofilin (n-cofilin), an actin-depolymerizing factor. N-cofilin, encoded by the CFL1 gene, is essential for cytokinesis, endocytosis, and plays a critical role in the development of all embryonic tissues. Inactivation of the Cfl1 gene in mice results in embryolethality and failure of neural tube closure by E10.5 . It has been suggested that the neural tube closure defects are due to compromised delamination and migration of neural crest cells in these animals. In vitro migration assays performed on neural crest cells from these knockout embryos demonstrated limited traveling distance, failure of cell polarization, and a lack of F-actin structures such as fibers, bundles and cortical F-actins . These findings suggest that n-cofilin regulates cytoskeletal dynamics during neural crest migration.
The human CFL1 gene (NM_005507), which maps to chromosome 11q13.1 , contains four exons and encodes an 18.5 kDa phosphoprotein. Human n-cofilin protein shares 98.8% homology with mouse protein, and the human CFL1 gene is 92.9% homologous with the mouse gene. In this study, we re-sequenced the genomic region on human chromosome 11 which encompasses the CFL1 gene, and tested the hypothesis that genetic polymorphisms in human CFL1 gene may modify human spina bifida risk. This hypothesis was evaluated in a population-based case-control study of infants with spina bifida.
Epidemiological data and biological specimens were derived from the California Birth Defects Monitoring Program, a population-based active surveillance system for collecting information on infants and fetuses with congenital malformations . Program staff collected diagnostic and demographic information from multiple sources of medical records for all live-born and stillborn fetuses (defined as 20 weeks gestation) and pregnancies electively or spontaneously terminated. Nearly all structural anomalies diagnosed within one year of delivery were ascertained. Overall ascertainment has been estimated as being 97% complete .
246 cases (infants with spina bifida) and 336 controls (non-malformed infants) were included in this study. Among the 246 cases, 86 (35.0%) were non-Hispanic white, 128 (52.0%) were Hispanic white, and 32 (13.0%) were of other ethnicities (African American, Asian, etc.). Among the 336 controls, 154 (45.8%) were non-Hispanic white, 113 (33.6%) were Hispanic white, and 69 (20.5%) were of other ethnicity (African American, Asian, etc.). These cases and controls were derived from 1983–86 and 1994–95 birth cohorts in selected California counties. Each case and control infant was linked to its newborn bloodspot, which served as the source of DNA in our genotyping analyses. All samples were obtained with approval from the State of California Health and Welfare Agency Committee for the Protection of Human Subjects. Genomic DNA was extracted from dried newborn screening bloodspots using the Puregene DNA Extraction Kit (Gentra, Minneapolis, MN).
Re-sequencing of CFL1 gene
Primer sequences for re-sequencing of CFL1 gene
Primer set name
Size of Amplicon (base pair)
Five common single nucleotide polymorphisms (SNPs) were identified from the re-sequenced CFL1 gene. SNPs rs652021, rs665306, and rs667555 were genotyped using TaqMan SNP assays (Assay-on-Demand, Applied Biosystems, Foster City, CA) on an ABI PRISM 7900 Sequence Detection System (Applied Biosystems, Foster City, CA) following the manufacturer's instructions. SNPs rs4621 and rs11227332 were genotyped using direct sequencing analyses. Genotyping assays were performed by laboratory technicians who were blinded to the case-control status of the samples. Five percent of the samples were duplicated in order to assess possible genotyping error.
Haploview software version 3.3.2 (Daly Lab at the Broad Institute, Cambridge, MA) was used to perform haplotype analysis . Haploview software predicts haplotype based on EM algorithm. Pair-wise Linkage Disequilibrium (LD) between markers was measured by D' (defined as the linkage disequilibrium measure, D, divided by the theoretical maximum for the observed allele frequencies) and the correlation coefficient r2. Spina bifida risk was measured in three ways: allelic association, haplotype association, and genotype association. Single SNP allele association and haplotype association analyses were performed using Haploview software. For each SNP, minor allele frequency (MAF) and heterozygosity (HET) were computed, and deviation from Hardy-Weinberg Equilibrium (HWE) was tested among control infants. For each haplotype, counts for case and control association tests were obtained by summing the fractional likelihoods of each individual. For example, if an individual was determined by the EM algorithm to have a 40% likelihood of haplotype A and 60% likelihood of haplotype B, 0.4 and 0.6 would be added to the counts for A and B respectively. 95% confidence intervals for allele frequencies and haplotype frequencies were calculated using an online program, VassarStats . Spina bifida risks were measured by odds ratios (ORs). Haplotype-specific odds ratios were computed taking the most common haplotype as reference. For genotype analysis, ORs were computed for each SNP by logistic regression models utilizing SAS software (version 9.1). Some models were adjusted for race/ethnicity (defined as non-Hispanic white and Hispanic white).
Characteristics of CFL1 SNP markers and allelic association analysis for spina bifida, non-Hispanic White vs. Hispanic White
Physical Location (NCBI Build 36.1)
Location in Gene
SNP Function Prediction§
HWE p value§§
Minor Allele Frequency
(95% Confidence Interval)
Allelic Association P value
no known function
splicing regulation (protein domain abolished)
0.253 (0.192, 0.326)
0.399 (0.328, 0.474)
0.386 (0.315, 0.461)
Pairwise LD of CFL1 SNP markers in different ethnicities
Associations between haplotypes and spina bifida risk, non-Hispanic White vs. Hispanic White
Freq. (95% CI) §
Estimated Haplotype Counts
Association between CFL1 genotypes and spina bifida risk: Odds Ratio estimation by Logistic Regression
No. of Control (%)
No. of Spina Bifida (%)
Crude OR (95% CI)
Non-Hispanic White OR (95% CI)
Hispanic White OR (95%CI)
1.5 (0.9, 2.2)
1.1 (0.6, 2.0)
1.8 (0.9, 3.7)
1.6 (1.0, 2.6)
2.8 (1.2, 6.5)
1.7 (0.8, 3.6)
1.4 (0.9, 2.2)
1.1 (0.6, 2.0)
1.7 (0.9, 3.5)
1.6 (1.0, 2.6)
2.8 (1.2, 6.5)
1.6 (0.8, 3.3)
1.3 (0.9, 2.0)
1.3 (0.7, 2.4)
1.6 (0.8, 3.0)
2.8 (1.1, 7.0)
5.1 (1.3, 20.2) §
2.8 (0.6, 14.4)
1.3 (0.9, 2.0)
1.0 (0.6, 1.9)
1.9 (0.9, 3.8)
1.5 (0.9, 2.4)
2.3 (1.0, 5.2)
1.7 (0.8, 3.5)
1.4 (0.9, 2.2)
1.1 (0.6, 2.0)
1.9 (0.9, 4.0)
1.4 (0.9, 2.3)
2.2 (0.9, 5.5)
1.6 (0.7, 3.4)
We performed functional prediction using an online program, FASTSNP, which predicts functional impact of SNPs according to current knowledge. FASTSNP provides a "risk score" for each SNP based on its putative biological function. Among the five SNPs we studied, rs4621, which is a synonymous change located in exon 1, has the highest "risk score". SNP rs652021 has a risk score of zero, which means that it has essentially no known functional impact. Strong linkage disequilibrium (LD) was observed among the five SNPs in both sub-populations (non-Hispanic white and Hispanic white) in our study. Therefore, the observed associations of these SNPs with spina bifida could be due to the strong LD to the "real" functional variation.
Discussion and conclusion
To our knowledge, this is the first investigation of human CFL1 gene polymorphisms as risk factors of spina bifida risk. Infants with the CGGTT haplotype had an increased spina bifida risk in the two studied ethnic groups. The observed increases in spina bifida risks indicated a possible role of actin depolymerizing factor in human neural tube morphogenesis.
One possible mechanism underlying the observed results that can be suggested follows. During early embryonic development, genes that regulate various morphogenetic activities must function cooperatively for the neural tube to close properly. Neural crest cell migration is one of these critical activities . Neural crest cells migrate from the neural tube to the periphery and give rise to a variety of cell types. Therefore faulty neural crest cell migration might interfere with normal neural tube closure. This suggestion is bolstered by observations derived from experimental models. Mammals have at least three highly conserved genes encoding F-actin polymerizing factors: n-cofilin , m-cofilin (muscle cofilin) , and ADF (actin depolymerizing factor) . Actin depolymerization is one of the key activities required for actin-driven motility . Recently, detailed analysis of n-cofilin function in mouse embryonic development using gene knockout technology suggested that n-cofilin is essential for neural crest cell migration . Our study, from a population perspective, revealed that genetic variations of actin depolymerizing factor n-cofilin among infants may contribute to the risk of spina bifida.
The strengths of this study include its population-based ascertainment of cases and controls and its evaluation of race/ethnicity as a potentially important modifier of risk in the presence of variant genotypes and haplotypes. Population variation needs to be evaluated in population-based case control studies of gene-risk association. Haplotype association analysis in both sub-populations suggested that haplotype CGGTT, which is comprised of minor alleles of all five SNPs, conferred an increased risk of spina bifida. Although not unexpected, it is noteworthy that the distribution of haplotypes differed between the two ethnic groups. For example, the frequency of the high risk haplotype CGGTT in non-Hispanic whites (0.161) was higher than that in Hispanic whites (0.108). Ordinal logistic regression analysis revealed that the homozygous status of minor alleles of all five SNPs were each associated with a greater than two-fold increase in spina bifida risk in the non-Hispanic white population. These increases were also seen in Hispanic whites, even though they did not reach statistical significance. There is currently no experimental data on the functional consequences of these variants. We performed functional prediction using an online program, FASTSNP, which predicts functional impact of SNPs according to current knowledge. FASTSNP provides a "risk score" for each SNP based on its putative biological function. Among the five SNPs we studied, rs4621, which is a synonymous change located in exon 1, has the highest "risk score". SNP rs652021 has a risk score of zero, which means that it has essentially no known functional impact. Strong linkage disequilibrium (LD) was observed among the five SNPs in both sub-populations (non-Hispanic white and Hispanic white) in our study. Therefore, the observed associations of these SNPs with spina bifida could be due to the strong LD to the "real" functional variation.
Our results, although based on multiple tests on small sample sizes and therefore lowered statistical power, represent a preliminary step in elucidating the association between CFL1 gene variations and spina bifida risk. The haplotype CGGTT, which is composed of minor alleles of all five SNPs and was commonly found in both major race/ethnic groups studied, conferred an increased risk of spina bifida. Combined with existing knowledge about NTDs, these polymorphisms could clarify the complex mechanism of genetic susceptibility underlying the risk for NTDs. Further studies of these variations using cellular experiments could potentially provide an in-depth view of the mechanism of our observations.
The authors are indebted to Drs. George Cunningham and Fred Lorey for making it possible to access newborn blood specimens. We also appreciate the technical support including DNA sample preparation from Ms. Consuelo Valdes. This work was supported by funds from NIH/NINDS RO1NS050249.
- Detrait ER, George TM, Etchevers HC, Gilbert JR, Vekemans M, Speer MC: Human neural tube defects: developmental biology, epidemiology, and genetics. Neurotoxicol Teratol. 2005, 27 (3): 515-524. 10.1016/j.ntt.2004.12.007.View ArticlePubMedPubMed CentralGoogle Scholar
- Blom HJ, Shaw GM, den Heijer M, Finnell RH: Neural tube defects and folate: case far from closed. Nat Rev Neurosci. 2006, 7 (9): 724-731. 10.1038/nrn1986.View ArticlePubMedPubMed CentralGoogle Scholar
- Berry RJ, Li Z, Erickson JD, Li S, Moore CA, Wang H, Mulinare J, Zhao P, Wong LY, Gindler J, Hong SX, Correa A: Prevention of neural-tube defects with folic acid in China. China-U.S. Collaborative Project for Neural Tube Defect Prevention. N Engl J Med. 1999, 341 (20): 1485-1490. 10.1056/NEJM199911113412001.View ArticlePubMedGoogle Scholar
- Czeizel AE, Dudas I: Prevention of the first occurrence of neural-tube defects by periconceptional vitamin supplementation. N Engl J Med. 1992, 327 (26): 1832-1835.View ArticlePubMedGoogle Scholar
- MRC: Prevention of neural tube defects: results of the Medical Research Council Vitamin Study. MRC Vitamin Study Research Group. Lancet. 1991, 338 (8760): 131-137. 10.1016/0140-6736(91)90133-A.View ArticleGoogle Scholar
- Barber R, Shalat S, Hendricks K, Joggerst B, Larsen R, Suarez L, Finnell R: Investigation of Folate Pathway Gene Polymorphisms and the Incidence of Neural Tube Defects in a Texas Hispanic Population. Molecular Genetics and Metabolism. 2000, 70 (1): 45-52. 10.1006/mgme.2000.2991.View ArticlePubMedGoogle Scholar
- Botto LD, Yang Q: 5,10-Methylenetetrahydrofolate reductase gene variants and congenital anomalies: a HuGE review. American Journal Of Epidemiology. 2000, 151 (9): 862-877.View ArticlePubMedGoogle Scholar
- Zhu H, Wicker NJ, Shaw GM, Lammer EJ, Hendricks K, Suarez L, Canfield M, Finnell RH: Homocysteine remethylation enzyme polymorphisms and increased risks for neural tube defects. Molecular Genetics and Metabolism. 2003, 78 (3): 216-221. 10.1016/S1096-7192(03)00008-8.View ArticlePubMedGoogle Scholar
- Shaw GM, Rozen R, Finnell RH, Wasserman CR, Lammer EJ: Maternal vitamin use, genetic variation of infant methylenetetrahydrofolate reductase, and risk for spina bifida. American Journal Of Epidemiology. 1998, 148 (1): 30-37.View ArticlePubMedGoogle Scholar
- Harris MJ, Juriloff DM: Mini-review: toward understanding mechanisms of genetic neural tube defects in mice. Teratology. 1999, 60 (5): 292-305. 10.1002/(SICI)1096-9926(199911)60:5<292::AID-TERA10>3.0.CO;2-6.View ArticlePubMedGoogle Scholar
- Gurniak CB, Perlas E, Witke W: The actin depolymerizing factor n-cofilin is essential for neural tube morphogenesis and neural crest cell migration. Dev Biol. 2005, 278 (1): 231-241. 10.1016/j.ydbio.2004.11.010.View ArticlePubMedGoogle Scholar
- Gillett GT, Fox MF, Rowe PS, Casimir CM, Povey S: Mapping of human non-muscle type cofilin (CFL1) to chromosome 11q13 and muscle-type cofilin (CFL2) to chromosome 14. Ann Hum Genet. 1996, 60 ( Pt 3): 201-211.View ArticleGoogle Scholar
- Croen LA, Shaw GM, Jensvold NG, Harris JA: Birth defects monitoring in California: a resource for epidemiological research. Paediatr Perinat Epidemiol. 1991, 5 (4): 423-427.View ArticlePubMedGoogle Scholar
- Schulman J, Edmonds LD, McClearn AB, Jensvold N, Shaw GM: Surveillance for and comparison of birth defect prevalences in two geographic areas--United States, 1983-88. MMWR CDC Surveill Summ. 1993, 42 (1): 1-7.PubMedGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, Higgins J, DeFelice M, Lochner A, Faggart M, Liu-Cordero SN, Rotimi C, Adeyemo A, Cooper R, Ward R, Lander ES, Daly MJ, Altshuler D: The structure of haplotype blocks in the human genome. Science. 2002, 296 (5576): 2225-2229. 10.1126/science.1069424.View ArticlePubMedGoogle Scholar
- Lowry R: VassarStats (http://faculty.vassar.edu/lowry/VassarStats.html).Google Scholar
- Yuan HY, Chiou JJ, Tseng WH, Liu CH, Liu CK, Lin YJ, Wang HH, Yao A, Chen YT, Hsu CN: FASTSNP: an always up-to-date and extendable service for SNP function analysis and prioritization. Nucleic Acids Res. 2006, 34 (Web Server issue): W635-41. 10.1093/nar/gkl236.View ArticlePubMedPubMed CentralGoogle Scholar
- Nishida E, Maekawa S, Sakai H: Cofilin, a protein in porcine brain that binds to actin filaments and inhibits their interactions with myosin and tropomyosin. Biochemistry. 1984, 23 (22): 5307-5313. 10.1021/bi00317a032.View ArticlePubMedGoogle Scholar
- Abe H, Ohshima S, Obinata T: A cofilin-like protein is involved in the regulation of actin assembly in developing skeletal muscle. J Biochem (Tokyo). 1989, 106 (4): 696-702.Google Scholar
- Moriyama K, Nishida E, Yonezawa N, Sakai H, Matsumoto S, Iida K, Yahara I: Destrin, a mammalian actin-depolymerizing protein, is closely related to cofilin. Cloning and expression of porcine brain destrin cDNA. J Biol Chem. 1990, 265 (10): 5768-5773.PubMedGoogle Scholar
- Loisel TP, Boujemaa R, Pantaloni D, Carlier MF: Reconstitution of actin-based motility of Listeria and Shigella using pure proteins. Nature. 1999, 401 (6753): 613-616. 10.1038/44183.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/8/12/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.