Skip to main content

SOX9 gene shows association with adolescent idiopathic scoliosis predisposition in Northwest Indians

Abstract

Background

Adolescent idiopathic scoliosis (AIS) is a common structural deformity of the spine affecting adolescent individuals globally. The disorder is polygenic and is accompanied by the association of various genetic loci. Genetic studies in Chinese and Japanese populations have shown the association of genetic variants of SOX9 with AIS curve severity. However, no genetic study evaluating the association of SRY-Box Transcription Factor 9 (SOX9) variants with AIS predisposition has been conducted in any Indian population. Thus, we aimed to investigate the association of the genetic variants of the SOX9 along with 0.88 Mb upstream region with AIS susceptibility in the population of Northwest India.

Methods

In total, 113 AIS cases and 500 non-AIS controls were recruited from the population of Northwest India in the study and screened for 155 genetic variants across the SOX9 gene and 0.88 Mb upstream region of the gene using Global Screening Array-24 v3.0 chip (Illumina). The statistical significance of the Bonferroni threshold was set at 0.000322.

Result

The results showed the association of 11 newly identified variants; rs9302936, rs7210997, rs77736349, rs12940821, rs9302937, rs77447012, rs8071904, rs74898711, rs9900249, rs2430514, and rs1042667 with the AIS susceptibility in the studied population. Only one variant, rs2430514, was inversely associated with AIS in the population, while the ten variants were associated with the AIS risk. Moreover, 47 variants clustered in the gene desert region of the SOX9 gene were associated at a p-value ≤ 0.05.

Conclusion

The present study is the first to demonstrate the association of SOX9 enhancer locus variants with AIS in any South Asian Indian population. The results are interesting as rs1042667, a 3' untranslated region (UTR) variant in the exon 3 and upstream variants of the SOX9 gene, were associated with AIS susceptibility in the Northwest Indian population. This provides evidence that the variants in the enhancer region of SOX9 might regulate its gene expression, thus leading to AIS pathology and might act as an important gene for AIS susceptibility.

Background

Numerous genetic loci have been identified to be associated with the adolescent idiopathic scoliosis (AIS) predisposition [1,2,3,4,5,6,7,8]. Several genetic studies have highlighted the importance of SOX9 in the pathogenesis of AIS [3, 9]. The SOX9 gene is a master transcription factor for chondrogenesis [10]; mutations in the gene were found to be associated with many disorders, including kyphoscoliosis, sex reversal, bending of long bones, and others [11]. The altered expression of the SOX9 gene in the mesenchymal stem cells of AIS patients has highlighted its role in the AIS pathogenesis [9]. It has been found that SOX9 is required for spine homeostasis in mice and drives the chondrogenesis in the spinal growth plate [12].

Moreover, a genetic study revealed that a genetic variant rs12946942, located near the SOX9 gene, was found to be associated with the severity of the AIS curve in Japanese and Chinese populations [4]. It has been established that the regulatory domain surrounding the SOX9 gene spans more than 2 Mb, consisting of multiple enhancers that alter the expression of the SOX9 gene during development [13], more specifically, a region located between 91 to 974 kb upstream of the SOX9 has shown enhancer marks [14]. Based on the rationale of the importance of SOX9 and its upstream region in the genetic predisposition of skeletal disorders and AIS, the screening of genetic variants across SOX9 and 0.88 Mb upstream region was performed using high-throughput genotyping. The present study primarily aims to evaluate the association of the SOX9 gene with AIS susceptibility in the Northwest Indian population in a case–control study design.

Materials and methods

Study design

The present study is a case–control population-based genetic epidemiology study.

Sample collection

Samples for the present study were collected from the school-based screening programme (the detailed characteristics of the population have been represented in our previous study) [15], followed by strict clinical evaluations after seeking informed consent from the individuals. The general characteristics of the cases and controls are provided in Table 1. Additionally, clinically diagnosed samples were recruited in the study from the All India Institute of Medical Sciences, New Delhi. All the cases were clinically diagnosed based on the Cobb angle measurement [16, 17]. The inclusion criteria for the cases were based on a Cobb angle measurement of 10° or greater. The cases and control samples recruited in the present study lie in the age group of 10–28 years. While cases were selected based on the age of onset of scoliosis, between 10 and 18 years. The control samples recruited in the present study were 21 years or older till 28 years identified during the population screening (the control selected in the present study were intentionally taken over 21 years to rule out any false negative case of late age onset of AIS, that would have remained undetected during adolescent but might have diagnosed few years later). The 2 millilitres of blood samples were collected from both the cases and controls in EDTA vials and stored at − 20 °C till further processing. A total sample size of 613 participants, comprising 113 cases and 500 non-AIS controls of Northwest Indian ethnicity (including various caste groups: Rajputs, Baniya, Brahmins, Kashmiri Pandits, Kashmiri Muslims, Sikhs, Buddhists, and many more) were recruited in the study. Of 113 AIS cases, 53 were males, and 60 were female.

