- Research article
- Open access
- Published:
Causal effects of skin microbiota on intervertebral disk degeneration, low back pain and sciatica: a two-sample Mendelian randomization study
Journal of Orthopaedic Surgery and Research volume 19, Article number: 476 (2024)
Abstract
Objective
The purpose of this study is to use two-sample Mendelian randomization (MR) to investigate the causal relationship between skin microbiota, especially Propionibacterium acnes, and intervertebral disc degeneration (IVDD), low back pain (LBP) and sciatica.
Methods
We conducted a two-sample MR using the aggregated data from the whole genome-wide association studies (GWAS). 150 skin microbiota were derived from the GWAS catalog and IVDD, LBP and sciatica were obtained from the IEU Open GWAS project. Inverse-variance weighted (IVW) was the primary research method, with MR-Egger and Weighted median as supplementary methods. Perform sensitivity analysis and reverse MR analysis on all MR results and use multivariate MR to adjust for confounding factors.
Results
MR revealed five skin microbiota associated with IVDD, four associated with LBP, and two with sciatica. Specifically, P.acnes in sebaceous skin environments were associated with reduced risk of IVDD; IVDD was found to increase the abundance of P.acnes in moist skin. Furthermore, ASV010 [Staphylococcus (unc.)] from dry skin was a risk factor for LBP and sciatica; ASV045 [Acinetobacter (unc.)] from dry skin and Genus Rothia from dry skin exhibited potential protective effects against LBP; ASV065 [Finegoldia (unc.)] from dry skin was a protective factor for IVDD and LBP. ASV054 [Enhydrobacter (unc.)] from moist skin, Genus Bacteroides from dry skin and Genus Kocuria from dry skin were identified as being associated with an increased risk of IVDD. Genus Streptococcus from moist skin was considered to be associated with an increased risk of sciatica.
Conclusions
This study identified a potential causal relationship between skin microbiota and IVDD, LBP, and sciatica. No evidence suggests skin-derived P.acnes is a risk factor for IVDD, LBP and sciatica. At the same time, IVDD can potentially cause an increase in P.acnes abundance, which supports the contamination theory.
Introduction
Low back pain (LBP) is a widespread condition that carries a significant burden, causing considerable disability and frequent absenteeism from work globally [1, 2]. According to statistics, approximately 80% of all people experience back pain at some point [3]. Intervertebral disc degenerative IVDD is a common spinal degenerative disease characterized by a gradual decrease of proteoglycan content and water content in the nucleus pulposus (NP) [4]. This pathological process may ultimately lead to the rupture of the intervertebral disc, making the disc more prone to protrusion. When such protrusion compresses the spinal cord or nerve roots, it can cause symptoms of LBP [4]. Notably, when the nerve roots in the L4-S1 region are compressed, patients often exhibit sciatica symptoms. Therefore, IVDD is not only an essential pathological basis for lower back pain (LBP) but also a critical factor in causing sciatica [5].
Research has shown the presence of microorganisms in intervertebral discs, and the abundance of bacteria varies between healthy and degenerative intervertebral discs [6]. The presence of low-virulence bacteria may be caused by colonization and can cause subclinical intervertebral disc infections, mediating IVDD and LBP [7, 8]. Currently, the vast majority of studies indicate that the primary pathogen of this low-virulence bacterium is Propionibacterium acnes [9], a Gram-positive, oxygen-tolerant anaerobic bacterium and an average bacterial population distributed in the skin. It belongs to an opportunistic pathogen that can cause infections outside of the skin, such as in bones and joints [10]. In 2001, Stirling et al. [11]. first detected P.acnes in degenerated intervertebral disc tissue and proposed that P.acnes may be one of the causes of intervertebral disc degeneration. Subsequently, multiple studies have explored the connection between them. A meta-analysis study reviewed 4109 patients and found that the low-grade bacterial infection rate was 30% (24 -37%), with P.acnes being the most common, accounting for 66.7% of bacterial-positive intervertebral discs [12]. Therefore, many researchers support the correlation between P.acnes and IVDD, LBP, and sciatica. However, some researchers oppose the impact of low-virulence bacteria on intervertebral discs and propose contamination theories. In short, current surgical or radiological guidance cannot wholly separate intervertebral disc tissue [13]. P.acnes can be cultured not only in intervertebral disc tissue, but also in muscle tissue, ligamentum flavum, and even operating room air [14], resulting in symbiotic bacterial contamination in the skin or muscles. A meta-analysis suggests that there is currently insufficient evidence to attribute low-virulence organisms to IDD, LBP, and disability [15]. There is still controversy over the relevant theories at present.
P.acnes, as a symbiotic bacterium located on the skin or oral cavity [10], may invade the circulatory system and colonize intervertebral discs when the skin or mucosa ruptures [9, 16]. Similarly, other skin microbiota, such as coagulase-negative Staphylococcus (CoNS) [17], may colonize intervertebral discs or indirectly affect intervertebral disc degeneration through a similar process. Rajasekaran et al. [6]. proposed the concept of the skin-spine microbiota axis. They found 29 common bacteria between the skin and IVD microbiota, indicating that the skin microbiota may be the source of the intervertebral disc microbiota. However, the impact and potential mechanisms of these skin microbiota on IVDD, LBP and sciatica are still unclear, and the intermediate processes are complex and susceptible to contamination theory, confounding factors and reverse causality [18]. Mendelian randomization (MR) provides an effective method for inferring causal relationships by utilizing genetic variation as an instrumental variable (IV). It imitates the causal reasoning of randomized controlled trials, effectively minimizing the influence of confounding factors and reversing the causal bias often encountered in observational studies [19, 20]. Therefore, MR provides a precise tool to address the limitations above. Our study utilized dual sample MR to investigate the causal relationship between skin microbiota, particularly Propionibacterium acnes, and IVDD, LBP, and sciatica.
Methods
Study design
This study complies the STROBE-MR guidelines [21]. Figure 1 illustrates the study’s overall framework. Our analysis used two-sample MR to assess the potential causal relationships between skin microbiota and IVDD, LBP and sciatica. Mendelian randomization depends on three crucial assumptions: firstly, that the chosen instrumental variables exhibit a robust association with the exposures being studied; secondly, that these instrumental variables are not influenced by confounding factors that might impact both exposures and outcomes; and thirdly, that there is no direct link between the instrumental variables and the outcomes themselves [21].
Data sources
The summary statistics for 150 skin microbiota were derived from the GWAS catalog (GWAS ID from GCST90133164 to GCST90133310). This extensive study encompassed the sequencing data of the 16 S rRNA gene, collected from meta-analyses of two population-based GWAS in Germany, including 1656 skin samples taken from dry skin (dorsal and volar forearm), moist (antecubital fossa) and sebaceous (retroauricular folds and forehead) skin microenvironment. The skin microbiota encompassed three phyla, four classes, seven orders, seven families, fifteen genera and forty-three ASV [22]. The GWAS summary data for IVDD (ncase = 20,001, ncontrol = 164,682), LBP (ncase = 13,178, ncontrol = 164,682) and sciatica (ncase = 19,509, ncontrol = 199,283) were obtained from the IEU Open GWAS project. Table 1 shows detailed information on the data sources.
Instrumental variables selection and data harmonization
The following steps filtered the instrumental variables: (1) SNPs were selected based on their genome-wide significance thresholds significantly associated with skin microbiota (P < 5 × 10− 8). However, the number of instrumental variables for some skin microbiota was less than 3 and was no longer suitable for further MR analysis. Therefore, for this part of exposre, we used a liberal significance threshold (5 × 10− 6) [23]; (2) linkage disequilibrium was removed from the selected SNPs, with an r2 threshold set to < 0.001 and a clustering distance of 10000kb [24]; (3) palindromic SNPs were excluded through data harmonization [25]; (4) the F-statistic was computed for each instrumental variable (IV), and those having an F-statistic value below ten were excluded to reduce weak instrument bias. The calculation of the F-statistic was based on the formula: F = r2 (n-2)/(1-r2) [26]. (5) The LDtrait tool (https://ldlink.nih.gov/?tab=ldtrait) was used to examine each IV and excluded SNPs associated with confounding factors [27]. Finally, we obtained the required instrumental variables.
Mendelian randomization analysis
We used a two-sample MR to assess the effects of skin microbiota on IVDD, LBP, and sciatica. The primary research method used was inverse variance weighting (IVW). Without directional pleiotropy, the IVW method can provide a relatively stable and accurate causal assessment by combining the Wald estimates of each IV [28]. When the Cochran Q test suggests the existence of potential heterogeneity, a random-effects IVW model is employed. Conversely, in cases where significant heterogeneity is absent, a fixed effects IVW model is utilized [29]. In addition, MR-Egger and weighted median serve as supplementary methods for estimating causal effects. MR-Egger is based on the assumption of InSIDE. This method can provide relatively robust estimates and evaluate horizontal pleiotropy through intercept terms. Weighted median can provide a robust estimate of causal effects, even with up to 50% of invalid IVs [28]. When the direction of IVW is consistent with that of the supplementary method, the results are stable [21]. To investigate the direct effects of skin microbiota on IVDD, LBP and sciatica, we used multivariate Mendelian randomization (MVMR), an MR method that allows for the co-detection of multiple risk factors, to adjust for potential confounding effects of BMI and smoking status [30]. Given the intimate association between Propionibacterium acnes and IVDD, LBP, and sciatica, we further investigated the impact of outcomes on exposure.
Sensitivity analysis
We conducted a sensitivity analysis on all obtained MR results, including heterogeneity and horizontal pleiotropy analysis. Specifically, Cochran’s Q statistics for the IVW method were used to evaluate whether MR results had heterogeneity, and when Q_pval was less than 0.05, it indicated the existence of heterogeneity [29]. The MR Egger intercept test was employed to detect the presence of horizontal pleiotropy in the MR results, where P < 0.05 indicated horizontal pleiotropy [31]. Additionally, the MR-PRESSO test was utilized to identify and exclude outliers, an additional means to detect horizontal pleiotropy. After excluding outliers, MR analysis will be conducted again [32]. The Leave-one-out sensitivity analysis was to remove each SNP one by one and subsequently perform the MR analysis again to identify influential outliers [33]. The results of these sensitivity analyses were used to evaluate the reliability of MR results.
All MR results were obtained using the “TwosampleMR” package in the R studio software.
Result
Instrumental variables selection
After the above screening steps, 553 (genome-wide significance threshold, P < 5 × 10− 8) and 994 (liberal significance threshold, P < 5 × 10− 6) SNPs were selected as instrumental variables for skin microbiota, and all F-statistics were greater than 10 (Supplementary Tables S1).
At the genome-wide significance threshold (5 × 10− 8), Causal effects of skin microbiota on IVDD LBP and sciatica
Using IVW analysis as our primary approach, we identified 6 skin microbiota taxa (five Amplicon Sequence Variant and one genus) that have a causal relationship with IVDD, LBP, and sciatica at the genome-wide significance threshold (5 × 10− 8). Specifically, ASV065 [Finegoldia (unc.)] from dry skin (OR = 0.968, 95%CI = 0.945 to 0.992, P = 0.008) was associated with a reduced risk of IVDD. Besides, ASV045 [Acinetobacter (unc.)] from dry skin (OR = 0.977, 95%CI = 0.957 to 0.997, P = 0.023), ASV065 [Finegoldia (unc.)] from dry skin (OR = 0.976, 95%CI = 0.957 to 0.995, P = 0.013) and Genus Rothia from dry skin (OR = 0.966, 95%CI = 0.933 to 1.000, P = 0.049) were related to a reduced risk of LBP. ASV010 [Staphylococcus (unc.)] from dry skin (OR = 1.038, 95%CI = 1.015 to 1.062, P = 0.001) was thought to be linked to an elevated risk of LBP. Moreover, ASV010 [Staphylococcus (unc.)] from dry skin (OR = 1.024, 95%CI = 1.005 to 1.044, P = 0.014) were considered to be associated with an increased risk of sciatica (Fig. 2). The scatter plots illustrates the causal relationships between 6 skin microbiota and IVDD, LBP and sciatica (Figure S1). The reverse MR indicates no evidence to suggest a reverse causal relationship between the above MR result (Supplementary Tables S3).
At the liberal significance threshold (5 × 10− 6), Causal effects of skin microbiota on IVDD LBP and sciatica
IVW analysis indicated that 5 skin microbiota taxa (two Amplicon Sequence Variants and three Genus) were causally related to IVDD, LBP and sciatica at the liberal significance threshold (5 × 10− 6). ASV001 [P.acnes] from sebaceous skin environments (OR = 0.955, 95%CI = 0.914 to 0.997, P = 0.038) were associated with reduced risk of IVDD. ASV054 [Enhydrobacter (unc.)] from moist skin (OR = 1.032, 95%CI = 1.009 to 1.056, P = 0.007), Genus Bacteroides from dry skin (OR = 1.020, 95%CI = 1.010 to 1.031, P < 0.001) and Genus Kocuria from dry skin (OR = 1.009, 95%CI = 1.001 to 1.018, P = 0.034) were identified as being association with an increased risk of IVDD. Besides, Genus Streptococcus from moist skin (OR = 1.035, 95%CI = 1.004 to 1.067, P = 0.026) was considered to be associated with an increased risk of sciatica (Fig. 3). The scatter plots illustrate the causal relationships between 5 skin microbiota and IVDD, LBP and sciatica (Figure S2). The reverse MR indicates no evidence to suggest a reverse causal relationship between the above MR result. (Supplementary Tables S3)
Causal effects of Propionibacterium acnes on IVDD, LBP and Sciatica
Numerous studies have implicated Propionibacterium acnes in the pathogenesis of IVDD, LBP and sciatica, so we further evaluated their association. The IVW analysis showed that P.acnes in sebaceous skin environments (OR = 0.955, 95%CI = 0.914 to 0.997, P = 0.038) were associated with reduced risk of IVDD. In addition, other MR results were not significant. Therefore, there was no evidence to suggest skin-derived P.acnes was a risk factor for IVDD, LBP, and sciatica (Fig. 4). However, when exploring the causal effects of IVDD, LBP and sciatica on P.acnes, we found that IVDD increased the abundance of P.acnes derived from moist skin (OR = 1.800, 95%CI = 1.131 to 2.866, P = 0.013, Fig. 5).
MVMR
To determine whether the statistically significant skin microbiota we had observed directly or indirectly affected the risk of IVDD, LBP and sciatica through confounding factors, we had conducted additional MVMR analysis to adjust for the effects of BMI and smoking. We found that after adjusting for BMI and smoking, the impact of Genus Bacteroides and Kocuria from dry skin on IVDD were still significant. In contrast, the impact of ASV065 [Finegoldia (unc.)] from dry skin on IVDD and various skin microbiota on LBP and sciatica were no longer significant. In addition, the impact of P.acnes in sebaceous skin environments on IVDD was attenuated with adjustment of BMI (Fig. 6).
Sensitivity analysis
We conducted a sensitivity analysis on all MR results. Specifically, Cochran’s Q test using the IVW method revealed heterogeneity in the MR results of ASV065 [Finegoldia (unc.)] from dry skin on IVDD (Q_pval = 0.038). Furthermore, the MR-Egger intercept and MR-PRESSO test demonstrated no horizontal pleiotropy in all MR results (Supplementary Tables S2). In addition, the MR leave-one-out sensitivity analysis revealed that all MR results remained stable after excluding IVs one by one (Figure S3-S4). The above sensitivity analysis results indicate that the MR results are stable and reliable.
Discussion
In this study, we first reported the relationship between skin microbiota, particularly P.acnes, IVDD, LBP, and sciatica, using two-sample Mendelian randomization. Our research indicated five skin microbiota were associated with IVDD, four with LBP, and two with sciatica. Specifically, IVDD was found to increase the abundance of P.acnes in moist skin. The sensitivity analysis results indicate that our MR results are stable and reliable, and the reverse MR results indicate that there is no reverse causal relationship between our MR results.
Propionibacterium acnes plays a pivotal role in maintaining skin health and homeostasis. However, they are also opportunistic pathogens closely associated with acne vulgaris [34]. When it infiltrates intervertebral discs, it can lead to IVDD, LBP and sciatica [10]. Senker et al. [35] found that 38.78% of samples cultured for intervertebral disc bacteria were positive, with P.acnes (18.05%) being one of the most commonly detected microorganisms. Albert et al. [36]. found that among 61 patients with lumbar disc herniation, 46% of them found microorganisms in the IVD tissue, with the most common being P.acnes. In addition, a randomized, controlled, double-blind trial was conducted on 162 patients with discogenic back pain, and the relief of back pain in the antibiotic treatment group was significantly higher than that in the placebo group [37]. Injecting p.acnes into the intervertebral discs of mice and rabbits can cause intervertebral disc modic changes [38, 39]. Therefore, P.acnes is likely associated with intervertebral disc degeneration.
Propionibacterium acnes may enter the intervertebral disc through multiple pathways. One possible method is for P.acnes to enter through the circulatory system [36]. P.acnes have the opportunity to invade the bloodstream when the skin or mucous membranes are damaged [9, 16]. Due to the aerobic environment of human blood, P.acnes has not shown any harm to the human body at this time. When the annulus fibrosus breaks, neovascularisation will provide a potential pathway for P.acnes to enter the nucleus pulposus [40]. The lack of blood vessels, sufficient water, low oxygen pressure, and pH-neutral “biofilm-like” structure in the intervertebral disc provides a suitable environment for the survival of P. acnes [17]. Another possible pathway is through the action of macrophages. Macrophages may engulf Propionibacterium acnes and accumulate around the ruptured annulus fibrosus [41, 42]. Finally, P.acnes entering the intervertebral disc may also originate from iatrogenic infections. Surgery or injections may introduce bacteria into the intervertebral disc [43].
Although there are multiple possible pathways, further research is needed to clarify the exact mechanism by which P.acnes enters intervertebral discs. When the abundance of P.acnes distributed on the skin increases, regardless of the potential mechanism through which P.acnes enters the intervertebral disc. It is possible to increase the risk of its colonization in the intervertebral disc. However, we concluded that P.acnes derived from the skin is not a risk factor for intervertebral disc degeneration. In other words, skin-derived P.acnes may be challenging to colonize intervertebral discs or may not significantly promote intervertebral disc degeneration. However, we cannot completely deny the impact of P.acnes on intervertebral disc degeneration. This is because P.acnes is widely distributed and may have other sources or enter the intervertebral disc through other means. For example, P.acnes is a resident bacterial community in the oral cavity [10]. Brushing teeth can cause P.acnes to enter the bloodstream, forming bacteremia [44]. The nucleus pulposus is in an anaerobic environment, and when the intervertebral disc degenerates, protrudes, or the annulus fibrosus ruptures, neovascularisation can be formed locally, becoming the entrance for P.acnes to enter the degenerated intervertebral disc [40]. P.acnes are classified into types I, II, and III, with type I further divided into four subtypes: IA1, IA2, IB, and IC [45]. Rollason et al. [46]. found a positive rate of 38% for P.acnes in 64 intervertebral disc specimens, with type II P.acnes accounting for the majority. However, the P.acnes that colonizes the skin are mainly type I, further indicating that these P.acnes may not originate from the skin.
Our research indicated that P.acnes derived from sebaceous glands can have a protective effect on IVDD. It may cause adaptive immunity, making it more difficult for P.acnes from other sources or other pathogens to evade the immune system and colonize the IVD [47]. In addition, studies have shown that the relative abundance of P.acnes in the skin of patients with acne vulgaris is comparable to that of healthy individuals. However, there are differences in classification [48]. Type IA1 is the primary strain type in the hair follicles of acne patients [49]. Different types of P.ances play various roles in IVDD. However, the GWAS summary data of this study did not strictly classify P.ances, which may be another reason why P.ances was not related to IVDD, LBP, and sciatica in this study.
Furthermore, numerous studies have proposed that the detection of P. acnes in intervertebral discs may be attributed to skin contaminants [50]. Carricajo et al. [14]. found only 2 cases of P.acnes infection in 54 lumbar disc tissue testing cases. Still, the detection rate of Propionibacterium acnes was as high as 29.6% in the control group samples of operating room air, yellow ligament, and muscle. Out of 120 samples obtained under strict aseptic conditions, 116 were sterile, 4 tested positive for coagulase-negative Staphylococcus aureus, and no samples tested positive for P.ances, supporting the contamination theory [51]. A recent prospective cohort study showed that less than 2% of patients with intervertebral disc herniation can detect the presence of infectious bacteria. P.acnes was not detected in any intervertebral disc tissue, and there was no significant correlation between positive tissue infection and any clinical prognostic factors [52]. Using anterior surgery with relatively less contamination, 379 intervertebral disc samples were cultured, and only 2 cases were cultured for P.acnes, while 4 cases were cultured for other bacteria, suggesting contamination [53]. Our research adds to this discourse, revealing that IVDD promotes an increased abundance of Propionibacterium acnes on the skin through specific undisclosed mechanisms. This finding supports the contamination theory, suggesting that the high abundance of P.acnes detected in cultures may originate from the skin.
Coagulase-negative staphylococci, such as S.epidermidis, S.capitis and S.hominis [54, 55], are the second most commonly detected group of bacteria after P.acnes [17]. Arndt et al. [56]. found that 16 out of 40 intervertebral discs tested positive for coagulase-negative staphylococci. Due to the similar characteristics of CoNS and P.acnes, they may implant into intervertebral discs resembling biofilms through a similar process [17]. In our study, ASV010 [Staphylococcus (unc.)] was identified as a risk factor for lower back pain (LBP) and sciatica. This bacterial community may belong to CoNS [22], which confirms previous research findings. S.epidermidis is the most common CoNS, which appropriately regulates protective immune responses when encountering pathogens [55]. The protease Esp produced by it can interfere with forming Staphylococcus aureus biofilm but can also degrade human complement component C5 and fibrinogen [57]. S.epidermidis can cause host diseases through virulence factors such as proteases (serine and cysteine proteases), lipases, and hemolysins [54]. After injecting S.epidermidis into the cervical intervertebral disc of mice, the intervertebral discs and subchondral bone structures are destroyed, and bone fusion occurs after discitis [58]. S.hominis is the second most common CoNS in healthy human skin [59] and is considered as a beneficial commensal bacterium that plays a vital role in maintaining skin health [60]. S.capitis infection has been described as associated with prosthetic joint infection [61]. However, although these bacterial communities belong to CoNS, we found no association between them and IVDD, LBP or sciatica in our study. There are numerous types of CoNS, and future research should attempt to distinguish the effects of different bacterial genera in CoNS on IVDD, as they may play various roles.
Furthermore, a study revealed the presence of Enhydroxide, Acinetobacter, Kocuria, Rothia and Streptococcus in IVD tissues [6]. The abundance of Enhydroactor is associated with psoriasis and diabetic foot ulcer [62, 63]. Our research indicates that Enhydroxybacter has the potential to increase the risk of IVDD. Unclassified Acinetobacter from dry skin environments exhibit potential protective effects against LBP and sciatica. Notably, Acinetobacter infrequently causes human infections and belongs to the aerobic bacterial group [64], so its colonization in IVD may be difficult. Interestingly, Acinetobacter johnsonii negatively correlates with the secretion of IL-4 and TNF-α [65], suggesting that certain skin-derived Acinetobacter species may possess anti-inflammatory properties. Kocuria species are catalase-positive and coagulase-negative Gram-positive coccoid bacteria [66]. Although it is relatively rare, it can cause serious infection [66]. We found that genus Kocuria is a risk factor for intervertebral disc degeneration. However, identifying Kocuria is complex, and there is a high possibility of misidentifying Kocuria as a coagulase-negative Staphylococcus based on Gram staining, catalase, and coagulase negativity [67]. In future research, attention should be paid to distinguishing between the two and being alert to their potential harm to intervertebral disc degeneration. Rothia spp. have been shown to interact with other microbiota through nitrate production, cross-feeding, and release of antimicrobial secondary metabolites to reduce the relative abundance of pathogens and maintain human health [68]. This may explain why it has a protective effect on LBP. Streptococcus alactolyticus and Streptococcus anginosus were found to only exist in degenerated intervertebral discs and not in healthy intervertebral discs [6], and their colonization in intervertebral discs may be harmful. The genus Bacteroides, characterized as Gram-negative, non-spore-forming, and strictly anaerobic bacteria [69], poses a potential risk for colonizing the IVD, thus explaining its association as a risk factor for IVDD. Finegoldia, an anaerobic bacterium considered an opportunistic pathogen [70] associated with physiological joint infections [71], surprisingly demonstrates protective potential against IVDD and LBP in our findings. From the MVMR results, we can conclude that the influence of skin microbiota on LBP and sciatica is easily influenced by BMI and smoking. Genius Bacteroides and Kocuria are not affected by the confounding effects of BMI and smoking, and P.acnes may be confounded by BMI.
This study also has limitations. Firstly, due to data limitations, we can only establish the association between skin microbiota, rather than the microbiota within IVD, and IVDD, LBP, and sciatica through MR. Therefore, we cannot completely deny the role of P.acnes in intervertebral disc degeneration, and previous studies have found that P.acnes in IVD may have other sources. Secondly, strict genome-wide significance thresholds (P < 5 × 10− 8) for some exposed data limit the number of instrumental variables suitable for analysis. Therefore, the threshold must be relaxed to P < 5 × 10− 6. However, this may reduce the accuracy of these results. Furthermore, our results only apply to the European population, and their applicability to other races is uncertain. In addition, the exposed sample only contains 597 cases, and the small sample size may affect the accuracy of the results. Furthermore, although we controlled for potential confounding factors for each IV, the results of MVMR suggest that BMI and smoking status may still contribute to confounding our findings. Therefore, future research must carefully consider and control these potential confounding factors to obtain more accurate and reliable results.
Conclusion
This study identified a potential causal relationship between skin microbiota and IVDD, LBP and sciatica. No evidence suggests skin-derived P.acnes is a risk factor for IVDD, LBP and sciatica. At the same time, IVDD can potentially cause an increase in P.acnes abundance, which supports the contamination theory.
Data availability
The GWAS data used in this study are publicly available and can be obtained from GWAS Catalog (https://www.ebi.ac.uk/gwas/); IEU Open GWAS (https://gwas.mrcieu.ac.uk/datasets/); The summary data analyzed in this study can be obtained from articles and supplementary materials.
Abbreviations
- GWAS:
-
Genome-wide association studies
- MR:
-
Mendelian randomization
- IVs:
-
Instrumental variables
- IVW:
-
Inverse-variance weighted
- SNP:
-
Single nucleotide polymorphism
- MVMR:
-
Multivariate Mendelian randomization
- MR-PRESSO:
-
MR pleiotropy residual sum and outlier
- OR:
-
Odds ratio
- CI:
-
Confidence interval
- ASVs:
-
Amplicon sequence variants
- IVDD:
-
Intervertebral disc degeneration
- LBP:
-
Low back pain
- CoNS:
-
Coagulase-negative Staphylococcus
References
Dagenais S, Caro J, Haldeman S. A systematic review of low back pain cost of illness studies in the United States and internationally. Spine J. 2008;8:8–20. https://doi.org/10.1016/j.spinee.2007.10.005.
Younes M, Béjia I, Aguir Z, Letaief M, Hassen-Zrour S, Touzi M, Bergaoui N. Prevalence and risk factors of disk-related sciatica in an urban population in Tunisia. Joint Bone Spine. 2006;73(5):538–42. https://doi.org/10.1016/j.jbspin.2005.10.022.
Andersson GB. Epidemiological features of chronic low-back pain. Lancet. 1999;354(9178):581–5. https://doi.org/10.1016/S0140-6736(99)01312-4.
Raj PP. Intervertebral disc: anatomy-physiology-pathophysiology-treatment. Pain Pract. 2008;8(1):18–44. https://doi.org/10.1111/j.1533-2500.2007.00171.x.
Jensen RK, Kongsted A, Kjaer P, Koes B. Diagnosis and treatment of sciatica. BMJ. 2019;367:l6273. https://doi.org/10.1136/bmj.l6273.
Rajasekaran S, Soundararajan DCR, Tangavel C, Muthurajan R, Sri Vijay Anand KS, Matchado MS, Nayagam SM, Shetty AP, Kanna RM, Dharmalingam K. Human intervertebral discs harbour a unique microbiome and dysbiosis determines health and disease. Eur Spine J. 2020;29(7):1621–40. https://doi.org/10.1007/s00586-020-06446-z.
Rajasekaran S, Tangavel C, Vasudevan G, Easwaran M, Muthurajan R, K SS, Murugan C, Nayagam SM, Kanna RM, Shetty AP. Bacteria in human lumbar discs - subclinical infection or contamination? Metabolomic evidence for colonization, multiplication, and cell-cell cross-talk of bacteria. Spine J. 2023;23(1):163–77. https://doi.org/10.1016/j.spinee.2022.05.001.
Urquhart DM, Zheng Y, Cheng AC, Rosenfeld JV, Chan P, Liew S, Hussain SM, Cicuttini FM. Could low grade bacterial infection contribute to low back pain? A systematic review. BMC Med. 2015;13:13. https://doi.org/10.1186/s12916-015-0267-x.
Chen Z, Cao P, Zhou Z, Yuan Y, Jiao Y, Zheng Y. Overview: the role of Propionibacterium acnes in nonpyogenic intervertebral discs. Int Orthop. 2016;40(6):1291–8. https://doi.org/10.1007/s00264-016-3115-5.
Perry A, Lambert P. Propionibacterium acnes: infection beyond the skin. Expert Rev Anti Infect Ther. 2011;9(12):1149–56. https://doi.org/10.1586/eri.11.137.
Stirling A, Worthington T, Rafiq M, Lambert PA, Elliott TS. Association between Sciatica and Propionibacterium acnes. Lancet. 2001;357(9273):2024–5. https://doi.org/10.1016/S0140-6736(00)05109-6.
Zhang M, Jia J, Deng L, Cao Z, Hu X, Lei S, Zhang G, Zhu D, Duan Y, Kang X. Risk factors associated with low-grade virulent infection in intervertebral disc degeneration: a systematic review and meta-analysis. Spine J. 2024;24(6):1034–45. https://doi.org/10.1016/j.spinee.2024.02.001.
Rao PJ, Maharaj M, Chau C, Taylor P, Phan K, Choy WJ, Scherman D, Mews P, Scholsem M, Coughlan M, Mobbs R. Degenerate-disc infection study with contaminant control (DISC): a multicenter prospective case-control trial. Spine J. 2020;20(10):1544–53. https://doi.org/10.1016/j.spinee.2020.03.013.
Carricajo A, Nuti C, Aubert E, Hatem O, Fonsale N, Mallaval FO, Vautrin AC, Brunon J, Aubert G. Propionibacterium acnes contamination in lumbar disc surgery. J Hosp Infect. 2007;66(3):275–7. https://doi.org/10.1016/j.jhin.2007.04.007.
Chen X, Lu S, Chopra N, Cui P, Zhang S, Narulla R, Diwan AD. The association between low virulence organisms in different locations and intervertebral disc structural failure: a meta-analysis and systematic review. JOR Spine. 2023;6(2):e1244. https://doi.org/10.1002/jsp2.1244.
Dokoshi T, Chen Y, Cavagnero KJ, Rahman G, Hakim D, Brinton S, Schwarz H, Brown EA, O’Neill A, Nakamura Y, et al. Dermal injury drives a skin to gut axis that disrupts the intestinal microbiome and intestinal immune homeostasis in mice. Nat Commun. 2024;15(1):3009. https://doi.org/10.1038/s41467-024-47072-3.
Coscia MF, Denys GA, Wack MF. Propionibacterium acnes, coagulase-negative Staphylococcus, and the Biofilm-like intervertebral disc. Spine (Phila Pa 1976). 2016;41(24):1860–5. https://doi.org/10.1097/BRS.0000000000001909.
Sekula P, Del Greco MF, Pattaro C, Köttgen A. Mendelian randomization as an Approach to assess causality using Observational Data. J Am Soc Nephrol. 2016;27(11):3253–65. https://doi.org/10.1681/asn.2016010098.
Davey Smith G, Hemani G. Mendelian randomization: genetic anchors for causal inference in epidemiological studies. Hum Mol Genet. 2014;23(R1):R89–98. https://doi.org/10.1093/hmg/ddu328.
Jansen H, Samani NJ, Schunkert H. Mendelian randomization studies in coronary artery disease. Eur Heart J. 2014;35(29):1917–24. https://doi.org/10.1093/eurheartj/ehu208.
Skrivankova VW, Richmond RC, Woolf BAR, Yarmolinsky J, Davies NM, Swanson SA, VanderWeele TJ, Higgins JPT, Timpson NJ, Dimou N, et al. Strengthening the reporting of Observational studies in Epidemiology using mendelian randomization: the STROBE-MR Statement. JAMA. 2021;326(16):1614–21. https://doi.org/10.1001/jama.2021.18236.
Moitinho-Silva L, Degenhardt F, Rodriguez E, Emmert H, Juzenas S, Möbus L, Uellendahl-Werth F, Sander N, Baurecht H, Tittmann L, et al. Host genetic factors related to innate immunity, environmental sensing and cellular functions are associated with human skin microbiota. Nat Commun. 2022;13(1):6204. https://doi.org/10.1038/s41467-022-33906-5.
Yu M, Li Y, Li B, Ge Q. Inflammatory biomarkers and delirium: a mendelian randomization study. Front Aging Neurosci. 2023;15:1221272. https://doi.org/10.3389/fnagi.2023.1221272.
Wang Y, Li T, Fu L, Yang S, Hu YQ. A Novel Method for mendelian randomization analyses with pleiotropy and linkage disequilibrium in genetic variants from Individual Data. Front Genet. 2021;12:634394. https://doi.org/10.3389/fgene.2021.634394.
Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, Laurin C, Burgess S, Bowden J, Langdon R, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408. https://doi.org/10.7554/eLife.34408.
Burgess S, Thompson SG. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol. 2011;40(3):755–64. https://doi.org/10.1093/ije/dyr036.
Fu H, Jiang S, Song S, Zhang C, Xie Q. Causal associations between chronic viral hepatitis and psychiatric disorders: a mendelian randomization study. Front Psychiatry. 2024;15:1359080. https://doi.org/10.3389/fpsyt.2024.1359080.
Burgess S, Davey Smith G, Davies NM, Dudbridge F, Gill D, Glymour MM, Hartwig FP, Kutalik Z, Holmes MV, Minelli C, et al. Guidelines for performing mendelian randomization investigations: update for summer 2023. Wellcome Open Res. 2019;4:186. https://doi.org/10.12688/wellcomeopenres.15555.3.
Bowden J, Del Greco MF, Minelli C, Zhao Q, Lawlor DA, Sheehan NA, Thompson J, Davey Smith G. Improving the accuracy of two-sample summary-data mendelian randomization: moving beyond the NOME assumption. Int J Epidemiol. 2019;48(3):728–42. https://doi.org/10.1093/ije/dyy258.
Burgess S, Thompson SG. Multivariable mendelian randomization: the use of pleiotropic genetic variants to estimate causal effects. Am J Epidemiol. 2015;181(4):251–60. https://doi.org/10.1093/aje/kwu283.
Bowden J, Davey Smith G, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44(2):512–25. https://doi.org/10.1093/ije/dyv080.
Verbanck M, Chen CY, Neale B, Do R. Detection of widespread horizontal pleiotropy in causal relationships inferred from mendelian randomization between complex traits and diseases. Nat Genet. 2018;50(5):693–8. https://doi.org/10.1038/s41588-018-0099-7.
Zheng J, Baird D, Borges MC, Bowden J, Hemani G, Haycock P, Evans DM, Smith GD. Recent developments in mendelian randomization studies. Curr Epidemiol Rep. 2017;4(4):330–45. https://doi.org/10.1007/s40471-017-0128-6.
Dréno B, Pécastaings S, Corvec S, Veraldi S, Khammari A, Roques C. Cutibacterium acnes (Propionibacterium acnes) and acne vulgaris: a brief look at the latest updates. J Eur Acad Dermatol Venereol. 2018;32(Suppl 2):5–14. https://doi.org/10.1111/jdv.15043.
Senker W, Aspalter S, Radl C, Pichler J, Doppler S, Weis S, Webersinke C, Wagner H, Hermann P, Aichholzer M, et al. Frequency and characteristics of bacterial and viral low-grade infections of the intervertebral discs: a prospective, observational study. J Orthop Traumatol. 2022;23(1):15. https://doi.org/10.1186/s10195-022-00633-y.
Albert HB, Lambert P, Rollason J, Sorensen JS, Worthington T, Pedersen MB, Nørgaard HS, Vernallis A, Busch F, Manniche C, Elliott T. Does nuclear tissue infected with bacteria following disc herniations lead to Modic changes in the adjacent vertebrae? Eur Spine J. 2013;22(4):690–6. https://doi.org/10.1007/s00586-013-2674-z.
Albert HB, Sorensen JS, Christensen BS, Manniche C. Antibiotic treatment in patients with chronic low back pain and vertebral bone edema (Modic type 1 changes): a double-blind randomized clinical controlled trial of efficacy. Eur Spine J. 2013;22(4):697–707. https://doi.org/10.1007/s00586-013-2675-y.
Dudli S, Liebenberg E, Magnitsky S, Miller S, Demir-Deviren S, Lotz JC. Propionibacterium acnes infected intervertebral discs cause vertebral bone marrow lesions consistent with Modic changes. J Orthop Res. 2016;34(8):1447–55. https://doi.org/10.1002/jor.23265.
Chen Z, Zheng Y, Yuan Y, Jiao Y, Xiao J, Zhou Z, Cao P. Modic changes and disc degeneration caused by Inoculation of Propionibacterium acnes inside intervertebral discs of rabbits: a pilot study. Biomed Res Int. 2016;2016:9612437. https://doi.org/10.1155/2016/9612437.
Albert HB, Kjaer P, Jensen TS, Sorensen JS, Bendix T, Manniche C. Modic changes, possible causes and relation to low back pain. Med Hypotheses. 2008;70(2):361–8. https://doi.org/10.1016/j.mehy.2007.05.014.
Fischer N, Mak TN, Shinohara DB, Sfanos KS, Meyer TF, Brüggemann H. Deciphering the intracellular fate of Propionibacterium acnes in macrophages. Biomed Res Int. 2013;2013:603046. https://doi.org/10.1155/2013/603046.
Feng P, Che Y, Gao C, Zhu L, Gao J, Vo NV. Immune exposure: how macrophages interact with the nucleus pulposus. Front Immunol. 2023;14:1155746. https://doi.org/10.3389/fimmu.2023.1155746.
Hernández-Palazón J, Puertas-García JP, Martínez-Lage JF, Tortosa JA. Lumbar spondylodiscitis caused by Propionibacterium acnes after epidural obstetric analgesia. Anesth Analg. 2003;96(5):1486–8. https://doi.org/10.1213/01.Ane.0000055819.35383.D3.
Bhanji S, Williams B, Sheller B, Elwood T, Mancl L. Transient bacteremia induced by toothbrushing a comparison of the Sonicare toothbrush with a conventional toothbrush. Pediatr Dent. 2002;24(4):295–9.
McDowell A, Barnard E, Liu J, Li H, Patrick S. Emendation of Propionibacterium acnes subsp. acnes (Deiko et al. 2015) and proposal of Propionibacterium acnes type II as Propionibacterium acnes subsp. defendens subsp. Nov. Int J Syst Evol Microbiol. 2016;66(12):5358–65. https://doi.org/10.1099/ijsem.0.001521.
Rollason J, McDowell A, Albert HB, Barnard E, Worthington T, Hilton AC, Vernallis A, Patrick S, Elliott T, Lambert P. Genotypic and antimicrobial characterisation of Propionibacterium acnes isolates from surgically excised lumbar disc herniations. Biomed Res Int. 2013;2013:530382. https://doi.org/10.1155/2013/530382.
O’Neill AM, Gallo RL. Host-microbiome interactions and recent progress into understanding the biology of acne vulgaris. Microbiome. 2018;6(1):177. https://doi.org/10.1186/s40168-018-0558-5.
Fitz-Gibbon S, Tomida S, Chiu BH, Nguyen L, Du C, Liu M, Elashoff D, Erfe MC, Loncaric A, Kim J, et al. Propionibacterium acnes strain populations in the human skin microbiome associated with acne. J Invest Dermatol. 2013;133(9):2152–60. https://doi.org/10.1038/jid.2013.21.
Cavallo I, Sivori F, Truglio M, De Maio F, Lucantoni F, Cardinali G, Pontone M, Bernardi T, Sanguinetti M, Capitanio B, et al. Skin dysbiosis and Cutibacterium acnes biofilm in inflammatory acne lesions of adolescents. Sci Rep. 2022;12(1):21104. https://doi.org/10.1038/s41598-022-25436-3.
Victoria MG, Gimeno-Gascón A, Ventero MP, Climent-Barberá JM, Cholbi-Llobell F, Rodríguez JC, Baño-Ruiz E, Caminero-Canas MA, Reus-Bañuls S. New contributions on chronic low back pain: disc infection or contaminated cultures? Apmis. 2023;131(6):277–83. https://doi.org/10.1111/apm.13307.
Ben-Galim P, Rand N, Giladi M, Schwartz D, Ashkenazi E, Millgram M, Dekel S, Floman Y. Association between Sciatica and microbial infection: true infection or culture contamination? Spine (Phila Pa 1976). 2006;31(21):2507–9. https://doi.org/10.1097/01.brs.0000238657.13263.b2.
Astur N, Martins DE, Kanas M, Doi AM, Martino MDV, Filho ENK, Wajchenberg M, Lenza M. Bacterial identification in herniated intervertebral discs: a prospective cohort study. Spine J. 2024;S1529–9430(24):00270–5. https://doi.org/10.1016/j.spinee.2024.05.015.
Rigal J, Thelen T, Byrne F, Cogniet A, Boissière L, Aunoble S, Le Huec JC. Prospective study using anterior approach did not show association between Modic 1 changes and low grade infection in lumbar spine. Eur Spine J. 2016;25(4):1000–5. https://doi.org/10.1007/s00586-016-4396-5.
Argemi X, Hansmann Y, Prola K, Prévost G. Coagulase-Negative Staphylococci Pathogenomics. Int J Mol Sci. 2019;20(5):1215. https://doi.org/10.3390/ijms20051215.
Alonzo F 3. Toward uncovering the complexities of bacterial Interspecies Communication and Competition on the skin. mBio. 2022;13(4):e0132022. https://doi.org/10.1128/mbio.01320-22.
Arndt J, Charles YP, Koebel C, Bogorin I, Steib JP. Bacteriology of degenerated lumbar intervertebral disks. J Spinal Disord Tech. 2012;25(7):E211–6. https://doi.org/10.1097/BSD.0b013e318269851a.
Flowers L, Grice EA. The skin microbiota: balancing risk and reward. Cell Host Microbe. 2020;28(2):190–200. https://doi.org/10.1016/j.chom.2020.06.017.
Li J, Chen Y, Wu H, Shan Z, Bei D, Gan K, Liu J, Zhang X, Chen B, Chen J, Zhao FD. Different responses of cervical intervertebral disc caused by low and high virulence bacterial infection: a comparative study in rats. Bioengineered. 2022;13(5):12446–61. https://doi.org/10.1080/21655979.2022.2075305.
Severn MM, Williams MR, Shahbandi A, Bunch ZL, Lyon LM, Nguyen A, Zaramela LS, Todd DA, Zengler K, Cech NB, et al. The ubiquitous human skin Commensal Staphylococcus hominis protects against opportunistic pathogens. mBio. 2022;13(3):e0093022. https://doi.org/10.1128/mbio.00930-22.
Ohshima H, Kurosumi M, Kanto H. New solution of beauty problem by Staphylococcus hominis: relevance between skin microbiome and skin condition in healthy subject. Skin Res Technol. 2021;27(5):692–700. https://doi.org/10.1111/srt.13001.
Tevell S, Hellmark B, Nilsdotter-Augustinsson Å, Söderquist B. Staphylococcus capitis isolated from prosthetic joint infections. Eur J Clin Microbiol Infect Dis. 2017;36(1):115–22. https://doi.org/10.1007/s10096-016-2777-7.
Hu J, Lu W, Li X, Yang J, Tan M, Hu K, Wang Q, Deng S, Liu Y, Chen J, et al. Microbiota differences of skin and pharyngeal microbiota between patients with plaque and guttate psoriasis in China. Front Microbiol. 2022;13:937666. https://doi.org/10.3389/fmicb.2022.937666.
Zhang XN, Wu CY, Wu ZW, Xu LX, Jiang FT, Chen HW. Association between the Diabetic Foot Ulcer and the bacterial colony of the skin based on 16S rRNA gene sequencing: an observational study. Clin Cosmet Investig Dermatol. 2023;16:2801–12. https://doi.org/10.2147/ccid.S425922.
Montaña S, Schramm ST, Traglia GM, Chiem K, Di Parmeciano G, Almuzara M, Barberis C, Vay C, Quiroga C, Tolmasky ME, et al. The genetic analysis of an Acinetobacter johnsonii clinical strain evidenced the Presence of Horizontal genetic transfer. PLoS ONE. 2016;11(8):e0161528. https://doi.org/10.1371/journal.pone.0161528.
Qi L, Peng J, Huang X, Zhou T, Tan G, Li F. Longitudinal dynamics of gut microbiota in the pathogenesis of acute graft-versus-host disease. Cancer Med. 2023;12(24):21567–78. https://doi.org/10.1002/cam4.6557.
Ziogou A, Giannakodimos I, Giannakodimos A, Baliou S, Ioannou P. Kocuria Species Infections in Humans-A narrative review. Microorganisms. 2023;11(9):2362. https://doi.org/10.3390/microorganisms11092362.
Kandi V, Palange P, Vaish R, Bhatti AB, Kale V, Kandi MR, Bhoomagiri MR. Emerging bacterial infection: identification and clinical significance of Kocuria Species. Cureus. 2016;8(8):e731. https://doi.org/10.7759/cureus.731.
West SR, Suddaby AB, Lewin GR, Ibberson CB, Rothia. Trends Microbiol. 2024;32(7):720–1. https://doi.org/10.1016/j.tim.2024.03.009.
Shin JH, Tillotson G, MacKenzie TN, Warren CA, Wexler HM, Goldstein EJC. Bacteroides and related species: the keystone taxa of the human gut microbiota. Anaerobe. 2024;85:102819. https://doi.org/10.1016/j.anaerobe.2024.102819.
Rosenthal ME, Rojtman AD, Frank E. Finegoldia magna (formerly Peptostreptococcus magnus): an overlooked etiology for toxic shock syndrome? Med Hypotheses. 2012;79(2):138–40. https://doi.org/10.1016/j.mehy.2012.04.013.
Levy PY, Fenollar F, Stein A, Borrione F, Raoult D. Finegoldia magna: a forgotten pathogen in prosthetic joint infection rediscovered by molecular biology. Clin Infect Dis. 2009;49(8):1244–7. https://doi.org/10.1086/605672.
Acknowledgements
Not applicable.
Funding
Tianjin Municipal Health Commission (Grant TJWJ2023XK023 to Baoshan Xu).
Author information
Authors and Affiliations
Contributions
YJ: designed the study; YJ, HC, SH and ZH: writing - original draft preparation; YJ and HC: formal analysis and investigation; All authors: writing - review and editing; BX: supervision. All authors agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
Corresponding author
Ethics declarations
Ethical approval and consent to participate
We used publicly available GWAS summary data in this study. therefore, this analysis does not require ethical approval and consent to participate.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Jia, Y., Chen, H., Huang, S. et al. Causal effects of skin microbiota on intervertebral disk degeneration, low back pain and sciatica: a two-sample Mendelian randomization study. J Orthop Surg Res 19, 476 (2024). https://doi.org/10.1186/s13018-024-04980-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13018-024-04980-w