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
. 2022 Feb 2;13(1):634.
doi: 10.1038/s41467-022-28167-1.

Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology

Collaborators, Affiliations

Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology

Gyda Bjornsdottir et al. Nat Commun. .

Erratum in

  • Author Correction: Rare SLC13A1 variants associate with intervertebral disc disorder highlighting role of sulfate in disc pathology.
    Bjornsdottir G, Stefansdottir L, Thorleifsson G, Sulem P, Norland K, Ferkingstad E, Oddsson A, Zink F, Lund SH, Nawaz MS, Bragi Walters G, Skuladottir AT, Gudjonsson SA, Einarsson G, Halldorsson GH, Bjarnadottir V, Sveinbjornsson G, Helgadottir A, Styrkarsdottir U, Gudmundsson LJ, Pedersen OB, Hansen TF, Werge T, Banasik K, Troelsen A, Skou ST, Thørner LW, Erikstrup C, Nielsen KR, Mikkelsen S; DBDS Genetic Consortium; GO Consortium; Jonsdottir I, Bjornsson A, Olafsson IH, Ulfarsson E, Blondal J, Vikingsson A, Brunak S, Ostrowski SR, Ullum H, Thorsteinsdottir U, Stefansson H, Gudbjartsson DF, Thorgeirsson TE, Stefansson K. Bjornsdottir G, et al. Nat Commun. 2022 Apr 27;13(1):2419. doi: 10.1038/s41467-022-30129-6. Nat Commun. 2022. PMID: 35477715 Free PMC article. No abstract available.

Abstract

Back pain is a common and debilitating disorder with largely unknown underlying biology. Here we report a genome-wide association study of back pain using diagnoses assigned in clinical practice; dorsalgia (119,100 cases, 909,847 controls) and intervertebral disc disorder (IDD) (58,854 cases, 922,958 controls). We identify 41 variants at 33 loci. The most significant association (ORIDD = 0.92, P = 1.6 × 10-39; ORdorsalgia = 0.92, P = 7.2 × 10-15) is with a 3'UTR variant (rs1871452-T) in CHST3, encoding a sulfotransferase enzyme expressed in intervertebral discs. The largest effects on IDD are conferred by rare (MAF = 0.07 - 0.32%) loss-of-function (LoF) variants in SLC13A1, encoding a sodium-sulfate co-transporter (LoF burden OR = 1.44, P = 3.1 × 10-11); variants that also associate with reduced serum sulfate. Genes implicated by this study are involved in cartilage and bone biology, as well as neurological and inflammatory processes.

PubMed Disclaimer

Conflict of interest statement

G.B., L.S., G.T., P.S., K.N., E.F., A.O., F.Z., S.H.L., M.S.N., G.B.W., A.T.S., S.A.G., G.E., G.H.H., G.S., A.H., U.S., L.J.G., I.J., U.T., H.S., D.F.G., T.E.T., and K.S. are employees of deCODE Genetics/Amgen Inc. The remaining authors declare no competing interests.

Figures

Fig. 1
Fig. 1. Manhattan plots showing results for meta-analyses of Intervertebral disc disorders (M51), dorsalgia (M54) and lumbar discectomy (LDHsurg).
The P values (−log 10) from meta-analyses of the studied phenotypes are plotted (y-axis) against their respective positions on each chromosome (x-axis). a Intervertebral disc disorders IDD (M51), additive model (four cohorts; 58,854 cases, 922,958 controls), (b) Dorsalgia (M54), additive model (four cohorts, 119,110 cases, 909,847 controls), and (c) severe lumbar IDD defined by surgery (LDHsurg) (three cohorts; 9188 cases, 780,233 controls). P values are two sided and derived from a likelihood-ratio test. The gray and black dots represent SNPs not reaching genome-wide significance threshold weighted for variant impact. The yellow dots represent genome-wide significant SNPs and the red dots represent genome-wide significant SNPs with moderate or high impact (Methods).
Fig. 2
Fig. 2. Genes likely to associate with (a) IDD and (b) Dorsalgia.
Sequence variants associated with (a) IDD and (b) Dorsalgia for which functional evidence supports implication of genes in back pain. The variants listed are either protein-coding variants or affect mRNA expression (top cis-eQTL) as depicted by gray boxes (lof loss-of-function) (Supplementary Data 5–8). *Variants also associated in cis with mRNA of other genes (Supplementary Data 7,8). The meta-analyses were performed using logistic regression, the risk (odds ratio OR in yellow) of (a) IDD and (b) Dorsalgia are here shown for the risk-increasing allele and significance in blue. COL11A1 and GSDMC are not included in the figure as evidence for their association with back pain was derived differently as described in results.
Fig. 3
Fig. 3. Mendelian randomization (MR) analyses of the genetic relationship between IDD and dorsalgia in terms of causality.
a shows effects of variants associating with Dorsalgia at genome-wide significance, on IDD and dorsalgia. b shows the effects of variants associating with IDD at genome-wide significance, on IDD and dorsalgia. Effects are expressed as logarithms of odds ratios (log(OR)) and black crosses indicate 95% confidence intervals (CI) around effects. To avoid sample overlap, exposure effects are from the cohorts from Iceland (IS), Denmark (DK) and Finland (FIN), while outcome effects are from UK-Biobank. The dashed blue lines show the linear regression fit through the origin, weighting variants according to the square of the standard error of their effect estimates (also known as inverse-variance weighted, IVW)a)). The IVW-MR method is a multiplicative random effects model, where the test statistic is from a t-distribution, the test is two sided. No multiple comparison adjustments were made. The dashed red lines show the weighted linear regression fit not constrained to go through the origin (also known as MR Eggerb)). For the IDD variants, the slopes of both regression lines are different from zero. For the dorsalgia variants, the slope of the regression line (IVW) through the origin is different from zero, but not the slope of the unconstrained regression line (MR Egger). Further, the effects of the dorsalgia variants deviate substantially more from the regression lines than the IDD variants. These results are not sensitive to outlier removal (Methods, Supplementary Fig. 3).

Similar articles

Cited by

References

    1. Global Diseases and Injuries. Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019. Lancet. 2020;396:1204–1222. - PMC - PubMed
    1. Hoy D, et al. A systematic review of the global prevalence of low back pain. Arthritis Rheum. 2012;64:2028–2037. - PubMed
    1. Wu A, et al. Global low back pain prevalence and years lived with disability from 1990 to 2017: estimates from the Global Burden of Disease Study 2017. Ann. Transl. Med. 2020;8:299. - PMC - PubMed
    1. Hoy D, Brooks P, Blyth F, Buchbinder R. The epidemiology of low back pain. Best Pract. Res. Clin. Rheumatol. 2010;24:769–781. - PubMed
    1. Fatoye F, Gebrye T, Odeyemi I. Real-world incidence and prevalence of low back pain using routinely collected data. Rheumatol. Int. 2019;39:619–626. - PubMed

Publication types

MeSH terms

Supplementary concepts