Evidence that TGFA influences risk to cleft lip with/without cleft palate through unconventional genetic mechanisms

This study examined the association between markers in transforming growth factor alpha (TGFA) and isolated, non-syndromic cleft lip with/without palate (CL/P) using a case–parent trio design, considering parent-of-origin effects. We also tested for gene–environmental interaction with common maternal exposures, and for gene–gene interaction using markers in TGFA and another recognized causal gene, IRF6. CL/P case–parent trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 17 single nucleotide polymorphisms (SNPs) in TGFA. The transmission disequilibrium test was used to test individual SNPs, and the parent-of-origin likelihood ratio test (PO-LRT) was used to assess parent-of-origin effects. We also screened for possible gene–environment interaction using PBAT, and tested for gene–gene interaction using conditional logistic regression models. When all trios were combined, four SNPs showed significant excess maternal transmission, two of which gave significant PO-LRT values [rs3821261: P = 0.004 and OR(imprinting) = 4.17; and rs3771475: P = 0.027 and OR(imprinting) = 2.44]. Haplotype analysis of these two SNPS also supported excess maternal transmission. We saw intriguing but suggestive evidence of G × E interaction for several SNPs in TGFA when either individual SNPs or haplotypes of adjacent SNPs were considered. Thus, TGFA appears to influence risk of CL/P through unconventional means with an apparent parent-of-origin effect (excess maternal transmission) and possible interaction with maternal exposures.     

Evidence, from family studies, for linkage disequilibrium between TGFA and a gene for nonsyndromic cleft lip with or without cleft palate.

The inheritance of alleles of the transforming growth factor alpha (TGFA) locus has been studied in families affected with cleft lip with or without cleft palate (CL/P), by using the transmission/disequilibrium test described by Spielman and colleagues. Only heterozygous parents with an affected child can be included in this test, but within such families a significantly greater frequency of C2 alleles were transmitted to affected children than would be expected by chance. There was no evidence that the total number of C2 alleles transmitted to affected and unaffected children differed significantly from random segregation. These data provide evidence from within families that a gene for susceptibility to CL/P is in significant linkage disequilibrium with the C2 allele of the TGFA locus.     

Few associations of candidate genes with nonsyndromic orofacial clefts in the population of Lithuania

Nonsyndromic orofacial clefting (NS-OFC) is a common complex multifactorial trait with a considerable genetic component and a number of candidate genes suggested by various approaches. Twenty biallelic and microsatellite DNA markers in the strong candidate loci TGFA, TGFB3, GABRB3, RARA, and BCL3 were analysed for allelic association with the NS-OFC phenotype in 112 nuclear families (proband + both parents) from Lithuania by using the transmission disequilibrium test (TDT). Associations were found between the TGFA gene marker rs2166975 and nonsyndromic cleft palate only (CPO) phenotype (p = 0.045, df 1) as well as between the D2S292 marker and the cleft lip with or without cleft palate (CL/CP) phenotype in allele-wise TDT (P = 0.005, df 9) and genotype-wise TDT (P = 0.021, df 24). A weak association (P = 0.085, df 3) of the BCL3 marker (BCL3 gene) with the risk of CPO was also found. Thus our initial results support the contribution of allelic variation in the TGFA locus to the aetiology of CL/CP in the population of Lithuania but they do not point to TGFA as a major causal gene. Different roles of the TGFA and BCL3 genes in the susceptibility to NS-OFC phenotypes are suggested.     

Confirming genes influencing risk to cleft lip with/without cleft palate in a case–parent trio study

A collection of 1,108 case–parent trios ascertained through an isolated, nonsyndromic cleft lip with or without cleft palate (CL/P) was used to replicate the findings from a genome-wide association study (GWAS) conducted by Beaty et al. (Nat Genet 42:525–529, 2010), where four different genes/regions were identified as influencing risk to CL/P. Tagging SNPs for 33 different genes were genotyped (1,269 SNPs). All four of the genes originally identified as showing genome-wide significance (IRF6, ABCA4 and MAF, plus the 8q24 region) were confirmed in this independent sample of trios (who were primarily of European and Southeast Asian ancestry). In addition, eight genes classified as ‘second tier’ hits in the original study (PAX7, THADA, COL8A1/FILIP1L, DCAF4L2, GADD45G, NTN1, RBFOX3 and FOXE1) showed evidence of linkage and association in this replication sample. Meta-analysis between the original GWAS trios and these replication trios showed PAX7, COL8A1/FILIP1L and NTN1 achieved genome-wide significance. Tests for gene–environment interaction between these 33 genes and maternal smoking found evidence for interaction with two additional genes: GRID2 and ELAVL2 among European mothers (who had a higher rate of smoking than Asian mothers). Formal tests for gene–gene interaction (epistasis) failed to show evidence of statistical interaction in any simple fashion. This study confirms that many different genes influence risk to CL/P.     

