Plastic Surgery InternationalVolume 2012 (2012), Article ID 247104, 5 pagesdoi:10.1155/2012/247104

Research Article

Preliminary Analysis of the Nonsynonymous Polymorphism rs17563 in BMP4 Gene in Brazilian Population Suggests Protection for Nonsyndromic Cleft Lip and Palate

Tânia Kawasaki Araújo,¹ Milena Simioni,¹ Têmis Maria Félix,² Liliane Todeschini de Souza,² Marshall Ítalo Barros Fontes,³ Isabella Lopes Monlleó,3,4 Josiane Souza,⁵ Agnes Cristina Fett-Conte,⁶ Rodrigo Secolin,¹ Iscia Lopes-Cendes,¹ Cláudia Vianna Maurer-Morelli,¹ and Vera Lúcia Gil-da-Silva-Lopes¹

¹Department of Medical Genetics, Faculty of Medical Sciences, University of Campinas (UNICAMP), 13083-887 Campinas, SP, Brazil²Medical Genetics Service, Hospital de Clinicas de Porto Alegre, 90035-903 Porto Alegre, RS, Brazil³Medical Genetics Sector, State University of Alagoas, 57010-300 Maceio, AL, Brazil⁴Clinical Genetics Service, University Hospital, Federal University of Alagoas, 57072-970 Maceio, AL, Brazil⁵Centro de Atendimento Integral ao Fissurado Labiopalatal (CAIF), 81050-00 Curitiba, PR, Brazil⁶Department of Molecular Biology, Medicine School of São José do Rio Preto (FAMERP FUNFARME), 15090-000 São José do Rio Preto, SP, Brazil

Received 27 July 2012; Accepted 3 November 2012

Academic Editor: Renato Da Silva Freitas

Copyright © 2012 Tânia Kawasaki Araújo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Cleft lip with or without palate (CL±P) is common congenital anomalies in humans. Experimental evidence has demonstrated that bone morphogenetic protein 4 gene (Bmp4) is involved in the etiology of CL±P in animal models. The nonsynonymous polymorphism rs17563 T>C (p.V152A) in the BMP4 gene has been associated to the risk of nonsyndromic CL±P in Chinese population and microforms from different ethnic backgrounds. The aim of this study was to investigate the role of BMP4 gene in CL±P in Brazilian sample using genetic association approach. Our sample was composed by 123 patients with nonsyndromic CL±P and 246 controls, in which absence of CL±P was confirmed in 3 generations. The rs17563 polymorphism was genotyped by PCR-RFLP technique. Logistic regression was performed to evaluate allele and genotype association. Our data showed statistical power to detect association (86.83%) in this sample. Logistic regression results showed significant association between C allele and CL±P (, OR , and 95% CI = 0.25–0.65), as well as CC genotype and CL±P (, OR , and 95% CI = 0.19–0.66). So, there is a strong association between nonsyndromic CL±P and BMP4 rs17563 polymorphism in our sample and the C allele had a protective effect against the occurrence of nonsyndromic CL±P.

1. Introduction

Nonsyndromic cleft lip with or without palate (CLP) is a congenital defect with multifatorial transmission. It is a birth defect that occurs with a prevalence of 1 in 500 to 1 in 2500 live births in different populations, varying with geographic location, ethnic group, and socioeconomic conditions [1–3]. Populations of either American native descent or Asian descent have the highest birth prevalence [4]. In Latin America, according to the ECLAMC (Latin American Collaborative Study of Congenital Malformations), the incidence of this malformation is 1 in 850 births [5]. CLP has a complex etiology that includes strong genetic and environmental factors.

Clinically, CLP are classified as nonsyndromic (isolated) or syndromic based on the presence of other congenital anomalies or mental disabilities [1, 6]. Around 400 recognizable syndromes that have CLP as part of the phenotype are described [7, 8], however, 70% of CLP cases are nonsyndromic. Approximately 12–25% of the genetic variations associated with nonsyndromic CLP have been identified [1, 6]. Although genetic studies have identified a number of candidate genes and chromosomal regions associated with CLP, findings from different studies have been inconsistent [9, 10].

