Association of a BMP9 Haplotype with Ossification of the Posterior Longitudinal Ligament (OPLL) in a Chinese Population

Abstract

Direct or ex vivo BMP9 adenoviral gene therapy can induce massive bone formation at the injection sites and clearly promote spinal fusion. A comprehensive analysis of the osteogenic activity indicated that BMP9 was one of the most potent inducers of osteogenic differentiation both in vitro and in vivo among 14 types of human BMPs. However, genetic variations and whether they correlated with OPLL were not considered. We have sequenced the complete BMP9 gene in 450 patients with OPLL and in 550 matched controls. Analyses were performed on single markers and haplotypes. Single marker tests identified 6 SNPs, among which the minor alleles of rs7923671 (T>C; P = 0.0026; OR: 1.33, CI: 1.10–1.60), rs75024165 (C>T, Thr304Met; P<0.001; OR: 1.76, CI: 1.47–2.12) and rs34379100 (A>C; P<0.001; OR: 1.52, CI: 1.27–1.82) were associated with OPLL. Logistic regression analysis showed that the additive model of rs75024165 (TT vs. CT vs. CC; P<0.001; OR: 1.74) and rs34379100 (CC vs. AC vs. AA; P = 0.003; OR: 1.95) retained statistical significance when adjusted for clinical and demographic characteristics. Linkage disequilibrium (LD) analysis identified one 3 kb block of intense LD in BMP9 and one specific haplotype, CTCA (P<0.001; OR: 2.37), that contained the OPLL-associated risk alleles and was a risk factor for OPLL. This haplotype is associated with increased severity of OPLL, as shown by the distribution of ossified vertebrae in patients with OPLL (P = 0.001). In summary, in the Chinese population studied, SNPs in the BMP9 gene appear to contribute to the risk of OPLL in association with certain clinical and demographic characteristics. The severity of OPLL seems to be mediated predominantly by genetic variations in a 3kb BMP9 locus with the specific haplotype CTCA.

About the Authors:

Yuan Ren

Affiliation: Department of Orthopaedics, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Zhi-zhong Liu

Affiliation: Department of Clinic Laboratory, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Jie Feng

Affiliation: Department of Institute of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Hong Wan

Affiliation: Department of Institute of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Jun-hua Li

Affiliation: Department of Institute of Neurosurgery, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Hao Wang

Affiliation: Department of Orthopaedics, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Xin Lin

* E-mail: [email protected]

Affiliation: Department of Orthopaedics, Beijing Tiantan Hospital, Capital Medical University, Chongwen District, Beijing, China

Introduction

Ossification of the posterior longitudinal ligament (OPLL) is characterized by ectopic new bone formation replacement of ligamentous tissues of the spine and often leads to various degrees of myeloradiculopathy as a result of compression of the spinal cord and nerve roots [1], [2]. OPLL of the spine has been predominantly reported for Asian populations with a prevalence of 1.9%–4.3% in Japanese (aged >30 years) and 0.44%–8.92% in Chinese [3], [4], [5], [6], [7]. The etiology of OPLL is not clear up to date, although high body mass index (BMI), high peripheral bone mineral density (BMD), dietary habits, aging and glucose intolerance, have been proposed causes of OPLL [8], [9], [10], [11], [12]. On the basis of nationwide pedigree surveys, and surveys of twins and analysis of human leukocyte antigen haplotypes, genetic factors have also been suggested to be involved in the development of OPLL [13], [14], [15], [16]. Previous studies have shown that single nucleotide polymorphisms (SNPs) in the transforming growth factor-β, nucleotide pyrophosphatase (NPPS), leptin receptor, collagen 6A1(COL6A1), bone morphogenetic protein 2 (BMP2) and bone morphogenetic protein 4 (BMP4) as well as estrogen receptor (ER) and interleukin-1 (IL-1) genes are associated with the development of OPLL [17], [18], [19], [20], [21]. However, only BMPs were shown to induce ectopic ossification when these proteins are implanted subcutaneously [22]. Our research group has previously shown that BMP2 and BMP4 synonymous and non-synonymous SNPs are associated with the occurrence and severity of OPLL in a small Chinese population, although little evaluation of the clinical and demographic backgrounds of the participants was carried out [18], [36]. In this study we focused on the role of Bone morphogenetic protein 9 (BMP9), also known as growth and differentiation factor–2 (GDF-2) on the development of OPLL. BMP9 was originally identified in fetal mouse liver cDNA libraries, but it is a relatively uncharacterized member of the BMP family [23]. It has been reported that BMP9 may play a role in regulating glucose and iron homeostasis in liver [24], [25]. Interestingly, a comprehensive analysis of the osteogenic activities of the 14 types of human BMPs demonstrated that BMP9 was one of the most potent inducers of osteogenic differentiation both in vitro and in vivo [26], [27]. BMP9 was also shown to be a potent modulator of cartilage development in vitro and may play a role in gene therapy for bony defects, craniofacial bone repair, and nonunion bone fracture repair [28], [29], [30]. Additional studies have demonstrated that human mesenchymal stem cells (hMSCs) transfected with hBMP9 could lead to ectopic bone formation in rabbit and induce spinal fusion in rodents [31], [32], [33]. Most recently, Leblanc et al. found that injecting BMP9 protein into damaged muscle can directly induce muscle heterotopic ossification [34]. Up to know only one study of BMP9 response to novel activin receptor-like kinase 1 (ALK1) mutants derived from non synonymous SNPs in hereditary hemorrhagic telangiectasia type 2 (HHT2) patients has been reported [35] and no previous studies have assessed BMP9 variations. In the current study we hypothesized that BMP9 might be involved in OPLL, and we evaluated the overall polymorphisms of the BMP9 gene in association with certain clinical and demographic characteristics in a large group of Chinese patients with OPLL.