Evidence of gene–environment interaction for the IRF6 gene and maternal multivitamin supplementation in controlling the risk of cleft lip with/without cleft palate

Although multiple genes have been identified as genetic risk factors for isolated, non-syndromic cleft lip with/without cleft palate (CL/P), a complex and heterogeneous birth defect, interferon regulatory factor 6 gene (IRF6) is one of the best documented genetic risk factors. In this study, we tested for association between markers in IRF6 and CL/P in 326 Chinese case–parent trios, considering gene–environment interaction for two common maternal exposures, and parent-of-origin effects. CL/P case–parent trios from three sites in mainland China and Taiwan were genotyped for 22 single nucleotide polymorphisms (SNPs) in IRF6. The transmission disequilibrium test was used to test for marginal effects of individual SNPs. We used PBAT to screen the SNPs and haplotypes for gene–environment (G × E) interaction and conditional logistic regression models to quantify effect sizes for SNP–environment interaction. After Bonferroni correction, 14 SNPs showed statistically significant association with CL/P. Evidence of G × E interaction was found for both maternal exposures, multivitamin supplementation and environmental tobacco smoke (ETS). Two SNPs showed evidence of interaction with multivitamin supplementation in conditional logistic regression models (rs2076153 nominal P = 0.019, rs17015218 nominal P = 0.012). In addition, rs1044516 yielded evidence for interaction with maternal ETS (nominal P = 0.041). Haplotype analysis using PBAT also suggested interaction between SNPs in IRF6 and both multivitamin supplementation and ETS. However, no evidence for maternal genotypic effects or significant parent-of-origin effects was seen in these data. These results suggest IRF6 gene may influence risk of CL/P through interaction with multivitamin supplementation and ETS in the Chinese population.     

Cleft lip with or without cleft palate: associations with transforming growth factor alpha and retinoic acid receptor loci.

The first association study of cleft lip with or without cleft palate (CL/P), with candidate genes, found an association with the transforming growth-factor alpha (TGFA) locus. This finding has since been replicated, in whole or in part, in three independent studies. Here we extend our original analysis of the TGFA TaqI RFLP to two other TGFA RFLPs and seven other RFLPs at five candidate genes in 117 nonsyndromic cases of CL/P and 113 controls. The other candidate genes were the retinoic acid receptor (RARA), the bcl-2 oncogene, and the homeobox genes 2F, 2G, and EN2. Significant associations with the TGFA TaqI and BamHI RFLPs were confirmed, although associations of clefting with previously reported haplotypes did not reach significance. Of particular interest, in view of the known teratogenic role of retinoic acid, was a significant association with the RARA PstI RFLP (P = .016; not corrected for multiple testing). The effect on risk of the A2 allele appears to be additive, and although the A2A2 homozygote only has an odds ratio of about 2 and recurrence risk to first-degree relatives (lambda 1) of 1.06, because it is so common it may account for as much as a third of the attributable risk of clefting. There is no evidence of interaction between the TGFA and RARA polymorphisms on risk, and jointly they appear to account for almost half the attributable risk of clefting.     

Strong evidence that KIAA0319 on chromosome 6p is a susceptibility gene for developmental dyslexia

Linkage between developmental dyslexia (DD) and chromosome 6p has been replicated in a number of independent samples. Recent attempts to identify the gene responsible for the linkage have produced inconsistent evidence for association of DD with a number of genes in a 575-kb region of chromosome 6p22.2, including VMP, DCDC2, KIAA0319, TTRAP, and THEM2. We aimed to identify the specific gene or genes involved by performing a systematic, high-density (approximately 2-3-kb intervals) linkage disequilibrium screen of these genes in an independent sample, incorporating family-based and case-control designs in which dyslexia was defined as an extreme representation of reading disability. Using DNA pooling, we first observed evidence for association with 17 single-nucleotide polymorphisms (SNPs), 13 of which were located in the KIAA0319 gene (P<.01-.003). After redundant SNPs were excluded, 10 SNPs were individually genotyped in 223 subjects with DD and 273 controls. Those SNPs that were significant at P相似文献   

