Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
J Prev Med Public Health. Author manuscript; available in PMC 2010 June 9.
Published in final edited form as:
PMCID: PMC2882880

Differential Parental Transmission of Markers in BCL3 among Korean Cleft Case-parent Trios



Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of approximately 1 in 700 live births. The B-Cell Leukemia/lymphoma 3 (BCL3) gene has been suggested as a candidate gene for CL/P based on association and linkage studies in some populations. This study tests for an association between markers in BCL3 and isolated, non-syndromic CL/P using a case-parent trio design, while considering parent-of-origin effects.


Forty case-parent trios were genotyped for two single nucleotide polymorphisms (SNPs) in the BCL3 gene. We performed a transmission disequilibrium test (TDT) on individual SNPs, and the FAMHAP package was used to estimate haplotype frequencies and to test for excess transmission of multi-SNP haplotypes.


The odds ratio for transmission of the minor allele, OR (transmission), was significant for SNP rs8100239 (OR=3.50, p=0.004) and rs2965169 (OR=2.08, p=0.027) when parent-of-origin was not considered. Parent-specific TDT revealed that SNP rs8100239 showed excess maternal transmission. Analysis of haplotypes of rs2965169 and rs8100239 also suggested excess maternal transmission.


BCL3 appears to influence risk of CL/P through a parent-of-origin effect with excess maternal transmission.

Keywords: BCL3, Oral cleft, Maternal transmission effects, Parent-of-origin


Oral clefts are one of the most common birth defects in humans, and represent a significant public health problem both in terms of the medical and economic burden for affected individuals and their families. Non-syndromic cleft lip with or without palate (CL/P) is ‘complex’ or ‘multifactorial’ in its etiology, in that both genes and environmental risk factors determine risk [1,2]. Although several candidate genes have been extensively studied in different populations (tranforming growth factor alpha, interferon regulatory factor 6, retinoic acid receptor alpha, etc), only a few genes have been shown to contain mutations that appear causal (msh homeobox 1, poliovirus receptor-related 1, etc.), and these are rare and often show incomplete penetrance [3-5].

The B-cell leukemia/lymphoma 3 (BCL3) gene is located on chromosome 19q13, where studies of multiplex families have yielded evidence for linkage to nonsyndromic orofacial clefts [6,7]. Several other studies also observed an association between markers in the BCL3 gene and CL/P [6,8]. Although these studies suggested candidate genes, there have not been many studies on whether the BCL3 gene is a risk factor for CL/P in Asian populations.

It is important to consider parent-of-origin effects when studying birth defects because maternal genotype controls the in utero environment of the developing fetus, and separating maternal genotypic effects from imprinting effects remains an important question [9,10]. Maternal parent-of-origin effects have been suggested for several genes associated with non-syndromic CL/P [11-14]. In this paper, we tested for an association between markers in BCL3 and the risk of CL/P in 40 Korean case-parent trios, with specific consideration of parent-of-origin effects.


I. Sample Description

As part of an international study of oral clefts, we enrolled 40 unrelated Korean patients aged 6 months to 19 years and their parents through the department of plastic surgery, Yonsei university medical center (Seoul, Korea) from January 2003 to March 2004. Parents of the cases were interviewed regarding family history, medical history, and exposure to suspected risk factors. The patient and his/her medical records were examined to confirm the classification of non-syndromic CL/P. There were 22 male cases and 18 female cases. The mother’s and father’s mean age at proband’s birth were 30.5 and 33.4. The institutional review boards of Yonsei university and the Johns Hopkins Bloomberg school of public health approved this study. All parents received adequate information about this study and gave written informed consent.

II. SNP Selection, DNA, & Genotyping

Single nucleotide polymorphisms (SNPs) were selected in a region surrounding BCL3 on chromosome 19q13, with a goal of identifying one SNP per 5 kb of physical distance. Variants with “SNP scores” (an assessment of design quality of the Illumina assay based on a proprietary algorithm) above 0.6, high validation levels in dbSNP (this included validation levels where the submitter had validated the SNP on multiple platforms), and high heterozygosity levels (particularly in multiple populations) were given higher priority during the selection process. From seven selected SNPs, two SNPs were found to be polymorphic in the Korean population (Table 1).

Table 1
SNP minor allele frequencies among parents of CL/P cases in Korea

