Association of a BMP9 haplotype with ossification of the posterior longitudinal ligament (OPLL) in a Chinese population

Direct or ex vivo BMP9 adenoviral gene therapy can induce massive bone formation at the injection sites and clearly promote spinal fusion. A comprehensive analysis of the osteogenic activity indicated that BMP9 was one of the most potent inducers of osteogenic differentiation both in vitro and in vivo among 14 types of human BMPs. However, genetic variations and whether they correlated with OPLL were not considered. We have sequenced the complete BMP9 gene in 450 patients with OPLL and in 550 matched controls. Analyses were performed on single markers and haplotypes. Single marker tests identified 6 SNPs, among which the minor alleles of rs7923671 (T>C; P=0.0026; OR: 1.33, CI: 1.10-1.60), rs75024165 (C>T, Thr304Met; P<0.001; OR: 1.76, CI: 1.47-2.12) and rs34379100 (A>C; P<0.001; OR: 1.52, CI: 1.27-1.82) were associated with OPLL. Logistic regression analysis showed that the additive model of rs75024165 (TT vs. CT vs. CC; P<0.001; OR: 1.74) and rs34379100 (CC vs. AC vs. AA; P=0.003; OR: 1.95) retained statistical significance when adjusted for clinical and demographic characteristics. Linkage disequilibrium (LD) analysis identified one 3 kb block of intense LD in BMP9 and one specific haplotype, CTCA (P<0.001; OR: 2.37), that contained the OPLL-associated risk alleles and was a risk factor for OPLL. This haplotype is associated with increased severity of OPLL, as shown by the distribution of ossified vertebrae in patients with OPLL (P=0.001). In summary, in the Chinese population studied, SNPs in the BMP9 gene appear to contribute to the risk of OPLL in association with certain clinical and demographic characteristics. The severity of OPLL seems to be mediated predominantly by genetic variations in a 3kb BMP9 locus with the specific haplotype CTCA.     

Multiple common susceptibility variants near BMP pathway loci GREM1, BMP4, and BMP2 explain part of the missing heritability of colorectal cancer

Genome-wide association studies (GWAS) have identified 14 tagging single nucleotide polymorphisms (tagSNPs) that are associated with the risk of colorectal cancer (CRC), and several of these tagSNPs are near bone morphogenetic protein (BMP) pathway loci. The penalty of multiple testing implicit in GWAS increases the attraction of complementary approaches for disease gene discovery, including candidate gene- or pathway-based analyses. The strongest candidate loci for additional predisposition SNPs are arguably those already known both to have functional relevance and to be involved in disease risk. To investigate this proposition, we searched for novel CRC susceptibility variants close to the BMP pathway genes GREM1 (15q13.3), BMP4 (14q22.2), and BMP2 (20p12.3) using sample sets totalling 24,910 CRC cases and 26,275 controls. We identified new, independent CRC predisposition SNPs close to BMP4 (rs1957636, P = 3.93×10(-10)) and BMP2 (rs4813802, P = 4.65×10(-11)). Near GREM1, we found using fine-mapping that the previously-identified association between tagSNP rs4779584 and CRC actually resulted from two independent signals represented by rs16969681 (P = 5.33×10(-8)) and rs11632715 (P = 2.30×10(-10)). As low-penetrance predisposition variants become harder to identify-owing to small effect sizes and/or low risk allele frequencies-approaches based on informed candidate gene selection may become increasingly attractive. Our data emphasise that genetic fine-mapping studies can deconvolute associations that have arisen owing to independent correlation of a tagSNP with more than one functional SNP, thus explaining some of the apparently missing heritability of common diseases.     

BMP2 variants in the risk of ankylosing spondylitis