Evidence for an association between nonsyndromic cleft lip with or without cleft palate and a gene located on the long arm of chromosome 4.

Recent studies suggest that the familial aggregation of nonsyndromic cleft lip with or without cleft palate (CL +/- P) is likely to be attributable to the effects of several susceptibility loci, acting in a multiplicative fashion. Two potential CL +/- P susceptibility loci (CSL), transforming growth factor alpha (TGFA) and retinoic acid receptor (RARA), have been identified through association studies. In addition, recent evidence of linkage between CL +/- P and two markers (D4S175 and D4S192) in the region 4q25-4q31.3 raised the possibility that a CSL, with a larger effect than either TGFA or RARA, may reside within this region of the human genome. The present analyses were undertaken to determine whether D4S175 or D4S192 is significantly associated with CL +/- P in a sample of unrelated patients that have previously provided evidence of associations between CL +/- P and both TGFA and RARA. The results of these analyses provide further, tentative, evidence for the presence of a CSL locus on the long arm of chromosome 4 and help to refine the location of this locus in the region of D4S175 and D4S192.     

Nonsyndromic cleft lip with or without cleft palate: evidence of linkage to BCL3 in 17 multigenerational families.

Nonsyndromic cleft lip with or without cleft palate (CL/P) is a common craniofacial developmental defect. Recent segregation analyses have suggested that major genes play a role in the etiology of CL/P. Linkage to 22 candidate genes was tested in 11 multigenerational families with CL/P, and 21 of these candidates were excluded. APOC2, 19q13.1, which is linked to the proto-oncogene BCL3, gave suggestive evidence for linkage to CL/P. The study was expanded to include a total of 39 multigenerational CL/P families. Linkage was tested in all families, using an anonymous marker, D19S178, and intragenic markers in BCL3 and APOC2. Linkage was tested under two models, autosomal dominant with reduced penetrance and affecteds only. Homogeneity testing on the two-point data gave evidence of heterogeneity at APOC2 under the affecteds-only model. Both models showed evidence of heterogeneity, with 43% of families linked at zero recombination to BCL3 when marker data from BCL3 and APOC2 were included. A maximum multipoint LOD score of 7.00 at BCL3 was found among the 17 families that had posterior probabilities > = 50% in favor of linkage. The transmission disequilibrium test provided additional evidence for linkage with the 3 allele of BCL3 more often transmitted to affected children. These results suggest that BCL3, or a nearby gene, plays a role in the etiology of CL/P in some families.     

Association of MSX1 and TGFB3 with nonsyndromic clefting in humans.

Nonsyndromic cleft lip with or without cleft palate (CL/P) and nonsyndromic cleft palate only (CPO) are common congenital anomalies with significant medical, psychological, social, and economic ramifications. Both CL/P and CPO are examples of complex genetic traits. There exists sufficient evidence to hypothesize that disease loci for CL/P and CPO can be identified by a candidate-gene linkage-disequilibrium (LD) strategy. Candidate genes for clefting, including TGFA, BCL3, DLX2, MSX1, and TGFB3, were screened for LD with either CL/P or CPO in a predominantly Caucasian population, with both case-control- and nuclear-family-based approaches. Previously reported LD for TGFA with both CL/P and CPO could not be confirmed, except in CL/P patients with a positive family history. Also, in contrast to previous studies, no LD was found between BCL3 and either CL/P or CPO. Significant LD was found between CL/P and both MSX1 and TGFB3 and between CPO and MSX1, suggesting that these genes are involved in the pathogenesis of clefting. In addition, a mutation search in the genes DLX2, MSX1, and TGFB3 was performed in 69 CPO patients and in a subset of the CL/P patients. No common mutations were found in the coding regions of these genes; however, several rare variants of MSX1 and TGFB3 were found that may alter the latters' normal function. These results form the basis for future research, including (a) mutation searches in the MSX1 and TGFB3 genes in Caucasian CL/P patients and (b) extension of the search for MSX1 mutations in CPO patients to the noncoding regions.     

