Journal of Peking University(Health Sciences) ›› 2019, Vol. 51 ›› Issue (3): 564-570. doi: 10.19723/j.issn.1671-167X.2019.03.028

Previous Articles     Next Articles

Exploring the association between SPRY gene family and non-syndromic oral clefts among Chinese populations using data of a next-generation sequencing study

Ren ZHOU1,Hong-chen ZHENG1,Wen-yong LI1,Meng-ying WANG1,Si-yue WANG1,Nan LI2,Jing LI3,Zhi-bo ZHOU2,Tao WU1,Hong-ping ZHU2△()   

  1. 1. Department of Epidemiology and Biostatistics, Peking University School of Public Health, Beijing 100191, China
    2. Department of Oral and Maxillofacial Surgery, Peking University School and Hospital of Stomatology & National Clinical Research Center for Oral Diseases & National Engineering Laboratory for Digital and Material Technology of Stomatology & Beijing Key Laboratory of Digital Stomatology, Beijing 100081, China
    3. Department of Pediatric Dentistry, Peking University School and Hospital of Stomatology, Beijing 100081, China
  • Received:2019-03-18 Online:2019-06-18 Published:2019-06-26
  • Supported by:
    Supported by the National Natural Science Foundation of China (81573225) and the Fundamental Research Funds for the Central Universities: Peking University Medicine Seed Fund for Interdisciplinary Research (BMU2017MX018)

Abstract: Objective: To explore the association between SPRY gene family and the risk of non-syndromic oral clefts among Chinese populations, in respect of single nucleotide polymorphisms (SNPs) association and parent-of-origin effects.Methods: Based on case-parent design, this study used the data of SPRY gene family in a next generation sequencing study of 183 non-syndromic cleft lip with or without cleft palate (NSCL/P) case-parent trios (549 participants) recruited from 2016 to 2018, to analyze the effects of SNP association and parent-of-origin. The sequencing study adopted a two-stage design. In the first stage, whole exome sequencing was conducted among 24 NSCL/P trios with family history to explore potential signals. Then in the second stage, another 159 NSCL/P trios were used as validation samples to verify the signals found in the first stage. The data of general information, disease features and parental environmental exposures for participants were collected through questionnaires. Blood samples were collected from each participant for DNA extraction and sequencing. Transmission disequilibrium tests (TDT) were conducted to test for the association between SNPs and NSCL/P, while Z score tests were applied to analyze parent-of-origin effects. The analyses were performed using Plink (v1.07). TRIO package in R (v3.5.1). Besides, famSKAT analyses were conducted in the first stage to combine the effect of SNPs located on the same gene, using famSKAT package in R(V3.5.1). Bonferroni method was adopted to correct multiple tests in the second stage. Results: Twenty-two SNPs in SPRY gene family were included for analyses after the quality control process in the first stage. Based on the variants annotation, functional prediction and statistical analysis, rs1298215244 (SPRY1) and rs504122 (SPRY2) were included in the second verification stage. TDTs in the verification stage revealed that rs1298215244: T>C, rs504122: G>C and rs504122: G>T were associated with the risk of NSCL/P after Bonferroni corrections, where rs504122: G>T was a rare variation. Although the test for parent-of-origin effect of rs1298215244: T>C reached nominal significance level, no SNP showed significant association with NSCL/P through parent-of-origin effect after Bonferroni corrections.Conclusion: This study found that SNPs (including both common and rare variants) among the SPRY gene family were associated with the risk of NSCL/P among Chinese populations. This study failed to detect parent-of-origin effects among the SPRY gene family.

Key words: Non-syndromic oral clefts, Association study, Case-parent trios, Next generation sequencing, SPRY gene family

CLC Number: 

  • R782.2

Table 1

The sex distribution of 183 NSCL/P cases in Chinese populations"

Stage Male Female Total
Stage 1 15 9 24
Stage 2 99 60 159
Total 114 69 183

Table 2

Transmission disequilibrium tests of SNPs in the SPRY gene family among 24 Chinese NSCL/P trios"

SNP Gene Chromosome Position Number of transmittedallele Number of non-transmitted allele P value
rs300574 SPRY1 4 124323738 18 9 0.083 3
rs4912844 SPRY4 5 141691371 1 5 0.102 5
rs4728 SPRY2 13 80910786 4 8 0.248 2
chr4: 124322600 SPRY1 4 124322600 1 0 0.317 3
rs1440299404 SPRY1 4 124322722 1 0 0.317 3

Table 3

