The risk of nonsyndromic cleft lip with or without cleft palate and Vax1 rs7078160 polymorphisms in southern Han Chinese

Abstract

Introduction:

Non-syndromic cleft lip with or without cleft palate is a common worldwide birth defect due to a combination of environmental and genetic factors. Genome-wide association studies reported the rs7078160 of Vax1 is closely related to non-syndromic cleft lip with or without cleft palate in European populations. The following studies showed the same results in Mongolian, Japanese, Filipino, Vietnamese populations etc. However, conflicting research had been reported in Chinese population,

Objective:

The aim of this study was to investigate the association between the rs7078160 polymorphism and non-syndromic cleft lip with or without cleft palate in Southern Chinese patients.

Methods:

In this study, we investigated the polymorphism distribution of rs7078160 in 100 complete patient trios (39 patients with non-syndromic cleft lip and palate; 36 patients with non-syndromic cleft lip only; 25 had non-syndromic cleft palate only; and their parents) from Southern ethnic Han Chinese. 60 healthy trios were selected as control. Polymerase chain reaction and Sanger sequencing were used to genotype rs7078160 in Vax1; both case-control and family-based associations were analyzed.

Results:

The case-control analyses revealed the rs7078160 polymorphism was significant, associated with non-syndromic cleft lip with or without cleft palate (p = 0.04) and non-syndromic cleft lip and palate (p = 0.01), but not associated with non-syndromic cleft lip only and nonsyndromic cleft palate only patients. The genotype composition of rs7078160 comprises mutated homozygous AA, heterozygous AG and wild homozygous GG. Cases with AG + AA genotypes compared with GG homozygotes showed an increased risk of non-syndromic cleft lip with or without cleft palate (p = 0.04, OR = 2.05, 95% CI: 1.01-4.16) and non-syndromic cleft lip and palate (p = 0.01, OR = 3.94, 95% CI: 1.34-11.54). In addition, we did not detect any transmissiondisequilibrium in rs7078160 (p = 0.68).

Conclusion:

This study suggests that rs7078160 polymorphism is a risk factor of non-syndromic cleft lip with or without cleft palate, and Vax1 is strongly associated with non-syndromic cleft lip with or without cleft palate in Southern Chinese Han populations.

KEYWORDS
Cleft lip; Cleft palate; Single nucleotide polymorphisms; Vax1 ; rs7078160

Resumo

Introdução:

A fenda labial não sindrômica, com ou sem fenda palatina, é um defeito congênito comum em todo o mundo, devido a uma combinação de fatores ambientais e genéticos. O genome-wide association studies relatou que o polimorfismo rs7078160 do Vax1 está intimamente relacionado à fenda labial não sindrômica, com ou sem fenda palatina em populações europeias. Estudos subsequentes mostraram os mesmos resultados nas populações mongol, japonesa, filipina e vietnamita etc. No entanto, pesquisas conflitantes foram relatadas na população chinesa.

Objetivo:

Investigar a associação entre o polimorfismo rs7078160 e fenda labial não sindrômica, com ou sem fenda palatina, em pacientes do sul da China.

Método:

Tentamos investigar a distribuição do polimorfismo rs7078160 em 100 trios completos de pacientes (39 pacientes com fenda labial e palatina não sindrômica; 36 pacientes com fenda labial somente, não sindrômica; 25 com fenda palatina somente, não sindrômica e seus pais), da etnia Han do sul da China, e em 60 trios saudáveis selecionados como controle. Reação de polimerase em cadeia e o sequenciamento de Sanger foram uszados para genotipar o polimorfismo rs7078160 do Vax1 e tanto os casos-controle quanto as associações baseadas na família foram analisadas.

Resultados:

As análises de caso-controle revelaram que o polimorfismo rs7078160 estava significativamente associado a fenda labial não sindrômica, com ou sem fenda palatina (p = 0,04) e fenda labial e palatina não sindrômica (p = 0,01), mas não estava associado a pacientes com fenda labial somente não sindrômica e fenda palatina somente não sindrômica. A composição do genótipo de rs7078160 compreende AA homozigoto mutado, AG heterozigoto e GG homozigoto selvagem. Casos com genótipos AG + AA comparados com GG homozigotos mostraram um risco aumentado de fenda labial não sindrômica, com ou sem fenda palatina (p = 0,04, OR = 2,05, IC de 95%: 1,01 ± 4,16) e fenda labial e palatina não sindrômica (p = 0,01, OR = 3,94, IC 95%: 1,34-11,54). Além disso, não detectamos desequilíbrio de transmissão em rs7078160 (p = 0,68). Conclusão: Este estudo sugere que o polimorfismo rs7078160 foi um fator de risco para fenda labial não sindrômica, com ou sem fenda palatina, e o gene Vax1 está fortemente associado com fenda labial não sindrômica, com ou sem fenda palatina em populações da etnia Han do sul da China.

PALAVRAS-CHAVE
Fenda labial; Fenda palatina; Polimorfismos de nucleotídeo único; Vax1 ; rs7078160

Introduction

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is one of the most common worldwide birth defect, with an incidence of 1.2 in 1000 live births.¹ In China, the prevalence of NSCL/P is higher than that in other racial and ethnic populations (1.48-3.27 in 1000 live births).²,³ Currently, with the development of plastic surgery, more patients with cleft lip or palate can be repaired by surgery⁴; however, they may still suffer from a handful of problems such as craniofacial deformation and ambiguous speech. Not only that, although international charities such as the Smile Train have provided support for many patients, the therapy of NSCL/P may be a heavy economic burden for some families, especially in the developing countries. Hence, clarifying the etiologies of NSCL/P is very important for disease diagnosis and prevention.

The etiologies and pathogenesis of NSCL/P are complex, including both genetic and environmental factors or their combination effort.³ In recent years, with the development of high-throughput sequencing technology and statistical methodology, Genome-Wide Association Studies (GWAS) have been made to identify candidate genes involved in the development of NSCL/P.⁵ GWAS evidences have identified dozens of potential causative loci and genes of NSCL/P, such as Irf6, Mafb, Grhl3, Abca4, Nog, Spry2, Tpm1, 8q21.3, Crebbp and Vax1.⁶-¹⁰ However, in the subsequent replications, few of those loci/genes showed consistently positive correlations across all studies. In fact, for more candidate genes, inconsistent results were reported in populations from different races and regions. Vax1 is one of such genes. A GWAS conducted by Mangold et al. identified a new loci rs7078160 associated with NSCL/P at 10q25 in European population.¹¹ Subsequently, Beaty et al. conducted another GWAS study in European and Asian groups, their founding confirmed rs7078160 of Vax1 was correlated with NSCL/P.⁶ Following studies in Mongolian, Japanese, Filipinos and Vietnamese populations supported above results.¹²-¹⁴ However, conflicting research has been reported in the Chinese population. A study by Li et al. found rs7078160 of Vax1 was not associated with NSCL/P in northern Chinese Han population.¹⁵ In this study, we first investigated the correlation between Vax1 rs7078160 polymorphisms and the risk of NSCL/P in the Southern Han Chinese population. Considering that complete trio samples could provide more reliable data than case-control samples, case/control-parents design was used in this study, so we could detect the genotypic and allelic distributions of target single nucleotide polymorphism (SNP) while analyzing the parental origin effects.

Methods

Subjects

The study was conducted in the Department of Oral and Maxillofacial Surgery, affiliated with the Stomatological Hospital of Zunyi Medical University (Zunyi, Guizhou, China) from 2016 to 2019. One hundred NSCL/P patients (39 patients with non-syndromic cleft lip and palate, NSCLP; 36 patients with non-syndromic cleft ip only, NSCLO; and 25 with non-syndromic cleft palate only, NSCPO) and their parents were recruited as subjects. Sixty healthy matched trios were selected as the control. The protocol of the study was approved by the Institutional Ethics Committee of Zunyi Medical University (ZMU-HG1/Sept/25/16 date: 11/01/2016). Informed consents were collected from the parents or their legal guardians.