Genomic DNA samples were prepared from peripheral blood using the previously described protein precipitation method [15]. DNA concentration was determined using the PicoGreen® dsDNA Quantitation Kit (Molecular Probes Inc., Eugene, OR, USA), and all DNA samples were stored at −20°C. A 4 μg aliquot of each genomic DNA sample was dispensed into a bar-coded 96-well microtiter plate at a concentration of 100 ng/μl, and was subsequently genotyped for SNP markers using the Illumina Golden-Gate™ chemistry with Sentrix® Array Matrices (Illumina, San Diego, USA) [16] at the SNP center of the genetic resources core facility (GRCF), a part of the McKusick-Nathans institute of genetic medicine, Johns Hopkins school of medicine. Two duplicates and four centre d’etude du polymorphisme humain (CEPH) controls were included on each plate to evaluate genotyping consistency within and between plates, and to insure correct orientation. Genotypes were generated on a BeadLab 1000 system (Illumina, San Diego, USA) [17].

III. Statistical Analysis

The minor allele frequency (MAF) was computed among parents, and pairwise linkage disequilibrium (LD) was computed as the R-square value for all SNPs using the Haploview program (Broad institute, Cambridge, USA) [18]. The standard transmission disequilibrium test (TDT) described by Speilman et al. [19] was used to test for excess transmission of individual alleles. Parent-of-origin effects were examined using Clayton’s extension of the TDT incorporated into STATA 8.2 (Stata Corporation, College station, USA), which stratifies the standard TDT into separate allelic tests for fathers and mothers [20].

The FAMHAP package (IMBIE, Bonn, Germany) was used to estimate haplotype frequencies and to test for excess transmission of multi-SNP haplotypes [6]. The FAMHAP package calculates maximum likelihood estimates (MLEs) of haplotype frequencies (for up to 20 SNPs) from nuclear families with varying numbers of children via the expectation-maximization algorithm, and is robust when handling missing SNPs [21]. This program provides a haplotype-based test for nuclear family data. This test statistic is based on Monte-Carlo simulations, in which the set of transmitted and non-transmitted genotypes/haplotypes is randomly permuted for each replicate [22,23]. In this analysis, the chi-square statistic for marker combinations is replaced by the maximum chi-square over single haplotypes (maximum TDT statistic). The program gives an empiric p-value, corrected for the multiple haplotypes being considered. This haplotype analysis was also carried out separately for maternal and paternal transmission.


Five of the seven SNPs were monomorphic, leaving only two SNPs with reasonable heterozygosity (Table 1). The R-square value between SNP rs2965169 and rs8100239 was 0.38. Only trios with complete data were used for the TDT. When all markers were screened using the TDT without considering the parent of origin, the odds ratio of transmission for the minor allele, OR (transmission), was significant for both SNP rs2965169 (OR=2.08, p=0.027) and SNP rs8100239 (OR=3.50, p=0.004)(Table 2).

Table 2
Number of Transmitted (T) or Non-Transmitted (NT) minor alleles in 40 CLP cases for TDT and estimated odds ratios of transmission OR* (transmission) ignoring parent-of-origin

Parent-of-origin effects were investigated by stratifying informative transmissions (T) and non-transmissions (NT) by parental source for these two SNPs (Table 3). This analysis revealed that SNP rs8100239 showed excess maternal transmission, significant at the p=0.004 level (OR=11.0).

Table 3
Number of Transmitted (T) or Non-Transmitted (NT) minor alleles to 40 CLP cases from TDT and estimated odds ratio considering parent-of-origin

Table 4 shows the results of haplotypes analysis for rs2965169 and rs8100239. In these Korean trios, haplotypes showed evidence of excess transmission of the 2-1 haplotype to CL/P children (p=0.018 for overall transmission). This can be largely attributed to excess maternal transmission (p=0.038).

Table 4
Analysis of haplotypes using rs2965169 and rs8100239 in BCL3 with analyses using the program FAMHAP


Our study of CL/P case-parent trios showed significant evidence of linkage and disequilibrium for SNP rs8100239 in BCL3. In screening for parent-of-origin effects, we found suggestive evidence of excess maternal transmission of this SNP. Haplotypes of rs8100239 and rs2965169 also showed significant deviation from the expected levels when transmitted from mothers, but not from fathers.

BCL3, a proto-oncogene that encodes a transcription factor involved in cell cycle regulation, has been suggested as a candidate gene for CL/P [6,24]. The BCL3 gene has been associated with oral clefts in some association studies [6,8], but not others [25,26]. A possible reason for these conflicting results is that the susceptibility loci may have different contributions in different populations [6]. The present study suggests that BCL3 is associated with oral clefts in Koreans.

