High levels of the tumor necrosis factor alpha (TNF-α) induce apoptosis and pro-inflammatory effects for primary degeneration of tendon and development of tendinopathy. The aim of this study was to investigate the association between the TNF-α polymorphisms and tendinopathy in athletes.
Two hundred and seventy athletes (135 tendinopathy cases and 135 controls) were included and genotyped (TNF-α -1031T > C; -857 C > T; -308G > A) using TaqMan validated assays. The association of the polymorphisms with tendinopathy was evaluated by a multivariate logistic regression model, using odds ratios (OR) and 95 % confidence intervals (CI).
The variant allele − 308 A was significantly associated with patellar (OR: 1.9; 95 % CI: 1.01–3.6) or Achilles tendinopathies (OR: 2.7; 95 % CI: 1.1–6.7). No significant differences were found in allele or genotype distributions of the − 1031T > C and − 857 C > T polymorphisms between cases and controls. TNF-α TCA haplotype was associated with increased tendinopathies risk, either considering all cases (OR: 2.6, 95 % CI: 1.3–5.3), patellar (OR: 3.3, 95 % CI: 1.5–7.3), rotator cuff (OR: 3.1, 95 % CI: 1.4–7.2) or Achilles tendinopathies (OR: 3.8, 95 % CI: 1.1–12.7).
These results suggest that the TNF-α polymorphisms could influence the susceptibility to developing tendinopathy among athletes. Knowledge of the TNF-α polymorphisms associated to tendinopathy in athletes can further understanding of the inflammatory role in the early stages of the disease and contribute for sports injury surveillance programmes, in which athletes with TNF-α TCA haplotype could be early subjected to cryotherapy after training and competition to avoid tendinopathy development.
Tendinopathy is characterized by pain, swelling, structural change and functional limitation of the tendon due to overuse [1, 2]. It is the main reason for clinical musculoskeletal complaint in athletes (15–50 %) , which can lead to reduced level of performance or end of one’s sport career . The commonly identified risk factor associated with tendinopathy in athletes are age, sex, metabolic and hormonal concentrations, and high physical load during training and matches according to each sport modality [4, 5].
Single nucleotide polymorphisms (SNPs), a variation of the nucleotide at a single position in DNA sequence, involved with inflammatory process were associated as non-modifiable factors for developing tendinopathies [6,7,8]. Pro-inflammatory cytokines, such as tumor necrosis factor alpha (TNF-α), and growth factors have been implicated as a mechanism in early stages of tendinopathies [1, 9]. Macrophages, mast cells, fibroblasts and endothelial cells synthesized and released TNF-α cytokine, which is the main chemokine inducer when tendons are mechanically overloaded . TNF-α signaling is mediated by two functionally distinct receptors: TNF-α receptor-1 (TNFR1) and TNF-α receptor-2 (TNFR2). The ligand receptor interaction between TNF-α–TNFR1 is responsible for induce apoptosis and proinflammatory effects, while the interaction between TNF-α–TNFR2 regulates tissue growth and repair . An experimental tendinopathy model produced by overuse shown that TNF-α mRNA was increased 11-fold in torn supraspinatus tendon compared to controls . In addition, TNF-α and its receptors were expressed in peritendinous tissue , and in rounded/enlarged nucleus human tenocytes, a typical characteristic of tendinopathy .
TNF-α is encoded by the gene of the same name, located on chromosome 6p21.3 between the human leukocyte antigen-B (HLA-B) and human leukocyte antigen-DR (HLA-DR) genes at the major histocompatibility complex class III region . The SNPs in the promoter region of TNF-α gene, such as -1031T > C (rs1799964), -857 C > T (rs1799724) and − 308G > A (rs1800629), have shown potential to alters the binding of transcription factors in the DNA, regulating TNF-a expression [15, 16]. The choice of TNF-α SNPs was due to their biological relevance in altering the expression of the gene and for have relatively high frequency in different populations . In addition, previous studies associated these TNF-α SNPs with some diseases, such as inflammatory bowel diseases , Congenital Zika syndrome , and cystic fibrosis .
We hypothesized that the TNF-α SNPs could be associated with a risk of developing tendinopathy in athletes; since high TNF-α expression could be modulated by polymorphisms and TNF-α induce apoptosis and pro-inflammatory effects for primary degeneration of tendon. As far as we know, there are no studies evaluating the influence of TNF-α SNPs as possible risk factors involved in the inflammatory molecular mechanism leading to tendinopathy. Thus, this study aimed to investigate the association between the TNF-α polymorphisms and tendinopathy in athletes.
Study design and population
This case-control study was approved by the Human Ethics Committee of the Instituto Nacional de Traumatologia e Ortopedia Jamil Haddad (protocol number 2.455.630/2017). All participating athletes provided written informed consent and answered a questionnaire about their epidemiological, clinical, sport and training characteristics, as well as tendon injury history and their specific information such as type, location and number of tendinopathy episodes, as previously described . At the end of data collection, a trained observer checked the questionnaire with each athlete, and the database was double-checked by different trained researchers.
The inclusion criteria were Brazilian competitive levels athletes aged 18–45 years old who were recruited between March 2018 and September 2019 at different sports training centres and competitions.
One hundred thirty-five athletes had tendinopathy clinically diagnosed by medical practitioners and confirmed with magnetic resonance image examination (MRI). All tendinopathy diagnoses were confirmed by two blinded specialized orthopaedic surgeons, as described in previous studies [6, 21]. The control group (N = 135) consisted of athletes without previous imaging diagnosis of tendinopathy and who were matched with tendinopathy cases for age (difference of ± 2 years), sex and sport modality. The sample size was calculated using Epi Info 7, version 7.1.3. (http://wwwn.cdc.gov/epiinfo/html/downloads.htm) to detect a difference between case and control groups, assuming an odds ratio (OR) of 2.0 with a power of 0.8 and 5 % type I error. The OR was based on previous evidence [22,23,24], and at least 128 athletes per group was necessary.
Genotyping of polymorphisms
Genomic DNA was obtained from oral mucosa collected from each athlete by swab. The TNF-α -1031T > C (rs1799964), -857 C > T (rs1799724) and − 308G > A (rs1800629) polymorphisms were genotyped using a TaqMan allelic discrimination assay obtained from Applied Biosystems (C___7514871_10, C__11918223_10 and C___7514879_10, respectively). For all polymorphisms real-time polymerase chain reaction (PCR) reactions were performed on a 7500 Real-Time System (Applied Biosystems, Foster City, CA, USA), and the genotypes were then determined directly. To assure genotyping quality, in each reaction two standardized positive controls of each polymorphism genotype were used.
The normally distribution of studied population was determined by the Shapiro-Wilk test. Comparisons of continuous variables between tendinopathy cases and controls groups were performed using the Student’s t test, and data were presented as mean ± standard deviation (SD). According to distribution and clinical significance the continuous variables (height, age at the beginning of sport practice, years of training and weekly training hours) were divided into quartiles. Categorical data were shown in proportions and differences between the two groups were evaluated using the Chi-squared (χ2) statistic test or Fischer exact test, when applicable.
Deviations from Hardy–Weinberg equilibrium (HWE) were assessed by the goodness-of-fit χ2 test. TNF-α (-1031T > C, -857 C > T, -308G > A) allele frequency and genotype distribution were derived by gene counting and frequencies between the two groups were compared using the χ2 test or, when appropriate, the Fisher’s exact test. The haplotype patterns and linkage disequilibrium coefficients (D’ is degree of imbalance in module and R2 is degree of correlation) were inferred using Haploview, as previously described .
Multivariate logistic regression analyses model were performed to evaluate the possible associations between epidemiological, clinic, sport and training characteristics as much as of the polymorphisms with tendinopathy, which was estimated by the OR with a 95 % confidence interval (95 % CI). As a final regression model used to control possible confounding factors, each variable was introduced considering the biological and statistical significance of the univariate analysis, which a input significance level less than 0.25 (P ≤ 0.25) and output significance was 0.05 (P ≤ 0.05) at the regression model. The difference was statistically significant when P < 0.05. All analyses were performed using the Statistical Package for Social Sciences (SPSS Inc., Chicago, IL, USA, version 20.0).
Of the 135 tendinopathy cases, 24 athletes (17.8 %) reported more than one diseased tendon. The cases reported tendinopathies of the patellar (N = 62, 45.9 %), rotator cuff (N = 50, 37.0 %), Achilles (N = 22, 16.3 %), wrist (N = 15, 11.1 %) and elbow tendinopathy (N = 9, 6.7 %) (Fig. 1).
Age, sex, sport modality of the controls and all cases, as well as cases divided into affected tendon groups, is summarized in Fig. 2. There was no significant difference of the age, sex and sports modality between tendinopathies subgroups (patellar, rotator cuff and Achilles) and control; however, these variables entered the multivariate model for stratified association analyzes, according to the biological importance for the tendinopathy development.
The demographic, clinical, sport and training characteristics variables of all tendinopathy cases and controls were presented in Table 1. In summary, all variables were analyzed to identify possible confounding variables of the true association between SNPs and tendinopathy. Initially, the variables BMI (P = 0.09), alcohol consumption (P = 0.25), nutritional follow-up by a specialist during a sports career (P = 0.002), declared preference member (P = 0.21) and weekly training hours (P = 0.15) were inserted in the logistic regression model. After multivariate analysis, only BMI and nutritional follow-up remained in this model.
The distribution of TNF-α (-1031T > C, -857 C > T and − 308G > A) SNPs was in Hardy–Weinberg equilibrium. The minor allele frequencies of the TNF-α SNPs in the study population are shown in Fig. 3. After adjustment by co-factors of the logistic regression model (age, sex, sport modality, BMI and nutritional follow-up) the TNF-α-308 A allele was significantly associated with patellar and Achilles tendinopathies. Moreover, the TNF-α -308AA genotype was only present in the tendinopathy cases, either considering all cases, patellar, rotator cuff or Achilles tendinopathies. Considering the recessive co-dominance model (TNF-α -308GG + GA versus AA) the TNF-α -308AA genotype was significantly associated with tendinopathy cases, either considering all cases, patellar and Achilles tendinopathies (Table 2). Despite of the TNF-α -308AA genotype suggests a more likely of having tendinopathy nis more than one tendon, when to compared only tendinopathy cases group (one versus 2 or more affected tendons) there was not statistical power due to the decrease of the sample size (data not shown). In addition, no significant differences were found in allele or genotype distributions of the TNF-α -1031T > C and TNF-α -857 C > T polymorphisms between tendinopathy cases and controls (data not shown).
Seven haplotypes of the TNF-α (-1031T > C, -857 C > T, -308G > A) SNPs were inferred, which account 100 % of the study population. The TCG haplotype was considered wild-type/reference haplotype due to present the highest frequency in the study population (N = 320, 59.2 %). After adjusting for confounding variables (age, sex, sport modality, BMI and nutritional follow-up), the TNF-α TCA haplotype was associated with increased tendinopathies risk, either considering all cases (P = 0.006), patellar (P = 0.004), rotator cuff (P = 0.008) or Achilles tendinopathies (P = 0.03) (Table 3).
Tendinopathy is a serious public health care problem and the knowledge of molecular mechanisms involved in its etiology remains an active area of ongoing research . Any tendon can undergo a tendinopathy process and some modifiable and non-modifiable risk factors are common for the disease in different affected tendon . Recent studies have challenged the “degenerative process” paradigm, suggesting that tendon overload is linked to a complex role of inflammation on tendon homeostasis dysregulation [6, 9, 26, 28]. Although the role of the inflammatory process is not clear, the dysregulation of the proinflammatory cytokines expression and release may contribute to chronic inflammatory responses [9, 26].
Metabolic diseases related to increased adiposity has been identified as important potentially modifiable risk factor for the onset and progression of a variety of tendinopathies . Adipose tissue is tightly associated with tendon inflammation and early tissue degeneration . High BMI and nutritional follow-up by a specialist during a sports career were non-modifiable risk factors for tendinopathy in our athletes. The increased BMI in athletes can result in nutritional monitoring for muscle mass gain, which optimizes the athlete’s performance and physical ability ; however, nutritional supplements may be a key component in the etiology of various diseases [31, 32], and diet can contribute negatively with tendon homeostasis .
Despite the different risk factors different types of tendinopathies, overloading and mechanical stress may induce the secretion of TNF-α by tenocytes and cause change cellular proliferation, onset of pain and ECM degradation . Under normal physiological conditions TNF-α is not detectable in tendon; however, TNF-α was detect in human tenocytes of Achilles tendinopathy samples, suggesting association with onset tissue apoptosis and in mechanotransduction failure to adapt tendon load . The variation in TNF-α cytokines production is tightly regulated by genetic variants . The present results indicate a positive association between TNF-α TCA haplotype and the risk of developing tendinopathy (2-4-fold), which is observed when analyzing only the patellar, rotator cuff or Achilles subgroups. The TNF-α TCA (-1031T > C, -857 C > T and − 308G > A) haplotype characterized by the presence of the variant allele of TNF-α – 308 A, which promotes loss of transcription factors like activator protein-2 binding, increasing the level of gene transcription . The TNF-α SNPs in the promoter region site are strongly in linkage disequilibrium and creates established haplotypes that affect differently gene expression and activity than those of each SNP evaluated separately [33, 34]. This may explain the increased level of TNF-α mRNA found in the degenerate tendon [11, 13]; and consequently, contribute to inter-individual variation in tendinopathy development.
Within the Brazilian population the TNF-α-308AA genotype is rarer (approximately 0–2 %) [23, 35, 36], and was only observed in the tendinopathy cases (~ 2 %). The total sample size was adequate to detect significant associations with 80 % statistical power; the small number of athletes with different locations of tendinopathies was the main limitation of this study. However, the strength of this study included the control group was matched with all tendinopathy case for age, sex, and sport modality to minimize the influence of the confounding factors. The results can be used to build a database from different populations to identify modifiable and non-modifiable risk factors associated with tendinopathy development in athletes.
It is essential to understand the molecular mechanism involved in the etiology of the disease and for control mechanical stress on the tendon of athletes most likely to develop overuse injuries . The changes in the cytokine production due different genotypes can have significant influence in the tendinopathy, which can impair early tissue regeneration. Identifying genetic changes may improve the prognosis of the disease and clarify new therapeutic targets or personalized training for the athlete, avoiding movement limitations, loss of physical performance and sports ability. Athletes with TNF-α TCA haplotype could be early subjected to cryotherapy after training and competition to avoid tendinopathy development. Whole-body cryotherapy decreased serum TNF-α (around 60 %) 24 h following exercise . Thus, this finding can be used in future studies to better understand the influence of genetic factors in the tendinopathy susceptibility and contribute to create sports injury surveillance programmes using genetic information aim reduce cases of the illness in athletes.
The TNF-α -308G > A SNP was potential non-modifiable risk associated with development of disease.
Availability of data and materials
Original data are available as Supplementary file 1.
Docking SI, Rio E, Cook J, et al. The prevalence of Achilles and patellar tendon injuries in Australian football players beyond a time-loss definition. Scand J Med Sci Sports. 2018;28(9):2016–22. https://doi.org/10.1111/sms.13086.
Salles JI, Lopes LR, Duarte MEL, et al. Fc receptor-like 3 (-169T > C) polymorphism increases the risk of tendinopathy in volleyball athletes: a case control study. BMC Med Genet. 2018;19(1):119. https://doi.org/10.1186/s12881-018-0633-6.
Diniz-Fernandes T, Godoy-Santos AL, Santos MC, et al. Matrix metalloproteinase-1 (MMP-1) and (MMP-8) gene polymorphisms promote increase and remodeling of the collagen III and V in posterior tibial tendinopathy. Histol Histopathol. 2018;33(9):929–36. https://doi.org/10.14670/HH-11-982.
Posthumus M, Collins M, Cook J, et al. Components of the transforming growth factor-beta family and the pathogenesis of human Achilles tendon pathology–a genetic association study. Rheumatology. 2010;49(11):2090–7. https://doi.org/10.1093/rheumatology/keq072.
Spang C, Renström L, Alfredson H, et al. Marked expression of TNF receptors in human peritendinous tissues including in nerve fascicles with axonal damage - Studies on tendinopathy and tennis elbow. J Musculoskelet Neuronal Interact. 2017;17(3):226–36.
Gaida JE, Bagge J, Purdam C, et al. Evidence of the TNF-α system in the human Achilles tendon: expression of TNF-α and TNF receptor at both protein and mRNA levels in the tenocytes. Cells Tissues Organs. 2012;196(4):339–52. https://doi.org/10.1159/000335475.
Nourian M, Chaleshi V, Pishkar L, et al. Evaluation of tumor necrosis factor (TNF)-α mRNA expression level and the rs1799964 polymorphism of the TNF-α gene in peripheral mononuclear cells of patients with inflammatory bowel diseases. Biomed Rep. 2017;6(6):698–702. https://doi.org/10.3892/br.2017.908.
Santos CNO, Ribeiro DR, Cardoso Alves J, et al. Association Between Zika Virus Microcephaly in Newborns With the rs3775291 Variant in Toll-Like Receptor 3 and rs1799964 Variant at Tumor Necrosis Factor-α Gene. J Infect Dis. 2019;220(11):1797–801. https://doi.org/10.1093/infdis/jiz392.
Hassanzad M, Farnia P, Ghanavi J, et al. TNFα -857 C/T and TNFR2 + 587 T/G polymorphisms are associated with cystic fibrosis in Iranian patients. Eur J Med Genet. 2019;62(11):103584. https://doi.org/10.1016/j.ejmg.2018.11.018.
Goes RA, Lopes LR, Cossich VRA, et al. Musculoskeletal injuries in athletes from five modalities: a cross-sectional study. BMC Musculoskelet Disord. 2020;21(1):122. https://doi.org/10.1186/s12891-020-3141-8.
Salles JI, Duarte ME, Guimarães JM, et al. Vascular Endothelial Growth Factor Receptor-2 Polymorphisms Have Protective Effect against the Development of Tendinopathy in Volleyball Athletes. PLoS One. 2016;11(12):e0167717. https://doi.org/10.1371/journal.pone.0167717.
Lomelí-Nieto JA, Muñoz-Valle JF, Baños-Hernández CJ, et al. TNFA – 308G > A and – 238G > A polymorphisms and risk to systemic sclerosis: impact on TNF-α serum levels, TNFA mRNA expression, and autoantibodies. Clin Exp Med. 2019;19(4):439–47. https://doi.org/10.1007/s10238-019-00569-4.
Rocha Loures MA, Macedo LC, Reis DM, et al. Influence of TNF and IL17 Gene Polymorphisms on the Spondyloarthritis Immunopathogenesis, Regardless of HLA-B27, in a Brazilian Population. Mediators Inflamm. 2018;2018:1395823. https://doi.org/10.1155/2018/1395823.
Perini JA, Cardoso JV, Berardo PT, et al. Role of vascular endothelial growth factor polymorphisms (-2578 C > A, -460 T > C, -1154G > A, + 405G > C and + 936 C > T) in endometriosis: a case-control study with Brazilians. BMC Womens Health. 2014;14:117. https://doi.org/10.1186/1472-6874-14-117.
Loiacono C, Palermi S, Massa B, et al. Tendinopathy: Pathophysiology, Therapeutic Options, and Role of Nutraceutics. A Narrative Literature Review. Medicina. 2019;55(8):447. https://doi.org/10.3390/medicina55080447.
Hulens M, Vansant G, Lysens R, et al. Study of differences in peripheral muscle strength of lean versus obese women: an allometric approach. Int J Obes Relat Metab Disord. 2001;25(5):676–81. https://doi.org/10.1038/sj.ijo.0801560.
DePhillipo NN, Aman ZS, Kennedy MI, et al. Efficacy of Vitamin C Supplementation on Collagen Synthesis and Oxidative Stress After Musculoskeletal Injuries: A Systematic Review. Orthop J Sports Med. 2018;6(10):2325967118804544. https://doi.org/10.1177/2325967118804544.
Mandal RK, Khan MA, Hussain A, et al. A trial sequential meta-analysis of TNF-α -308G > A (rs800629) gene polymorphism and susceptibility to colorectal cancer. Biosci Rep. 2019;39(1):BSR20181052. https://doi.org/10.1042/BSR20181052.
Karimi M, Goldie LC, Cruickshank MN, et al. A critical assessment of the factors affecting reporter gene assays for promoter SNP function: a reassessment of -308 TNF polymorphism function using a novel integrated reporter system. Eur J Hum Genet. 2009;17(11):1454–62. https://doi.org/10.1038/ejhg.2009.80.
Oki E, Norde MN, Carioca AAF, et al. Polymorphisms of the TNF-α gene interact with plasma fatty acids on inflammatory biomarker profile: a population-based, cross-sectional study in São Paulo, Brazil. Br J Nutr. 2017;117(12):1663–73. https://doi.org/10.1017/S0007114517001416.
Patente TA, Monteiro MB, Vieira SM, et al. Linkage disequilibrium with HLA-DRB1-DQB1 haplotypes explains the association of TNF-308G > A variant with type 1 diabetes in a Brazilian cohort. Gene. 2015;568(1):50–4. https://doi.org/10.1016/j.gene.2015.05.017.
The authors thank the support of Sports Trauma Center of National Institute of Traumatology and Orthopaedics (INTO), Sports Training Center and National or State Sports Championships for an opportunity to recruit athletes. And the technical assistance of Ana Carolina Leocadio de Souza, Camili Gomes Pereira and Jade Pires do Nascimento from Laboratório de Pesquisa de Ciências Farmacêuticas, Centro Universitário Estadual da Zona Oeste (UEZO) and Research Division of INTO. This study was partially supported by the Coordination for the Improvement of Higher Education Personnel (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES) - Finance Code 001.
This study was supported by the Brazilian agency Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro - FAPERJ, Brazil. Funding body contributed to acquisition of research inputs.
Authors and Affiliations
Divisão de Pesquisa, Instituto Nacional de Traumatologia e Ortopedia (INTO), Avenida Brasil, 500, RJ, zip code 20940-070, Rio de Janeiro, Brazil
Lucas Rafael Lopes, João Antonio Matheus Guimarães, Victor Soares Wainchtock & Jamila Alessandra Perini
Laboratório de Pesquisa de Ciências Farmacêuticas, Centro Universitário Estadual da Zona Oeste (UEZO), Rio de Janeiro, Brazil
Lucas Rafael Lopes, Victor Soares Wainchtock & Jamila Alessandra Perini
Programa de Pós-graduação em Saúde Pública e Meio Ambiente, Escola Nacional de Saúde Pública, Fundação Oswaldo Cruz (Fiocruz), Rio de Janeiro, Brazil
Lucas Rafael Lopes & Jamila Alessandra Perini
Centro de Trauma do Esporte, Instituto Nacional de Traumatologia e Ortopedia (INTO), Rio de Janeiro, Brazil
Vitor Almeida Ribeiro de Miranda, Gabriel Garcez de Araujo Souza, João Alves Grangeiro Neto & Rodrigo de Araújo Goes
JAP participated in conception and design of study. LRL, VARM, GGAS, VWS, RAG and JAP collated the data and developed the database. LRL, VSW, and JAP helped to experiments. LRL performed the statistical analysis. LRL, VARM, GGAS, JAMG, JAGN and JAP analysis and interpretation of data. LRL, VARM and JAP wrote the manuscript. JAMG, JAGN, RAG and JAP critical revision of the manuscript for important intellectual content. All authors read and approved the final manuscript.
This study was approved by the Human Research Ethics Committee of the Instituto Nacional de Traumatologia e Ortopedia, Rio de Janeiro, Brazil (protocol number 2.455.630/2017). All participating provided written informed consent.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
Lopes, L.R., de Miranda, V.A.R., Guimarães, J.A.M. et al. Association of TNF-α -308G > A polymorphism with susceptibility to tendinopathy in athletes: a case–control study.
BMC Sports Sci Med Rehabil13, 51 (2021). https://doi.org/10.1186/s13102-021-00276-2