Clinical assessment was conducted by experienced oral and maxillofacial surgeons, to exclude the syndromic CL/P patients associated with any other congenital disabilities. All NSCLP, NSCLO and NSCPO cases were unilateral and none of the patients had affected relatives; patients with family history of genetic disease were also excluded. Healthy volunteers without a family history of cleft lip and/or palate or other abnormality were chosen as controls. All cases and controls were Chinese Han population, with the same age and gender distribution, and lived in Southern China for at least three generations. Three milliliters peripheral venous blood was collected from each subject and stored at -80^◦C.

Figure 1
DNA sequencing results of rs7078160. (A) Mutated homozygous genotype AA. (B) Heterozygous genotype AG. (C) Wild homozygous genotype GG.

DNA extraction, PCR amplification and sequencing

Genomic DNA was extracted from blood sample by TIANamp Blood DNA kit (Tiangen Biotech, Beijing, China). The specific rs7078160 primers were designed according to Vax1 reference sequences in NCBI (National Center for Biotechnology Information, Bethesda, MD, USA) by Primer 5.0 software (Premier Biosoft International, Palo Alto, CA, USA). The sequence of forward primer was 5^, TGGGAAGTGGGTGAGATGGA 3^,, reverse primer was 5^, ATTGGGCGGACCCAGTAAAG 3^,. All PCR amplifications were performed on a Gene Amp PCR system 9700 thermal cycler (Applied Biosystems, Foster City, CA, USA). The reaction mixture consisted of 10-30 ng of genomic DNA, 1 µ,L of each primer (10 µ,moL L^-1), 4 µ,L dNTPs (2.5 mmoL L^-1 of each nucleobase), 5 µ,L 10 Buffer (with MgCl₂) and 1.0 U of Taq polymerase (TaKaRa, Dalian, China), add ddH₂O to 50 µ,L volume. The PCR amplification program was: (1) denaturing (94^◦C, 4 min); (2) 34 cycles including denaturing (94^◦C, 30s), annealing (58^◦C, 30 s), extension (72^◦C, 40 s), (3) final extension (72^◦C, 5 min). PCR product was detected by agarose gel electrophoresis and purified. Sequencing was conducted by BGI (Beijing Genomics Institute, Beijing, China). We aligned and edited sequences using DNA sequencing software Chromas 2.6.6 (Technelysium, Queensland, Australia) and MEGA X.¹⁶ The DNA sequencing results of rs7078160 were shown in Figure 1.

Statistical methods

The genotype and allele frequencies were summarized in the cases (total group, NSCLP, NSCLO and NSCPO subgroup) and the control. The Hardy--Weinberg equilibrium (HWE) states that both allele and genotype frequencies in a random mating population remain constant, in order to demonstrate the genetic stability of our subjects, HWE was assessed according to the genotype frequency. Results were tested with a Chi-square test in SPSS 16.0 (SPSS Inc., Chicago, IL, USA). In case-parents’ trios, the transmission disequilibrium test (TDT) was conducted by statistical analysis for genetic epidemiology (S.A.G.E.) software package. TDT is a family-based test which offers a powerful way to test the transmission of target alleles from heterozygous parents to affected offspring. The odds ratio (OR) and 95% confidence intervals (95% CI) were also calculated, OR > 1 signified an increased risk of NSCL/P. All data were tested by two-sided p-value, p < 0.05 were accepted as significant.

Results

