Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Birth Defects Res A Clin Mol Teratol. Author manuscript; available in PMC 2013 March 26.
Published in final edited form as:
PMCID: PMC3608124

Evidence of Gene-Environment Interaction for the RUNX2 Gene and Environmental Tobacco Smoke in Controlling the Risk of Cleft Lip with/without Cleft Palate


This study examined the association between 49 markers in the Runt-related transcription factor 2 (RUNX2) gene and nonsyndromic cleft lip with/without cleft palate (CL/P) among 326 Chinese case-parent trios, while considering gene-environment (GxE) interaction and parent-of-origin effects. Five single-nucleotide polymorphisms (SNPs) showed significant evidence of linkage and association with CL/P and these results were replicated in an independent European sample of 825 case-parent trios. We also report compelling evidence for interaction between markers in RUNX2 and environmental tobacco smoke (ETS). Although most marginal SNP effects (i.e., ignoring maternal exposures) were not statistically significant, eight SNPs were significant when considering possible interaction with ETS when testing for gene (G) and GxE interaction simultaneously or when considering GxE alone. Independent samples from European populations showed consistent evidence of significant GxETS interaction at two SNPs (rs6904353 and rs7748231). Our results suggest genetic variation in RUNX2 may influence susceptibility to CL/P through interacting with ETS.

Keywords: RUNX2, oral clefts, gene-environment interaction, parent-of-origin effects, imprinting


Nonsyndromic cleft lip with or without cleft palate (CL/P), which includes cases with cleft lip alone plus cases with cleft lip and cleft palate is one of the most common human birth defects (Mossey, 2007). Chinese people have higher prevalence rates of CL/P among liveborn infants compared to other racial groups (Cooper et al., 2006). CL/P is considered to be a complex disease because both genetic and environmental risk factors contribute to its etiology (Mossey et al., 2009). The Runt-related transcription factor 2 (RUNX2) gene is a member of the RUNX family of transcription factors, and its gene product is involved in osteoblast differentiation and bone morphogenesis (Komori et al., 1997). Mutations in RUNX2 cause cleidocranial dysplasia, a rare autosomal dominant inherited skeleton condition where clefts of the palate or submucous palate can be present (Lee et al., 1997; Mundlos et al., 1997; Otto et al., 1997; Cooper et al., 2001). Previous studies of CL/P multiplex families yielded some evidence of linkage on chromosome 6p, where this gene is located (Eiberg et al., 1987). An association analyses of case-parent trios from four populations also suggested RUNX2 may influence risk of CL/P through imprinting effects (Sull et al., 2008). In this study, case-parent trios from Maryland, Taiwan, Singapore, and Korea were genotyped for the RUNX2 gene and three single nucleotide polymorphisms (SNPs) showed significant excess maternal transmission. Gene-environment (GxE) interaction has been suggested for several genes associated with CL/P (Beaty et al., 2002; Mossey, 2007; Shi et al., 2008; Sull et al., 2009; Zhu et al., 2009), although to our knowledge no studies have investigated this question for RUNX2.

In this study, we tested for association between SNPs in RUNX2 and CL/P using 326 Chinese case-parent trios, which are independent of the previous report by Sull et al. (2008), which considered maternal genotypic effects. Here we explicitly test for imprinting effects and GxE interaction with common maternal exposures, including environmental tobacco smoke (ETS) and multivitamin supplementation.


Sample description

The “International Genetic Epidemiology of Oral Clefts” study is a multicenter, international family-based study initiated in 2003 to investigate the genetic etiology of oral clefts. As part of this study, case-parent trios were recruited from three sites in mainland China (Weifang, Shandong Province; Wuhan, Hubei Province; Chengdu, Sichuan Province) and Taiwan. Research protocols were reviewed and approved by institutional review boards at each institution. Informed consent was obtained from parents. The majority of cases were infants seen during a surgical or postsurgical visit. All probands were examined for other congenital anomalies or major developmental delays, and were classified as having an isolated nonsyndromic CL/P. Ethnicity and other demographic data were obtained through structured interviews. Maternal exposure information, including cigarette smoking, ETS, multivitamin supplementation, and alcohol consumption was collected through direct interview of the mothers. Environmental exposures were defined as being exposed from 3 months before pregnancy through the first trimester, except for ETS whereas exposure was defined as being exposed at any time from 3 months before pregnancy through the entire pregnancy. The proportion of infants exposed to maternal cigarette smoking and alcohol consumption was very low (around 1%), so only maternal ETS and multivitamin supplementation were analyzed. Table 1 presents information on gender, family history, and maternal exposures among CL/P probands. None of the parents of these CL/P cases were themselves affected. The proportion of exposure to ETS and multivitamin supplementation were 40.5 and 8.6%, respectively.

