Genetic Analysis of the Single-Nucleotide Polymorphisms rs880810, rs545793, rs80094639, and rs13251901 in Nonsyndromic Oral Clefts: A Case–Parent Trio Study

Abstract

Oral clefts, including cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP), are the most common types of congenital anomalies of the human face. Various genetic and environmental factors play a role in developing oral clefts. Several studies have shown the association of the PAX7 gene and the 8q24 region with these oral clefts in different populations worldwide. However, there are no reported studies on the possible connection between the PAX7 gene and the 8q24 region nucleotide variants and the risk of developing nonsyndromic oral clefts (NSOC) in the Indian population. Hence, this study aimed to test the possible association between PAX7 gene single-nucleotide polymorphisms (SNPs) rs880810, rs545793,rs80094639, and rs13251901 of the 8q24 region using a case-parent trio design. Forty case-parent trios were selected from the CLP center. Genomic DNA was isolated from the cases and their parents. The rs880810, rs545793, rs80094639, and rs13251901 were genotyped by the MassARRAY technique. PLINK software was used for statistical analysis. All the SNPs were tested for Hardy-Weinberg equilibrium. No statistical significance was found with any SNPs, as none of the genotyped SNPs showed a p-value of less than 0.05. Hence, the rs880810, rs545793, and rs80094639 of the PAX7 gene, and rs13251901 of the 8q24 region are not associated with NSOC in the Indian population.

Keywords

oral clefts - nonsyndromic - PAX7 - gene - genotyping

Introduction

Oral clefts, including cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP), are the most common congenital maxillofacial deformities[1] and form a major public health burden both in terms of medical expenses and disturbing to the patients and their families.[2] These oral clefts not only result in maxillofacial deformities but also occur with complications, such as regurgitation, middle ear infections, speech disorders, and other psychosocial problems.[3] [4]

The expression of oral clefts is usually unique and their occurrence varies by ethnicity, geography, and socioeconomic status.[5] Asians show higher prevalence rates, Caucasians show intermediate, and Africans show lower rates.[6] Literature suggests that 70% of oral cleft cases are nonsyndromic, whereas 30% are syndromic.[7] [8]

The paired box 7 (PAX7) gene has been indicated as a candidate gene for oral clefts. It also plays a role in neural crest development and maxillary growth. In addition, during craniofacial development, the PAX7 involves in the regulation of morphogenesis, survival, patterning, and specification of the frontonasal structures. So, any embryological disturbance in the neural crest development may lead to the development of oral clefts such as cleft lip and cleft palate.[9] A small number of human studies are available in the literature on this gene.[10] [11] [12] Traditionally, case–control studies have been used to study the etiology of oral clefts. Alternate to the traditional case–control study is the case-parent trio design in which cleft patient and their parents are analyzed for the genotypes.[13]