Meta-analysis of 13 genome scans reveals multiple cleft lip/palate genes with novel loci on 9q21 and 2q32-35

Isolated or nonsyndromic cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex etiology. A 10-cM genome scan of 388 extended multiplex families with CL/P from seven diverse populations (2,551 genotyped individuals) revealed CL/P genes in six chromosomal regions, including a novel region at 9q21 (heterogeneity LOD score [HLOD]=6.6). In addition, meta-analyses with the addition of results from 186 more families (six populations; 1,033 genotyped individuals) showed genomewide significance for 10 more regions, including another novel region at 2q32-35 (P=.0004). These are the first genomewide significant linkage results ever reported for CL/P, and they represent an unprecedented demonstration of the power of linkage analysis to detect multiple genes simultaneously for a complex disorder.     

Linkage disequilibrium between GABRB3 gene and nonsyndromic familial cleft lip with or without cleft palate

The malformation of nonsyndromic cleft lip with or without cleft palate (CL/P) is a common congenital disease that affects approximately 1/1000 newborns in Caucasian populations. Genetic studies indicate that CL/P has the characteristics of a complex genetic trait. Linkage analysis and mouse-model knockout studies have suggested several candidate genes mapping in different chromosome regions for CL/P malformation. On these grounds, we have investigated, by linkage disequilibrium (LD) and parametric and nonparametric linkage analyses, five different candidate genes, including those for the beta3 subunit of the gamma-aminobutyric acid receptor (GABRB3), glutamic acid decarboxylase 1 (GAD1), retinoic acid receptor alpha (RARA), transforming growth factor beta3 (TGFB3), and msh ( Drosophila) homeobox homolog 1 (MSX1). Interestingly, a significant LD between GABRB3 and CL/P was obtained ( P-value=0.008 in the allele-wise analysis for multiallelic markers), suggesting that the GABRB3 gene is involved in this congenital disease. This new finding in humans is in agreement with previously reported data obtained with the murine model. Indeed, mouse studies indicate a role for gamma-aminobutyric acid (GABA) and its receptor in normal palate development. Exclusion of the GAD1 gene, which encodes the GABA-producing enzyme, in CL/P pathogenesis was obtained in our study. Moreover, we were unable to confirm the involvement of the MSX1 gene in nonsyndromic CL/P. Modest evidence of LD between marker alleles and CL/P was found at the RARA and TGFB3 loci suggesting a minor role for these genes in our family set of nonsyndromic CL/P.     

Pleiotropy method reveals genetic overlap between orofacial clefts at multiple novel loci from GWAS of multi-ethnic trios

Based on epidemiologic and embryologic patterns, nonsyndromic orofacial clefts– the most common craniofacial birth defects in humans– are commonly categorized into cleft lip with or without cleft palate (CL/P) and cleft palate alone (CP), which are traditionally considered to be etiologically distinct. However, some evidence of shared genetic risk in IRF6, GRHL3 and ARHGAP29 regions exists; only FOXE1 has been recognized as significantly associated with both CL/P and CP in genome-wide association studies (GWAS). We used a new statistical approach, PLACO (pleiotropic analysis under composite null), on a combined multi-ethnic GWAS of 2,771 CL/P and 611 CP case-parent trios. At the genome-wide significance threshold of 5 × 10⁻⁸, PLACO identified 1 locus in 1q32.2 (IRF6) that appears to increase risk for one OFC subgroup but decrease risk for the other. At a suggestive significance threshold of 10⁻⁶, we found 5 more loci with compelling candidate genes having opposite effects on CL/P and CP: 1p36.13 (PAX7), 3q29 (DLG1), 4p13 (LIMCH1), 4q21.1 (SHROOM3) and 17q22 (NOG). Additionally, we replicated the recognized shared locus 9q22.33 (FOXE1), and identified 2 loci in 19p13.12 (RAB8A) and 20q12 (MAFB) that appear to influence risk of both CL/P and CP in the same direction. We found locus-specific effects may vary by racial/ethnic group at these regions of genetic overlap, and failed to find evidence of sex-specific differences. We confirmed shared etiology of the two OFC subtypes comprising CL/P, and additionally found suggestive evidence of differences in their pathogenesis at 2 loci of genetic overlap. Our novel findings include 6 new loci of genetic overlap between CL/P and CP; 3 new loci between pairwise OFC subtypes; and 4 loci not previously implicated in OFCs. Our in-silico validation showed PLACO is robust to subtype-specific effects, and can achieve massive power gains over existing approaches for identifying genetic overlap between disease subtypes. In summary, we found suggestive evidence for new genetic regions and confirmed some recognized OFC genes either exerting shared risk or with opposite effects on risk to OFC subtypes.     

