|Year : 2017 | Volume
| Issue : 1 | Page : 9-14
Linkage evidence and methylation at 2q region in nonsyndromic cleft lip and/or palate of Malay population
Nurul Syazana Mohamad Shah1, Sarina Sulong2, Wan Azman Wan Sulaiman1, Ahmad Sukari Halim1
1 Reconstructive Science Unit, School of Medical Sciences, Health Campus, Universiti Sains , 16150 Kubang Kerian, Kelantan, Malaysia
2 Human Genome Centre, School of Medical Sciences, Health Campus, Universiti Sains , 16150 Kubang Kerian, Kelantan, Malaysia
|Date of Web Publication||2-May-2017|
Wan Azman Wan Sulaiman
Reconstructive Science Unit, School of Medical Sciences, Universiti Sains Malaysia, Health Campus, 16150 Kubang Kerian, Kelantan
Source of Support: None, Conflict of Interest: None
Introduction: Nonsyndromic cleft lip and/or palate (NSCLP) occurs as a result of multifactorial determinants, involving both genetic and environmental factors. Several candidate genes associated with NSCLP have been discovered through genetic approach, but there is paucity of studies focusing on epigenetic determinants in NSCLP. We are interested to reveal linkage evidence of SATB2 at 2q region in large-extended NSCLP families of Malay population and its methylation activity in causing cleft formation. Materials and Methods: Eight large-extended families were included in this study. Microarray analysis was carried out and genome-wide linkage was determined using GeneHunter Multipoint Linkage Analysis v2.1r5. SATB2 methylation was tested on 100 NSCLP patients by DNA sequencing. Results: Genome-wide linkage analysis has revealed significant nonparametric linkage score and suggestive logarithm of the odds (LOD) score at 2q region in family 50 and family 100. Genome-wide heterogeneity LOD score of 2.63 and α =0.122 were found in total families at 2q33.1-q35 region. Significant copy number loss (P < 0.05) in NSCLP family compared with the normal control supports the linkage evidence of SATB2 in those families with positive linkage. Epigenetic testing found SATB2 unmethylation at DNA promoter region. Discussion: Linkage evidence and significant low copy number of SATB2 in NSCLP family of Malay population confirmed that genetic factors play a major role in causing cleft defects. SATB2 unmethylation could not support the epigenetic occurrence in causing craniofacial deformities. Conclusions: Linkage evidence and significant low copy number of SATB2 in NSCLP family of Malay population confirmed that genetic factors play a major role in causing cleft defects. SATB2 unmethylation could not support the epigenetic occurrence in causing craniofacial deformities.
Keywords: Epigenetic, linkage, microarray, nonsyndromic cleft lip and/or palate
|How to cite this article:|
Shah NS, Sulong S, Wan Sulaiman WA, Halim AS. Linkage evidence and methylation at 2q region in nonsyndromic cleft lip and/or palate of Malay population. J Cleft Lip Palate Craniofac Anomal 2017;4:9-14
|How to cite this URL:|
Shah NS, Sulong S, Wan Sulaiman WA, Halim AS. Linkage evidence and methylation at 2q region in nonsyndromic cleft lip and/or palate of Malay population. J Cleft Lip Palate Craniofac Anomal [serial online] 2017 [cited 2017 Dec 15];4:9-14. Available from: http://www.jclpca.org/text.asp?2017/4/1/9/205413
| Introduction|| |
Orofacial clefts (OFCs) are congenital malformations that comprise a large fraction of human congenital disabilities that affect the lip and/or palate. Cleft lip and/or palate (CLP) are immediately recognizable disruptions of normal facial structure. Several studies have reported that Asian and Native North American descent populations with prevalence of 2/1000 births have highest birth prevalence compared to the other populations., Several candidate genes have been reported to play a role in OFC including MSX1, SATB2, FGF, and IRF6.,,, Main focus of this study is to discover evidence of linkage at chromosome 2q, particularly SATB2 in multiplex families of nonsyndromic CLP (NSCLP), and to investigate its methylation activity to be associated with craniofacial defects among the Malays.
Although SATB2 is known as one of the candidate genes to NSCLP, none of the genes seem to play a major role to this deformity due to population variation. The first linkage study was done by Beiraghi et al. that reported the possible gene in 4q region tested on a single five-generation family of cleft lip/palate (CL/P) but later reported to have no significant evidence for linkage in 56 families tested., Linkage studies revealed several candidate regions that have been studied in different populations.,,
As OFC occurs as a result of multifactorial determinants, epigenetic is one of the environmental effects that could also contribute to the cleft formation. In eukaryotes, epigenetic occurs by silencing the mechanism through reversible biochemical modifications of DNA and associated proteins without altering DNA sequence, whereas DNA methyltransferase catalyzes the reaction of cytosine to form 5-methylcytosine by an addition of a methyl group.,, Long-term gene silencing occurred due to DNA methylation; hence, it correlates the DNA methylation to the inhibition of gene expression. Since both genetic and environmental factors play a role in craniofacial deformities, we are interested to discover the functional role of SATB2 at 2q region from genetic and epigenetic findings among the Malays. Therefore, we could clarify SATB2 mutation presence as a result of either one or both factors.
| Materials and Methods|| |
Eight large-extended families with a history of NSCLP were included while 100 NSCLP patients were included for epigenetic analysis. All the subjects were recruited from Hospital Universiti Sains Malaysia and blood withdrawal was undertaken. Healthy controls with no family history of OFC were included as a normal control. This study was approved by the Research Ethics Committee (Human) of the Universiti Sains Malaysia, and written informed consent was obtained from all the subjects.
Inclusion and exclusion criteria
All family members either affected with cleft lip, CL/P, or cleft palate (CP), or healthy individuals were included in this study. Patients were first screened by specialists of plastic surgery. Any major abnormalities or syndromes detected during the screening were excluded. Patients with one minor anomaly such as low-set ears, hypertelorism, clinodactyly, and single palmar crease were included.
Genome-wide linkage analysis
Genomic DNA was isolated from blood samples using the QiaAmp DNA Blood Mini Kit (USA), following the manufacturer's instructions. Extracted DNA of the family members of NSCLP was genotyped using Illumina Infinium Human Linkage-24 Beadchip (Illumina, California, USA). Following DNA denaturation, the samples were hybridized, fragmented, and precipitated before the staining process. Data were analyzed using Illumina's GenomeStudio Genotyping Module. Markers with a minor allele frequency >5% and call rate higher than 95% were used in the analyses.
Confirmation by copy numzber variation assay real-time polymerase chain reaction
Copy number variation (CNV) was carried out on the two selected families (family 50 and family 100) showed suggestive or significant linkage, healthy control as a normal and commercialized human control as a calibrator. Specific SATB2 primer-probe pair was designed; Homo sapiens SATB2 National Center for Biotechnology Information (NCBI) location: Chr.2:200134223-200335989; cytoband: 2q33.1d (Hs07541174_cn). CNV was performed using TaqMan ® Genotyping Master Mix for absolute quantitation of copy number using real-time polymerase chain reaction (qPCR). RNase P was used as an endogenous control and NTC (reaction mixture without DNA template) as a negative control. qPCR was run using ABI 7500 system with the following parameter; hold at 95°C for 10 min, followed by 95°C for 15 seconds and 60°C for 1 min for 40 cycles.
Epigenetic status of SATB2 by sequencing
Bisulfite-treated DNA and polymerase chain reaction
Bisulfite conversion of DNA was done using Methyledge ™ Bisulfite Conversion System (Promega, USA), following the manufacturer's instructions. Specific primers were designed using MethPrimer program (www.urogene.org/methprimer/) to produce a set of primer pair to detect methylated and unmethylated DNA sequences for SATB2. PCR amplification was performed using GoTaq ® G2 Hot Start Green Master Mix (Promega, USA) kit with a total volume of 25 μl. PCR conditions for both methylated and unmethylated primers were as follows: initial denaturation at 95°C for 2 min, denaturation at 95°C for 15 seconds, annealing at 61°C for 1 min, extension at 72°C for 1 min for 40 cycles, and final extension 72°C for 5 min.
DNA sequencing analysis
Purified PCR products were sequenced using BigDye Terminator V3.1 Sequencing Standard Kit (Applied Biosystems, Foster City, USA) as recommended by the manufacturer. The cycle sequencing was carried out in a total volume of 20 μl, containing around 40–90 ng of purified PCR products, 2 μl BigDye reaction premix, 4 μl BigDye sequencing buffer, 1 μl primer (10 μM), and distilled water. The cycle sequencing was run using ABI GeneAmp 9700 PCR System with initial denaturation at 96°C for 1 min, denaturation at 96°C for 10 seconds, annealing at 50°C for 5 seconds for 25 cycles, and extension at 60°C for 4 min. DNA sequences were analyzed using the Bioedit Sequence Alignment Editor Software v5.0.9 and were blasted with reference sequence.
For linkage analysis, easyLinkage Plus v5.08 software was used to analyze a large-scale single nucleotide polymorphism data. It would check for Mendelian/non-Mendelian genotyping errors and Hardy–Weinberg equilibrium. GeneHunter Multipoint Linkage Analysis v2.1r5 was carried out on both parametric and nonparametric analysis with logarithm of the odds (LOD) and nonparametric linkage (NPL) score. The thresholds for genome-wide linkage analysis were categorized based on the criteria proposed by Lander and Lander and Kruglyak (1995). LOD score was indicated as if 1.9 ≤ LOD <3.3 is a suggestive linkage and LOD ≥3.3 is a significant linkage. The allele sharing NPL score was categorized into three stages; score in between 2.2 and 3.5 is a suggestive linkage, 3.6–5.3 is a significant linkage, and ≥5.4 is highly significant linkage. Data output were then imported to NCBI web-based followed by GeneDistiller database to sort the genes according to the information obtained from the linkage data.
Triplicates of two independent experiments were carried out, and data were expressed as mean ± standard deviation. Differences between the groups (NSCLP and normal) were evaluated using Mann–Whitney U-test, and P< 0.05 was considered statistically significant. These values were calculated using the SPSS Statistics (SPSS 22, SPSS, Inc., New York, USA).
| Results|| |
Linkage evidence in nonsyndromic cleft lip and/or palate families
Suggestive linkage peaks have been attained in family 100 with parametric LOD score of 2.12 and 2.18 at chromosomal region 2q31.2-q33.3, respectively [Table 1]. Genome-wide heterogeneity LOD (HLOD) score of 2.63 and α = 0.122 was also found in total families at similar 2q33.1-q35 region [Table 1]. The score for HLOD ≥1 and α ≤ 0.71 indicate that the proportion of families linked and showed the presence of genetic heterogeneity in this population. Nonparametric analysis attained suggestive and significant NPL score (2.49 and 3.69) at 2q34-q36.3 and 2q31.1-q35 in family 50 and family 100, respectively [Table 1].
|Table 1: Linkage intervals with a heterogeneity logarithm of the odds score .1 was detected in total families, logarithm of the odds score >2 at chromosomal region 2 detected in family 100 and nonparametric linkage score >2 in family 50 and family 100|
Click here to view
Copy number variation analysis
Copy number calculated (CNC) value for the calibrator was 2.06 and most normal controls (5 out of 6) had CNC value 2.49 ≥ CNC ≥1.50 [Figure 1]. From nine members of family 100, only one nonaffected member (100-4) showed CNC <1.50 (CNC = 1.29) while the other four nonaffected members including four affected members; 100-P1, 100-P2, 100-3, and 100-14 had 2.49 ≥ CNC ≥1.50. Similarly, all the family members (n = 8) in family 50 including the two affected individuals (50-9 and 50-12) had 2.49 ≥ CNC ≥1.50 [Figure 1].
|Figure 1: Bars indicated copy number calculated of SATB2 in family 100 (purple bars), family 50 (green bars), healthy control (orange bars), and calibrator (blue bar). Error bars represent the standard deviation of uncertainty in each sample tested|
Click here to view
Scatter plot data showed a distinct distribution of CNC between the target NSCLP families and normal controls [Figure 2]. The straight line depicted the normal copy number fell at two. SATB2 copy number for family 50 and family 100 was found scattered <2 meanwhile SATB2 copy number in normal controls were dispersed to >2. By comparing the CNC of SATB2 between family 50 (1.90 ± 0.12, P = 0.039) and family 100 (1.87 ± 0.24, P = 0.045) to the normal control group (2.18 ± 0.30), a significant decrease of SATB2 copy number was found in both families in comparison to normal, indicating that SATB2 copy loss was present among the NSCLP families.
|Figure 2: Scatter plot indicated copy number calculated value in family 50 and family 100 in comparison to the normal controls|
Click here to view
With DNA sequencing, target sequence that was aligned together with the reference sequence revealed high repetition of SATB2 unmethylation for all the target and normal controls. C to T (C > T) conversion was frequently observed on the sequence alignment in both groups as indicated in blue circles [Figure 3]. This sequence confirmed that no SATB2 methylation occurred.
|Figure 3: Methylation analysis on SATB2. Mapping of reference (upper) and target (lower) sequence in SATB2 for comparison. Unmethylation status was identified as repetitive C > T conversions detected, indicating in blue circles|
Click here to view
| Discussion|| |
Linkage evidence at 2q31.1-q35 region
A suggestive LOD score in family 100 and an HLOD score for total family (HLOD = 2.63, α = 0.124) was attained on chromosome 2q region. In addition, highest significant NPL score in family 100 and suggestive NPL score in family 50 were also attained in chromosomal region 2q31.1-q36.3. GeneDistiller database extracted several genes fell within this chromosomal region including ZNF533, FSIP2, HOXD4, SUMO1, MSTN, NRP2, ABCA12, and ERBB4. The most common one that plays a role in craniofacial development was SATB2. Although SATB2 has been widely reported to be associated with clefting, the results have not been consistent among populations. Therefore, it was our interest to study SATB2 since it has never been reported in Malay population linkage studies to date.
Several genome-wide linkage studies conducted using multiplex families from various ethnic origins showed evidence of linkage at chromosome 2q.,, The genome scan meta-analysis from 13 populations of NSCLP conducted by Marazita et al. yielded highly significant result in 2q32-q35 region that contains the gene for DNA-binding protein SATB2. In addition, linkage evidence and linkage heterogeneity have been attained in two regions of chromosome 2 and 2q region, respectively, in 10 multiplex families at three different regions. Nevertheless, a previous study on case-parent trios in three distinct populations has found SATB2 as a candidate gene for NSCLP. Therefore, positive linkage found in the 2q region in NSCLP families of Malay population has provided an evidence of gene susceptibility on the loci.
Copy number loss of SATB2
CNVs is able in altering gene expression or gene dosage and disturbing gene sequence that finally lead to genetic defects and phenotypic variation., As anticipated, absolute copy number loss of SATB2 was significantly detected among the affected families compared to the healthy controls. This finding plausibly explained the significant role of SATB2 in normal craniofacial development, whereas any loss of it may cause cleft defects. Previous finding has reported that SATB2 lied within the deleted region, with low copy number of SATB2 being detected in an NSCLP offspring rather than in unaffected parents and healthy controls.
SATB2 copy number loss detected among the NSCLP patients could also be caused by a de novo change. Previous studies had detected a de novo mutation through array comparative genomic hybridization method. A patient with CP had submicroscopic copy number alterations whereas de novo 4.5 Mb deletion on loci 2q33.1 region has been detected. However, here, we did not further investigate whether the copy number loss occurs due to a de novo process or inherited from an affected/unaffected parent. In addition, previous report on decreased SATB2 copy number caused increased of cell death at the developing jaw primordial including palate and hindered regional development and therefore caused the craniofacial defects. The loss of SATB2 copy number in the family of affected cleft increased the risk and susceptibility of SATB2 to cause OFC.
Epigenetic alteration works by triggered changes in gene expression or inherited changes of phenotype through DNA methylation without altering the DNA sequences. Aberrant DNA methylation of susceptible genes could result in gene silencing and spontaneously affect the normal biochemical pathway and palatal developmental process. It has been said that methylation often detected in a promoter region of genes. Meanwhile, two separate studies on human methylome and murine secondary palate methylome found different findings. Methylated regions of interest were found in intragenic regions so-called gene bodies which include exons, introns and 3'- and 5'-untranslated regions rather than in a promoter region., Interestingly, methylated promoter regions were found associated with decreased gene expression, but intragenic methylation correlated with increased gene transcription.
SATB2 is a homeobox transcription factor that is known to play a role in nonsyndromic CP. It was first reported whereas translocation in 2q32-q33 region disturbed the gene regulation. Two similar findings have been reported with regard to SATB2 mutation in both syndromic (patient with CP, osteoporosis, and mental retardation) and nonsyndromic CP cases., Several studies have reported deletion or balanced translocations at 2q32-q34 region with commonly shared features such as developmental delay, cleft/high palate, defects in jaw formation, feeding difficulties, and abnormal dentition, which may relate to NSCLP as well.,, Recently, there is no report regarding SATB2 methylation in association with NSCLP formation. SATB2 unmethylation in all the targets and normal controls found in our study confirm the epigenetic event did not occur and hence has no role in craniofacial development.
| Conclusion|| |
Finding more samples from families with strong NSCLP history and collecting the samples from all the affected and unaffected members were the limitations of our study. The limitation was obvious whereas the linkage evidence of SATB2 at 2q region was only attained in two families, indicating that different family had different causative gene. Significant linkage and low copy number of SATB2 in NSCLP family compared to the normal indicated that genetic factor is the causal factor for SATB2 in NSCLP formation. This study revealed SATB2 methylation did not affect craniofacial development and could not support the role of epigenetic process in association with OFC.
We are grateful to all the patients, parents, and families for their cooperation. We thank the staff of Reconstructive Science Unit, Universiti Sains Malaysia. The study was supported by Research University grant 1001/PPSP/812083. We also thank the scholarship SLAB/SLAI under Ministry of Education Malaysia and USM for supporting the student.
Financial support and sponsorship
The study was supported by Research University Grant (1001/PPSP/812083).
Conflicts of interest
There are no conflicts of interest.
| References|| |
Jugessur A, Murray JC. Orofacial clefting: Recent insights into a complex trait. Curr Opin Genet Dev 2005;15:270-8.
Dixon MJ, Marazita ML, Beaty TH, Murray JC. Cleft lip and palate: Understanding genetic and environmental influences. Nat Rev Genet 2011;12:167-78.
Cooper ME, Ratay JS, Marazita ML. Asian oral-facial cleft birth prevalence. Cleft Palate Craniofac J 2006;43:580-9.
Vieira AR, McHenry TG, Daack-Hirsch S, Murray JC, Marazita ML. A genome wide linkage scan for cleft lip and palate and dental anomalies. J Med Genet A 2008;146A:1406-13.
Lace B, Kempa I, Piekuse L, Grinfelde I, Klovins J, Pliss L, et al.
Association studies of candidate genes and cleft lip and palate taking into consideration geographical origin. Eur J Oral Sci 2011;119:413-7.
Zeiger JS, Hetmanski JB, Beaty TH, VanderKolk CA, Wyszynski DF, Bailey-Wilson JE, et al.
Evidence for linkage of nonsyndromic cleft lip with or without cleft palate to a region on chromosome 2. Eur J Hum Genet 2003;11:835-9.
Riley BM, Mansilla MA, Ma J, Daack-Hirsch S, Maher BS, Raffensperger LM, et al.
Impaired FGF signaling contributes to cleft lip and palate. Proc Natl Acad Sci U S A 2007;104:4512-7.
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:769-80.
Vieira AR, Avila JR, Daack-Hirsch S, Dragan E, Félix TM, Rahimov F, et al.
Medical sequencing of candidate genes for nonsyndromic cleft lip and palate. PLoS Genet 2005;1:e64.
Blanco R, Suazo J, Santos JL, Carreño H, Palomino H, Jara L. 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 2005;42:267-71.
Beiraghi S, Foroud T, Diouhy S, Bixler D, Conneally PM, Delozier-Blanchet D, et al.
Possible localization of a major gene for cleft lip and palate to 4q. Clin Genet 1994;46:255-6.
Yildirim Y, Kerem M, Köroglu Ç, Tolun A. A homozygous 237-kb deletion at 1p31 identified as the locus for midline cleft of the upper and lower lip in a consanguineous family. Eur J Hum Genet 2014;22:333-7.
Carinci F, Pezzetti F, Scapoli L, Martinelli M, Carinci P, Tognon M. Genetics of nonsyndromic cleft lip and palate: A review of international studies and data regarding the Italian population. Cleft Palate Craniofac J 2000;37:33-40.
Martinelli M, Scapoli L, Pezzetti F, Carinci F, Francioso F, Baciliero U, et al.
Linkage analysis of three candidate regions of chromosome 1 in nonsyndromic familial orofacial cleft. Ann Hum Genet 2001;65(Pt 5):465-71.
Jaenisch R, Bird A. Epigenetic regulation of gene expression: How the genome integrates intrinsic and environmental signals. Nat Genet 2003;33 Suppl: 245-54.
Ehrlich M, Wang RY. 5-Methylcytosine in eukaryotic DNA. Science 1981;212:1350-7.
He XJ, Chen T, Zhu JK. Regulation and function of DNA methylation in plants and animals. Cell Res 2011;21:442-65.
Mohamad Shah NS, Salahshourifar I, Sulong S, Wan Sulaiman WA, Halim AS. Discovery of candidate genes for nonsyndromic cleft lip palate through genome-wide linkage analysis of large extended families in the Malay population. BMC Genet 2016;17:39.
Lander ES, Schork NJ. Genetic dissection of complex traits. Science 1994;265:2037-48.
Beaty TH, Hetmanski JB, Fallin MD, Park JW, Sull JW, McIntosh I, et al.
Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations. Hum Genet 2006;120:501-18.
Marazita ML, Murray JC, Lidral AC, Arcos-Burgos M, Cooper ME, Goldstein T, et al.
Meta-analysis of 13 genome scans reveals multiple cleft lip/palate genes with novel loci on 9q21 and 2q32-35. Am J Hum Genet 2004;75:161-73.
Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, et al.
Global variation in copy number in the human genome. Nature 2006;444:444-54.
Simioni M, Araujo TK, Monlleo IL, Maurer-Morelli CV, Gil-da-Silva-Lopes VL. Investigation of genetic factors underlying typical orofacial clefts: Mutational screening and copy number variation. J Hum Genet 2015;60:17-25.
Urquhart J, Black GC, Clayton-Smith J. 4.5 Mb microdeletion in chromosome band 2q33.1 associated with learning disability and cleft palate. Eur J Med Genet 2009;52:454-7.
Britanova O, Depew MJ, Schwark M, Thomas BL, Miletich I, Sharpe P, et al.
Satb2 haploinsufficiency phenocopies 2q32-q33 deletions, whereas loss suggests a fundamental role in the coordination of jaw development. Am J Hum Genet 2006;79:668-78.
Seelan RS, Mukhopadhyay P, Pisano MM, Greene RM. Developmental epigenetics of the murine secondary palate. ILAR J 2012;53:240-52.
Rauch TA, Wu X, Zhong X, Riggs AD, Pfeifer GP. A human B cell methylome at 100-base pair resolution. Proc Natl Acad Sci U S A 2009;106:671-8.
FitzPatrick DR, Carr IM, McLaren L, Leek JP, Wightman P, Williamson K, et al.
Identification of SATB2 as the cleft palate gene on 2q32-q33. Hum Mol Genet 2003;12:2491-501.
Leoyklang P, Suphapeetiporn K, Siriwan P, Desudchit T, Chaowanapanja P, Gahl WA, et al.
Heterozygous nonsense mutation SATB2 associated with cleft palate, osteoporosis, and cognitive defects. Hum Mutat 2007;28:732-8.
Van Buggenhout G, Van Ravenswaaij-Arts C, Mc Maas N, Thoelen R, Vogels A, Smeets D, et al.
The del(2)(q32.2q33) deletion syndrome defined by clinical and molecular characterization of four patients. Eur J Med Genet 2005;48:276-89.
Brewer CM, Leek JP, Green AJ, Holloway S, Bonthron DT, Markham AF, et al.
A locus for isolated cleft palate, located on human chromosome 2q32. Am J Hum Genet 1999;65:387-96.
[Figure 1], [Figure 2], [Figure 3]