The purpose of the study was to explore the genetic effects of bone morphogenetic protein (BMP2) polymorphisms on the susceptibility to ankylosing spondylitis (AS) in Chinese Han population. The case-control study included 120 AS cases and 110 healthy controls. Hardy-Weinberg equilibrium test was performed in control group. BMP2 rs235768 and rs3178250 polymorphisms were analyzed by polymerase chain reaction and direct sequencing. Additionally, the χ² test was used to estimate association strength between BMP2 genetic polymorphisms and AS susceptibility, and the results were assessed via odds ratio (OR) with the corresponding 95% confidence interval (95%CI). Results adjustment was performed using logistic regression analysis. AA, AT, TT genotype and A, T allele frequencies of BMP2 rs235768 polymorphism presented no significant differences between case and control groups (P > .05 for all). TC genotype of rs3178250 polymorphism showed significantly higher in case group than that in control group (P = .048). After adjusting, TC genotype was a risk factor for AS (OR = 2.095; 95%CI = 1.086-4.038; P = .027). BMP2 rs3178250 polymorphism may increase individual susceptibility to AS in Chinese Han population.     

The association of polymorphisms in BMP2/MYO1H and skeletal Class II div.1 maxillary and mandibular dimensions. A preliminary ‘report

IntroductionThe genetic impact directly or indirectly predefines maxillofacial dimensions, potentially leading to an inappropriate relationship of the jaws and subsequently skeletal malocclusion. Previous studies focused mainly on genetic polymorphisms and class III malocclusion. This study was set out to investigate the association between genetic polymorphisms in two genes BMP2 (rs235768) and MYO1H (rs11066446) with Class II division 1 malocclusion, skeletal variation in vertical plane, and maxillary and mandibular jaws length.Subjects and methodsSixty patients classified as Skeletal Class I (n = 30) and Class II division 1 (n = 30) were recruited. DNA was extracted from saliva and analyzed by Sanger sequencing. Lateral cephalometric radiographs were measured for the anterio-posterior relationship of maxillary and mandibular arch using digital tracing. Hardy-Weinberg equilibrium analysis of genotype frequencies was performed using Chi-square test to compare genotype distribution among groups and multiple logistic regression analysis adjusted by gender was also performed.ResultsThe rs235768 polymorphism in BMP2 was associated with hypodivergent face, increased maxillary length, and decreased mandibular length. Meanwhile, the rs11066446 polymorphism in MYO1H was associated with decreased mandibular length. New polymorphism was identified in MYO1H (rs10850090) in association with decreased mandibular length.ConclusionA potential association between polymorpisms in BMP2 rs235768 and MOY1H rs11066446 and rs10850090 and Class II division 1 skeletal malocclusion related phenotypes exists, however, the degree of it has to be further investigated and yet to be discovered.     

Association between BH3 interacting domain death agonist (BID) gene polymorphism and ossification of the posterior longitudinal ligament in Korean population

The purpose of this study was to investigate single nucleotide polymorphisms (SNPs) of the BH3 interacting domain death agonist (BID) gene as a risk factor in Korean patients with ossification of the posterior longitudinal ligament (OPLL). To investigate the genetic association, two coding SNPs (rs8190315, Ser10Gly; rs2072392, Asp60Asp) of BID were genotyped in 157 OPLL patients and 209 control subjects. SNPStats, SNPAnalyzer Pro, Helixtree, and Haploview 4.2 programs were used for association analysis. Multiple logistic regression models (codominant, dominant, and recessive) were calculated for the odds ratios (ORs), 95 % confidence intervals (CIs), and corresponding P values. For multiple testing, Bonferroni correction was performed. After Bonferroni correction, genotype analysis of both rs8190315 and rs2072392 showed association between the OPLL group and the control group in the codominant model (P = 0.042, OR 1.86, 95 % CI 1.10–3.15). A complete linkage disequilibrium block was estimated between the two SNPs. Both of the G allele of rs8190315 and C allele of rs2072392 were strongly associated with an increased risk in the development of OPLL (P = 0.0052, OR 2.66, 95 % CI 1.51–4.68). These results suggest that BID is associated with OPLL, and both the G allele of a missense SNP (rs8190315, Ser10Gly) and C allele of a synonymous SNP (rs2072392, Asp60Asp) are risk factors for the development of OPLL in Korean population.     

