Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2012 Dec;44(12):1360-4.
doi: 10.1038/ng.2463. Epub 2012 Nov 18.

A genome-wide association study identifies susceptibility loci for nonsyndromic sagittal craniosynostosis near BMP2 and within BBS9

Cristina M Justice  1 , Garima Yagnik, Yoonhee Kim, Inga Peter, Ethylin Wang Jabs, Monica Erazo, Xiaoqian Ye, Edmond Ainehsazan, Lisong Shi, Michael L Cunningham, Virginia Kimonis, Tony Roscioli, Steven A Wall, Andrew O M Wilkie, Joan Stoler, Joan T Richtsmeier, Yann Heuzé, Pedro A Sanchez-Lara, Michael F Buckley, Charlotte M Druschel, James L Mills, Michele Caggana, Paul A Romitti, Denise M Kay, Craig Senders, Peter J Taub, Ophir D Klein, James Boggan, Marike Zwienenberg-Lee, Cyrill Naydenov, Jinoh Kim, Alexander F Wilson, Simeon A Boyadjiev
Affiliations

A genome-wide association study identifies susceptibility loci for nonsyndromic sagittal craniosynostosis near BMP2 and within BBS9

Cristina M Justice et al. Nat Genet. 2012 Dec.

Abstract

Sagittal craniosynostosis is the most common form of craniosynostosis, affecting approximately one in 5,000 newborns. We conducted, to our knowledge, the first genome-wide association study for nonsyndromic sagittal craniosynostosis (sNSC) using 130 non-Hispanic case-parent trios of European ancestry (NHW). We found robust associations in a 120-kb region downstream of BMP2 flanked by rs1884302 (P = 1.13 × 10(-14), odds ratio (OR) = 4.58) and rs6140226 (P = 3.40 × 10(-11), OR = 0.24) and within a 167-kb region of BBS9 between rs10262453 (P = 1.61 × 10(-10), OR = 0.19) and rs17724206 (P = 1.50 × 10(-8), OR = 0.22). We replicated the associations to both loci (rs1884302, P = 4.39 × 10(-31) and rs10262453, P = 3.50 × 10(-14)) in an independent NHW population of 172 unrelated probands with sNSC and 548 controls. Both BMP2 and BBS9 are genes with roles in skeletal development that warrant functional studies to further understand the etiology of sNSC.

PubMed Disclaimer

Figures

Figure 1
Figure 1. Manhattan plot of P-values obtained from the genome-wide TDT of 130 trios (N=914,402)
The x-axis corresponds to the genomic position of the autosomes, while the y-axis is the −log10 of the P-value. Horizontal dashed line corresponds to the genome-wide significance threshold of P ≤ 5 × 10−8.
Figure 2
Figure 2. Regional association plots for associations with sNSC at genome-wide significance
Upper panel: P-values (-log10) of the GWAS (solid circle), replication (triangle) and meta-analysis (asterisk) on the y-axis are plotted against the genomic positions of each SNP on the x-axis on (a) chromosome 20 and on (b) chromosome 7. Horizontal dashed line corresponds to the genome-wide significance threshold of P ≤ 5 × 10−8. Genes in region are depicted below (not to scale). The LD (r2) between the lead SNP and other SNPs are indicated in different colors. The blue lines indicate the recombination rates in cM per Mb using HapMap controls. Lower panel: LD plots (based on r2 value) using genotyped subjects from the GWAS study.
Figure 2
Figure 2. Regional association plots for associations with sNSC at genome-wide significance
Upper panel: P-values (-log10) of the GWAS (solid circle), replication (triangle) and meta-analysis (asterisk) on the y-axis are plotted against the genomic positions of each SNP on the x-axis on (a) chromosome 20 and on (b) chromosome 7. Horizontal dashed line corresponds to the genome-wide significance threshold of P ≤ 5 × 10−8. Genes in region are depicted below (not to scale). The LD (r2) between the lead SNP and other SNPs are indicated in different colors. The blue lines indicate the recombination rates in cM per Mb using HapMap controls. Lower panel: LD plots (based on r2 value) using genotyped subjects from the GWAS study.
See this image and copyright information in PMC

References

    1. Cohen MM. Craniosynostosis: Diagnosis, Evaluation, and Management. Oxford University Press; New York: 2000.
    1. Johnson D, et al. A novel mutation, Ala315Ser, in FGFR2: a gene-environment interaction leading to craniosynostosis? Eur J Hum Genet. 2000;8(8):571–577. - PubMed
    1. Seto ML, et al. Isolated sagittal and coronal craniosynostosis associated with TWIST box mutations. Am J Med Genet A. 2007;143(7):678–86. - PubMed
    1. Weber I, et al. Molecular analyisis of 74 patients with craniosynostosis. Eur J Hum Genet. 2001;9(Sup.1):P0409, 179.
    1. Merrill AE, et al. Cell mixing at a neural crest-mesoderm boundary and deficient ephrin-Eph signaling in the pathogenesis of craniosynostosis. Hum Mol Genet. 2006;15(8):1319–28. - PubMed

Publication types