Aryl hydrocarbon receptor nuclear translocator 2 (ARNT2): structure,gene mapping,polymorphisms, and candidate evaluation for human orofacial clefts

BACKGROUND: Nonsyndromic orofacial clefts have an estimated incidence of 1/1000 live births. Population genetic and embryologic studies suggest that cleft palate only (CPO) may be a distinct clinical entity from cleft lip with or without cleft palate (CL/P). Both CPO and CL/P are thought to be multifactorial in etiology, with evidence indicating that genetic, environmental, and developmental determinants may all play a role. The ARNT2 gene localizes to a conserved linkage group on mouse chromosome 7 that is syntenic with human chromosome 15q23-25. This chromosomal region was previously identified as a teratogen-induced clefting susceptibility locus in a genome-wide scan of AXB and BXA recombinant inbred mice. Arnt2 is expressed in the first branchial arch in mice. The teratogen 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) acts through the aryl hydrocarbon receptor (Ahr) pathway to produce dose-dependent CPO and thymic wasting in mice exposed in utero. Arnt2 and Ahr proteins dimerize in vitro. TCDD exposure is also associated with orofacial clefting in children of parents involved in agricultural work. METHODS: To determine whether ARNT2 influences human craniofacial development, we identified the human ARNT2 gene and conducted genomic structural analysis. Mutational screening was performed in infants with nonsyndromic CPO or CL/P who were identified by the Iowa Birth Defects Registry. RESULTS: A common amino acid polymorphism was detected but, no obvious disease-causing mutations were detected by SSCP analysis. The microsatellite marker, GATA89D04 (D15S823) was identified within intron 11 of the human ARNT2 gene, and linkage disequilibrium of nonsyndromic CPO and CL/P parent-infant trios was conducted. CONCLUSIONS: No association was demonstrated with CPO (n = 45) and CL/P (n = 37). Teratology 66:85-90, 2002.     

Joint Testing of Genotypic and Gene-Environment Interaction Identified Novel Association for BMP4 with Non-Syndromic CL/P in an Asian Population Using Data from an International Cleft Consortium

Background

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common disorder with complex etiology. The Bone Morphogenetic Protein 4 gene (BMP4) has been considered a prime candidate gene with evidence accumulated from animal experimental studies, human linkage studies, as well as candidate gene association studies. The aim of the current study is to test for linkage and association between BMP4 and NSCL/P that could be missed in genome-wide association studies (GWAS) when genotypic (G) main effects alone were considered.

Methodology/Principal Findings

We performed the analysis considering G and interactions with multiple maternal environmental exposures using additive conditional logistic regression models in 895 Asian and 681 European complete NSCL/P trios. Single nucleotide polymorphisms (SNPs) that passed the quality control criteria among 122 genotyped and 25 imputed single nucleotide variants in and around the gene were used in analysis. Selected maternal environmental exposures during 3 months prior to and through the first trimester of pregnancy included any personal tobacco smoking, any environmental tobacco smoke in home, work place or any nearby places, any alcohol consumption and any use of multivitamin supplements. A novel significant association held for rs7156227 among Asian NSCL/P and non-syndromic cleft lip and palate (NSCLP) trios after Bonferroni correction which was not seen when G main effects alone were considered in either allelic or genotypic transmission disequilibrium tests. Odds ratios for carrying one copy of the minor allele without maternal exposure to any of the four environmental exposures were 0.58 (95%CI = 0.44, 0.75) and 0.54 (95%CI = 0.40, 0.73) for Asian NSCL/P and NSCLP trios, respectively. The Bonferroni P values corrected for the total number of 117 tested SNPs were 0.0051 (asymptotic P = 4.39*10⁻⁵) and 0.0065 (asymptotic P = 5.54*10⁻⁵), accordingly. In European trios, no significant association was seen for any SNPs after Bonferroni corrections for the total number of 120 tested SNPs.