High-risk single-nucleotide polymorphisms (SNPs rs880810, rs545793, rs80094639 of the PAX7 gene, and rs13251901 of the 8q24 region are involved in the etiology of nonsyndromic oral clefts (NSOC) have been evaluated in different populations worldwide. A literature review showed the SNP rs13251901 belongs 8q24 region to have gene–gene (GxG)interaction with the PAX7 gene in cleft case-parent trios. Furthermore, several genomic studies have successfully replicated the associations between PAX7 and NSOC in other populations like Singaporean, Korean, Taiwanese, Philippines, Japanese, and Chinese.[10] [11] [12] However, the current literature search shows no studies on the association between these SNPs with NSOC in the Indian population. Hence, in this study, we tested for the possible association between PAX7 gene SNPsrs880810, rs545793, rs80094639, and rs13251901 of the 8q24 region in NSOC using a case-parent trio design. The advantages of the case-parent trio design include robustness in sample collection, testing maternal versus paternal effects, studying the parent-of-origin effects, and obtaining the accuracy of the genotyping.

Materials and Methods

Study Samples and Ethics Statement

Forty case-parent trios were selected from the CLP clinic of the Mallige Medical Centre, Bangalore. The sample size was estimated with a desired statistical power of 80% and a significance level of 0.05. Therefore, it was determined to include 40 case-parent trios. Geneticists examined each case to exclude the syndromic forms of oral clefts. The present research was approved by the Institutional Review Board (IRB) of DAPM RV Dental College, Bangalore (IRB no. 230/Vol-2/2017), and the research complied with the Helsinki declaration on medical research and ethics.

Inclusion and Exclusion Criteria

Oral clefts, either having CL or CP, or CLP only, were included in this study. Geneticists screened cleft children for any syndromes or abnormalities. Patients with a congenital anomaly and any syndromes associated with cleft were excluded from the study. The mean age and standard deviation of the included study subjects are cleft patients (6.62 ± 5.30 years), father (36.07 ± 6.42 years), and mother (30.37 ± 5.45 years). Written informed consent was obtained from all the cases and their parents.

SNP Selection, DNA Extraction, and Genotyping

The SNP srs880810, rs545793, rs80094639, and rs13251901 were selected from the previous genome-wide association studies, NCBI database( www.ncbi.nlh.nih.gov/snp/ ) and Ensembl genome database ( http://asia.ensembl.org/Homo_sapiens/Info/Index ). Genomic DNA was extracted from 3 mL of peripheral blood using QIAamp DNA Mini Kit (Qiagen Inc, California, United States) following the manufactures instructions. All the SNPs were genotyped using the Agena Bio MassARRAY (Agena Bioscience, Inc., San Diego, California, United States) platform using matrix-assisted laser desorption/ionization-time of flight mass spectrometry. The genomic data acquired from the genotyping were evaluated using statistical tests.

Statistical Data Analyses

All the case-parent trios underwent Hardy-Weinberg equilibrium (HWE) test, minor allele frequency (MAF) analysis, and the standard transmission disequilibrium test (TDT). The HWE, MAF, allelic TDT, and parent-of-origin effects in all the case-parent trios were calculated using PLINK software.[14] The 95% confidence interval for the assessed odds ratios was calculated, and a p-value less than 0.05 was set as statistically significant.

Results

Forty trios were selected for this study, and [Table 1] presents the demographic characteristics of study participants. The SNPs rs880810, rs545793, rs80094639, and rs13251901 ([Table 2]) followed the HWE test. The family-based TDT results ([Table 3]) showed no significant association with NSOC in the allelic frequencies of the Indian case-parent trios, as none of the genotyped SNPs showed a p-value less than 0.05.

Table 1 Summary of demographic variables

Sl. no.	Variables		Sample distribution n = 40 (%)
1	Age in years (mean ± standard deviation [SD])	Cleft patient	6.62 ± 5.30
		Father	36.07 ± 6.42
		Mother	30.37 ± 5.45
2	Gender of cleft patient	Male	15 (28%)
		Female	25 (72%)
3	Socioeconomic status	Lower middle class (III)	31 (72%)
		Lower class (V)	9 (28%)
4	Oral clefts	Cleft lip (CL)	13 (22%)
		Cleft palate (CP)	6 (16%)
		Cleft lip and palate (CLP)	21 (62%)

Table 2 List of the SNPs examined

Gene	SNP	Position	Call rate (%)	Alleles	Ancestral allele	Global MAF
PAX7	rs880810	1:18668060	100	C/G	C	0.11
PAX7	rs545793	1:18690028	100	G/A/C/T	G	0.34
PAX7	rs80094639	1:18723759	100	A/G	A	0.31
8q24	rs13251901	8:128287015	100	A/G	A	0.41

Abbreviations: A, adenine; G, guanine; C, cytosine; MAF, minor allele frequency (cleft cases); PAX7, paired box protein 7; SNP, single nucleotide polymorphism; T, thymine.

Table 3 Association between the PAX7 gene SNPs with NSOC

Gene	SNP	A1	A2	OR	CHISQ	p-Value
PAX7	rs880810	G	C	3	1	0.31
PAX7	rs545793	G	C	1.5	3	0.08
PAX7	rs80094639	G	A	2	4	0.24
8q24	rs13251901	G	A	0.68	0.92	0.33

Abbreviations: A, adenine; A1, major allele (wild allele); A2, minor allele (variant); C, cytosine; CHISQ, chi-squared test; G, guanine; NSOC, nonsyndromic oral clefts; OR, odds ratio; PAX7, paired box protein 7; SNP, single nucleotide polymorphism; p-Value < 0.05 is significant.

Discussion

Several genes and environmental factors play a significant role in the etiology of oral clefts. Advances in genetics and molecular biology techniques have discovered the genetic basis of these craniofacial defects. Numerous genes associated with the etiology of oral clefts have been discovered and registered in the Online Mendelian Inheritance in Man database. Different genetic approaches have been used to identify these genes causing oral clefts, including animal models,[15] linkage studies,[16] candidate gene analyses,[17] and genome-wide association studies.[18]

Case-parent trio design studies are commonly used in genetics as an alternative to case–control studies. In these studies, the affected cleft children and their parents (father and mother) are genotyped.[19] Therefore, it does not require control samples because the untransmitted parental genotypes are used as “controls” for the transmitted alleles or genotypes.[20] In addition, the case–parent trio design has the advantage of testing maternal versus paternal effects. It separates these from the effects of the fetal genotype versus parental origin robustly.[21]

The PAX7 gene encodes specific DNA-binding transcription factors and is involved in neural crest development which is closely associated with cranial and maxillofacial development. In turn, any defect in neural crest development may cause the development of oral clefts.[22] The PAX7 gene, along with PAX3, plays an essential role during craniofacial development in the regulation of morphogenesis, survival, patterning, and specification of the frontonasal structures. PAX3/PAX7 double mutant mice showed impaired function against environment-related teratogenesis and exhibited more severe facial morphogenesis defects. The transcription factors of this PAX7 and PAX3 help to maintain the proliferative cells during the development of embryonic and fetal muscles of the trunk and limbs.

Several studies have documented their findings regarding the 1p36.13 chromosomal region and shown the association of the PAX7 gene in the development of oral clefts in different populations worldwide. Furthermore, considering the G × G interaction, the polymorphisms of the PAX7 showed a significant association with NSOC in some populations. But most of the studies were largely focused on the Western population,[10] [11] and there are no reported case-parent trio studies of the PAX7 gene in the Indian population.

Hence, the SNP srs880810, rs545793, rs80094639, and rs13251901 of the 8q24 region were analyzed in the NSOC patients in a case-parent trio design; there is no reported trio study on the role of these SNPs in the Indian population. Forty case-parent trios were selected from the cleft center with cleft samples consisting of males (15) and females (25). The gender difference in our study samples was more compared to the global prevalence of NSOC in males and females. This study results suggest that rs880810, rs545793, rs80094639, and rs13251901 are not significantly associated with NSOC.

In a study by Sull et al,[10] a total of 297 CLP trios from four different populations (Singaporean, Korean, Taiwanese, and Maryland) showed that the rs880810 and rs545793 of PAX7were not associated with NSOC.

Also, in multiplex families of CLP, the rs80094639 showed no significant association.[23] A study showed the SNP rs13251901 belongs 8q24 region to have GxG interaction with the PAX7 gene in cleft case-parent trios and showed no significance with NSOC.[24]

The difference or varying results related to these SNPs in the genetic etiology of oral cleft might be due to the different ethnicity, epigenetic factors, and environmental factors.

Limitations

The limitations of our study are a relatively smaller sample size and the analysis of only four polymorphisms. Further studies should be conducted with more SNPs of the PAX7 gene and the 8q24 region with a bigger sample size for the Indian population to substantiate their role in the etiology of NSOC.

Conclusion

This study revealed that the SNPs rs880810, rs545793, rs80094639 of the PAX7 gene and rs13251901 of the 8q24 region are not associated with the risk of NSOC in the Indian population.

Conflict of Interest

None declared.

Acknowledgments

We are thankful to all the patients and families who participated in this genetic study. We also thank Dr. D. V. S Sudhakar, Scientist C, ICMR-NIRRH, and Dr. Archana Krishnamurthy for the statistical help.

• References

• 1 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 (01) 1-11
• 2 Mitchell LE, Beaty TH, Lidral AC. et al; International Consortium for Oral Clefts Genetics. Guidelines for the design and analysis of studies on nonsyndromic cleft lip and cleft palate in humans: summary report from a Workshop of the International Consortium for Oral Clefts Genetics. Cleft Palate Craniofac J 2002; 39 (01) 93-100
• 3 Kahraman A, Yüce S, Koçak ÖF, Canbaz Y, Işik D. Congenital lateral cleft palate of unknown etiology. J Craniofac Surg 2015; 26 (04) 1332-1333
• 4 Indencleef K, Hoskens H, Lee MK. et al. The intersection of the genetic architectures of orofacial clefts and normal facial variation. Front Genet 2021; 12: 626403
• 5 Vanderas AP. Incidence of cleft lip, cleft palate, and cleft lip and palate among races: a review. Cleft Palate J 1987; 24 (03) 216-225
• 6 Tolarová MM, Cervenka J. Classification and birth prevalence of orofacial clefts. Am J Med Genet 1998; 75 (02) 126-137
• 7 Neela PK, Gosla SR, Husain A, Mohan V. CRISPLD2 gene polymorphisms with nonsyndromic cleft lip palate in Indian population. Glob Med Genet 2020; 7 (01) 22-25
• 8 Huang L, Jia Z, Shi Y. et al. Genetic factors define CPO and CLO subtypes of nonsyndromic orofacial cleft. PLoS Genet 2019; 15 (10) e1008357
• 9 Khan AN, Prashanth CS, Srinath N. Genetic etiology of cleft lip and palate. AIMS Mol Sci 2020; 7 (04) 328-348
• 10 Sull JW, Liang KY, Hetmanski JB. et al. Maternal transmission effects of the PAX genes among cleft case-parent trios from four populations. Eur J Hum Genet 2009; 17 (06) 831-839
• 11 Guo Q, Li D, Meng X. et al. Association between PAX7 and NTN1 gene polymorphisms and nonsyndromic orofacial clefts in a northern Chinese population. Medicine (Baltimore) 2017; 96 (19) e6724
• 12 Butali A, Suzuki S, Cooper ME. et al. Replication of genome wide association identified candidate genes confirm the role of common and rare variants in PAX7 and VAX1 in the etiology of nonsyndromic CL(P). Am J Med Genet A 2013; 161A [phone]
• 13 Khan MI, Cs P. Case-parent trio studies in cleft lip and palate. Glob Med Genet 2020; 7 (03) 75-79
• 14 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 (03) 559-575
• 15 Mansouri A, Stoykova A, Torres M, Gruss P. Dysgenesis of cephalic neural crest derivatives in Pax7-/- mutant mice. Development 1996; 122 (03) 831-838
• 16 Scapoli C, Collins A, Martinelli M, Pezzetti F, Scapoli L, Tognon M. Combined segregation and linkage analysis of nonsyndromic orofacial cleft in two candidate regions. Ann Hum Genet 1999; 63 (Pt 1): 17-25
• 17 Watanabe A, Akita S, Tin NT. et al. A mutation in RYK is a genetic factor for nonsyndromic cleft lip and palate. Cleft Palate Craniofac J 2006; 43 (03) 310-316
• 18 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 (06) 525-529
• 19 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 (06) 505-512
• 20 Schwender H, Li Q, Neumann C. et al. Detecting disease variants in case-parent trio studies using the ioconductor software package trio. Genet Epidemiol 2014; 38 (06) 516-522
• 21 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 (03) 294-303
• 22 Neela PK, Gosla SR, Husain A, Mohan V, Thumoju S, Bv R. Association of MAPK4 and SOX1-OT gene polymorphisms with cleft lip palate in multiplex families: a genetic study. J Dent Res Dent Clin Dent Prospect 2020; 14 (02) 93-96
• 23 Neela PK, Reddy GS, Husain A, Mohan V, Thumoju S, Bv R. Association of single nucleotide polymorphisms on locus 18q21.1 in the etiology of nonsyndromic cleft lip palate (NSCLP) in Indian multiplex families. Glob Med Genet 2021; 8 (01) 24-31
• 24 Xiao Y, Taub MA, Ruczinski I. et al. Evidence for SNP-SNP interaction identified through targeted sequencing of cleft case-parent trios. Genet Epidemiol 2017; 41 (03) 244-250

Publikationsverlauf

Artikel online veröffentlicht:
28. März 2023