Materials and Methods

Participants

A case group that consisted of 450 Chinese patients with the common disease type of cervical OPLL and a control group of 550 healthy participants without cervical OPLL were enrolled in this study. All participants were recruited from Grade III, class A hospitals in mainland China. The diagnosis of OPLL was based on imaging findings, which included radiographs, computed tomography (CT) scans, and magnetic resonance imaging (MRI) of the cervical spine, in accordance with the criteria reported by Tsuyama [37]. The severity of OPLL was determined by the number of ossified cervical vertebrae observed using lateral radiographs and CT images. Patients with ankylosing spondylitis and metabolic diseases associated with OPLL, such as hypophosphatemic rickets/osteomalacia, osteoporosis, diffuse idiopathic skeletal hyperostosis (DISH), and hyperparathyroidism, were excluded from the study on the basis of radiographic and biochemical examinations. Participants who had taken drugs such as estrogen, progesterone, glucocorticoids, bisphosphonates, alfacalcidol, and calcitriol were also excluded.

Blood samples were collected from all participants. The study protocol was approved by the ethical committee of Beijing Tiantan Hospital, Capital Medical University with the following reference number: 200824 and written informed consent was obtained from all participants before the study. To evaluate the potential risk factors of OPLL, a standardized questionnaire [11] was used and related anthropometric and lifestyle determination methods are shown in Material S1.

Genotyping of BMP9

Genomic DNA was extracted from peripheral blood leukocytes using a Wizard Genomic DNA Purification Kit (Promega, Madison, WI, USA). The BMP9 gene was amplified by the polymerase chain reaction (PCR) with a standard protocol [38] in a total of 5 overlapping fragments and sequenced. The reactions were performed. Different thermocycling parameters were used for different primer pairs (Table S1). Sequencing reactions were performed using the Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA, USA ), and the extension products were analyzed on an ABI 3730XL POP7 DNA sequencing analysis 5.2 system (Applied Biosystems). All described data about polymorphisms were discovered in our study, but before our experiments, all these rs-ID (rs3758496, rs12252199, rs7923671, rs75024165, rs34379100, rs9421799) have also been reported in GeneBank.

Statistical Analysis

Hardy–Weinberg equilibrium and the genotypic and allelic distributions were evaluated using χ2 tests. The statistical power was calculated using Quanto v1.2.3 software. To control for potential confounding variables, we constructed a logistic regression model with the occurrence of OPLL as an outcome variable. The predictor variables included age, gender, sleeping habits, BMI, FBG, BMD, exercise, smoking, alcohol consumption, and the additive model for each SNP in BMP9. These clinical and demographic factors, which have been associated with the occurrence of OPLL in previous studies of other Asiatic populations, were re-evaluated in our study [8], [9], [10], [11], [12]. A backwards stepwise regression analysis (binary logistic regression) was performed to assess the independent effects of these covariates in the logistic regression model for OPLL and estimate the odds ratios (ORs) and 95% confidence intervals (95% CIs) for the two groups with regard to the risk of developing OPLL. A nonparametric test (Mann–Whitney exact test) was employed to compare the number of ossified cervical vertebrae between the two groups. The data are expressed as the means ± SD. SPSS 13 software (SPSS, Chicago, IL, USA) was used to analyze the data and calculate P values, odds ratios, and 95% CIs. Haplotype-based association studies were performed using both the Haploview V.3.32 and the SHESIS software platforms. Unless indicated otherwise, a P value less than 0.05 was considered statistically significant.

Results

Allele and Genotype Findings

