Introduction
============

Cleft palate (CP) is a common birth defect, which has a lower birth prevalence compared to cleft lip with/without cleft palate (CL/P): 1/2500 live births vs. 1/700; but CP shows less variability in birth prevalence across populations compared to CL/P ([@B1]).

In nonsyndromic cleft palate only (NSCPO), affected individuals have no other physical or developmental anomalies. Most studies suggest that about 50% of CP is nonsyndromic ([@B2]). Both population studies and family studies suggested that genetic factors played a critical role in the etiology of NSCPO ([@B3],[@B4]). Among first degree relatives, the relative risks of recurrence risks were 56 for cleft palate only and 32 for any cleft lip when compared to the general population in Norway ([@B3]); several genes have been identified for syndromic forms of CP, few have been identified as influencing risk to NSCPO. The etiology of this complex trait has been widely studied in order to search for the risk factors and to design strategies for prevention.

Genome wide association studies entail the matching of a given human genome sequence with an annotated, high-resolution map of common genetic variation; They are contributing a lot to our understanding of diseases to which there is a genetic predisposition ([@B5]). Genome wide association studies (GWAS) already have identified tens of susceptible loci for cleft lip with or without cleft palate (CL/P).

However, whether these loci associated with nonsyndromic cleft palate only for Han Chinese remains unknown. In this study, we replicated 38 SNPs from 19 genes/ regions of 11 chromosomes from previous GWAS ([@B6]-[@B10]) and other studies ([@B11]-[@B13]) with prior compelling evidence contributing to NSCL/P to investigate their roles in Han Chinese population.

Material and Methods
====================

\- Samples

Our samples consisted of 144 complete case-parent trios with nonsydromic cleft palate only (NSCPO), 59 males, 82 females and 3 unknown gender of the probands. All subjects were self-identified as Western Han Chinese, they were recruited between 2008 and 2013 from the Cleft Surgery Department of West China Hospital of Stomatology, Sichuan University. Written informed consent was obtained from parents on behalf of the children and all affected individuals old enough to give their own consent in this study. The consent procedure and this study approved by the Hospital Ethics Committee (HEC) of West China Hospital of Stomatology, Sichuan University.

\- Genotyping