Table 1 General characteristics of AIS cases and controls of the Northwest Indian population

DNA isolation and genotyping

The genomic DNA of each participant was extracted from the peripheral blood sample utilising the XpressDNA Blood Kit (MagGenome®). The DNA was further processed for qualitative and quantitative analysis through the utilisation of Nanodrop 2000/2000c (ThermoFisher, USA) and Agarose gel electrophoresis. All the samples were genotyped using an Infinium Global Screening Array (GSA)-24 v3.0 chip (Illumina).

Statistical analysis

The genotyping data were analysed using PLINK v1.07 and v1.09 [18]. A comprehensive analysis of the SOX9 gene and its 0.88 Mb upstream region revealed a total of 155 variations. The variants were obtained through the genotyping data using the PLINK flag “–extract” and subsequently subjected to filtration based on genotype and minor allele frequency, with thresholds set at “–mind 0.05” and “–maf 0.05”, respectively, as well as a genotyping rate threshold of “–geno 0.1”. Furthermore, the study utilised a predetermined threshold of 0.0005 to assess the Hardy–Weinberg equilibrium (HWE). Ten cases and five controls were excluded from the study due to a low call rate of less than 95% or being identified as ethnic outliers. This determination was made through principal component analysis after merging data from the UK Biobank and 1000 Genomes Phase 3 with our population data. The variants that passed the quality control assessment were subjected to an association test utilising the parameters –assoc –ci 0.95. The statistical significance of the association with AIS was evaluated at a Bonferroni correction of p-value = 0.000322. The regional association plot of filtered variants was generated using the LDassoc tool [19] using South Asians (SAS) as the reference population. The linkage disequilibrium (LD) of all variants was assessed using the Haploview tool [20]. The associated variants were further annotated using the variant effect predictor (VEP) tool from Ensembl [21, 22] and the UCSC genome browser [23, 24].

Results

The Cobb angle in the cases ranges from 18.9° to 101° bending of the spine. The average Cobb angle of the AIS cases recruited in the study was 57.3°, with thoracic curvature being the prevalent type.

In the present study, out of 155 variants, 135 variants passed the quality test after filtering the genotyping data (HWE in controls ≥ 0.0005 and minor allele frequency in controls ≥ 0.05). From 135 variants, 11 variants were found to be significantly associated with the AIS genetic predisposition in the population of Northwest India (Table 2). Moreover, 47 variants were able to cross the p-value threshold of ≤ 0.05 in the present study, indicating that the association of the variants of the SOX9 gene does not appear by chance in the studied population group. All the variants other than the 11 significantly associated variants are tabulated in Additional file 2: Table S1.

Table 2 List of variants significantly associated with the AIS susceptibility in the Northwest Indian population at Bonferroni p-value threshold of 0.000322

The associated variants were further annotated using the Ensembl VEP tool and the UCSC genome browser. It was observed that variant rs1042667 was present at the 3' UTR region in exon 3 of the SOX9 gene; thus, it might play a role in the gene expression of the SOX9 gene. In contrast, all other associated variants were intergenic (Additional file 2: Table S1). The regional association of 135 variants that passed after filtering the genotyping data is represented in Fig. 1. Furthermore, the linkage disequilibrium plot of all 155 variants of SOX9 was generated using the r-squared values. The dark black colour of the boxes represents variants in high LD in the studied population (Additional file 1: Figure S1).

Fig. 1
figure 1

Regional association plot representing 135 variants that passed the quality test in the SOX9 and 0.88 Mb upstream of the gene

Discussion

Studies have extensively investigated the functional role of the SOX9 gene. 517–595 Kb gene desert upstream of the SOX9 gene promoter is related to the regulation of the gene in gonads [25]. The role of variants in and near the SOX9 gene in the pathogenesis of various skeletal disorders has been noted. The reciprocal translocation in the breakpoint position located 203 kb upstream of the SOX9 gene was found to be associated with acampomelic campomelic dysplasia [26]. Similarly, microdeletion on either side of the SOX9 was associated with Pierre Robin sequence malformation [27, 28]. Genetic variations mapped to 1.3 Mb upstream and downstream of SOX9 have been associated with several skeletal abnormalities and other disorders [29]. Numerous studies have indicated the importance of SOX9 expression in the chondrocyte and bone development [10, 30,31,32,33].

