Print version ISSN 1415-4757
Genet. Mol. Biol. vol.31 no.3 São Paulo 2008
HUMAN AND MEDICAL GENETICS
José SuazoI; José Luis SantosII; Lilian JaraI; Rafael BlancoI
IHuman Genetics Program, Institute of Biomedical Sciences, School of Medicine, University of Chile, Santiago, Chile
IIDepartment of Epidemiology, Institute for Technological Nutrition of Foods, University of Chile, Santiago, Chile
Nonsyndromic cleft lip/palate (NSCLP) is a congenital malformation with features of a complex genetic trait. Several studies have reported positive association and linkage between NSCLP and microsatellite markers in the 4q28-4q33 region particularly with the D4S192 (4q31) marker. We hypothesized that the candidate genes SMAD1 and HHIP (4q31) could be involved in the etiology of NSCLP based on previous positive linkage results and their important role in maxillofacial development. We evaluated the possible association between microsatellite markers located at less than 1 cM from these genes and NSCLP using a sample of 58 Chilean case-parent trios. Microsatellite markers were analyzed using the polymerase chain reaction (PCR) with fluorescent labeled primers. Electrophoresis of the PCR products was performed on a laser-fluorescent automatic DNA sequencer. The extended transmission disequilibrium test (ETDT) was used to analyze allelic transmissions from the parents to their affected progeny. No significant association due to linkage disequilibrium was detected between both markers and NSCLP.
Key words: candidate genes, case-parents trio design, linkage disequilibrium, nonsyndromic cleft lip/palate.
Nonsyndromic cleft lip/palate (NSCLP), is a common birth defect with features of a genetically complex trait. Attempts to localize NSCLP loci in the human genome have generated considerable but sometimes discordant information (Murray, 1995). Thus, the nature of the genetic contribution to the etiology of NSCLP is still being studied and remains unresolved. Several candidate genes and loci mapping in various chromosome regions have been claimed to be involved in cleft determination by parametric and non-parametric linkage and association analysis (Carinci et al. 2000; Murray, 2002). In this context, some studies have reported positive and negative findings between NSCLP and the 4q28-4q33 region. Beiraghi et al. (1994) reported the first study in a single five-generation family where two markers (D4S175 and D4S192) were informative for linkage. Mitchell et al. (1995) and Paredes et al. (1999) found significant association between D4S192 and NSCLP in case-control studies. Additionally, Marazita et al. (2002) and Wyszynski et al. (2003) in genome wide scan studies also proposed the involvement of a cleft susceptibility locus in 4q region. Notwithstanding, Blanton et al. (1996), Prescott et al. (2000), and Wong et al. (2000) report exclusion of linkage between NSCLP and this chromosomal region. Recently, our group screened the 4q24-q33 region using five short tandem repeat (STR) microsatellite markers but found no evidence for linkage and association between them and NSCLP in the Chilean population using a case-parents trio design (Blanco et al., 2005).
In a previous study in the Chilean population (Paredes et al., 1999), we reported a positive association between D4S192 and NSCLP. Based on this finding, we searched for possible candidate genes within the 4q31 region using as a reference point the aforementioned STR. A search in genomic databases, showed that there were two possible candidate genes both centromeric and telomeric to D4S192. Moreover, experimental evidence has shown that these two genes are involved in craniofacial development. The SMAD1 gene is a homolog of the drosophila gene denominated mothers against decapentaplegic 1 (MAD1) and encodes a protein which participates in a signaling cascade triggered by transforming growth factor betas (TGFbs) and bone morphogenetic proteins (BMPs) (Cohen, 2003). Nonaka et al. (1999) reported that SMAD1 serves as a point of convergence for the integration of two different growth factor signaling pathways, BMP4 and epidermal growth factor (EGF), during chondrogenesis. Hatakeyama et al. (2003) showed that in chondroprogenitor cells, BMP stimulates differentiation through mechanisms mediated by SMAD1. On the other hand, the hedgehog-interacting protein (HHIP) gene encodes a protein which is a member of the sonic hedgehog (SHH) signaling pathway (Jiang et al., 2006). Also Cobourne and Sharpe (2002) have postulated that a role for HHIP may be to allow differential responses to the SHH signal within different regions of odontogenic-specific cells in regions of the mandibular process. Rice et al. (2006) studying the expression patterns of hedgehog signaling pathway members during mouse palatal development reported that at 13 dpc, HHIP was expressed at a low level throughout the palatal mesenchyme and that the highest levels of transcripts localize adjacent to the thickened palatal oral epithelium in the anterior sections. At 14.5 dpc HHIP was expressed in the developing palatine bones.
The purpose of the study described in this paper was to use a case-parents trio design to assess a Chilean population for a possible association due to linkage disequilibrium between STR D4S1586 located at 0.3 cM telomeric from SMAD1 and STR D4S2998 located at 0.1 cM centromeric from HHIP and NSCLP. These markers had not been included in previous association and linkage studies of Chilean populations (Paredes et al.,1999; Blanco et al. 2005).
Our sample consisted of 58 case-parents trios selected because they had an offspring affected with NSCLP. The unrelated NSCLP probands were identified and interviewed during the course of clinical examinations at the Cleft Lip/Palate Clinic of the School of Dentistry of the University of Chile, at the Dr. Alfredo Gantz Foundation located in the city of Santiago, Chile. One of the authors (RB) identified the affected individuals in the course of examinations conducted between 2004 and 2005. In-depth interviews of at least three family members were conducted to provide detailed information for pedigree construction. All the families enrolled were of Chilean ancestry including subjects presenting NSCLP as the unique familiar disease. Families using clefting drugs, such as phenytoin, warfarin and ethanol were excluded from the study. The pedigree history corresponds to individuals belonging to low to middle low socioeconomic strata given the genetic composition of the Chilean population which presents a relationship between ethnicity, Amerindian admixture, genetic markers, socioeconomic strata and prevalence of NSCLP (Valenzuela, 1988; Palomino et al., 1997). The Institutional Review Board of the School of Medicine of the University of Chile approved the study and all participants gave their informed consent. Genomic DNA was extracted from peripheral blood cells (Poncz et al., 1982). The microsatellite markers were D4S1586 located at 0.3 cM telomeric from the SMAD1 gene and D4S2998 located at 0.1 cM centromeric from HHIP. The polymerase chain reaction (PCR) was carried out according to a standard amplification protocol (Suazo et al., 2004) using the primers described by Dib (1996). Annealing temperatures were 62 °C for D4S1586 and 52 °C for D4S2998. The amplification reaction was performed using a fluorescent-dye-labeled forward primer. The products were analyzed in an ABI PRISM 310 genetic analyzer (Applied Biosystems). The electrophoretic results were processed by GENESCAN 3.1.2 software, and allele assignation was carried out using Genotyper software, version 2.5. The Extended Transmission Disequilibrium Test (ETDT) for multiple alleles was carried out to assess the differential pattern of excess transmission of alleles from heterozygous parents to diseased children (Sham and Curtis, 1995). By sampling case-parent trios through an affected child, the association between alleles of genetic markers and the disease would cause transmission to appear different from the expected probability of 0.5. When cases are unrelated probands, ETDT represents a valid test of association even if population stratification is present. Given the relatively reduced sample size of our study and the low a priori statistical power to detect weak associations we also computed p-values based on simulations (10,000 per marker) made with the Monte Carlo ETDT program (MCETDT) (Zhao et al., 1999), which avoids the problems associated with p-values based on chi-square distributions applied to sparse transmission tables. Additionally, individual transmission of each allele with respect to the rest of the alleles was evaluated (this procedure would imply applying a correction for multiple comparisons). Exact p-values were evaluated to assess significance of individual alleles (Cleves et al., 1997) using a method implemented in the statistical package STATA 8.2 (Stata Statistical Software, 2004).
Multiple-allele ETDT analyses for marker-disease association based on case-parents trios yielded a p-value of 0.13 for the D4S1586 marker and 0.10 for the D4S2998 marker. Individual allele variants also showed no preferential transmissions. The Monte-Carlo ETDT yielded p-values of 0.17 for D4S1586 and 0.20 for the D4S2998 (Tables 1 and 2).
The genes SMAD1 and HHIP are located in 4q31.21 at a distance of 3.7 and 2.6 cM respectively from D4S192 (Figure 1). Several authors have reported positive findings in this chromosomal region for NSCLP. The results of experimental evidence cited in the preceding paragraphs support the role of these genes in craniofacial development but, nevertheless, the results of our study did not show positive evidence of the involvement of these genes in NSCLP in the Chilean population studied. It must also be taken into consideration the low a priori statistical power to detect weak associations and therefore we cannot exclude the possibility that these genes may have a role in NSCLP either as low susceptibility or modifier genes.
We are grateful to the family members for their participation in this study and also to the staff members of the Dr. Alfredo Gantz Foundation; the Cleft Lip Palate Clinic; School of Dentitstry, University of Chile; and the Cleft Children Center for their help in recruiting the families for this study. This work was supported by a grant DI SAL 03/13-2 from the University of Chile.
Alvarez Martínez C, Binato R, González S, Pereira M, Robert B and Abdelhay E (2001) Characterization of a smad motif similar to Drosophila mad in the mouse Msx1 promoter. Biochem Biophys Res Commun 291:655-662. [ Links ]
Beiraghi S, Foroud T, Diouhy S, Bixler D, Conneally PM, Delozier-Blanchet D and Hodes MS (1994) Possible localization of a major gene for cleft lip and palate to 4q. Clin Genet 46:255-256. [ Links ]
Blanco R, Suazo J, Santos JL, Carreño H, Palomino H and Jara L (2005) No evidence for linkage and association between 4q microsatellite markers and nonsyndromic cleft lip and palate in Chilean case-parents trios. Cleft Palate Craniofac J 42:267-271. [ Links ]
Blanton SH, Crowder E, Malcolm S, Winter R, Gasser DL, Stal S, Mulliken J and Hecht JJ (1996) Exclusion of linkage between cleft lip with or without cleft palate and markers on chromosome 4 and 6. Am J Hum Genet 58:241-243. [ Links ]
Carinci F, Pezzetti F, Scapoli L, Martinelli M, Carinci P and Tognon M (2000) Genetics of nonsyndromic cleft lip and palate: A review of international studies and data regarding the Italian population. Cleft Palate Craniofac J 37:33-40. [ Links ]
Cleves MA, Olson JM and Jacobs KB (1997) Exact transmission disequilibrium test with multiallelic markers. Genet Epidemiol 14:337-347. [ Links ]
Cobourne M and Sharpe P (2002) Expression and regulation of hedgehog-interacting protein during early tooth development. Connect Tissue Res 43:143-147. [ Links ]
Cohen M (2003) TGFβ/smad signaling system and its pathologic correlates. Am J Med Genet 116A:1-10. [ Links ]
Dib C (1996) A comprehensive genetic map of the human genome based on 5,264 microsatellites. Nature 380:152-154. [ Links ]
Hatakeyama B, Nguyen J, Wang X, Nuckolls G and Shum L (2003) Smad signaling in mesenchymal and chondroprogenitor cells. J Bone Joint Surg 85-A Suppl 3:13-18. [ Links ]
Jiang R, Busch J and Lidral A (2006) Development of the upper lip: Morphogenetic and molecular mechanisms. Dev Dyn 235:1152-1166. [ Links ]
Marazita ML, Field L, Cooper M, Tobias R, Maher B, Peanchitlertkajorn S and Liu Y (2002) Genome scan for loci envolved in cleft lip with or without cleft palate, in Chinese multiplex families. Am J Hum Genet 71:349-364. [ Links ]
Mitchell LE, Healey SC and Chenevix-Trench G (1995) Evidence for an association between nonsyndromic cleft lip with or without cleft palate and a gene located on the long arm of chromosome 4. Am J Hum Genet 57:1130-1136. [ Links ]
Murray JC (1995) Face facts: Genes, environment and clefts. Am J Hum Genet 57:227-232. [ Links ]
Murray JC (2002) Gene/environment causes of cleft lip and/or palate. Clin Genet 61:248-256. [ Links ]
Nonaka Z, Shum L, Takahashi I, Takahashi K, Ikura T, Dashner R, Nucolls G and Slavkin H (1999) Convergence of the BMP and EGF signaling pathway on Samd1 in the regulation of chondrogenesis. Int J Ved Biol 43:795-807. [ Links ]
Palomino HM, Palomino H, Cauvi D, Barton SA and Chackraborty R (1997) Facial clefting and Amerindian admixture in populations of Chile. Am J Hum Biol 9:225-232. [ Links ]
Paredes M, Carreño H, Solá JA, Segú J, Palomino H and Blanco R (1999) Asociación entre el fenotipo fisura labiopalatina no sindrómica y marcadores de microsatélite ubicados en 4q. Rev Med Chile 127:1431-1438. [ Links ]
Poncz M, Solowiejcyk D, Harvel B, Mory Y, Schwartz E and Surrey S (1982) Construction of human gene libraries from small amounts of peripheral blood: Analysis of beta-like globin genes. Hemoglobin 6:27-36. [ Links ]
Prescott N, Lees M, Winter R and Malcolm S (2000) Identification of susceptibility loci for nonsyndromic cleft lip with or without cleft palate in a two stage genome scan of affected sib-pairs. Hum Genet 106:345-350. [ Links ]
Rice R, Connor E and Rice D (2006) Expression patterns of hedgehog signalling pathway members during mouse palate development. Gene Expr Pattern 6:206-212. [ Links ]
Sham PC and Curtis D (1995) An extended transmission disequilibrium test for multi-allele marker loci. Ann Hum Genet 59:323-336. [ Links ]
Suazo J, Santos JL, Carreño H, Jara L and Blanco R (2004) Linkage disequilibrium between MSX1 gene and nonsyndromic cleft lip palate in the Chilean population. J Dent Res 83:782-785. [ Links ]
Valenzuela C (1988) On sociogenetic clines. Ethol Sociobiol 9:259-269. [ Links ]
Wong FK, Hagberg C, Karsten A, Larson O, Gustavsson M, Huggare J, Larsson C, Teh BT and Linder-Aronson S (2000) Linkage analysis of candidate regions in Swedish nonsyndromic cleft li with or without cleft palate families. Cleft Palate Craniofac J 37:357-362. [ Links ]
Wyszynski D, Albach-Hejazi H, Aldirani M, Hammod M, Shkair H, Karam A, Alashkar J, Holmes T, Pugh E, Doheny K, et al. (2003) A genome-wide scan for loci predisposing to non-syndromic cleft lip with or without cleft palate in two large Syrian families. Am J Med Genet 123:140-147. [ Links ]
Zhao JH, Sham PC and Curtis D (1999) A program for the Monte Carlo evaluation of significance of the extended transmission disequilibrium test. Am J Hum Genet 64:1484-1485. [ Links ]
Send correspondence to:
Human Genetics Program
Institute of Biomedical Sciences, School of Medicine, University of Chile
Av. Independencia 1027, Casilla
70061, Santiago, Chile
Received: June 21, 2007; Accepted: November 22, 2007.
Associate Editor: Paulo A. Otto