Venous blood samples were drawn from participants and DNA was extracted by phenol chloroform extraction protocol. The SNP genotyping work was performed using a custom-by-design 2x48-Plex SNPscanTM Kit (Cat\#:G0104, Genesky Biotechnologies Inc., Shanghai, China).

\- Statistical analysis

The unaffected parents were underwent Hardy-Weinberg Equilibrium (HWE) analysis and minor allele frequency (MAF) determination. HWE, MAF, allelic TDT and parent-of-origin effects, were calculated using PLINK ([@B14]). Pairwise LD as both their D' and r2 were computed for all the SNPs using the haploview program (http://www.broad.mit.edu/haploview/haploview). Genotypic TDT and Likelihood ratio test for epistatic interactions based on genotypic TDTs were determined with R Package trio (v1.4.23) ([@B15]), all two-way interactions comprised a matrix in genotype format were tested using the function colGxG, without specifying the genes, and the interactions between all 38 (38-1)/2 pairs of the SNPs in a matrix were tested. We used a Bonferonni correction for 38 tests to determine a threshold for formal significance of *p*=0.0013.

Results
=======

\- Transmission Disequilibrium Test

All SNPs were passed the HWE test (*p* \>0.05) ([Table 1](#T1){ref-type="table"}). Allelic TDT results showed that T allele at rs742071 (PAX7) (*p*=0.025, ORtransmission=3.0, 95%CI: 1.09-8.25) and G allele at rs2485893 (10kb 3' of SYT14) were associated with NSCPO (*p*= 0.0036, ORtransmission= 0.60, 95%CI: 0.42-0.85) ([Table 2](#T2){ref-type="table"}). Genotypic TDT based on 3 pseudo controls further confirmed that rs742071 (*p*-value=0.03, ORtransmission=3.0, 95%CI: 1.09-8.25) and rs2485893 are associated with NSCPO under additive model (*p*-value=0.02, ORtransmission=0.66, 95%CI: 0.47-0.92) ([Table 3](#T3){ref-type="table"}).

Table 1Minor allele frequency and Hardy-Weinberg Equilibrium Test of the SNPs for NSCPO.

Table 2Allelic TDT results of the SNPs for NSCPO.

Table 3Genotypic TDT Based on 3 Pseudo Controls for NSCPO under different genetic models.

\- Parent-of-Origin Effects

There was no significant difference of minor allele transmissions between the maternal and paternal for all SNPs (data not shown). However, we found an excess of maternal transmission of the allele G at rs2485893 (*p*=0.026), allele A at rs8001641 (*p*=0.018) and allele G at rs17563 (*p*=0.045) compared with the paternal ([Table 4](#T4){ref-type="table"}), which might warrant future investigations.

Table 4Parent-of-origin effect of effect of the SNPs.

\- Gene by Gene Interactions

The Genotypic TDT for epistatic interactions showed that rs4844913 (43kb 3' of DIEXF) interacted with rs11119388 (SYT14) (*p*=1.80E-08) and rs6072081 (53kb 3' of MAFB) interacted with rs6102085 (33kb 3' of MAFB) (*p*= 3.60E-04) for NSCPO (Fig. [1](#F1){ref-type="fig"}).

Figure 1Genotypic TDT for Epistatic Interactions of Pair-wise SNPs.

\- Pair-wise Linkage Disequilibrium and Haplotype Analysis

We calculated the pair-wise linkage disequilibrium (LD) of SNPs on chromosome 1 based the association results above. There were strong LD between six pairs of SNPs (rs4920552-rs766325, rs126280-rs642961, rs2064163-12063989, rs4844913-rs9429830, rs4844913-rs227178 and rs9429830-227178) with D'\>0.93 and r2\>0.80, which distributed on two haplotypes (Fig. [2](#F2){ref-type="fig"}). Based on this, we tried to test if these pairs of SNPs segregate together among NSCPO by carrying out the haplotype analysis. The results did not show any significance (data not shown).

Figure 2Pair-wise Linkage Disequilibrium (A) and Haplotype (B) of SNPs at Chromosome 1.

Discussion
==========

Rs742071 is located in the intron of PAX7 which is involved in neural crest induction and is expressed in cranial neural crest cells, and mice lacking Pax7 have malformations of the nasal and maxillary structures ([@B16]). The PAX7 was a second tier GWAS hit ([@B8]), later it was confirmed by replication among European and Southeast Asian ([@B17]) and GWAS meta-analysis ([@B10]). Recently, Leslie *et al.* 2015 performed targeted sequencing of 13 regions from GWASs and other studies in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses, the results indicated that de novo mutation p.Ala259Val disrupted PAX7 function and might contribute to CLP pathogenesis in this individual ([@B18]). Although intronic SNPs do not typically alter protein structure, associations with intronic variants have been reported for a number of complex diseases. In this study, we found rs742071 (PAX7) was associated with NSCPO and had a larger genetic effects compared with the associations with NSCL/P from previous GWASs ([@B8],[@B10]); Motif analyses by HaploReg indicated that the T allele of rs742071 could greatly alter the affinity of Sin3Ak-20_disc4 (score: 2.1-14.1).

Rs742071 had the minor allele frequency as 0.04 among Han Chinese population, indicating it may need larger sample size to validate its significance; and with the limited sample size, rs742071 did not pass the threshold of the Bonferonni correction p value in the current study. We will add more samples to study it and other variants at PAX7 gene among NSCPO among Han Chinese population.

Rs2485893 was associated with CL/P among Asian ancestry with *p* value 7.86E-07 by GWAS ([@B8]). In this study, we found rs2485893 (10kb 3' of SYT14) was also found to be associated with NSCPO. Marked parent of origin effects were seen for rs2485893 alleles, over-transmission was seen preferentially from mothers compared with fathers ([Table 4](#T4){ref-type="table"}). Motif analyses by HaploReg indicated that the G allele of rs2485893 could greatly alter the affinity of AFP1 (score: 2.7-12.3).

Numerous studies have shown that highly conserved non-coding elements act as developmental enhancers *in vivo* ([@B19]-[@B21]). Non-coding conserved elements around rs742071 and rs2485893 therefore might represent putative regulatory elements for PAX7 and SYT14, we will perform functional studies to elucidate their roles in human craniofacial development.

Gene-gene interactions have been proposed as a potential source of the remaining heritability. Genotypic TDT for epistatic interactions showed that rs4844913 interacts with rs11119388 (SYT14) and rs6072081 interacts with rs6102085 for NSCPO, which provided new insights for the previous GWASs.

In summary, we replicated 38 SNPs contributing to NSCL/P to investigate their roles in NSCPO among Han Chinese population. In this study, we found that rs742071 and rs2485893 were associated NSCPO from Han Chinese population; also, interactions of rs4844913:rs11119388 and rs6072081:rs6102085 for NSCPO were identified, which may provide new insights for the previous GWASs.

We are particularly thankful to the patients and their families who participated in this study. This project was supported by the National Natural Science Foundation of China (No. 81600849 and 30530730) and Scientific Research Funds for Young Teachers of Sichuan University (No. 2015SCU11999). Shi-Jun Duan, Ning Huang, Zhong-Lin Jia, these authors contributed equally to this work.

[^1]: **Conflict of interest statement:**The authors have declared that no conflict of interest exist.