The clinical and demographic characteristics of OPLL patients and the control group are shown in Table 1; the two groups did not differ significantly in each item, implying that they are comparable. The BMP9 genes of all participants were sequenced, and 6 SNPs were identified and genotyped (minor allele frequency in cases >0.05). The characteristics of the SNPs are given in Table S2. Allele and genotype frequencies of all SNPs were in Hardy–Weinberg equilibrium in both groups. Single-marker tests, which were based on an allelic model, revealed associations between an increased risk of OPLL and the minor alleles of rs7923671 (T>C; P = 0.0026; OR: 1.33,CI: 1.10–1.60), rs75024165 (C>T, Thr304Met; P<0.001; OR: 1.76, CI: 1.47–2.12), and rs34379100 (A>C; P<0.001; OR: 1.52,CI: 1.27–1.82).

[Figure omitted. See PDF.]

Table 1. Clinical and demographic characteristics.

https://doi.org/10.1371/journal.pone.0040587.t001

Table S2 and Table S3 show that the power of the minor allele frequencies was more than 90% for all SNPs with P<0.001. The power of the ossified vertebrae averages was more than 90% for related clinical and demographic characteristics with P<0.05 (except regular sleeping habit: 79.5%, normal exercise: 79.5%, and smoker: 87.9%).

The backward regression analysis, which was based on an additive model and adjusted for clinical and demographic characteristics, showed significant associations between the occurrence of OPLL and the additive model of rs75024165 (TT vs. CT vs. CC; P<0.001; OR: 1.82,CI: 1.42–2.32) and rs34379100 (CC vs. AC vs. AA; P = 0.003; OR:1.95,CI:1.15–1.96). Gender (male vs. female; P = 0.043; OR: 1.32, CI: 1.01–1.73) and FBG (≥7.0 vs. <7.0, mmol/L; P<0.001; OR: 5.42, CI: 3.60–8.15) showed statistically significant differences between the case and control groups (Table 2).

[Figure omitted. See PDF.]

Table 2. Analysis of the genotype distribution in BMP9 between the OPLL and control groups.

https://doi.org/10.1371/journal.pone.0040587.t002

Association between BMP9 Haplotypes and Occurrence of OPLL

We identified one 3 kb block of intense LD. This block encompassed SNPs rs9421799–rs7923671, and the linkage equilibrium coefficient D’ reduced gradually from the middle of the block toward rs9421799 and rs7923671 at the boundaries (Fig. 1). In this block, we identified several common single SNP haplotypes that covered together more than 60% of the individuals within both groups. Analysis of haplotype-based associations between the case and control groups found that the haplotype CTCA, which included the risk (*) alleles for two SNPs (rs7923671 C, rs75024165 *T, rs34379100 *C, and rs9421799 A) was associated strongly with OPLL (P<0.01; OR: 2.37, CI: 1.93–2.90) (Table 3).

[Figure omitted. See PDF.]

Figure 1. Genomic organization and linkage disequilibrium (LD) mapping of SNPs in a 6 kb genomic region within the BMP9 gene.

The left black bar indicates the exon 2 and the right black bar indicates the exon 1 coding area. Each box represents the D’ value of the LD (range from 0 to 1) between pairs of SNPs. Red, strong LD; grey, high D’/low LOD value; white, low D’.

https://doi.org/10.1371/journal.pone.0040587.g001

[Figure omitted. See PDF.]

Table 3. Haplotypes of the BMP9 SNPs in cases and controls.

https://doi.org/10.1371/journal.pone.0040587.t003

Association between BMP9 Haplotypes and Extent of OPLL Related to Clinical and Demographic Characteristics

The number of ossified vertebrae was analyzed in relation to the haplotype distribution in patients with OPLL who had been stratified by clinical and demographic characteristics. The results are shown in Fig. 2 and Table S3. Figure 2 shows two-dimensional computed tomography (2D-CT) scan pictures of OPLL patients with others (A) and the CTCA haplotype (B–H). As visible the degrees of ossifications in the CTCA haplotype group (B–H) is varying and the lower panel shows a scheme of the occurrence of ossification severities in the haplotype groups related to demographic and clinical characteristics. Patients with haplotype CTCA (71/450) developed a significantly (P = 0.001) greater number of ossified cervical vertebrae (3.41±1.58) than those who did not carry this haplotype (2.76±1.44). This effect was more obvious in patients with certain characteristics. The characteristics of smoker (P = 0.003), normal exercise (P = 0.007) and regular sleeping habit (P = 0.007) were associated with increased severity of OPLL in patients with the CTCA haplotype. In addition, normal FBG (P<0.001), middle age (P<0.001), regular drink (P = 0.002), high BMD (P<0.001), and being over weighted (P = 0.002) increased the extent of OPLL in patients with the CTCA haplotype significantly. The characteristics of male sex and high FBG might affect the initiation of OPLL significantly, in accordance with the findings of the logistic regression analysis, but only male sex (P = 0.001) was associated with an increase in the extent of heterotopic ossification in patients with the CTCA haplotype. Patients with the CTCA haplotype may suffer OPLL more severely as other haplotype patients, because they develop multiple cervical ossifications, but the severity is also related to other clinical and demographic characteristics. In addition, high FBG or low BMD seem to protect even CTCA haplotype patients against ossification. In summary, the CTCA haplotype is a major factor for ossification, but other factors as shown in Fig. 2 and Table S3 are contributors to the ossification severity.