Table 1
Characteristics of 326 Chinese Nonsyndromic Cleft Lip with or without Cleft Palate (CL/P) Cases

Single-Nucleotide Polymorphism Selection and Genotyping

SNPs in and around RUNX2 on chromosome 6p21 were identified from the literature, dbSNP (, and based on results of our previous study (Sull et al., 2008). A set of 69 SNPs was chosen based on the criteria of high “design scores” as provided by Illumina (San Diego, CA), high heterozygosity, and HapMap validation ( Genomic DNA samples were prepared from peripheral blood lymphocytes by the protein precipitation method (Bellus et al., 1995). Primers for each SNP were synthesized using the Oligator technology by Illumina as part of an oligo pool for the BeadLab 1000 system (San Diego, California). A 4 μg aliquot of each genomic DNA sample was dispensed into a barcoded 96-well microtiter plate at a concentration of 100 ng/μl and genotyped using the Illumina Golden-Gate chemistry at the SNP Center of the Genetic Resources Core Facility, a part of the McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins School of Medicine. The average distance between neighboring markers was 4.814 kb (based on Build 36.1 of dbSNP). From these 69 selected SNPs, 20 SNPs were monomorphic in this sample and were dropped. No Mendelian inconsistencies were found for the remaining 49 SNPs. Two duplicates and four Centre d’Etude du Polymorphisme Humain control DNA samples were included to evaluate genotyping consistency within and between plates and to insure correct plate orientation. The minor allele frequencies of the 49 SNPs ranged from 5.0% to 47.6% (Supplementary Table S1).

Statistical analysis

Minor allele frequencies were computed using the parents. Pairwise linkage disequilibrium (LD) was measured as r2 for all SNPs using the Haploview program, and was used to identify LD blocks (Barrett et al., 2005). Nine independent SNPs and eight blocks of LD were identified (Fig. 1). Clayton’s extension of the Transmission Disequilibrium Test (TDT) incorporated into STATA 10 was used on individual SNPs to test for evidence of linkage or LD (Spielman et al., 1993; Cordell et al., 2004). To assess the significance of TDT findings while considering multiple tests, the permutation test in PLINK was used to get corrected p values (PLINK,; Purcell et al., 2007). A Bonferroni test would be too conservative because there is strong LD among these markers in RUNX2. In the permutation test, the LD structure was preserved and the permuted p values represent the chance of getting a test statistic as large as, or larger than the observed one over 10,000 tests. The family-based association test program (FBAT, was used for haplotype analysis (Rabinowitz and Laird, 2000). FBAT uses a score test to compare the observed-to-expected genotypes among offspring under the null hypothesis of no association between observed genotype counts and the phenotype. This approach has been extended to consider haplotypes of several SNPs (Laird et al., 2000; Rabinowitz and Laird, 2000; Horvath et al., 2004). We used a sliding window approach for haplotypes to systematically analyze all adjacent 2-SNP and 3-SNP combinations. Plots of −log10 (p value) for individual SNPs and haplotypes were generated using the R package snp.plotter (Luna and Nicodemus, 2007).

Figure 1
Significance of individual single-nucleotide polymorphisms (SNPs) and sliding window haplotypes for the RUNX2 gene in 326 nonsyndromic cleft lip with or without cleft palate trios. The −log10 (p value) for the overall test for an individual SNP ...

We screened for possible GxE interaction between all SNPs in RUNX2 and two common maternal exposures: ETS and multivitamin supplementation. We used conditional logistic regression under additive models to estimate the relative risk (RR) for being a case with or without environmental exposures for each SNP (Vansteelandt et al., 2008). Estimates of RR (CL/P|G without E) and RR (CL/P|G with E) provide an estimate of the allelic impact of an SNP on the risk of having CL/P. To ensure validly testing for GxE interaction, we also used conditional logistic regression to test for a log-additive GxE interaction while saturating the genetic main effects (by allowing two RR parameters in the model, one for heterozygotes and one for homozygotes). We also used the general strategy proposed by Vansteelandt et al. (2008) to estimate the haplotype-environment interaction, where family-based association tests are evaluated while allowing for a potential GxE interaction in a 2 degree of freedom (df) test, followed by a separate 1 df test for GxE interaction alone (Schaid, 1996; Kraft et al., 2007). This approach is implemented in the PBAT package (v 3.6; The 2 df test examines the genetic effect of a haplotype while taking into account possible effects of GxE interaction, while the 1 df test investigates the effect of GxE interaction alone.

Parent-of-origin effects were also tested using several methods. First, the transmission asymmetry test (TAT) was used to check for potential parent-of-origin effects (Weinberg et al., 1998). A significant TAT could reflect an imprinting effect or a maternal genotype effect. Therefore, we used the max_Z2 test implemented in the TRIMM package to test for maternal genotypic effects separately from effects of the child’s genotype (Shi et al., 2007). This approach maintains LD patterns among SNPs, and generates empiric p values through permutation. Finally, we used the parent-of-origin likelihood ratio test (PO-LRT) developed by Weinberg (1999) to test for imprinting effects as implemented in the program LEM (van Den Oord and Vermunt, 2000).

Replication Study

We carried out an in silico replication study using samples of European and Asian ancestry drawn from a genome-wide association study (GWAS) of nonsyndromic CL/P using the case-parent design (Beaty et al., 2010). The Chinese samples presented here were part of this International Cleft Consortium GWAS, but were dropped for this replication study to insure an independent sample.

Seventy-eight SNPs in RUNX2 were genotyped in the independent replication samples from the International Cleft Consortium with 2491 individuals from 825 European CL/P families and 2151 individuals from 717 Asian CL/P trios (Supplementary Table S2, online only). Samples were genotyped using the Illumina Human610-Quad v.1_B BeadChip (Illumina). Supplementary Table S3, online only, presents detailed information on these 78 SNPs. TDT analysis in PLINK, tests of GxE interaction using conditional logistic regression models, and PBAT were performed to confirm the findings seen among these Chinese trios.


No SNPs showed evidence of deviating from Hardy–Weinberg equilibrium (all were p > 0.10). Figure 1 shows results of the FBAT analyses on individual SNPs and haplotypes. Among 49 SNPs in RUNX2, five SNPs showed a statistically significant association with CL/P at a nominal value of p < 0.05 (Table 2). SNP rs545289 remained significant after permutation test with permuted p value = 0.035. Sliding windows of haplotypes consisting of two, three, four, and five SNPs were also tested. All haplotypes, including rs545289, also showed nominally significant p values (0.035 < p < 0.05 in Fig. 1), possibly reflecting effects of this single SNP. In the replication sample, five SNPs (rs2153277, rs1200426, rs1200425, rs9395112, and rs2782660) yielded significant evidence of association in the combined European samples (Supplementary Table S4, online only). SNP rs1200425 was identified in both original Chinese samples and independent European replication samples, and it was in high LD with rs545289, the most significant SNP in the Chinese data, based on LD information of Hap-Map CEU samples (D’ = 1, r2 = 0.924). However, none of these SNPs showed significant association in TDT among the remaining Asian samples in the International Cleft Consortium.

Table 2
TDT analysis of 5 significant SNPs in RUNX2 gene among 326 Chinese non-syndromic CL/P case parent trios

Conditional logistic regression under an additive model was used to test for possible GxE interaction and estimate exposure-specific RRs. Twelve SNPs yielded nominal significance in the likelihood ratio test (LRT) examining G effects and GxETS interaction with 2 df, although eight of these did not show significant SNP effects in the TDT alone. Figure 2A shows estimated RRs and their 95% confidence intervals (CIs) in exposed and unexposed groups for these SNPs. Among these 12 SNPs, 8 SNPs (rs910586, rs765724, rs2790103, rs9472494, rs2396442, rs6904353, rs7748231, and rs10948237) also showed significant p values in the 1 df test for GxETS alone. For SNP rs765724, being a heterozygous child of an exposed mother was associated with a slightly increased risk (RR = 5.00; 95% CI = 1.45-17.27), but not among children of unexposed mothers (RR = 0.57; 95% CI = 0.31-1.06). The LRT for GxETS interaction in this conditional logistic regression model was also statistically significant (p = 0.00052). We also used a saturated conditional logistic regression model with two genotypic RR parameters in the model to validate these results, but again saw no substantial difference (data not shown).

Figure 2
Testing for gene-environment (GxE) interaction for common maternal exposures in 326 cleft palate (CL/P) case-parent trios from Chinese populations. (A) Single-nucleotide polymorphisms (SNPs) with nominal significant p values in the likelihood ratio test ...

Figure 2B shows results of screening for GxETS interaction between SNPs in RUNX2 using PBAT. Three 2-SNP haplotypes and six 3-SNP haplotypes yielded p < 0.01 in the 2 df test for G and GxETS interaction, although these haplotypes did not include any SNP showing significant individual effects.

SNPs rs6930053, rs2772395, and rs485817 showed nominally significant p values for Gx(vitamin) interaction using the LRT (0.027 < p < 0.05). Haplotypes involving these three SNPs also yielded significant evidence of interaction with exposure. Like ETS, SNPs with no significant marginal effects also gave intriguing evidence of GxE interaction. For example, two 2-SNP haplotypes and three 3-SNP haplotypes (none of which contained SNPs which were individually significant) gave nominally significant p values in both the 2 df and 1 df tests for Gx(vitamin) interaction (0.003 < p < 0.05).

We replicated the intriguing results of GxE interaction among the International Cleft Consortium samples. Five SNPs in RUNX2 yielded significant p values in the LRT with 2 df for effects of G and GxETS jointly among European trios, although none of them showed significant gene effect ignoring exposures (Supplementary Table S5, online only). Among these markers, four SNPs also showed significant p values in the 1 df test for GxE effect alone. In the test of Gx(vitamin), two SNPs in the European trios and one SNP in Asian trios were statistically significant in the 2 df LRT. PBAT was used to test for haplotype-environment interactions. In the test of Gx(ETS) interaction, two 3-SNP haplotypes among the European samples, and two 2-SNP haplotypes and one 3-SNP haplotype among the Asian samples yielded p < 0.01 in the 2 df test for G and GxE interaction, although these haplotypes did not include any SNPs showing significant individual effects. The results of tests for haplotype-vitamin interaction also confirmed these findings from the Chinese samples. Three 2-SNP and two 3-SNP haplotypes without SNPs showing any individual gene effects gave significant p values < 0.01 in the 1 df test for GxE interaction alone among European trios.

Parent-of-origin effects were first investigated by stratifying informative transmissions (T) and non-transmissions (NT) by parental source in the TAT. Among SNPs not in perfect LD (i.e., those with r2 < 0.8), three SNPs showed excess maternal transmission at p < 0.05 level (Table 3). TRIMM was used to test for maternal genotype effects (Shi et al., 2007). However, none showed any evidence of maternal genotype effect. The PO-LRT was used to test for imprinting effects controlling for maternal genotype effect. SNP rs675613 yielded a marginally significant p value (RR = 2.2; p = 0.046), suggesting the maternally derived copy of the minor allele at rs675613 was associated with a two-fold increase in risk compared to receiving a paternally derived copy of this allele (Table 3).

Table 3
Parental transmission and maternal genotypic effect of three SNPs in RUNX2 in 326 CL/P triosa


Our results showed consistent evidence of linkage in the presence of association (LD) for five SNPs in the RUNX2 gene (nominal p < 0.05) and one SNP rs545289 remained significant in permutation tests (empiric p = 0.035). Single SNP analysis using conditional logistic regression models and haplotype analysis using PBAT suggested the RUNX2 gene may also influence the risk of CL/P through interaction with ETS and multivitamin supplementation. We also found one SNP (rs675613) showed marginally significant maternal over-transmission after considering maternal genotypic effects (nominal p = 0.046).

The RUNX2 gene is located on chromosome 6p21, and encodes a transcription factor critical for osteoblastic differentiation (Zelzer and Olsen, 2003; Schroeder et al., 2005; Komori, 2010). It has been suggested that RUNX2 may be under the regulation of estrogen in both osteoblasts and osteoclastogenesis (Mundlos et al., 1997; McCarthy et al., 2003). Mutations in RUNX2 lead to cleidocranial dysplasia, a heritable congenital skeletal disorder characterized by partly or completely missing clavicles, dental anomalies, and craniofacial abnormalities including submucous CL/P and cleft lip (Lee et al., 1997; Otto et al., 1997; Cooper et al., 2001; Yamachika et al., 2001; Aberg et al., 2004). Animal models have suggested mutations or expression levels of RUNX2 may be involved in the development of several craniofacial defects (Aberg et al., 2004; Eswarakumar et al., 2004).

Our study found five SNPs in the RUNX2 gene that showed significant evidence of linkage and association with nonsyndromic CL/P among Chinese case-parent trios. In our study of 326 Chinese CL/P case-parent trios, rs545289 yielded marginal association possibly due to LD with untyped variant influencing risk of nonsyndromic CL/P (empiric p = 0.035), and this finding was confirmed among independent Europeans trios. Five SNPs in RUNX2 yielded statistical significance in TDT analysis among these trios of European ancestry. SNP rs1200425 was statistically significant in both the Chinese samples and the independent replication samples, and was in high LD with the most significant SNP (rs545289) among Chinese trios. Under-transmission of the minor allele (T allele) at rs1200425 was associated with decreased risk of CL/P among the Chinese and the replication European trios. However, the association was not replicated in the Asian subjects. The combined Asian replication samples were from different regions (Singapore, Korea, and Philippines). Principal component analysis showed genetic distances between Singapore and Philippines were relatively large compared to the distances among other East Asian samples, which might explain why none of the SNPs yielded significant evidence of association among Asians drawn from the International Cleft Consortium. Most Chinese trios of current study were also genotyped for 78 SNPs in RUNX2 in the GWAS. TDT analysis (Supplementary Table S6, online only) showed rs16873396 yielded significant association with CL/P (nominal p = 0.008) in the GWAS and this further confirmed the finding of the original analysis (nominal p = 0.012; Table 2).

The interaction between susceptibility genes and common environment factors can also influence the etiology of nonsyndromic oral clefts (Graham and Shaw, 2005; Young et al., 2007). Considering GxE interaction will identify more risk loci, lead to better understanding of underlying biologic mechanisms, and provide information for designing effective preventive strategies. However, there is little consistent evidence of interactions between genes and environment risk factors of nonsyndromic CL/P (Mossey et al., 2009; Shi et al., 2008). Moreover, to our knowledge, no study to date has focused on whether GxE interaction influences the association between RUNX2 and the risk of isolated nonsyndromic CL/P. Our study found intriguing evidence of GxETS and Gx(vitamin) interaction among both the Chinese trios and independent replication samples.

Several studies have shown ETS is not only an independent risk factor of CL/P, but may interact with genetic variants (MSX1, CYP1A1, BMP4T538C, and MYH9) to increase risk (Chevrier et al., 2008; Jia et al., 2010; Jianyan et al., 2010). Our study suggests an interaction between ETS and markers in RUNX2. Both single marker and haplotype analysis yielded consistent results of Gx(ETS) interaction among Chinese trios and replication trios. The 2 df test identified SNPs influencing risk of disease when considering maternal exposure. Eight SNPs in RUNX2 showing no significant individual effects yielded nominal significance in both the 2 df and 1 df tests when exposure to ETS was included in the model. This result illustrates the importance of considering possible GxE interaction in the etiology of CL/P.

Maternal vitamin use has also been reported to be inversely associated with the risk of CL/P, and potential interaction between candidate genes and vitamin supplementation has been suggested (Shaw et al., 1998; Johnson and Little, 2008). We also found evidence of GxE interaction for maternal multivitamin supplementation at three SNPs in RUNX2 and these findings were replicated in single SNP and haplotype analysis among the replication samples. These results suggest being exposed to maternal multivitamin supplementation and carrying certain genotypes may lower the risk of the baby developing CL/P.

Maternal genotypic effects for nonsyndromic CL/P have also been reported for several candidate genes (RUNX2, TGFA, and CBS; Rubini et al., 2005; Sull et al., 2008; Sull et al., 2009). Several studies using cellular models or animal models suggested RUNX2 may function through an epigenetic mechanism or by regulating expression of a maternally imprinted gene (Prescott et al., 2001; Flores et al., 2004; Weinstein et al., 2004; Bertaux et al., 2006). Sull et al. (2008) found several SNPs in RUNX2 suggesting possible imprinting effects in a 46.7 kb region (from rs910586 to rs2396442). In our study, we found suggestive evidence of excess maternal transmission for three SNPs (rs16873396, rs16873401, and rs675613) using the TAT, but only rs675613 showed significant imprinting effects (with an estimated RR = 2.2 by the PO-LRT). This SNP is located about 67 kb downstream of the region previously reported by Sull et al. (2008). We did not find evidence for imprinting using the most significant SNPs in Sull’s study and the effect size for rs675613 was smaller in these data. Subjects in Sull’s study were drawn from four populations including Maryland, Singapore, Taiwan, and Korea, which might account for this difference. The PO-LRT estimates an imprinting effect while controlling for possible maternal genotypic effects. By conditioning on the child’s genotype, it also removes dependence on genotype effects of the cases, and remains valid when applied to markers in LD with an unobserved causal variant (Shi et al., 2007). These results suggest that RUNX2 may influence the risk of CL/P through possible maternal-specific effects, although further studies are clearly needed to confirm this finding.

The case-parent trio design is robust to population stratification, and provides a unique opportunity to investigate parent-of-origin effects (Cordell et al., 2004; Starr et al., 2005). The present study found some evidence of potential maternal over-transmission for markers in RUNX2 and risk of nonsyndromic CL/P. Moreover, these results showed evidence for association in the context of GxE interaction between markers in RUNX2 and both ETS and maternal multivitamin supplementation. ETS, in particular, seems to increase the risk of nonsyndromic CL/P in cases carrying certain genotypes. The exposure rate of ETS in this sample was high (40.5%), reflecting the high prevalence of smoking among Chinese men (about 60%; Gu et al., 2004). In addition, mothers in mainland China had much lower use of multivitamin supplementation compared to those in Taiwan (6.8 vs. 14.1%). If this observation is confirmed, such a GxE interaction will create opportunities for effective intervention to reduce the risk of nonsyndromic CL/P.

Supplementary Material

supplemental materials


We thank all of the families at each recruitment site for participating in this study, and we gratefully acknowledge the invaluable assistance of clinical, field, and laboratory staff who contributed to making this work possible. This research was supported by R21-DE-013707 and R01-DE-014581 from the National Institute of Dental & Craniofacial Research. The International Cleft Consortium including genotyping and analysis was supported by the National Institute for Dental and Craniofacial Research through U01-DE-004425; “International Consortium to Identify Genes & Interactions Controlling Oral Clefts”, 2007 to 2009; T.H. Beaty, Principal Investigator. Tao Wu is supported by the International Collaborative Genetics Research Training Program (ICGRTP), NIH D43 TW06176. This research was also supported by the Intramural Research Program of the NIH, National Institute of Environmental Health Sciences. We thank the Smile Train Foundation for supporting cleft research in China. We are particularly grateful to Dr. Jae Woong Sull for his assistance in performing the analyses using LEM.

Supported by a grant from the National Institute of Dental & Craniofacial Research: R21-DE-013707 and R01-DE-014581


Additional Supporting Information may be found in the online version of this article.


  • Aberg T, Cavender A, Gaikwad JS, et al. Phenotypic changes in dentition of Runx2 homozygote-null mutant mice. J Histochem Cytochem. 2004;52:131–139. [PubMed]
  • Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21:263–265. [PubMed]
  • Beaty TH, Hetmanski JB, Zeiger JS, et al. Testing candidate genes for non-syndromic oral clefts using a case-parent trio design. Genet Epidemiol. 2002;22:1–11. [PubMed]
  • Beaty TH, Murray JC, Marazita ML, et al. A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4. Nat Genet. 2010;42:525–529. [PMC free article] [PubMed]
  • Bellus GA, Hefferon TW, Ortiz de Luna RI, et al. Achondroplasia is defined by recurrent G380R mutations of FGFR3. Am J Hum Genet. 1995;56:368–373. [PubMed]
  • Bertaux K, Broux O, Chauveau C, et al. Runx2 regulates the expression of GNAS on SaOs-2 cells. Bone. 2006;38:943–950. [PubMed]
  • Chevrier C, Bahuau M, Perret C, et al. Genetic susceptibilities in the association between maternal exposure to tobacco smoke and the risk of nonsyndromic oral cleft. Am J Med Genet A. 2008;146A:2396–2406. [PubMed]
  • Cooper ME, Ratay JS, Marazita ML. Asian oral-facial cleft birth prevalence. Cleft Palate Craniofac J. 2006;43:580–589. [PubMed]
  • Cooper SC, Flaitz CM, Johnston DA, et al. A natural history of cleidocranial dysplasia. Am J Med Genet. 2001;104:1–6. [PubMed]
  • Cordell HJ, Barratt BJ, Clayton DG. Case/pseudocontrol analysis in genetic association studies: a unified framework for detection of genotype and haplotype associations, gene-gene and gene-environment interactions, and parent-of-origin effects. Genet Epidemiol. 2004;26:167–185. [PubMed]
  • Eiberg H, Bixler D, Nielsen LS, et al. Suggestion of linkage of a major locus for nonsyndromic orofacial cleft with F13A and tentative assignment to chromosome 6. Clin Genet. 1987;32:129–132. [PubMed]
  • Eswarakumar VP, Horowitz MC, Locklin R, et al. A gain-of-function mutation of Fgfr2c demonstrates the roles of this receptor variant in osteogenesis. Proc Natl Acad Sci U S A. 2004;101:12555–12560. [PubMed]
  • Flores MV, Tsang VW, Hu W, et al. Duplicate zebrafish runx2 orthologues are expressed in developing skeletal elements. Gene Expr Patterns. 2004;4:573–581. [PubMed]
  • Graham JM, Jr, Shaw GM. Gene-environment interactions in rare diseases that include common birth defects. Birth Defects Res A Clin Mol Teratol. 2005;73:865–867. [PubMed]
  • Gu D, Wu X, Reynolds K, et al. Cigarette smoking and exposure to environmental tobacco smoke in China: the international collaborative study of cardiovascular disease in Asia. Am J Public Health. 2004;94:1972–1976. [PubMed]
  • Horvath S, Xu X, Lake SL, et al. Family-based tests for associating haplotypes with general phenotype data: application to asthma genetics. Genet Epidemiol. 2004;26:61–69. [PubMed]
  • Jia ZL, Li Y, Chen CH, et al. Association among polymorphisms at MYH9, environmental factors, and nonsyndromic orofacial clefts in western China. DNA Cell Biol. 2010;29:25–32. [PubMed]
  • Jianyan L, Zeqiang G, Yongjuan C, et al. Analysis of interactions between genetic variants of BMP4 and environmental factors with nonsyndromic cleft lip with or without cleft palate susceptibility. Int J Oral Maxillofac Surg. 2010;39:50–56. [PubMed]
  • Johnson CY, Little J. Folate intake, markers of folate status and oral clefts: is the evidence converging? Int J Epidemiol. 2008;37:1041–1058. [PubMed]
  • Komori T, Yagi H, Nomura S, et al. Targeted disruption of Cbfa1 results in a complete lack of bone formation owing to maturational arrest of osteoblasts. Cell. 1997;89:755–764. [PubMed]
  • Komori T. Regulation of bone development and extracellular matrix protein genes by RUNX2. Cell Tissue Res. 2010;339:189–195. [PubMed]
  • Kraft P, Yen YC, Stram DO, et al. Exploiting gene-environment interaction to detect genetic associations. Hum Hered. 2007;63:111–119. [PubMed]
  • Laird NM, Horvath S, Xu X. Implementing a unified approach to family-based tests of association. Genet Epidemiol. 2000;19(Suppl 1):S36–S42. [PubMed]
  • Lee B, Thirunavukkarasu K, Zhou L, et al. Missense mutations abolishing DNA binding of the osteoblast-specific transcription factor OSF2/CBFA1 in cleidocranial dysplasia. Nat Genet. 1997;16:307–310. [PubMed]
  • Luna A, Nicodemus KK. snp.plotter: an R-based SNP/haplotype association and linkage disequilibrium plotting package. Bioinformatics. 2007;23:774–776. [PubMed]
  • McCarthy TL, Chang WZ, Liu Y, Centrella M. Runx2 integrates estrogen activity in osteoblasts. J Biol Chem. 2003;278:43121–43129. [PubMed]
  • Mossey P. Epidemiology underpinning research in the aetiology of orofacial clefts. Orthod Craniofac Res. 2007;10:114–120. [PubMed]
  • Mossey PA, Little J, Munger RG, et al. Cleft lip and palate. Lancet. 2009;374:1773–1785. [PubMed]
  • Mundlos S, Otto F, Mundlos C, et al. Mutations involving the transcription factor CBFA1 cause cleidocranial dysplasia. Cell. 1997;89:773–779. [PubMed]
  • Otto F, Thornell AP, Crompton T, et al. Cbfa1, a candidate gene for cleidocranial dysplasia syndrome, is essential for osteoblast differentiation and bone development. Cell. 1997;89:765–771. [PubMed]
  • Prescott NJ, Winter RM, Malcolm S. Nonsyndromic cleft lip and palate: complex genetics and environmental effects. Ann Hum Genet. 2001;65(Pt 6):505–515. [PubMed]
  • Purcell S, Neale B, Todd–Brown K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]
  • Rabinowitz D, Laird N. A unified approach to adjusting association tests for population admixture with arbitrary pedigree structure and arbitrary missing marker information. Hum Hered. 2000;50:211–223. [PubMed]
  • Rubini M, Brusati R, Garattini G, et al. Cystathionine beta-synthase c.844ins68 gene variant and non-syndromic cleft lip and palate. Am J Med Genet A. 2005;136A:368–372. [PubMed]
  • Schaid DJ. General score tests for associations of genetic markers with disease using cases and their parents. Genet Epidemiol. 1996;13:423–449. [PubMed]
  • Schroeder TM, Jensen ED, Westendorf JJ. Runx2: a master organizer of gene transcription in developing and maturing osteoblasts. Birth Defects Res C Embryo Today. 2005;75:213–225. [PubMed]
  • Shaw GM, Wasserman CR, Murray JC, Lammer EJ. Infant TGF-alpha genotype, orofacial clefts, and maternal periconceptional multivitamin use. Cleft Palate Craniofac J. 1998;35:366–370. [PubMed]
  • Shi M, Umbach DM, Weinberg CR. Identification of risk-related haplotypes with the use of multiple SNPs from nuclear families. Am J Hum Genet. 2007;81:53–66. [PubMed]
  • Shi M, Wehby GL, Murray JC. Review on genetic variants and maternal smoking in the etiology of oral clefts and other birth defects. Birth Defects Res C Embryo Today. 2008;84:16–29. [PMC free article] [PubMed]
  • Spielman RS, McGinnis RE, Ewens WJ. Transmission test for linkage disequilibrium: the insulin gene region and insulin-dependent diabetes mellitus (IDDM) Am J Hum Genet. 1993;52:506–516. [PubMed]
  • Starr JR, Hsu L, Schwartz SM. Assessing maternal genetic associations: a comparison of the log-linear approach to case-parent triad data and a case-control approach. Epidemiology. 2005;16:294–303. [PubMed]
  • Sull JW, Liang KY, Hetmanski JB, et al. Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations. Genet Epidemiol. 2008;32:505–512. [PMC free article] [PubMed]
  • Sull JW, Liang KY, Hetmanski JB, et al. Evidence that TGFA influences risk to cleft lip with/without cleft palate through unconventional genetic mechanisms. Hum Genet. 2009;126:385–394. [PMC free article] [PubMed]
  • van Den Oord EJ, Vermunt JK. Testing for linkage disequilibrium, maternal effects, and imprinting with (in)complete case-parent triads, by use of the computer program LEM. Am J Hum Genet. 2000;66:335–338. [PubMed]
  • Vansteelandt S, Demeo DL, Lasky–Su J, et al. Testing and estimating gene-environment interactions in family-based association studies. Biometrics. 2008;64:458–467. [PubMed]
  • Weinberg CR, Wilcox AJ, Lie RT. A log-linear approach to case-parent-triad data: assessing effects of disease genes that act either directly or through maternal effects and that may be subject to parental imprinting. Am J Hum Genet. 1998;62:969–978. [PubMed]
  • Weinberg CR. Methods for detection of parent-of-origin effects in genetic studies of case–parents triads. Am J Hum Genet. 1999;65:229–235. [PubMed]
  • Weinstein LS, Liu J, Sakamoto A, et al. Minireview: GNAS: normal and abnormal functions. Endocrinology. 2004;145:5459–5464. [PubMed]
  • Yamachika E, Tsujigiwa H, Ishiwari Y, et al. Identification of a stop codon mutation in the CBFA1 runt domain from a patient with cleidocranial dysplasia and cleft lip. J Oral Pathol Med. 2001;30:381–383. [PubMed]
  • Young DW, Hassan MQ, Yang XQ, et al. Mitotic retention of gene expression patterns by the cell fate-determining transcription factor Runx2. Proc Natl Acad Sci U S A. 2007;104:3189–3194. [PubMed]
  • Zelzer E, Olsen BR. The genetic basis for skeletal diseases. Nature. 2003;423:343–348. [PubMed]
  • Zhu H, Kartiko S, Finnell RH. Importance of gene–environment interactions in the etiology of selected birth defects. Clin Genet. 2009;75:409–423. [PubMed]