RUNX2对BMP9诱导的间充质干细胞C3H10T1/2成骨分化的影响

目的:研究和确认RUNX2在骨形态发生蛋白9(BMP9)诱导的间充质干细胞C3H10T1/2成骨分化中的作用。方法:通过Western blot、RT-PCR、荧光素酶活性分析检测BMP9对RUNX2表达的影响;分别在过表达RUNX2和RNA干扰抑制RUNX2表达的情况下,利用碱性磷酸酶(ALP)活性测定和染色、钙盐沉积实验,免疫细胞化学和裸鼠皮下异位成骨实验分析RUNX2对于BMP9诱导的间充质干细胞成骨分化的影响。结果:BMP9可以促进RUNX2的表达;RUNX2体外可促进BMP9诱导的C3H10T1/2的ALP活性和钙盐沉积,却抑制了OCN表达,RUNX2还可促进BMP9诱导的裸鼠皮下异位成骨;而在降低RUNX2表达后,BMP9诱导的C3H10T1/2细胞的ALP活性、钙盐沉积、OCN表达和裸鼠皮下异位成骨均受到抑制。结论:RUNX2可以促进BMP9诱导的间充质干细胞C3H10T1/2细胞成骨分化。     

BMP2具有诱导C3H10T1/2细胞成脂肪分化的能力

探讨骨形态发生蛋白2（BMP2）诱导鼠胚胎间充质干细胞C3H10T1/2成脂肪分化能力,为临床脂肪代谢疾病的治疗提供理论基础．培养多潜能的间充质干细胞C3H10T1/2,用20 μg/ml BMP2对其诱导一定时间后,RT-PCR检测是否存在BMP信号通路中关键分子BMP受体BMPR I, BMPR Ⅱ及Smad 1/5/8的表达.Western印迹检测Smad 蛋白及MAPK 信号通路中p38磷酸化水平变化,QRT PCR检测成脂肪标志基因aP2以及成脂肪相关转录因子PPARγ,C/EBPα,C/EBPβ表达水平,同时用油红O染色,观测C3H10T1/2细胞成脂肪分化情况．经BMP2诱导后,C3H10T1/2细胞成脂肪分化标志(油红O染色)显著增加,Smad 蛋白及p38磷酸化水平有所上升,同时成脂肪标志基因aP2以及成脂肪相关转录因子PPARγ,C/EBPα,C/EBPβ表达水平各有一定程度提高．BMP2具有诱导C3H10T1/2细胞成脂肪分化能力,其成脂肪分化呈现对BMP2作用的时间依赖性．     

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.     

Effect of ANXA2 gene single nucleotide polymorphism (SNP) on the development of osteonecrosis in Indian sickle cell patient: a PCR-RFLP approach

Osteonecrosis is a serious complication in sickle cell patients. The common sites of the necrosis are femoral head, head of the humerus and acetabulam. Annexin A2 (ANXA2) protein mainly functions in bone formation and bone resorption. Alteration of ANXA2 gene may affect the manifestations of osteonecrosis in the patients. PCR-RFLP is a common applicable technique for the detection of known mutation/polymorphisms. Here we are presenting application of the PCR-RFLP technique for determination of the ANXA2 gene single nucleotide polymorphism frequency and their clinical association among Indian sickle cell patients. Five known SNPs of ANXA2 gene (rs7170178, rs73435133, rs73418020, rs72746635 and rs73418025) were determined using the HpyCH4V, DdeI, HpyCH4III and Sau 961 restriction enzyme respectively. Restriction enzyme DdeI was common for rs73435133 and rs72746635 SNP. Only the rs7170178 SNP was detected among patient and control and the other four SNPs were absent in the studied groups. The frequency of ANXA2 gene rs7170178 SNP (A/G, G/G) was comparatively higher in sickle cell patients than controls and it was clinically associated with sickle cell osteonecrosis. The P value of heterozygotes (A/G) and homozygotes (G/G) genotypes were <0.001 and 0.001 respectively, which were highly significant. This study established the application of PCR-RFLP in detection of ANXA2 SNPs in sickle cell patients.     