[Figure omitted. See PDF.]

Figure 2. Association between BMP9 Haplotypes and Extent of OPLL related to clinical and demographic characteristics: The upper panel shows Two-dimensional computed tomography (2D-CT) pictures of OPLL severities in cervical spines (C2–7).

Patients A: other haplotypes, Patients B–H: haplotype CTCA. The Table shows ossification severity (B–H) in the CTCA and other haplotype groups related to demographic and clinical characteristics.

https://doi.org/10.1371/journal.pone.0040587.g002

Discussion

To confirm that the BMP9 locus is associated with the occurrence of OPLL, we sequenced the complete 6 kb BMP9 genomic region and identified 6 SNPs. A LD test showed one block of intense LD, which encompassed SNPs rs9421799–rs7923671 (Fig. 1). Haplotype block analysis identified the haplotype CTCA as a risk factor for OPLL; this haplotype included the risk alleles of rs75024165 (T) and rs34379100(C) (Table 3). Finally, the results of our study suggested that the CTCA haplotype was associated with increased severity of OPLL, especially in patients with certain clinical and demographic characteristics in a Chinese population (Fig. 2 and Table S3). As mentioned above, we identified a 3 kb block of intense LD that included the translated region of exon 2 of the BMP9 gene and contained four SNPs. Among the four SNPs in this region, the additive models of rs75024165 (TT vs. CT vs. CC) and rs34379100 (CC vs. AC vs. AA) were associated with the occurrence of OPLL. rs34379100 is located in the 3′ untranslated region (3′-UTR) of BMP9. We suggest here for the first time that this synonymous SNP might affect the degree of post-transcriptional regulation in OPLL patients, with a resulting protein over expression. rs75024165 C>T is located in the translated region of BMP9 exon 2. This non synonymous SNP codes an amino acid change: the T allele forms part of a methionine codon, whereas the C allele is part of a codon for threonine what might affect the BMP9 protein function. Radiographic and histological evaluations have demonstrated that overproduction of BMP6 and BMP9 protein induces a more robust and larger ossification mass than other BMPs in an orthotopic ossification animal model [26]. Dumont and Helm found that direct or ex vivo BMP9 adenoviral gene therapy could induce massive bone formation at the injection sites and clearly promote spinal fusion on CT films [32], [33]. However, additional studies are needed to understand the functional effects of the CTCA haplotype in vivo. In addition the study on the effect of BMP9 polymorphism should be extended to the ossification of the thoracic and lumbar spine. We did not find OPLL patients in the CTCA haplotype who suffered from high FBG in the case group. Chen et al. indicated that human BMP9 proteins reduced glycemia to near-normal levels after a subcutaneous injection in diabetic mice. Overproduction of BMP9 might cause the same effect in regulating blood glucose concentration in OPLL patients with the haplotype CTCA [24], but a further analysis of this negative correlation is necessary.

We conclude that the BMP9 gene variants which are associated with susceptibility to OPLL in our Chinese population are located between rs9421799 and rs7923671. Among the four SNPs in this region, we obtained supporting evidence for an association between the occurrence of OPLL and rs75024165 and rs34379100. Moreover patients who carried the CTCA haplotype had a higher number of ossified vertebrae than those who did not (Fig. 2), but after analyzing the distribution of the number of ossified vertebrae in relation to clinical and demographic characteristics, smoking, regular drinking, normal FBG, middle age, male, regular sleeping habit, high BMD, normal exercise, and being over weighted were additional significant factors for the ossification severity in CTCA haplotype patients.

In summary, this study provides evidence for an association between a specific polymorphism in the BMP9 gene and the occurrence of OPLL. The CTCA haplotype is not only associated with the occurrence of OPLL but also with increased severity of OPLL in Chinese patients, when it does occur in combination with other clinical and demographic characteristics. The results of the current study may contribute to the understanding of the background mechanisms and molecular etiology of OPLL, which should lead to further biomedical or functional studies and recommendations for a healthy lifestyle that will reduce the incidence of this disease.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Primer pairs and their annealing temperatures, which have been used to amplify 5 BMP9 genomic areas spanning SNPs rs3758496, rs12252199, rs7923671, rs75024165, rs34379100 and rs9421799.