The Hardy--Weinberg equilibrium analysis revealed all the subjects showed a genetic equilibrium at rs7078160 (Table 1). The distribution of the rs7078160 variant genotypes and alleles in both case groups and controls are presented in Table 1. The proportions of genotypes were 21% AA, 57% AG, 22% GG in cases and 11.7% AA, 51.7% AG, 36.7% GG in controls. The G allele frequency was 50.5% in cases and 62.5% in controls. The allele frequencies of rs7078160 in the total case group (p = 0.04) and NSCLP subgroup (p = 0.01) were substantially different when contrasted to the control group, however, the allele frequencies in the NSCLO (p = 0.19) and NSCPO subgroup (p = 0.43) were not significant when compared with the control group (Table 1).

The OR and 95% CI between AG + AA and GG for rs7078160 were also significant difference in the total case group (p = 0.04, OR = 2.05, 95% CI: 1.01-4.16) and the NSCLP sub-group (p = 0.01, OR = 3.94 ± 95% CI: 1.34-11.54); the results were not significant in the NSCLO (p = 0.24) and NSCPO subgroup (p = 0.68) (Table 2). The transmission disequilibrium test was conducted for trios with heterozygote genetic structure, and no statistically significant transmitted disequilibrium was detected in rs7078160 (p = 0.68) (Table 3).

Discussion

The important etiological role of rs7078160 was first reported by Mangold et al. in a genome-wide association study of European NSCL/P populations.¹¹ The top susceptive NSCL/P SNP in 10q25 was rs7078160 in Vax1 (Ventral Anterior Homeobox 1). Vax1 was a crucial gene for the development of the basal forebrain and visual system in animal experiments, and it also participated in the rostral formation of vertebrates.¹⁷ Slavotinek et al. reported that Vax1 mutation is associated with human cleft lip and palate, absence of the pineal gland, corpus callosum agenesis, hippocampal malformations and small optic nerves.¹⁸

To date, several studies have confirmed the association between rs7078160 and NSCL/P in different racial and ethnic populations, including Japanese, Polish, central Africans and southeast Asian populations.¹²-¹⁴ However, inconsistent results were found in Chinese, Kenyan and another Polish cases.¹⁵,¹⁹,²⁰ In Chinese, research of Li et al. fail to identify the risk of Vax1 in NSCL/P in northern Chinese population.¹⁵ However, recently a study conducted by Zhang et al. found the Vax1 gene was a risk factor for NSCL/P in western Han Chinese population.⁵ Apart from the geographical difference, the different studies method might be an explanation; the former studies used case-control design, while Zhang et al. adopted a case-parent trio design.⁵ In this study, we chose the Southern Chinese Han population as the object. The allele frequencies of rs7078160 in the NSCLP group were substantially different when contrasted to the control group (p = 0.02, respectively) (Table 1). OR and 95% CI between AG + AA and GG showed a significant difference in the total case group (p = 0.04, OR = 2.05, 95% CI: 1.01-4.16) and NSCLP subgroup (p = 0.01, OR = 3.94, 95% CI: 1.34-11.54), suggested that allele A of rs7078160 may be a risk factor for NSCL/P, but not associated with NSCLO and NSCPO (Table 2). The TDT test did not detect significant transmitted disequilibrium in rs7078160 (p = 0.68) (Table 3). This result was conflicting with the previous NSCL/P studies by Li et al., in which they failed to find an association of Vax1 and NSCL/P in northern Chinese Han population.¹⁵ The inconsistency might be due to the geographical difference and different sampling design. The complete case-parent trio might provide more comprehensive genetic information compared to the case-control design.

Table 1
Genotypic and allelic distributions of rs7078160 in NSCL/P patients and controls (^ap < 0.05).

