Disc herniation is a complex spinal disorder associated with disability and high healthcare cost. Lumbar disc herniation is strongly associated with disc degeneration. Candidate genes of the aggrecan metabolic pathway may associate with the severity of lumbar disc herniation.

This study evaluated the association of single nucleotide variants (SNVs) of the candidate genes of the aggrecan metabolic pathway with the severity of lumbar disc herniation in patients with chronic mechanical low back pain. In addition, we assessed the in-silico functional analysis of the significant SNVs and association of their haplotypes with the severity of lumbar disc herniation.

A descriptive cross sectional study was carried out on 106 patients. Severity of disc herniation and disc degeneration were assessed on T2-weighted mid sagittal lumbar MRI scan. Sixty two exonic SNVs of ten candidate genes of aggrecan metabolic pathway (ACAN, IL1A, IL1B, IL6, MMP3, ADAMTS4, ADAMTS5, TIMP1, TIMP2 and TIMP3) were genotyped on a Sequenom MassARRAY iPLEX platform. Multivariable linear regression analysis was carried out using PLINK 1.9 software adjusting for age, gender, body mass index and severity of disc degeneration. Four online bioinformatics tools (Provean, SIFT, PolyPhen and Mutation Taster) were used for in-silico functional analysis.

Mean age was 52.42 ± 9.42 years and 69.8% were females. The mean severity of disc herniation was 2.81 ± 1.98. The rs2272023, rs35430524, rs2882676, rs2351491, rs938609, rs3825994, rs1042630, rs698621 and rs3817428 variants and their haplotypes of ACAN were associated with the severity of lumbar disc herniation. However, only the rs35430524, rs938609 and rs3817428 variants of ACAN were detected as pathogenic by in-silico functional analysis.

SNVs of ACAN and their haplotypes are associated with the severity of lumbar disc herniation. Functional genetic studies are necessary to identify the role of these significant SNVs in the pathogenesis of disc herniation.