https://doi.org/10.1371/journal.pone.0040587.s001

(TIF)

Table S2.

Map of the BMP9 gene region. The upper bar on the left picture indicates exon 1 and the lower bar exon 2. The locations and rs-IDs of the 6 SNPs are given relative to the Ensembl ENSG00000128802 chromosome 10 sequence, P values of the alleles in our cohort and the amino acid substitutions for non-synonymous SNPs are indicated.

https://doi.org/10.1371/journal.pone.0040587.s002

(TIF)

Table S3.

The distribution of ossified vertebrae in relation to the haplotype and clinical and demographic characteristics of the patients. Note: The data are expressed as [number of CTCA carriers vs. non carriers in case group (mean of ossified cervical vertebrae±SD)].

https://doi.org/10.1371/journal.pone.0040587.s003

(TIF)

Material S1.

https://doi.org/10.1371/journal.pone.0040587.s004

(DOC)

Acknowledgments

The authors thank the participants for making this study possible and the Institute of Neurosurgery of Beijing for technical support. Written consent for publication was obtained from the patient or their relative. All authors confirm that the content has not been published elsewhere and does not overlap with or duplicate their published work.

Author Contributions

Conceived and designed the experiments: XL. Performed the experiments: YR. Analyzed the data: YR Z-zL JF. Contributed reagents/materials/analysis tools: H. Wan H. Wang J-hL. Wrote the paper: YR.

Citation: Ren Y, Liu Z-z, Feng J, Wan H, Li J-h, Wang H, et al. (2012) Association of a BMP9 Haplotype with Ossification of the Posterior Longitudinal Ligament (OPLL) in a Chinese Population. PLoS ONE 7(7): e40587. https://doi.org/10.1371/journal.pone.0040587

References

1. Sakou T, Matsunaga S, Koga H (2000) Recent progress in the study of pathogenesis of ossification of the posterior longitudinal ligament. J Orthop Sci 5: 310–315.T. SakouS. MatsunagaH. Koga2000Recent progress in the study of pathogenesis of ossification of the posterior longitudinal ligament.J Orthop Sci5310315

2. Schmidt MH, Quinones-Hinojosa A, Rosenberg WS (2002) Cervical myelopathy associated with degenerative spine disease and ossification of the posterior longitudinal ligament. Semin Neurol 22: 143–148.MH SchmidtA. Quinones-HinojosaWS Rosenberg2002Cervical myelopathy associated with degenerative spine disease and ossification of the posterior longitudinal ligament.Semin Neurol22143148

3. Matsunaga S (2009) Updates on ossification of posterior longitudinal ligament. Epidemiology and pathogenesis of OPLL. Clin Calcium 19: 1415–1420.S. Matsunaga2009Updates on ossification of posterior longitudinal ligament. Epidemiology and pathogenesis of OPLL.Clin Calcium1914151420

4. Liu KC (1990) Epidemiological study on ossification of the posterior longitudinal ligament (OPLL) in the cervical spine–comparison of the prevalence between Japanese and Taiwanese. Nihon Seikeigeka Gakkai Zasshi 64: 401–408.KC Liu1990Epidemiological study on ossification of the posterior longitudinal ligament (OPLL) in the cervical spine–comparison of the prevalence between Japanese and Taiwanese.Nihon Seikeigeka Gakkai Zasshi64401408

5. ZGSC Li ZS (1999) The incidence of the ossification of the posterior longitudinal ligament in northern region of China. Chin J Spine Spinal Cord. pp. 285–286.ZS ZGSC Li1999The incidence of the ossification of the posterior longitudinal ligament in northern region of China.Chin J Spine Spinal Cord285286

6. Kim TJ, Bae KW, Uhm WS, Kim TH, Joo KB, et al. (2008) Prevalence of ossification of the posterior longitudinal ligament of the cervical spine. Joint Bone Spine 75: 471–474.TJ KimKW BaeWS UhmTH KimKB Joo2008Prevalence of ossification of the posterior longitudinal ligament of the cervical spine.Joint Bone Spine75471474

7. Inamasu J, Guiot BH, Sachs DC (2006) Ossification of the posterior longitudinal ligament: an update on its biology, epidemiology, and natural history. Neurosurgery 58: 1027–1039; discussion 1027–1039. J. InamasuBH GuiotDC Sachs2006Ossification of the posterior longitudinal ligament: an update on its biology, epidemiology, and natural history.Neurosurgery 58: 1027–1039; discussion 1027–1039

8. Shingyouchi Y, Nagahama A, Niida M (1996) Ligamentous ossification of the cervical spine in the late middle-aged Japanese men. Its relation to body mass index and glucose metabolism. Spine (Phila Pa 1976) 21: 2474–2478.Y. ShingyouchiA. NagahamaM. Niida1996Ligamentous ossification of the cervical spine in the late middle-aged Japanese men. Its relation to body mass index and glucose metabolism.Spine (Phila Pa 1976)2124742478