Excess maternal transmission could reflect genomic imprinting or maternal genotype effects. Maternal genotypic effects for non-syndromic cleft lip with/without palate (CL/P) have also been reported for several other candidate genes (5,10 methylenetetrahydrofolate reductase (MTHFR) and cystathionine beta synthase), but these have yet to be confirmed [3,12]. Recently, Reuter et al. [27] found parent-of-origin effects in the transforming growth factor beta 3 gene among central Europeans with nonsyndromic cleft lip and palate, and found a lower risk of maternal transmission compared to paternal transmission. Our results for rs8100239 in BCL3 (analyzed both alone and as a haplotype) showed evidence of excess maternal transmission, which could reflect an imprinting effect or a maternal genotype effect.

Several studies have also suggested that the BCL3 gene may be a modifier or may exert an additive effect [6,28]. In a family-based association study in a Brazilian population, Gaspar et al. [28] observed an interaction between the maternal MTHFR and offspring’s BCL3 genotypes.

The case-parent trio design offers the advantage of testing directly for maternal vs. paternal effects, and allows the separation of these effects from the effects of the fetal genotype vs. parental origin in a robust manner [20,29,30]. Another advantage of this study design is that it minimizes issues of confounding that plague traditional case-control designs. The small sample size of the present study may not provide a reliable result (e.g. OR=11.0 in Table 3), and further replication studies are necessary. Another limitation of this study was that only two SNPs of seven SNPs were polymorphic in this study because we selected SNPs from Caucasian samples. Further studies should be conducted with more polymorphic SNPs for the Korean population.

The present study suggests maternal transmission effects for markers in BCL3 and risk of non-syndromic CL/P. Further work will be needed to confirm this suggestion that maternal transmission of alleles in BCL3 influences the risk of CL/P, and to determine its ultimate impact on risk.


This research was supported by the Korean Research Foundation (2004-041-E00104).