BMP9经JNK激酶途径调控间充质干细胞C3H10T1/2成骨分化

为了证实JNK激酶在骨形态发生蛋白9(bone morphogenetic proteins 9,BMP9) 诱导间充质干细胞C3H10T1/2成骨分化中的作用,利用重组腺病毒将BMP9导入间充质干细胞C3H10T1/2. 通过碱性磷酸酶(ALP)活性测定、钙盐沉积实验、荧光素酶报告基因检测、Western印迹和组织化学染色等方法,检测BMP9是否可经JNK激酶途径调控间充质干细胞C3H10T1/2向成骨分化.动物实验验证在RNA沉默JNK蛋白激酶后,对BMP9诱导间充质干细胞C3H10T1/2向成骨分化的影响．结果发现,BMP9可以增强JNK 激酶的磷酸化;利用JNK抑制剂SP600125抑制JNK激酶活性后,BMP9诱导的间充质干细胞C3H10T1/2的早期成骨指标ALP活性和晚期指标钙盐沉积均受到抑制,而且经典SMAD信号的活化也相应受到抑制;RNA干扰沉默JNK基因表达后,同样也可抑制BMP9 诱导的C3H10T1/2细胞的ALP活性和裸鼠皮下异位成骨．因此表明,BMP9可活化JNK激酶途径从而诱导间充质干细胞C3H10T1/2向成骨分化．     

MiR-21协同BMP9促进间充质干细胞C3H10T1/2成骨分化

目的:探讨miR-21与BMP9之间的关系,明确miR-21在BMP9诱导间充质干细胞成骨分化中的作用。方法:(1)Ad-BMP9感染C3H10T1/2细胞,Real-time-PCR检测miR-21表达。RT-PCR检测ALP的表达。(2)MiR-21转染C3H10T1/2细胞,Real-time-PCR检测miR-21和BMP9表达。(3)MiR-21和BMP9-CM处理C3H10 T1/2细胞,ALP活性和染色实验检测C3H10 T1/2细胞早期成骨能力。茜素红S染色实验检测钙盐沉积情况。(4)MiR-21和BMP9-CM处理C3H10 T1/2细胞,Real-time-PCR检测成骨分化相关因子ALP,OCN的表达。(5)MiR-21和BMP9-CM处理C3H10T1/2细胞,Western blot检测p-Smad1/5蛋白水平的表达。结果:(1)BMP9暂时降低miR-21的表达。MiR-21也可以暂时降低BMP9的表达。(2)MiR-21可以协同BMP9增强ALP和钙盐沉积。(3)MiR-21协同BMP9增加了p-Smad1/5蛋白水平的表达。结论:MiR-21与BMP9存在相互关系,两者可以互相调节表达。MiR-21可以协同BMP9促进间充质干细胞C3H10T1/2细胞成骨分化,这一过程与增强BMP9/Smad信号的激活程度有关。     

CHOP T/C and C/T haplotypes contribute to early-onset type 2 diabetes in Italians