9. Kobashi G, Washio M, Okamoto K, Sasaki S, Yokoyama T, et al. (2004) High body mass index after age 20 and diabetes mellitus are independent risk factors for ossification of the posterior longitudinal ligament of the spine in Japanese subjects: a case-control study in multiple hospitals. Spine (Phila Pa 1976) 29: 1006–1010.G. KobashiM. WashioK. OkamotoS. SasakiT. Yokoyama2004High body mass index after age 20 and diabetes mellitus are independent risk factors for ossification of the posterior longitudinal ligament of the spine in Japanese subjects: a case-control study in multiple hospitals.Spine (Phila Pa 1976)2910061010

10. Okamoto K, Kobashi G, Washio M, Sasaki S, Yokoyama T, et al. (2004) Dietary habits and risk of ossification of the posterior longitudinal ligaments of the spine (OPLL); findings from a case-control study in Japan. J Bone Miner Metab 22: 612–617.K. OkamotoG. KobashiM. WashioS. SasakiT. Yokoyama2004Dietary habits and risk of ossification of the posterior longitudinal ligaments of the spine (OPLL); findings from a case-control study in Japan.J Bone Miner Metab22612617

11. Washio M, Kobashi G, Okamoto K, Sasaki S, Yokoyama T, et al. (2004) Sleeping habit and other life styles in the prime of life and risk for ossification of the posterior longitudinal ligament of the spine (OPLL): a case-control study in Japan. J Epidemiol 14: 168–173.M. WashioG. KobashiK. OkamotoS. SasakiT. Yokoyama2004Sleeping habit and other life styles in the prime of life and risk for ossification of the posterior longitudinal ligament of the spine (OPLL): a case-control study in Japan.J Epidemiol14168173

12. Yamauchi T, Taketomi E, Matsunaga S, Sakou T (1999) Bone mineral density in patients with ossification of the posterior longitudinal ligament in the cervical spine. J Bone Miner Metab 17: 296–300.T. YamauchiE. TaketomiS. MatsunagaT. Sakou1999Bone mineral density in patients with ossification of the posterior longitudinal ligament in the cervical spine.J Bone Miner Metab17296300

13. Matsunaga S, Yamaguchi M, Hayashi K, Sakou T (1999) Genetic analysis of ossification of the posterior longitudinal ligament. Spine (Phila Pa 1976) 24: 937–938; discussion 939. S. MatsunagaM. YamaguchiK. HayashiT. Sakou1999Genetic analysis of ossification of the posterior longitudinal ligament.Spine (Phila Pa 1976) 24: 937–938; discussion 939

14. Terayama K (1989) Genetic studies on ossification of the posterior longitudinal ligament of the spine. Spine (Phila Pa 1976) 14: 1184–1191.K. Terayama1989Genetic studies on ossification of the posterior longitudinal ligament of the spine.Spine (Phila Pa 1976)1411841191

15. Sakou T, Taketomi E, Matsunaga S, Yamaguchi M, Sonoda S, et al. (1991) Genetic study of ossification of the posterior longitudinal ligament in the cervical spine with human leukocyte antigen haplotype. Spine (Phila Pa 1976) 16: 1249–1252.T. SakouE. TaketomiS. MatsunagaM. YamaguchiS. Sonoda1991Genetic study of ossification of the posterior longitudinal ligament in the cervical spine with human leukocyte antigen haplotype.Spine (Phila Pa 1976)1612491252

16. Taketomi E, Sakou T, Matsunaga S, Yamaguchi M (1992) Family study of a twin with ossification of the posterior longitudinal ligament in the cervical spine. Spine (Phila Pa 1976) 17: S55–56.E. TaketomiT. SakouS. MatsunagaM. Yamaguchi1992Family study of a twin with ossification of the posterior longitudinal ligament in the cervical spine.Spine (Phila Pa 1976)17S5556

17. Kong Q, Ma X, Li F, Guo Z, Qi Q, et al. (2007) COL6A1 polymorphisms associated with ossification of the ligamentum flavum and ossification of the posterior longitudinal ligament. Spine (Phila Pa 1976) 32: 2834–2838.Q. KongX. MaF. LiZ. GuoQ. Qi2007COL6A1 polymorphisms associated with ossification of the ligamentum flavum and ossification of the posterior longitudinal ligament.Spine (Phila Pa 1976)3228342838