1. Cobourne MT. The complex genetics of cleft lip and palate. Eur J Orthod. 2004;26(1):7–16. [PubMed]
2. Wyszynski DF, Duffy DL, Beaty TH. Maternal cigarette smoking and oral clefts: A meta-analysis. Cleft Palate Craniofac J. 1997;34(3):206–210. [PubMed]
3. Carinci F, Scapoli L, Palmieri A, Zollino I, Pezzetti F. Human genetic factors in non-syndromic cleft lip and palate: An update. Int J Pediatr Otorhinolaryngol. 2007;71(10):1509–1519. [PubMed]
4. van den Boogaard MJ, Dorland M, Beemer FA, van Amstel HK. MSX1 mutation is associated with orofacial clefting and tooth agenesis in humans. Nat Genet. 2000;24(4):342–343. [PubMed]
5. Zucchero TM, Cooper ME, Maher BS, Daack-Hirsch S, Nepomuceno B, Ribeiro L, et al. Interferon regulatory factor 6 (IRF6) gene variants and the risk of isolated cleft lip or palate. N Engl J Med. 2004;351(8):769–780. [PubMed]
6. Gaspar DA, Matioli SR, Pavanello RC, Araújo BC, André M, Steman S, et al. Evidence that BCL3 plays a role in the etiology of nonsyndromic oral clefts in Brazilian families. Genet Epidemiol. 2002;23(4):364–374. [PubMed]
7. Martinelli M, Scapoli L, Pezzetti F, Carinci F, Carinci P, Baciliero U, et al. Suggestive linkage between markers on chromosome 19q13.2 and nonsyndromic orofacial cleft malformation. Genomics. 1998;51(2):177–181. [PubMed]
8. Beaty TH, Wang H, Hetmanski JB, Fan YT, Zeiger JS, Liang KY, et al. A case-control study of nonsyndromic oral clefts in Maryland. Ann Epidemiol. 2001;11(6):434–442. [PubMed]
9. Weinberg CR, Umbach DM. A hybrid design for studying genetic influences on risk of diseases with onset early in life. Am J Hum Genet. 2005;77(4):627–636. [PubMed]
10. Wilkins JF, Haig D. What good is genomic imprinting: The function of parent-specific gene expression. Nat Rev Genet. 2003;4(5):359–368. [PubMed]
11. Jugessur A, Wilcox AJ, Lie RT, Murray JC, Taylor JA, Ulvik A, et al. Exploring the effects of methylenetetrahydrofolate reductase gene variants C677T and A1298C on the risk of orofacial clefts in 261 Norwegian case-parent triads. Am J Epidemiol. 2003;157(12):1083–1091. [PubMed]
12. Rubini M, Brusati R, Garattini G, Magnani C, Liviero F, Bianchi F, et al. Cystathionine beta-synthase c.844ins68 gene variant and non-syndromic cleft lip and palate. Am J Med Genet A. 2005;136A(4):368–372. [PubMed]
13. Sull JW, Liang KY, Hetmanski JB, Fallin MD, Ingersoll RG, Park J, et al. Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations. Genet Epidemiol. 2008;32(6):505–512. [PMC free article] [PubMed]
14. van Rooij IA, Vermeij-Keers C, Kluijtmans LA, Ocke MC, Zielhuis GA, Goorhuis-Brouwer SM, et al. Does the interaction between maternal folate intake and the methylenetetrahydrofolate reductase polymorphisms affect the risk of cleft lip with or without cleft palate? Am J Epidemiol. 2003;157(7):583–591. [PubMed]
15. Bellus GA, Hefferon TW, de Luna RI Ortiz, Hecht JT, Horton WA, Machado M, et al. Achondroplasia is defined by recurrent G380R mutations of FGFR3. Am J Hum Genet. 1995;56(2):368–373. [PubMed]
16. Oliphant A, Barker DL, Stuelpnagel JR, Chee MS. BeadArray™ technology: Enabling an accurate, cost-efficient approach to high-throughput genotyping. Biotechniques. 2002;(Suppl):56–61. [PubMed]
17. Fan JB, Oliphant A, Shen R, Kermani BG, Garcia F, Gunderson KL, et al. Highly parallel SNP genotyping. Cold Spring Harb Symp Quant Biol. 2003;68:69–78. [PubMed]
18. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: Analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–265. [PubMed]
19. 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(3):506–516. [PubMed]
20. 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(3):167–185. [PubMed]
21. Becker T, Knapp M. Maximum-likelihood estimation of haplotype frequencies in nuclear families. Genet Epidemiol. 2004;27(1):21–32. [PubMed]
22. Knapp M, Becker T. Family-based association analysis with tightly linked markers. Hum Hered. 2003;56(1-3):2–9. [PubMed]
23. Zhao H, Zhang S, Merikangas KR, Trixler M, Wildenauer DB, Sun F, et al. Transmission/disequilibrium tests using multiple tightly linked markers. Am J Hum Genet. 2000;67(4):936–946. [PubMed]
24. Franzoso G, Bours V, Azarenko V, Park S, Tomita-Yamaguchi M, Kanno T, et al. The oncoprotein Bcl-3 can facilitate NF-kappa B-mediated transactivation by removing inhibiting p50 homodimers from select kappa B sites. EMBO J. 1993;12(10):3893–3901. [PubMed]
25. Fujita H, Nagata M, Ono K, Okubo H, Takagi R. Linkage analysis between BCL3 and nearby genes on 19q13.2 and non-syndromic cleft lip with or without cleft palate in multigenerational Japanese families. Oral Dis. 2004;10(6):353–359. [PubMed]
26. Suazo J, Santos JL, Silva V, Jara L, Palomino H, Blanco R. Possible association due to linkage disequilibrium of TGFA, RARA and BCL3 with nonsyndromic cleft lip with or without cleft palate in the Chilean population. Rev Med Chil. 2005;133(9):1051–1058. (Spanish) [PubMed]
27. Reutter H, Birnbaum S, Mende M, Lauster C, Schmidt G, Henschke H, et al. TGFB3 displays parent-of-origin effects among central Europeans with nonsyndromic cleft lip and palate. J Hum Genet. 2008;53(7):656–661. [PubMed]
28. Gaspar DA, Matioli SR, de Cássia Pavanello R, Araújo BC, Alonso N, Wyszynski D, et al. Maternal MTHFR interacts with the offspring’s BCL3 genotypes, but not with TGFA, in increasing risk to nonsyndromic cleft lip with or without cleft palate. Eur J Hum Genet. 2004;12(7):521–526. [PubMed]
29. Sinsheimer JS, Palmer CG, Woodward JA. Detecting genotype combinations that increase risk for disease: Maternal-fetal genotype incompatibility test. Genet Epidemiol. 2003;24(1):1–13. [PubMed]
30. 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(3):294–303. [PubMed]