Type 2 diabetes (T2D) is characterized by impaired insulin secretion, insulin insensitivity and decreased beta-cell mass. Multiple genes contribute to T2D. The chromosome 12q13.1 region is in linkage to T2D in different populations, including our Italian dataset. CHOP is a candidate gene for the linkage, as it is located in the chromosome 12q13.1 region, and may contribute to T2D by increasing beta-cell apoptosis susceptibility and by impairing insulin sensitivity. Our goal was to identify any potential CHOP gene variants contributing to T2D in our Italian early-onset T2D families, which show linkage to the CHOP region. We directly sequenced the CHOP gene in 28 Italian probands of the linked T2D families and in 115 control subjects. We performed genotype and haplotype association tests with T2D of the identified single nucleotide polymorphisms (SNPs). We performed model-free and parametric association haplotype tests with T2D. We identified three SNPs [5'UTR-c.279T > C, 5'UTR-c.120A > G and + nt30C > T (F10F)] in CHOP. These SNPs are in complete linkage disequilibrium. The genotype association test showed an association trend with T2D of TT (F10F) and AG (-c.120A > G). The haplotype association test provided significant results for the haplotypes T/C (frequency = 0.33) and C/T (frequency = 0.01) (at 5'UTR-c.279T > C and + nt30C > T, respectively) under non-parametric analysis (P-value = 0.0000), recessive model (P-value = 0.0000) and additive model (P-value = 0.0014). Our data show that CHOP described haplotypes T/C and C/T, as an additive and as a homozygous variant, contribute significantly to T2D in our Italian early-onset group. We conclude that the CHOP T/C and C/T haplotype contributes to our T2D linkage signal on chromosome 12q13.1.     

Polymorphisms of MRP2 (ABCC2) are associated with susceptibility to nonalcoholic fatty liver disease

AimsWe hypothesized that ABCC2 gene variants may contribute to susceptibility to nonalcoholic fatty liver disease (NAFLD). Additionally, we tested the hypothesis of a relation between gene variants and disease severity.Patients and MethodsThe study involved 167 individuals: 109 consecutively presenting unrelated patients with features of NAFLD and different stages of disease severity, and a group of 58 healthy individuals. Four tag single-nucleotide polymorphisms (SNPs; rs717620 A/G, rs2756105 C/T, rs2002042 C/T and rs3740066 A/G) representing 46 polymorphic sites (r²>.8) were genotyped. Furthermore, two additional SNPs (rs17222723 A/T and rs8187710 G/A) were included.ResultsOn univariate analysis, after multiple comparison correction by permutation tests, there were significant differences observed in the allele frequencies of rs17222723 and rs8187710 between healthy individuals and NAFLD patients (empirical P=.037 and .035, respectively). Allelic odds ratios [95% confidence interval] for rs17222723 and rs8187710 were 2.80 [1.11–7.04] and 2.80 [1.11–7.04], respectively. When we tested the hypothesis of a relation between gene variants and the clinical and histological spectra of NAFLD by multinomial regression analysis, a significant association was observed with the same markers: rs17222723 (P=.0029) and rs8187710 (P=.015).ConclusionsOur study suggests a potential role of ABCC2 in susceptibility to NAFLD and disease severity.     

Genetic Variants in CASP3, BMP5, and IRS2 Genes May Influence Survival in Prostate Cancer Patients Receiving Androgen-Deprivation Therapy

Several genome-wide association studies (GWAS) have been conducted to identify the common single nucleotide polymorphisms (SNPs) that influence the risk of prostate cancer. It was hypothesized that some prostate cancer-associated SNPs might relate to the clinical outcomes in patients treated for prostate cancer using androgen-deprivation therapy (ADT). A cohort of 601 patients who have received ADT for prostate cancer was genotyped for 29 SNPs that have been associated with prostate cancer in Cancer Genetic Markers of Susceptibility GWAS, and within the genes that have been implicated in cancer. Prognostic significance of these SNPs on the disease progression, prostate cancer-specific mortality (PCSM) and all-cause mortality (ACM) after ADT were assessed by Kaplan-Meier analysis and Cox regression model. Three SNPs, namely CASP3 rs4862396, BMP5 rs3734444 and IRS2 rs7986346, were found to be closely associated with the ACM (P≤0.042), and BMP5 rs3734444 and IRS2 rs7986346 were also noted to be significantly related to the PCSM (P≤0.032) after adjusting for the known clinicopathologic predictors. Moreover, patients carrying a greater number of unfavorable genotypes at the loci of interest had a shorter time to ACM and PCSM during ADT (P for trend <0.001). Our results suggest that CASP3 rs4862396, BMP5 rs3734444 and IRS2 rs7986346 may affect the survival in patients after ADT for prostate cancer, and the analysis of these SNPs can help identify patients at higher risk of poor outcome.     

孤儿核受体SHP敲减抑制BMP9诱导C3H10T1/2细胞的成骨分化

孤儿核受体SHP(small heterodimer partner)是核受体超家族中的一员,具有LXXLL模体及配体结合域,但无经典的DNA结合域.它可与多种转录因子结合,调节细胞的增殖、分化和代谢等生物学过程.但目前关于SHP在BMP9诱导成骨分化中的确切作用却尚不清楚.本研究证明,SHP参与BMP9诱导的C3H10T1/2细胞成骨分化. RT-PCR结合Western印迹方法检测蛋白揭示,异位表达BMP9上调了SHP在C3H10T1/2细胞中的表达. 小干扰RNA敲减SHP基因在C3H10T1/2细胞的表达下调了成骨相关基因Runx2、Id1、Id2及CTGF的表达,而过表达BMP9则可上调这些基因的表达.碱性磷酸酶(ALP)活性测定/染色及茜素红染色显示,敲减核受体SHP基因可抑制BMP9的成骨分化作用,而过表达BMP9可部分消除SHP 敲减导致的成骨抑制作用.上述结果提示,核受体SHP为BMP9诱导的C3H10T1/2细胞成骨分化所必需. 究竟BMP9如何上调SHP基因表达,以及SHP究竟通过何种机制上调BMP9下游成骨分化相关基因的表达尚待进一步研究.     

CREB/Wnt10b mediates the effect of COX-2 on promoting BMP9-induced osteogenic differentiation via reducing adipogenic differentiation in mesenchymal stem cells

Bone morphogenetic protein 9 (BMP9) is one of the most potent osteogenic factors, which may be a potential candidate for bone tissue engineering. However, the osteogenic capacity of BMP9 still need to be further enhanced. In this study, we determined the effect of Wnt10b on BMP9-induced osteogenic differentiation in mesenchymal stem cell (MSCs) and the possible mechanism underlying this process. We introduced the polymerase chain reaction (PCR), Western blot analysis, histochemical stain, ectopic bone formation, and microcomputed tomography analysis to evaluate the effect of Wnt10b on BMP9-induced osteogenic differentiation. Meanwhile, PCR, Western blot analysis, chromatin immunoprecipitation, and immunoprecipitation were used to analyze the possible relationship between BMP9 and Wnt10b. We found that BMP9 upregulates Wnt10b in C3H10T1/2 cells. Wnt10b increases the osteogenic markers and bone formation induced by BMP9 in C3H10T1/2 cells, and silencing Wnt10b decreases these effects of BMP9. Meanwhile, Wnt10b enhances the level of phosphorylated Smad1/5/8 (p-Smad1/5/8) induced by BMP9, which can be reduced by silencing Wnt10b. On the contrary, Wnt10b inhibits adipogenic markers induced by BMP9, which can be decreased by silencing Wnt10b. Further analysis indicated that BMP9 upregulates cyclooxygenase-2 (COX-2) and phosphorylation of cAMP-responsive element binding (p-CREB) simultaneously. COX-2 potentiates the effect of BMP9 on increasing p-CREB and Wnt10b, while silencing COX-2 decreases these effects. p-CREB interacts with p-Smad1/5/8 to bind the promoter of Wnt10b in C3H10T1/2 cells. Our findings suggested that Wnt10b can promote BMP9-induced osteogenic differentiation in MSCs, which may be mediated through enhancing BMP/Smad signal and reducing adipogenic differentiation; BMP9 may upregulate Wnt10b via the COX-2/p-CREB-dependent manner.     

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.