Genotype	Total	NSCLP	NSCLO	NSCPO	Controls
Test for HWE
Chi-square	1.96	3.41	0.47	0.02	0.63
p-value	0.16	0.06	0.49	0.90	0.43
Genotype distribution
AA (n, %)	21 (21)	9 (23.1)	7 (19.4)	5 (20.0)	7 (11.7)
AG (n, %)	57 (57)	25 (64.1)	20 (55.6)	12 (48.0)	31 (51.7)
GG (n, %)	22 (22)	5 (12.8)	9 (25.0)	8 (32.0)	22 (36.7)
Chi-square	4.99	7.48	1.95	1.02
p-value	0.08	0.02a	0.38	0.60
Allele frequency A (n, %)	99 (49.5)	43 (54.4)	34 (47.2)	22 (44.0)	45 (37.5)
G (n, %)	101 (50.5)	35 (44.3)	38 (52.8)	28 (56.0)	75 (62.5)
Chi-square	4.36	5.95	1.76	0.62
p-value	0.04a	0.01a	0.19	0.43

NSCLP, non-syndromic cleft lip and palate; NSCLO, non-syndromic cleft lip only; NSCPO, non-syndromic cleft palate only.

Table 2
Results of association tests with rs7078160 in non-syndromic cleft lip with or without cleft palate (NSCL/P) (^ap < 0.05).

				Variable
	AG + AA	GG	Chi-square	p-value	Odds ratio	95% CI
Controls	38	22
Total	78	22	4.05	0.04a	2.05	1.01-4.16
Non-syndromic cleft lip and palate (NSCLP)	34	5	6.78	0.01a	3.94	1.34-11.54
Non-syndromic cleft lip only (NSCLO)	27	9	1.40	0.24	1.74	0.69-4.35
Non-syndromic cleft palate only (NSCPO)	17	8	0.17	0.68	1.23	0.46-3.31

Table 3
Results of the transmission disequilibrium test.

rs7078160	Transmitted	Not transmitted	Total
	A	G
A	27	50	77
G	45	22	67
Total	72	72	144
Chi-square	0.26
p-value	0.68

Collectively, this study has demonstrated that rs7078160 is involved in the development of NSCL/P in Southern Han Chinese. Our study provided additional evidence for the etiology of Vax1 in NSCL/P of Chinese Han population. However, our study still has several potential limitations. First, we focused on only one SNP rs7078160; the study of more candidates with SNP could enhance the significance of our study. Second, due to the difficulties in collecting the complete trio samples, the number of cases and control trios were relatively limited. To increase the number of the subjects would strengthen our results. Thirdly, the environmental interactions and maternal factors during pregnancy were not evaluated in our study. Further efforts are still needed to clarify the deeper association between heredity and environment.

Conclusion

This study had revealed that rs7078160 was involved in the development of NSCL/P in Southern Chinese Han population, demonstrating Vax1 might be associated with NSCL/P. Those findings could be helpful in extending our knowledge about NSCL/P etiology.

Acknowledgments

We are grateful to all of the individuals who made this study possible, especially the patients, their parents, and the control subjects who voluntarily cooperated with this study. This research was supported by Top Science and Technology Talent Project in Higher Education Institution of Guizhou Province (Qian-Jiao-He KY Zi [2016]080), Outstanding Young Talent Project of Zunyi Medical University (17zy-002) and Science and Technology Project of Zunyi Science and Technology Bureau (ZSKHZC NS 2019-15).

References

• ¹
  Rahimov F, Jugessur A, Murray JC. Genetics of nonsyndromic orofacial clefts. Cleft Palate-Craniofac J. 2012;49:73-91.
• ²
  Chang WJ, See LC, Lo LJ. Time trend of incidence rates of cleft lip/palate in Taiwan from 1994 to 2013. Biomed J. 2016;39:150-4.
• ³
  Dixon MJ, Marazita ML, Beaty TH, Murray JC. Cleft lip and palate: synthesizing genetic and environmental influences. Nat Rev Genet. 2011;12:167-78.
• ⁴
  Fitzsimons KJ, Copley LP, Deacon SA, van der Meulen JH. Hospital care of children with a cleft in England. Arch Dis Child. 2013;98:970-4.
• ⁵
  Zhang BH, Shi JY, Lin YS, Shi B, Jia ZL. VAX1 gene associated non-syndromic cleft lip with or without palate in Western Han Chinese. Arch Oral Biol. 2018;95:40-3.