In the present study, screening of genetic variants in the upstream region and across the SOX9 gene was performed in the population of Northwest India. Interestingly, out of 135 variants, 11 newly identified variants (Table 2) were associated with AIS genetic predisposition in the studied population group. All the variants except rs1042667 were observed to be intergenic/intronic in nature. The annotation of the variants using VEP and UCSC showed no significant functional consequences. One of the associated variants rs1042667 [OR 1.776 (1.31–2.408, 95% CI), p-value 0.000195] that passed the Bonferroni correction threshold is from the 3' UTR region of SOX9 gene, and we believe the variant might have a regulatory role on the SOX9 gene expression and thus the association.

Interestingly, the earlier reported intergenic variant, rs12946942 in the 17q24.3 locus, upstream of the SOX9 gene, did not show an association in the present study. This variant has shown an association with the AIS curve severity in the Japanese [4] and Han Chinese [4] populations and was replicated in an international meta-analysis with AIS severity [34]. However, a study in Chinese ancestry has also shown a lack of association of the variant rs12946942 with AIS curve severity [35].

We hypothesised that the differences are potentially due to genetic heterogeneity and differences in LD structure in the region. We further evaluated the LD of rs12946942 with the associated variants that were found in the present study. As anticipated, the variants did not show strong LD with the previously reported variant rs12946942 (~ D´ = 0.5, and ~ R2 = 0.02) in the SAS population. This substantiates our hypothesis of the possibility of the LD structure differences and evidence of genetic heterogeneity in different populations for AIS predisposition.

Limitations of the study

This is the first genetic replication study of AIS susceptibility in any South Asian Indian (SAI) population group. The study still includes a small size of the population. It is pertinent to focus on the large sample size-based replication studies in the Indian population to replicate the association results.

Strength of the study

No dataset from the region is available for comparisons of the association of the SOX9 genetic variants. The present study is baseline data for SOX9 genetic variants association with AIS in the Indian population. Also, even though the incidence of AIS is low, the present study will be a reference for future studies of AIS in SAI populations.

Conclusion

The study hints towards the genetic heterogeneity in the Asian population. More replication studies are warranted in other population groups to identify the role of this region in AIS susceptibility to the respective population group and understand genetic heterogeneity in a better way. Also, functional studies are needed to validate the region, as an enhancer of the SOX9 gene, related to scoliosis and associated phenotypes.

Availability of data and materials

All the data related to the study have been provided in the manuscript.

Abbreviations

AIS:

Adolescent idiopathic scoliosis

SOX9:

SRY-Box Transcription Factor 9

UTR:

Untranslated region

HWE:

Hardy–Weinberg equilibrium

SAS:

South Asians

LD:

Linkage disequilibrium

VEP:

Variant effect predictor