The famSKAT analysis of SPRY genes among 24 Chinese NSCL/P trios"

phi Gene Number of alleles P value
0 SPRY1 5 1
SPRY2 5 0.533 2
SPRY4 12 0.754 2
0.2 SPRY1 5 1
SPRY2 5 0.587 5
SPRY4 12 0.918 3
0.5 SPRY1 5 1
SPRY2 5 0.662 8
SPRY4 12 0.964 0
0.8 SPRY1 5 1
SPRY2 5 0.676 8
SPRY4 12 0.962 2
1 SPRY1 5 1
SPRY2 5 1
SPRY4 12 0.943 6

Table 4

Transmission disequilibrium tests of potential signals in the SPRY gene family among 159 Chinese NSCL/P trios"

Allele Gene Variation MAF OR P value
rs1298215244 SPRY1 T>C 0.114 3.53 4.06×10-6
rs504122 SPRY2 G>C 0.116 2.45 4.81×10-4
rs504122 SPRY2 G>T 0.017 10.00 0.006 7

Table 5

Analyses of parent-of-origin effects for potential signals in the SPRY gene family among 159 Chinese NSCL/P trios"

Allele Gene Variation Paternal transmitted ∶untransmitted counts Maternal transmitted ∶untransmitted counts P value
rs1298215244 SPRY1 T>C 2.5 ∶25.5 12.5 ∶27.5 0.039 6
rs504122 SPRY2 G>C 11 ∶28 9 ∶21 0.870 6
rs504122 SPRY2 G>T 1 ∶3 0 ∶7 -
[1] Panamonta V, Pradubwong S, Panamonta M , et al. Global birth prevalence of orofacial clefts: a systematic review[J]. J Med Assoc Thai, 2015,98(Suppl 7):S11-S21.
[2] Dai L, Zhu J, Mao M , et al. Time trends in oral clefts in Chinese newborns: data from the Chinese National Birth Defects Monitoring Network[J]. Birth Defects Res A Clin Mol Teratol, 2010,88(1):41-47.
[3] Beaty TH, Murray JC, Marazita ML , et al. A genome-wide association study of cleft lip with andwithout cleft palate identifies risk variants near MAFB and ABCA4[J]. Nat Genet, 2010,42(6):525-529.
doi: 10.1038/ng.580
[4] Mangold E, Ludwig KU, Birnbaum S , et al. Genome-wide association study identifies two susceptibility loci for nonsyndromic cleft lip with or without cleft palate[J]. Nat Genet, 2010,42(1):24-26.
doi: 10.1038/ng.506
[5] Leslie EJ, Carlson JC, Shaffer JR , et al. A multi-ethnic genome-wide association study identifies novel loci for non-syndromic cleft lip with or without cleft palate on 2p24.2, 17q23 and 19q13[J]. Hum Mol Genet, 2016,25(13):2862-2872.
[6] Yu Y, Zuo X, He M , et al. Genome-wide analyses of non-syndromic cleft lip with palate identify 14 novel loci and genetic heterogeneity[J]. Nat Commun, 2017,8:14364-14374.
doi: 10.1038/ncomms14364
[7] Leslie EJ, Carlson JC, Shaffer JR , et al. Association studies of low-frequency coding variants in nonsyndromic cleft lip with or without cleft palate[J]. Am J Med Genet A, 2017,173(6):1531-1538.
doi: 10.1002/ajmg.a.38210
[8] Collins FS, Guyer MS, Charkravarti A . Variations on a theme: cataloging human DNA sequence variation[J]. Science, 1997,278(5343):1580-1581.
doi: 10.1126/science.278.5343.1580
[9] McCarthy MI, Abecasis GR, Cardon LR , et al. Genome-wide association studies for complex traits: consensus, uncertainty and challenges[J]. Nat Rev Genet, 2008,9(5):356-369.
doi: 10.1038/nrg2344
[10] Zuk O, Hechter E, Sunyaev SR , et al. The mystery of missing heritability: Genetic interactions create phantom heritability[J]. Proc Natl Acad Sci USA, 2012,109(4):1193-1198.
doi: 10.1073/pnas.1119675109
[11] Thisse B, Thisse C . Functions and regulations of fibroblast growth factor signaling during embryonic development[J]. Dev Biol, 2005,287(2):390-402.
doi: 10.1016/j.ydbio.2005.09.011
[12] Mason JM, Morrison DJ, Basson MA , et al. Sprouty proteins: multifaceted negative-feedback regulators of receptor tyrosine kinase signaling[J]. Trends Cell Biol, 2006,16(1):45-54.
doi: 10.1016/j.tcb.2005.11.004
[13] Stanier P, Pauws E . Development of the lip and palate: FGF signalling[J]. Front Oral Biol, 2012,16:71-80.
doi: 10.1159/000337618
[14] Ludwig KU, Mangold E, Herms S , et al. Genome-wide meta-analyses of nonsyndromic cleft lip with or without cleft palate identify six new risk loci[J]. Nat Genet, 2012,44(9):968-971.
doi: 10.1038/ng.2360
[15] Jia Z, Leslie EJ, Cooper ME , et al. Replication of 13q31.1 association in nonsyndromic cleft lip with cleft palate in Europeans[J]. Am J Med Genet A, 2015,167A(5):1054-1060.
[16] Moreno Uribe LM, Fomina T, Munger RG , et al. A population-based study of effects of genetic loci on orofacial clefts[J]. J Dent Res, 2017,96(11):1322-1329.
doi: 10.1177/0022034517716914
[17] Ward LD, Kellis M . HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants[J]. Nucleic Acids Res, 2012,40(Database issue):D930-D934.
doi: 10.1093/nar/gkr917
[18] Iyengar SFarnham PJ . KAP1 protein: an enigmatic master regulator of the genome[J]. J Biol Chem, 2011,286(30):26267-26276.
doi: 10.1074/jbc.R111.252569
[19] Hacohen N, Kramer S, Sutherland D , et al. sprouty encodes a novel antagonist of FGF signaling that patterns apical branching of the Drosophila airways[J]. Cell, 1998,92(2):253-263.
doi: 10.1016/S0092-8674(00)80919-8
[20] Gross I, Bassit B, Benezra M , et al. Mammalian sprouty proteins inhibit cell growth and differentiation by preventing ras activation[J]. J Biol Chem, 2001,276(49):46460-46468.
doi: 10.1074/jbc.M108234200
[21] Impagnatiello MA, Weitzer S, Gannon G , et al. Mammalian sprouty-1 and -2 are membrane-anchored phosphoprotein inhibitors of growth factor signaling in endothelial cells[J]. J Cell Biol, 2001,152(5):1087-1098.
doi: 10.1083/jcb.152.5.1087
[22] Yang X, Kilgallen S, Andreeva V , et al. Conditional expression of Spry1 in neural crest causes craniofacial and cardiac defects[J]. BMC Dev Biol, 2010,10:48-59.
doi: 10.1186/1471-213X-10-48
[23] Goodnough LH, Brugmann SA, Hu D , et al. Stage-dependent craniofacial defects resulting from Sprouty2 overexpression[J]. Dev Dyn, 2007,236(7):1918-1928.
doi: 10.1002/(ISSN)1097-0177
[24] Matsumura K, Taketomi T, Yoshizaki K , et al. Sprouty2 controls proliferation of palate mesenchymal cells via fibroblast growth factor signaling[J]. Biochem Biophys Res Commun, 2011,404(4):1076-1082.
doi: 10.1016/j.bbrc.2010.12.116
[25] Guilmatre A, Sharp AJ . Parent of origin effects[J]. Clin Genet, 2012,81(3):201-209.
doi: 10.1111/cge.2012.81.issue-3
[1] Meng-ying WANG,Wen-yong LI,Ren ZHOU,Si-yue WANG,Dong-jing LIU,Hong-chen ZHENG,Jing LI,Nan LI,Zhi-bo ZHOU,Hong-ping ZHU,Tao WU,Yong-hua HU. Evaluating the effect of WNT pathway genes considering interactions on the risk of non-syndromic oral clefts among Chinese populations [J]. Journal of Peking University (Health Sciences), 2020, 52(5): 815-820.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] Author. English Title Test[J]. Journal of Peking University(Health Sciences), 2010, 42(1): 1 -10 .
[2] . [J]. Journal of Peking University(Health Sciences), 2009, 41(2): 188 -191 .
[3] . [J]. Journal of Peking University(Health Sciences), 2009, 41(3): 376 -379 .
[4] . [J]. Journal of Peking University(Health Sciences), 2009, 41(4): 459 -462 .
[5] . [J]. Journal of Peking University(Health Sciences), 2010, 42(1): 82 -84 .
[6] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 319 -322 .
[7] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 333 -336 .
[8] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 337 -340 .
[9] . [J]. Journal of Peking University(Health Sciences), 2007, 39(3): 225 -328 .
[10] . [J]. Journal of Peking University(Health Sciences), 2007, 39(4): 346 -350 .