18. Wang H, Liu D, Yang Z, Tian B, Li J, et al. (2008) Association of bone morphogenetic protein-2 gene polymorphisms with susceptibility to ossification of the posterior longitudinal ligament of the spine and its severity in Chinese patients. Eur Spine J 17: 956–964.H. WangD. LiuZ. YangB. TianJ. Li2008Association of bone morphogenetic protein-2 gene polymorphisms with susceptibility to ossification of the posterior longitudinal ligament of the spine and its severity in Chinese patients.Eur Spine J17956964

19. Nakamura I, Ikegawa S, Okawa A, Okuda S, Koshizuka Y, et al. (1999) Association of the human NPPS gene with ossification of the posterior longitudinal ligament of the spine (OPLL). Hum Genet 104: 492–497.I. NakamuraS. IkegawaA. OkawaS. OkudaY. Koshizuka1999Association of the human NPPS gene with ossification of the posterior longitudinal ligament of the spine (OPLL).Hum Genet104492497

20. Kawaguchi Y, Furushima K, Sugimori K, Inoue I, Kimura T (2003) Association between polymorphism of the transforming growth factor-beta1 gene with the radiologic characteristic and ossification of the posterior longitudinal ligament. Spine (Phila Pa 1976) 28: 1424–1426.Y. KawaguchiK. FurushimaK. SugimoriI. InoueT. Kimura2003Association between polymorphism of the transforming growth factor-beta1 gene with the radiologic characteristic and ossification of the posterior longitudinal ligament.Spine (Phila Pa 1976)2814241426

21. Tahara M, Aiba A, Yamazaki M, Ikeda Y, Goto S, et al. (2005) The extent of ossification of posterior longitudinal ligament of the spine associated with nucleotide pyrophosphatase gene and leptin receptor gene polymorphisms. Spine (Phila Pa 1976) 30: 877–880; discussion 881. M. TaharaA. AibaM. YamazakiY. IkedaS. Goto2005The extent of ossification of posterior longitudinal ligament of the spine associated with nucleotide pyrophosphatase gene and leptin receptor gene polymorphisms.Spine (Phila Pa 1976) 30: 877–880; discussion 881

22. Wang EA, Rosen V, Cordes P, Hewick RM, Kriz MJ, et al. (1988) Purification and characterization of other distinct bone-inducing factors. Proc Natl Acad Sci U S A 85: 9484–9488.EA WangV. RosenP. CordesRM HewickMJ Kriz1988Purification and characterization of other distinct bone-inducing factors.Proc Natl Acad Sci U S A8594849488

23. Song JJ, Celeste AJ, Kong FM, Jirtle RL, Rosen V, et al. (1995) Bone morphogenetic protein-9 binds to liver cells and stimulates proliferation. Endocrinology 136: 4293–4297.JJ SongAJ CelesteFM KongRL JirtleV. Rosen1995Bone morphogenetic protein-9 binds to liver cells and stimulates proliferation.Endocrinology13642934297

24. Chen C, Grzegorzewski KJ, Barash S, Zhao Q, Schneider H, et al. (2003) An integrated functional genomics screening program reveals a role for BMP-9 in glucose homeostasis. Nat Biotechnol 21: 294–301.C. ChenKJ GrzegorzewskiS. BarashQ. ZhaoH. Schneider2003An integrated functional genomics screening program reveals a role for BMP-9 in glucose homeostasis.Nat Biotechnol21294301

25. Truksa J, Peng H, Lee P, Beutler E (2006) Bone morphogenetic proteins 2, 4, and 9 stimulate murine hepcidin 1 expression independently of Hfe, transferrin receptor 2 (Tfr2), and IL-6. Proc Natl Acad Sci U S A 103: 10289–10293.J. TruksaH. PengP. LeeE. Beutler2006Bone morphogenetic proteins 2, 4, and 9 stimulate murine hepcidin 1 expression independently of Hfe, transferrin receptor 2 (Tfr2), and IL-6.Proc Natl Acad Sci U S A1031028910293

26. Kang Q, Sun MH, Cheng H, Peng Y, Montag AG, et al. (2004) Characterization of the distinct orthotopic bone-forming activity of 14 BMPs using recombinant adenovirus-mediated gene delivery. Gene Ther 11: 1312–1320.Q. KangMH SunH. ChengY. PengAG Montag2004Characterization of the distinct orthotopic bone-forming activity of 14 BMPs using recombinant adenovirus-mediated gene delivery.Gene Ther1113121320

27. Luu HH, Song WX, Luo X, Manning D, Luo J, et al. (2007) Distinct roles of bone morphogenetic proteins in osteogenic differentiation of mesenchymal stem cells. J Orthop Res 25: 665–677.HH LuuWX SongX. LuoD. ManningJ. Luo2007Distinct roles of bone morphogenetic proteins in osteogenic differentiation of mesenchymal stem cells.J Orthop Res25665677