• ⁶
  Beaty TH, Murray JC, Marazita ML, Munger RG, Ruczinski I, Hetmanski JB, 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:525-9.
• ⁷
  Birnbaum S, Ludwig KU, Reutter H, Herms S, Steffens M, Rubini M, et al. Key susceptibility locus for nonsyndromic cleft lip with or without cleft palate on chromosome 8q24. Nat Genet. 2009;41:473-7.
• ⁸
  Grant SFA, Wang K, Zhang H, Glaberson W, Annaiah K, Kim CE, et al. A Genome-wide association study identifies a locus for nonsyndromic cleft lip with or without cleft palate on 8q24. J Pediatr. 2009;155:909-13.
• ⁹
  Ludwig KU, Mangold E, Herms S, Nowak S, Reutter H, Paul A, et al. Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci. Nat Genet. 2012;44:968-71.
• ¹⁰
  Sun YM, Huang YQ, Yin AH, Pan YH, Wang YR, Wang C, et al. Genome-wide association study identifies a new susceptibility locus for cleft lip with or without a cleft palate. Nat Commun. 2015;6:6414.
• ¹¹
  Mangold E, Ludwig KU, Birnbaum S, Baluardo C, Ferrian M, Herms S, et al. Genome-wide association study identifies two susceptibility loci for nonsyndromic cleft lip with or without cleft palate. Nat Genet. 2010;42:24-6.
• ¹²
  Butali A, Suzuki S, Cooper ME, Mansilla AM, Cuenco K, Leslie EJ, 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;161:965-72.
• ¹³
  Figueiredo JC, Stephanie LY, Raimondi H, Magee K, Baurley JW, Sanchez-Lara PA, et al. Genetic risk factors for orofacial clefts in Central Africans and Southeast Asians. Am J Med Genet A. 2014;164A:2572-80.
• ¹⁴
  Mostowska A, Hozyasz KK, Wojcicka K, Biedziak B, Jagodzinski PP. Polymorphic variants at 10q25.3 and 17q22 loci and the risk of non-syndromic cleft lip and palate in the polish population. Birth Defects Res A Clin Mol Teratol. 2012;94:42-6.
• ¹⁵
  Li DM, Liu TT, Meng XB, Guo Q, Shi JN, Hao YR, et al. Polymorphic variants in VAX1 and the risk of nonsyndromic cleft lip with or without cleft palate in a population from northern China. Medicine. 2017;96:e6550.
• ¹⁶
  Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547-9.
• ¹⁷
  Hallonet M, Hollemann T, Wehr R, Jenkins NA, Gruss P. Vax1 is a novel homeobox-containing gene expressed in the developing anterior ventral forebrain. Development. 1998;125:2599-610.
• ¹⁸
  Slavotinek AM, Chao R, Vacik T, Yahyavi M, Abouzeid H, Bardakjian T, et al. VAX1 mutation associated with microphthalmia, corpus callosum agenesis, and orofacial clefting: the first description of a VAX1 phenotype in humans. Hum Mutat. 2012;33:364-8.
• ¹⁹
  Weatherley-White RC, Ben S, Jin Y, Riccardi S, Spritz RA. Analysis of genomewide association signals for nonsyndromic cleft lip/palate in a Kenya African cohort. Am J Med Genet A. 2011;155A:2422-5.
• ²⁰
  Zawi´slak A, Wo´zniak K, Jakubowska A, Lubin´ski J, Znamirowska-Bajowska A. Polymorphic variants in vax1 gene (rs7078160) and bmp4 gene (rs762642) and the risk of non-syndromic orofacial clefts in the polish population. Med Wieku Rozwoj. 2014;18:16-22.

Publication Dates

• Publication in this collection
  06 Dec 2021
• Date of issue
  2021

History

• Received
  03 June 2020
• Accepted
  28 Aug 2020
• Published
  11 Oct 2020