References

  1. Khanshour AM, Kou I, Fan Y, Einarsdottir E, Makki N, Kidane YH, et al. Genome-wide meta-analysis and replication studies in multiple ethnicities identify novel adolescent idiopathic scoliosis susceptibility loci. Hum Mol Genet. 2018;27(22):3986–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Kou I, Takahashi Y, Johnson TA, Takahashi A, Guo L, Dai J, et al. Genetic variants in GPR126 are associated with adolescent idiopathic scoliosis. Nat Genet. 2013;45(6):676–9.

    Article  CAS  PubMed  Google Scholar 

  3. Londono D, Kou I, Johnson TA, Sharma S, Ogura Y, Tsunoda T, et al. A meta-analysis identifies adolescent idiopathic scoliosis association with LBX1 locus in multiple ethnic groups. J Med Genet. 2014;51(6):401–6.

    Article  CAS  PubMed  Google Scholar 

  4. Miyake A, Kou I, Takahashi Y, Johnson TA, Ogura Y, Dai J, et al. Identification of a susceptibility locus for severe adolescent idiopathic scoliosis on chromosome 17q24.3. PLoS ONE. 2013;8(9):e72802.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Ogura Y, Kou I, Miura S, Takahashi A, Xu L, Takeda K, et al. A functional SNP in BNC2 is associated with adolescent idiopathic scoliosis. Am J Hum Genet. 2015;97(2):337–42.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Sharma S, Londono D, Eckalbar WL, Gao X, Zhang D, Mauldin K, et al. A PAX1 enhancer locus is associated with susceptibility to idiopathic scoliosis in females. Nat Commun. 2015;6(1):6452.

    Article  CAS  PubMed  Google Scholar 

  7. Takahashi Y, Kou I, Takahashi A, Johnson TA, Kono K, Kawakami N, et al. A genome-wide association study identifies common variants near LBX1 associated with adolescent idiopathic scoliosis. Nat Genet. 2011;43(12):1237–40.

    Article  CAS  PubMed  Google Scholar 

  8. Zhu Z, Tang NL-S, Xu L, Qin X, Mao S, Song Y, et al. Genome-wide association study identifies new susceptibility loci for adolescent idiopathic scoliosis in Chinese girls. Nat Commun. 2015;6(1):8355.

    Article  CAS  PubMed  Google Scholar 

  9. Chen C, Xu C, Zhou T, Gao B, Zhou H, Chen C, et al. Abnormal osteogenic and chondrogenic differentiation of human mesenchymal stem cells from patients with adolescent idiopathic scoliosis in response to melatonin. Mol Med Rep. 2016;14(2):1201–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Dy P, Wang W, Bhattaram P, Wang Q, Wang L, Ballock RT, et al. Sox9 directs hypertrophic maturation and blocks osteoblast differentiation of growth plate chondrocytes. Dev Cell. 2012;22(3):597–609.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Lekovic GP, Rekate HL, Dickman CA, Pearson M. Congenital cervical instability in a patient with camptomelic dysplasia. Childs Nerv Syst. 2006;22:1212–4.

    Article  PubMed  Google Scholar 

  12. Wise CA, Sepich D, Ushiki A, Khanshour AM, Kidane YH, Makki N, et al. The cartilage matrisome in adolescent idiopathic scoliosis. Bone Research. 2020;8(1):13.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Lefebvre V, Angelozzi M, Haseeb A. SOX9 in cartilage development and disease. Curr Opin Cell Biol. 2019;61:39–47.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Cha S, Lim JE, Park AY, Do J-H, Lee SW, Shin C, et al. Identification of five novel genetic loci related to facial morphology by genome-wide association studies. BMC Genomics. 2018;19(1):481.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Singh H, Shipra SV, Sharma I, Sharma A, Modeel S, et al. The first study of epidemiology of adolescent idiopathic scoliosis shows lower prevalence in females of Jammu and Kashmir. India Am J Transl Res. 2022;14(2):1100–6.

    PubMed  Google Scholar 

  16. Wang J, Zhang J, Xu R, Chen TG, Zhou KS, Zhang HH. Measurement of scoliosis Cobb angle by end vertebra tilt angle method. J Orthop Surg Res. 2018;13(1):223.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Hurtado-Avilés J, Santonja-Medina F, León-Muñoz VJ, Sainz de Baranda P, Collazo-Diéguez M, Cabañero-Castillo M, et al. Validity and Absolute Reliability of the Cobb Angle in Idiopathic Scoliosis with TraumaMeter Software. Int J Environ Res Public Health. 2022;19(8):4655.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Machiela MJ, Chanock SJ. LDassoc: an online tool for interactively exploring genome-wide association study results and prioritizing variants for functional investigation. Bioinformatics. 2018;34(5):887–9.

    Article  CAS  PubMed  Google Scholar 

  20. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5.

    Article  CAS  PubMed  Google Scholar 

  21. McLaren W, Gil L, Hunt SE, Riat HS, Ritchie GR, Thormann A, et al. The ensembl variant effect predictor. Genome Biol. 2016;17(1):1–14.

    Article  Google Scholar 

  22. Cunningham F, Allen JE, Allen J, Alvarez-Jarreta J, Amode MR, Armean IM, et al. Ensembl 2022. Nucleic Acids Res. 2022;50(D1):D988–95.

    Article  CAS  PubMed  Google Scholar 

  23. Karolchik D, Baertsch R, Diekhans M, Furey TS, Hinrichs A, Lu Y, et al. The UCSC genome browser database. Nucleic Acids Res. 2003;31(1):51–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Nassar LR, Barber GP, Benet-Pagès A, Casper J, Clawson H, Diekhans M, et al. The UCSC genome browser database: 2023 update. Nucleic Acids Res. 2023;51(D1):D1188–95.

    Article  CAS  PubMed  Google Scholar 

  25. Benko S, Gordon CT, Mallet D, Sreenivasan R, Thauvin-Robinet C, Brendehaug A, et al. Disruption of a long distance regulatory region upstream of <em>SOX9</em> in isolated disorders of sex development. J Med Genet. 2011;48(12):825–30.

    Article  CAS  PubMed  Google Scholar 

  26. Takano T, Ota H, Ohishi H, Hata K, Furukawa R, Nakabayashi K. Adult acampomelic campomelic dysplasia and disorders of sex development due to a reciprocal translocation involving chromosome 17q24.3 upstream of the SOX9 gene. Eur J Med Genet. 2021;64(11):104332.

    Article  CAS  PubMed  Google Scholar 

  27. Amarillo IE, Dipple KM, Quintero-Rivera F. Familial microdeletion of 17q243 upstream of SOX9 is associated with isolated Pierre Robin sequence due to position effect. Am J Med Genet Part A. 2013;161(5):1167–72.

    Article  CAS  Google Scholar 

  28. Benko S, Fantes JA, Amiel J, Kleinjan D-J, Thomas S, Ramsay J, et al. Highly conserved non-coding elements on either side of SOX9 associated with Pierre Robin sequence. Nat Genet. 2009;41(3):359–64.

    Article  CAS  PubMed  Google Scholar 

  29. Smyk M, Akdemir KC, Stankiewicz P. SOX9 chromatin folding domains correlate with its real and putative distant cis-regulatory elements. Nucleus. 2017;8(2):182–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bi W, Deng JM, Zhang Z, Behringer RR, De Crombrugghe B. Sox9 is required for cartilage formation. Nat Genet. 1999;22(1):85–9.

    Article  CAS  PubMed  Google Scholar 

  31. Akiyama H. Control of chondrogenesis by the transcription factor Sox9. Mod Rheumatol. 2008;18(3):213–9.

    Article  CAS  PubMed  Google Scholar 

  32. Foster JW, Dominguez-Steglich MA, Guioli S, Kwok C, Weller PA, Stevanović M, et al. Campomelic dysplasia and autosomal sex reversal caused by mutations in an SRY-related gene. Nature. 1994;372(6506):525–30.

    Article  CAS  PubMed  Google Scholar 

  33. Seda M, Crespo B, Corcelli M, Osborn DP, Jenkins D. A CRISPR/Cas9-generated mutation in the zebrafish orthologue of PPP2R3B causes idiopathic scoliosis. Sci Rep. 2023;13(1):6783.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Takeda K, Kou I, Otomo N, Grauers A, Fan Y-H, Ogura Y, et al. A multiethnic meta-analysis defined the association of rs12946942 with severe adolescent idiopathic scoliosis. J Hum Genet. 2019;64(5):493–8.

    Article  CAS  PubMed  Google Scholar 

  35. Man GC-W, Tang NL-S, Chan TF, Lam TP, Li JW, Ng BK-W, et al. Replication study for the association of GWAS-associated loci with adolescent idiopathic scoliosis susceptibility and curve progression in a Chinese population. Spine. 2019;44(7):464.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