Candidate genes for CLP have been proposed as result of studies using knockout mice models. One interesting candidate gene is the Bone morphogenetic protein 4 (BMP4; MIM#112262). Experimental evidences have postulated that Bmp4 is involved in the etiology of CLP in mice. Bilateral cleft lip was presented in all knockout mice for Bmp4 at 12 days after conception, however after 14.2 days only 22% of the embryos showed cleft lip. Moreover, Bmp signaling is also required for cell proliferation at earlier stages, in the maxillary process mesenchyme [11].

BMP4 gene, located in 14q22-q23, is a member of the BMP family and a member of the transforming growth factor-b (TGF-b) of secretory signaling molecules that play essential roles in embryonic development [12]. This gene is a well-characterized mammalian growth factor. Expression of the Bmp4 gene occurs in a wide spectrum of normal tissue in mouse embryos, as cardiac [11], limb [13], and face tissue [11, 14].

In addition, possible association between BMP4 gene and nonsyndromic CL±P in humans was firstly suggested in a case-control study performed with Chinese children in 2008 [15]. Authors analyzed the single-nucleotide polymorphism (SNP) rs17563 (p.Val152Ala) and found an increased risk of nonsyndromic CL±P among carriers of the C allele.

A missense and nonsense mutations were found in BMP4 gene in a sample of patients from different genetic background (Mongolia, Philippines, United States of America, Colombia, Guatemala, and Europe) diagnosed with subepithelial, microform, and overt cleft lip. These mutations were not found in a control population [16].

Furthermore, Suazo et al. [17] reported a mutation screening analysis of BMP4 in a sample of 150 Chilean nonsyndromic CLP case-parent trios. Significant deviations from expected transmissions were observed for haplotypes conformed by rs1957860 and rs762642. These polymorphisms delimitate a genomic region where a promoter and an enhancer of BMP4 are located [18]. Then, Suazo et al. [19] searched for nonsyndromic CLP risk variants within the two BMP4 promoters by direct sequencing in a 167 Chilean nonsyndromic CLP cases and 336 controls. They found three novel variants considered as cleft risk factors, as they are not present in controls.

Lately, Chen et al. [20] provided further evidence of association between BMP4 gene and nonsyndromic CLP. They tested for possible association between markers in and around the BMP4 gene and nonsyndromic CLP in Asian and Maryland trios. Nominal significant evidence of linkage and association was observed for three SNPs (rs10130587, rs2738265, and rs2761887) in 221 Asian trios and for one SNP (rs762642) in 76 Maryland trios.

Brazilian population presents a heterogeneous ethnic origin which includes Native Americans, Europeans, Africans, and some degree of Asians. Relative contribution of each one of these ethnic groups is variable according to the geographic region of the country. This characteristic gives a unique and richly blended genetic and cultural background to the country [21–25].

Therefore, since Brazilian population presents a diverse genetic background, the purpose of this study was to investigate the role of the rs17563 polymorphism of BMP4 gene in CLP in a sample of Brazilian population.

2. Materials and Methods

2.1. Study Population

The sample consisted of 369 individuals enrolled in Craniofacial Brazil Project. All cases were evaluated by clinical geneticist using the same clinical protocol, including family history. Subjects were recruited between 2009 and 2011 from five centers of research involved in Craniofacial Brazil Project located in six cities from three different regions of Brazil (Southeast region: Campinas and São José do Rio Preto, South region: Porto Alegre and Curitiba, Northeast region: Maceió). Case group was compound of 123 individuals (62 males and 61 females) with nonsyndromic CLP (98 cleft lip and palate and 25 cleft lip only). Control group were recruited from Southeast region by universitary e-mail included 246 phenotypically normal unrelated individuals (102 males and 144 females) without family history of oral cleft in three generations and no Japanese or Chinese ancestry. Their family history was collected by a participant researcher.

The study was approved by the local Institutional Review Boards and the National Research Ethics Committee (CEP 059/2008). All participants provided informed consent. The informed consent for patients under 18 years old was obtained from their parents.

2.2. Molecular Analysis

Genomic DNA was extracted from peripheral blood lymphocytes according to the method described by De Araujo et al. [26].

SNP genotyping was performed using a polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay. The PCR primers were designed based on the SNP flanking sequence described in the Ensembl Genome Browser (http://www.ensembl.org/index.html).

Primer sequences for PCR reactions were 5′AGTTTGGCTGCTTCTCCC3′ (forward) and 5′AGTTTGGCTGCTTCTCCC3′ (reverse). The PCR reaction were performed in a total volume of 15 L containing 1.5 L KCl₂ (10X), 0.6 L MgCl₂ (25 M), 1.5 L dNTP (2 M), 0.35 L of each primer primer sense (5pmoles), 0.5 L DMSO₄, 0.2 L Taq DNA Polimerase (5U) (New England BioLabs, Ipswich, MA, USA), and 1 L DNA (200 ng/L) e 9 L H₂O. The PCR cycle conditions consisted of an initial denaturation step at 94°C for 5 min followed by 35 cycles of 1 min at 94°C, 1 min at 56°C, 1 min at 72°C, and a final elongation at 72°C for 7 min.

Adequate restriction endonuclease for SNP was selected using the Gene Runner software (Hastings Software Inc., Hastings, NY, USA; http://www.generunner.net/). Primers, restriction enzyme, and the length of digested fragments for the SNP are listed in Table 1.

Table 1: Primers, restriction enzyme, and fragment lengths for rs17563 of BMP4 gene.

SNP
Primer
Annealing temperature
Restriction enzyme
Fragment lengths

rs17563
Forward 5′CACCATTCATTGCCCAAC3′
56°C
HphI
T Allele: 195pb + 229pb

Reverse 5′AGTTTGGCTGCTTCTCCC3′
C Allele: 424pb

The digestion reaction were performed in a total volume of 20 L containing 15 L PCR, 0.5 L enzyme, and 2.0 L buffer e 2.5 L H₂O. The PCR products were digested overnight at 37°C with HphI enzyme (New England BioLabs, Ipswich, MA, USA). PCR products were electrophoresed through 1% agarose gel in the presence of ethidium bromide and visualized by fluorescence in UV light; the digested PCR products were resolved on 12% polyacrylamide gels and visualized by a silver staining protocol.

2.3. SNP Association Analysis

Gender distribution was evaluated using Fisher’s exact test in statistical environment [27]. Minor allele frequency (MAF) and Hardy-Weinberg equilibrium were estimated using Haploview software [28] (http://www.broadinstitute.org/haploview). Allele and genotype association analysis were performed by logistic regression using statistical environment [27]. We evaluated posthoc statistical power of our sample by GPOWER software [29], using the following parameters: logistic regression test; two-tail; ; statistical significance level ; .

2.4. In Silico Test

In silico tests were performed to analyze whether the substitution c.538T>C in the BMP4 gene (rs17563) modifies the structure or function of the encoded protein. The sequence of amino acids used in the tests was obtained from NCBI database (http://www.ncbi.nlm.nih.gov/).

Three in silico tests were used: Grantham scale [30], PolyPhen [31] (http://genetics.bwh.harvard.edu/pph/), and SNAP [32] (http://rostlab.org/services/snap/) for the analysis of rs17563.

3. Results

The gender distribution was similar between groups (Fisher’s exact test ). MAF was 0.478 and SNP rs17563 was in Hardy-Weinberg equilibrium () in the sample. There were significant differences in the genotypes and allele frequencies of the rs17563 between nonsyndromic CLP and control groups.

Logistic regression revealed association between C allele and CL±P (, , and 95% ). Genotype logistic regression analysis showed rs17563 CC genotype was associated with decreased of CLP susceptibility (, , and 95% –0.66) (Table 2). The sample analyzed showed statistical power to detect genetic association (86.83%).

Table 2: The genotype and allele distribution of rs17563 in nonsyndromic CL±P patients and controls and result of logistic regression.

SNP rs17563
Controls  (%)
Cases  (%)

OR
95% CI

Genotype
TT¹
52 (21,1)
49 (39,8)
—
—
1.00
—

TC
130 (52,8)
53 (43,1)
10.645
0.00110
0.43
0.26–0.71

CC
64 (26,1)
21 (17,1)
11.217
0.00081
0.35
0.19–0.66

Allele
T¹
234 (47,6)
151 (61,4)
—
—
1.00
—

C
258 (52,4)
95 (38,6)
13.992
0.00018
0.40
0.25–0.65

1Reference factor.

The results obtained of in silico tests were similar for the three. Grantham’s scale showed that SNP was considered “moderately conservative” (value 64). The PolyPhen algorithm aligned sequences “wild” and “altered” to 75 homologous sequences and estimated the change as “benign” (score 0416). In addition, SNAP demonstrated that p.Val152Ala is neutral (reliability index = 3, expected accuracy = 78%).

4. Discussion

Mutant mice for Bmp4 gene with CLP phenotype suggest a possible role for this gene in lip and palate development [11]. Moreover, have been reported few studies of association between BMP4 gene and CLP in population.

Lin et al. [15] found significant differences in genotype distribution and allele polymorphism rs17563 of BMP4 gene between 184 nonsyndromic CLP cases and 205 controls. In a study with 1614 individuals from different countries (Mongolia, Philippines, United States of America, Colombia, Guatemala, and Europe), missense and nonsense mutations in BMP4 gene were detected in 5 of 968 individuals with overt CLP, 1 of 30 cases with microforms and 2 of 87 cases with orbicular oris muscle. They also found a borderline difference in the frequency of SNP rs17563 in cases compared with controls [16]. Thus, the results of this study support the hypothesis that this polymorphism might be clinically important in the genesis of CLP.

Recently, Chen et al. [20] tested for possible association between markers in BMP4 gene and nonsyndromic CLP in 297 Asian and Maryland trios. Their results does not support the evidence of linkage and association for the SNP rs17563.

However, in the present study the polymorphism rs17563 in BMP4 gene was associated with decreased susceptibility for nonsyndromic CLP (, , and 95% –0.66). This result was the opposite found by Lin et al. [15] and Suzuki et al. [16]. The difference could be attributed to the clinical selection of individuals, the well-defined control group, and the statistics power herein used. As well, this study obtained a strong association between nonsyndromic CLP and BMP4 rs17563 polymorphism and suggested that C allele had a protective role in the occurrence of nonsyndromic CLP in the Brazilian sample.

The rs17563 polymorphism replaces the amino acid valine by alanine at position 152 of the protein. It would seem that the structure and function of the protein are not significantly affected by the substitution due to physicochemical similarities between the involved amino acids. In our study, results of the in silico test using Grantham’s scale, PolyPhen, and SNAP were similar and showed the moderately conservative, benign, and neutral character of the rs17563 polymorphism. These results corroborate the hypothesis of a protective effect of this polymorphism.

Although cleft lip with cleft palate and cleft lip only are usually grouped into the same category, Dixon et al. [33] suggested that they could be etiologically distinct. Thus, further genetic association studies designed to examine these two conditions separately would be helpful in elucidating the role BMP4 rs17563 in orofacial clefts.

5. Conclusions

This study found a strong association between nonsyndromic CLP and BMP4 rs17563 polymorphism and a possible protective effect of the C allele against nonsyndromic CLP in a Brazilian sample. In this regard it could be speculated that rs17563 polymorphism plays different roles in oral cleft pathogenesis in specific populations. It would be interesting to replicate this study in individuals with similar and different ethnicity. As well, this association reinforces the importance of this gene as candidate for oral cleft.

Acknowledgments

The authors would like to thank the subjects and their families for cooperation in this study and Renata Barbosa for technical assistance. This study was supported by FAPESP (Grant no. 2008/10596-0) and CNPq (Grants nos. 01300.000433/2008-00, 502438/2010-0, and 136414/2009-5).

References

B. C. Schutte and J. C. Murray, “The many faces and factors of orofacial clefts,” Human Molecular Genetics, vol. 8, no. 10, pp. 1853–1859, 1999.

P. L. Bender, “Genetics of cleft lip and palate,” Journal of Pediatric Nursing, vol. 15, no. 4, pp. 242–249, 2000.

J. C. Murray, “Gene/environment causes of cleft lip and/or palate,” Clinical Genetics, vol. 61, no. 4, pp. 248–256, 2002.

M. B. Forrester and R. D. Merz, “Descriptive epidemiology of oral clefts in a multiethnic population, Hawaii, 1986–2000,” Cleft Palate-Craniofacial Journal, vol. 41, no. 6, pp. 622–628, 2004.

“ECLAMC (Latin American Collaborative Study of Congenital Malformations) final document. Angra dos Reis: ECLAMC VI,” 2004, (Electronic version: February 2005, Rio de Janeiro)

M. T. Cobourne, “The complex genetics of cleft lip and palate,” European Journal of Orthodontics, vol. 26, no. 1, pp. 7–16, 2004.

A. C. Lidral and L. M. Moreno, “Progress toward discerning the genetics of cleft lip,” Current Opinion in Pediatrics, vol. 17, no. 6, pp. 731–739, 2005.

Å Sivertsen, A. J. Wilcox, R. Skjærven et al., “Familial risk of oral clefts by morphological type and severity: population based cohort study of first degree relatives,” British Medical Journal, vol. 336, no. 7641, pp. 432–434, 2008.

A. R. Vieira, J. R. Avila, S. Daack-Hirsch et al., “Medical sequencing of candidate genes for nonsyndromic cleft lip and palate.,” PLoS Genetics, vol. 1, no. 6, article e64, 2005.

A. Mostowska, K. K. Hozyasz, P. Wojcicki, B. Biedziak, P. Paradowska, and P. P. Jagodzinski, “Association between genetic variants of reported candidate genes or regions and risk of cleft lip with or without cleft palate in the polish population,” Birth Defects Research A, vol. 88, no. 7, pp. 538–545, 2010.

W. Liu, X. Sun, A. Braut et al., “Distinct functions for Bmp signaling in lip and palate fusion in mice,” Development, vol. 132, no. 6, pp. 1453–1461, 2005.

L. M. Leong and P. M. Brickell, “Bone morphogenetic protein-4,” International Journal of Biochemistry and Cell Biology, vol. 28, no. 12, pp. 1293–1296, 1996.

A. Bandyopadhyay, K. Tsuji, K. Cox, B. D. Harfe, V. Rosen, and C. J. Tabin, “Genetic analysis of the roles of BMP2, BMP4, and BMP7 in limb patterning and skeletogenesis.,” PLoS Genetics, vol. 2, no. 12, article e216, 2006.

S. G. Gong and C. Guo, “Bmp4 gene is expressed at the putative site of fusion in the midfacial region,” Differentiation, vol. 71, no. 3, pp. 228–236, 2003.

J. Y. Lin, Y. J. Chen, Y. L. Huang et al., “Association of bone morphogenetic protein 4 gene polymorphisms with nonsyndromic cleft lip with or without cleft palate in Chinese children,” DNA and Cell Biology, vol. 27, no. 11, pp. 601–605, 2008.

S. Suzuki, M. L. Marazita, M. E. Cooper et al., “Mutations in BMP4 are associated with Subepithelial, Microform, and Overt Cleft Lip,” American Journal of Human Genetics, vol. 84, no. 3, pp. 406–411, 2009.

J. Suazo, J. L. Santos, L. Jara, and R. Blanco, “Association between bone morphogenetic protein 4 gene polymorphisms with nonsyndromic cleft lip with or without cleft palate in a chilean population,” DNA and Cell Biology, vol. 29, no. 2, pp. 59–64, 2010.

S. Kawai and T. Sugiura, “Characterization of human bone morphogenetic protein (BMP)-4 and -7 gene promoters: activation of BMP promoters by Gli, a sonic hedgehog mediator,” Bone, vol. 29, no. 1, pp. 54–61, 2001.

J. Suazo, J. C. Tapia, J. L. Santos, V. G. Castro, A. Colombo, and R. Blanco, “Risk variants in BMP4 promoters for nonsyndromic cleft lip/palate in a Chilean population,” BMC Medical Genetics, vol. 12, article 163, 2011.

Q. Chen, H. Wang, J. B. Hetmanski et al., “BMP4 was associated with NSCL/P in an Asian population,” PLOS ONE, vol. 7, no. 4, article e35347, 2012.

J. Pellegrino Jr., L. Castilho, M. Rios, and C. A. De Souza, “Blood group genotyping in a population of highly diverse ancestry,” Journal of Clinical Laboratory Analysis, vol. 15, no. 1, pp. 8–13, 2001.

S. M. Callegari-Jacques, D. Grattapaglia, F. M. Salzano et al., “Historical genetics: spatiotemporal analysis of the formation of the Brazilian population,” American Journal of Human Biology, vol. 15, no. 6, pp. 824–834, 2003.

L. B. Ferreira, C. T. Mendes-Junior, C. E. V. Wiezel, M. R. Luizon, and A. L. Simões, “Genomic ancestry of a sample population from the State of São Paulo, Brazil,” American Journal of Human Biology, vol. 18, no. 5, pp. 702–705, 2006.

M. R. Luizon, C. T. Mendes, S. F. De Oliveira, and A. L. Simões, “Ancestry informative markers in Amerindians from Brazilian Amazon,” American Journal of Human Biology, vol. 20, no. 1, pp. 86–90, 2008.

Y. C. N. Muniz, L. B. Ferreira, C. T. Mendes-Junior, C. E. V. Wiezel, and A. L. Simões, “Genomic ancestry in urban Afro-Brazilians,” Annals of Human Biology, vol. 35, no. 1, pp. 104–111, 2008.

M. De-Araujo, M. R. Sanches, L. A. Suzuki, G. Guerra, S. B. Farah, and M. P. De-Mello, “Molecular analysis of CYP21 and C4 genes in Brazilian families with the classical form of steroid 21-hydroxylase deficiency,” Brazilian Journal of Medical and Biological Research, vol. 29, no. 1, pp. 1–13, 1996.

R Development Core Team, R: A Language and Environment For Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2.12.2 edition, 2010.

J. C. Barrett, B. Fry, J. Maller, and M. J. Daly, “Haploview: analysis and visualization of LD and haplotype maps,” Bioinformatics, vol. 21, no. 2, pp. 263–265, 2005.

F. Faul, E. Erdfelder, A. G. Lang, and A. Buchner, “G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences,” Behavior Research Methods, vol. 39, no. 2, pp. 175–191, 2007.

R. Grantham, “Amino acid difference formula to help explain protein evolution,” Science, vol. 185, no. 4154, pp. 862–864, 1974.

V. Ramensky, P. Bork, and S. Sunyaev, “Human non-synonymous SNPs: server and survey,” Nucleic Acids Research, vol. 30, no. 17, pp. 3894–3900, 2002.

J. Thusberg, A. Olatubosun, and M. Vihinen, “Performance of mutation pathogenicity prediction methods on missense variants,” Human Mutation, vol. 32, no. 4, pp. 358–368, 2011.

M. J. Dixon, M. L. Marazita, T. H. Beaty, and J. C. Murray, “Cleft lip and palate: understanding genetic and environmental influences,” Nature Reviews Genetics, vol. 12, no. 3, pp. 167–178, 2011.