Lumbar disc degeneration (LDD) is associated with both genetic and environmental factors and affects many people worldwide. A hallmark of LDD is loss of proteoglycan and water content in the nucleus pulposus of intervertebral discs. While some genetic determinants have been reported, the etiology of LDD is largely unknown. Here we report the findings from linkage and association studies on a total of 32,642 subjects consisting of 4,043 LDD cases and 28,599 control subjects. We identified carbohydrate sulfotransferase 3 (CHST3), an enzyme that catalyzes proteoglycan sulfation, as a susceptibility gene for LDD. The strongest genome-wide linkage peak encompassed CHST3 from a Southern Chinese family-based data set, while a genome-wide association was observed at rs4148941 in the gene in a meta-analysis using multiethnic population cohorts. rs4148941 lies within a potential microRNA-513a-5p (miR-513a-5p) binding site. Interaction between miR-513a-5p and mRNA transcribed from the susceptibility allele (A allele) of rs4148941 was enhanced in vitro compared with transcripts from other alleles. Additionally, expression of CHST3 mRNA was significantly reduced in the intervertebral disc cells of human subjects carrying the A allele of rs4148941. Together, our data provide new insights into the etiology of LDD, implicating an interplay between genetic risk factors and miRNA.

Lumbar-disc herniation (LDH), one of the most common musculoskeletal diseases, has strong genetic determinants. Recently, several genes that encode extracellular matrix (ECM) proteins in the intervertebral disc have been reported to associate with LDH. Thrombospondins (THBSs) 1 and 2 are good candidates for the LDH susceptibility gene: They are intervertebral disc ECM proteins that regulate the effective levels of matrix metalloproteinases (MMPs) 2 and 9, which are key effectors of ECM remodeling. Here, we report that THBS2 is associated with LDH in Japanese populations. An intronic SNP in THBS2 (IVS10-8C --> T; rs9406328) showed significant association (p = 0.0000028) with LDH in two independent Japanese populations. This SNP, located in a polypyrimidine tract upstream of the 3' splice site of intron 10, exerts allelic differences on exon 11 skipping rates in vivo, with the susceptibility allele showing increased skipping. Skipping of exon 11 results in decreased THBS2 interaction with MMP2 and MMP9. Further, a missense SNP in MMP9 (Q279R; rs17576) is also strongly associated with LDH in the Japanese population (p = 0.00049) and shows a combinatorial effect with THBS2 (odds ratio 3.03, 95% confidence interval 1.58-5.77). Thus, a splicing-affecting SNP in THBS2 and a missense SNP in MMP9 are associated with susceptibility to LDH. Our data indicate that regulation of intervertebral disc ECM metabolism by the THBS2-MMP system plays an essential role in the etiology and pathogenesis of LDH.

No abstract