Conclusions/Significance

Our findings add evidence from GWAS to support the role of BMP4 in susceptibility to NSCL/P originally identified in linkage and candidate gene association studies.     

Medical sequencing of candidate genes for nonsyndromic cleft lip and palate

Nonsyndromic or isolated cleft lip with or without cleft palate (CL/P) occurs in wide geographic distribution with an average birth prevalence of 1/700. We used direct sequencing as an approach to study candidate genes for CL/P. We report here the results of sequencing on 20 candidate genes for clefts in 184 cases with CL/P selected with an emphasis on severity and positive family history. Genes were selected based on expression patterns, animal models, and/or role in known human clefting syndromes. For seven genes with identified coding mutations that are potentially etiologic, we performed linkage disequilibrium studies as well in 501 family triads (affected child/mother/father). The recently reported MSX1 P147Q mutation was also studied in an additional 1,098 cleft cases. Selected missense mutations were screened in 1,064 controls from unrelated individuals on the Centre d'Étude du Polymorphisme Humain (CEPH) diversity cell line panel. Our aggregate data suggest that point mutations in these candidate genes are likely to contribute to 6% of isolated clefts, particularly those with more severe phenotypes (bilateral cleft of the lip with cleft palate). Additional cases, possibly due to microdeletions or isodisomy, were also detected and may contribute to clefts as well. Sequence analysis alone suggests that point mutations in FOXE1, GLI2, JAG2, LHX8, MSX1, MSX2, SATB2, SKI, SPRY2, and TBX10 may be rare causes of isolated cleft lip with or without cleft palate, and the linkage disequilibrium data support a larger, as yet unspecified, role for variants in or near MSX2, JAG2, and SKI. This study also illustrates the need to test large numbers of controls to distinguish rare polymorphic variants and prioritize functional studies for rare point mutations.     

Association of genetic variation of the transforming growth factor-alpha gene with cleft lip and palate.

Complex segregation analysis of pedigrees having nonsyndromic cleft lip with or without cleft palate (CL/P) (Chung et al. 1986; Marazita et al. 1986) has shown that a major-locus model best explains the observed recurrence of CL/P in Caucasian families. To identify this major gene, we compared the frequencies of 12 RFLPs at five loci-epidermal growth factor, transforming growth factor-alpha, epidermal growth factor receptor, glucocorticoid receptor, and estrogen receptor-in both a group of 80 subjects with nonsyndromic CL/P and 102 controls. These candidate genes were selected because studies in rodents had suggested their possible involvement in palatogenesis. A significant association was observed between two RFLPs at the transforming-growth-factor-alpha (TGFA) locus and the occurrence of clefting (P = .0047 and P = .0052). This suggests that either the TGFA gene itself or DNA sequences in an adjacent region contribute to the development of a portion of cases of CL/P in humans and provides an opportunity to begin to examine the molecular events underlying lip and palate formation.     

BMP4 was associated with NSCL/P in an Asian population

Background

The Bone Morphogenetic Protein 4 gene (BMP4) is located in chromosome 14q22-q23 which has shown evidence of linkage for isolated nonsyndromic cleft lip with or without cleft palate (NSCL/P) in a genome wide linkage analysis of human multiplex families. BMP4 has been shown to play crucial roles in lip and palatal development in animal models. Several candidate gene association analyses also supported its potential risk for NSCL/P, however, results across these association studies have been inconsistent. The aim of the current study was to test for possible association between markers in and around the BMP4 gene and NSCL/P in Asian and Maryland trios.

Methodology/Principal Findings

Family Based Association Test was used to test for deviation from Mendelian assortment for 12 SNPs in and around BMP4. Nominal significant evidence of linkage and association was seen for three SNPs (rs10130587, rs2738265 and rs2761887) in 221 Asian trios and for one SNP (rs762642) in 76 Maryland trios. Statistical significance still held for rs10130587 after Bonferroni correction (corrected p = 0.019) among the Asian group. Estimated odds ratio for carrying the apparent high risk allele at this SNP was 1.61 (95%CI = 1.20, 2.18).

Conclusions

Our results provided further evidence of association between BMP4 and NSCL/P.