28. Abdelaal MM, Tholpady SS, Kessler JD, Morgan RF, Ogle RC (2004) BMP-9-transduced prefabricated muscular flaps for the treatment of bony defects. J Craniofac Surg 15: 736–741; discussion 742–734. MM AbdelaalSS TholpadyJD KesslerRF MorganRC Ogle2004BMP-9-transduced prefabricated muscular flaps for the treatment of bony defects.J Craniofac Surg 15: 736–741; discussion 742–734

29. Alden TD, Beres EJ, Laurent JS, Engh JA, Das S, et al. (2000) The use of bone morphogenetic protein gene therapy in craniofacial bone repair. J Craniofac Surg 11: 24–30.TD AldenEJ BeresJS LaurentJA EnghS. Das2000The use of bone morphogenetic protein gene therapy in craniofacial bone repair.J Craniofac Surg112430

30. Kimelman-Bleich N, Pelled G, Zilberman Y, Kallai I, Mizrahi O, et al. Targeted gene-and-host progenitor cell therapy for nonunion bone fracture repair. Mol Ther 19: 53–59.Kimelman-Bleich N, Pelled G, Zilberman Y, Kallai I, Mizrahi O, et al. Targeted gene-and-host progenitor cell therapy for nonunion bone fracture repair.Mol Ther195359

31. Jiang B, Li M, Cao YJ (2008) Ectopic bone formation of bone marrow mesenchymal stem cells with gene transfer of human bone morphogenetic protein-9 gene in rabbit. Sichuan Da Xue Xue Bao Yi Xue Ban 39: 723–727.B. JiangM. LiYJ Cao2008Ectopic bone formation of bone marrow mesenchymal stem cells with gene transfer of human bone morphogenetic protein-9 gene in rabbit.Sichuan Da Xue Xue Bao Yi Xue Ban39723727

32. Dumont RJ, Dayoub H, Li JZ, Dumont AS, Kallmes DF, et al. (2002) Ex vivo bone morphogenetic protein-9 gene therapy using human mesenchymal stem cells induces spinal fusion in rodents. Neurosurgery 51: 1239–1244; discussion 1244–1235. RJ DumontH. DayoubJZ LiAS DumontDF Kallmes2002Ex vivo bone morphogenetic protein-9 gene therapy using human mesenchymal stem cells induces spinal fusion in rodents.Neurosurgery 51: 1239–1244; discussion 1244–1235

33. Helm GA, Alden TD, Beres EJ, Hudson SB, Das S, et al. (2000) Use of bone morphogenetic protein-9 gene therapy to induce spinal arthrodesis in the rodent. J Neurosurg 92: 191–196.GA HelmTD AldenEJ BeresSB HudsonS. Das2000Use of bone morphogenetic protein-9 gene therapy to induce spinal arthrodesis in the rodent.J Neurosurg92191196

34. Leblanc E, Trensz F, Haroun S, Drouin G, Bergeron E, et al. (2011) BMP-9-induced muscle heterotopic ossification requires changes to the skeletal muscle microenvironment. J Bone Miner Res 26: 1166–1177.E. LeblancF. TrenszS. HarounG. DrouinE. Bergeron2011BMP-9-induced muscle heterotopic ossification requires changes to the skeletal muscle microenvironment.J Bone Miner Res2611661177

35. Ricard N, Bidart M, Mallet C, Lesca G, Giraud S, et al. (2010) Functional analysis of the BMP9 response of ALK1 mutants from HHT2 patients: a diagnostic tool for novel ACVRL1 mutations. Blood 116: 1604–1612.N. RicardM. BidartC. MalletG. LescaS. Giraud2010Functional analysis of the BMP9 response of ALK1 mutants from HHT2 patients: a diagnostic tool for novel ACVRL1 mutations.Blood11616041612

36. Meng XL, Wang H, Yang H, Hai Y, Tian BP, et al. (2010) T allele at site 6007 of bone morphogenetic protein-4 gene increases genetic susceptibility to ossification of the posterior longitudinal ligament in male Chinese Han population. Chin Med J (Engl) 123: 2537–2542.XL MengH. WangH. YangY. HaiBP Tian2010T allele at site 6007 of bone morphogenetic protein-4 gene increases genetic susceptibility to ossification of the posterior longitudinal ligament in male Chinese Han population.Chin Med J (Engl)12325372542

37. Tsuyama N (1984) Ossification of the posterior longitudinal ligament of the spine. Clin Orthop Relat Res. pp. 71–84.N. Tsuyama1984Ossification of the posterior longitudinal ligament of the spine.Clin Orthop Relat Res7184

38. Oste C (1988) Polymerase chain reaction. Biotechniques 6: 162–167.C. Oste1988Polymerase chain reaction.Biotechniques6162167