All the authors acknowledge NSS Volunteers and M.Sc. Dissertation students at SMVDU, Mr. Hetender Singh for their assistance in population screening and sample collection. Authors acknowledge Central MassArray Facility, SoBT SMVDU for facility for processing of DNA samples. SS acknowledges specifically grant from National Geographic Society GR-000046670 for developing baseline dataset for population structure and appropriate stratification analyses for population based studies in Jammu and Kashmir.

Funding

No funding was received to conduct the study.

Author information

Authors and Affiliations

Authors

Contributions

SS, ER conceived the study, and VS assisted in study design and execution. HS, S, carried out population screening, sample collection, DNA isolation, and genotyping of the samples. HS performed statistical analyses of the data. HS wrote the manuscript. MG, NG, GG, AKP, RS, SP, BG are clinical collaborators who evaluated the cases and samples. HS, SS, ER and VS critically reviewed the manuscript.

Corresponding authors

Correspondence to Ekta Rai or Swarkar Sharma.

Ethics declarations

Ethics approval and consent to participate

Shri Mata Vaishno Devi University, Institutional Ethical Review Board (IERB) (IERB Serial No: SMVDU/IERB/18/69 and SMVDU/IERB/18/68, on 23rd April 2018) and Institute Ethics committee, All India Institute of Medical Sciences, New Delhi (Ref No.: IEC-20/08.01.2021, RP-25/2021 on 8th January 2021) approved the study.

Consent for publication

Written informed consent was taken from all the subjects participated in the study.

Competing interests

All 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.

Supplementary Information

Additional file 1: Figure S1.

Linkage disequilibrium plot of all the genetic variants screened across the locus, including the SOX9 gene in the Northwest Indian population.

Additional file 2: Table S1.

List of variants studied in the Northwest Indian population in addition to variations that passed the Bonferroni Correction Threshold.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, 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 changes were made. 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/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Singh, H., Shipra, Gupta, M. et al. SOX9 gene shows association with adolescent idiopathic scoliosis predisposition in Northwest Indians. Eur J Med Res 29, 66 (2024). https://0-doi-org.brum.beds.ac.uk/10.1186/s40001-024-01635-8

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s40001-024-01635-8

Keywords