Skip to main content
Download PDF

A thorough analysis of data on the correlation between COL9A1 polymorphisms and the susceptibility to congenital talipes equinovarus: a meta-analysis

Journal of Orthopaedic Surgery and Research volume 19, Article number: 345 (2024) Cite this article

Abstract

Background

Congenital talipes equinovarus (CTEV) is a prevalent pediatric deformity with a multifactorial etiology. The objective of this meta-analysis was to explore the association between genetic variations in COL9A1 and the susceptibility to CTEV.

Methods

A comprehensive analysis of pertinent literature released before November 15, 2023, in electronic bibliographic databases was carried out. The importance of the connection was clarified through odds ratios (ORs) with 95% confidence intervals (CIs), utilizing random or fixed-effects models depending on study heterogeneity. Statistical analysis was executed using Comprehensive Meta-Analysis software (Version 4.0).

Results

A total of eight case-control studies involving 833 CTEV patients and 1280 healthy individuals were included in the analysis. Among these, four studies investigated the rs1135056 variant, encompassing 432 CTEV cases and 603 controls; two studies examined the rs35470562 variant, with 189 CTEV cases and 378 controls; and two studies explored the rs592121 variant, including 212 CTEV cases and 299 controls. The results revealed a significant association between the rs1135056 and rs35470562 polymorphisms in the COL9A1 gene, suggesting an increased risk of CTEV in the overall population. Conversely, no such association was found for the rs592121 variant.

Conclusion

Our findings reveal a substantial association between the genetic variants COL9A1 rs1135056 and rs35470562 and susceptibility to CTEV. Conversely, the variant rs592121 did not exhibit any corresponding link. However, the limitations imposed by the small study population have compromised the statistical reliability and generalizability of the results.

Introduction

Congenital talipes equinovarus (CTEV), commonly referred to as clubfoot, is a prevalent musculoskeletal anomaly in the pediatric population, with an estimated birth prevalence of approximately 1 in 1,000 individuals [1, 2]. The clinical presentation of CTEV includes adduction of the forefoot, plantar flexion of the ankle, excessive inward rotation of the tibia, and overpronation of the calcaneus [3]. While CTEV is typically idiopathic in nature, it can also be associated with other congenital malformations, such as arthrogryposis, myelomeningocele, spina bifida, tethered cord syndrome, and amniotic band syndrome [4,5,6]. The global incidence of CTEV is influenced by the birth rate, with an estimated 100,000 cases reported annually worldwide. Furthermore, bilateral involvement is observed in approximately 30–50% of cases [7, 8]. It is noteworthy that a significant proportion of affected individuals reside in low- and middle-income countries, particularly in Southeast Asia and Africa, where access to appropriate treatment options may be limited [8,9,10]. The higher incidence of CTEV in these regions is postulated to have a genetic basis [11, 12]. Over the past few decades, non-surgical interventions, specifically the Ponseti technique, have become the preferred approach to managing CTEV. This technique involves a series of manipulations and castings, followed by Achilles tenotomy performed percutaneously, and the subsequent use of a foot abduction brace to minimize the need for additional surgical procedures [13, 14]. The occurrence of CTEV has been reported to vary across different geographical regions, with the lowest prevalence documented among individuals of Chinese descent and the highest among those of Polynesian ethnicity [11, 15, 16]. Additionally, the female-to-male ratio in CTEV cases remains consistent across diverse ethnic groups, ranging from approximately 1:2.0 to 1:2.5 [1, 17]. Moreover, developmental dysplasia of the hip (DDH) is frequently observed in syndromic malformations, such as clubfoot, due to its impact on connective tissue [18]. This condition is linked to various complications, including instability, dislocation, and the early onset of osteoarthritis in the hip joint [19]. Investigating the correlation between clubfoot and DDH is an area of ongoing exploration and clinical significance. Although these congenital conditions may seem unrelated, emerging data suggests a potential connection, with certain studies indicating a higher occurrence of DDH among individuals with clubfoot [18, 19]. The precise mechanisms underlying this correlation are not yet fully understood, but researchers have proposed several theories, including shared genetic and environmental factors, as well as possible biomechanical interactions during prenatal development.

The etiology of CTEV, although speculative in nature, is widely believed to be multifactorial in origin [17, 20]. While genetics undoubtedly play a significant role, various environmental factors also exert influence on the development of CTEV [21,22,23]. Notably, maternal smoking, advanced maternal age, and maternal nutrition have all been linked to an increased risk of CTEV. Moreover, it is worth noting that the interplay between genetic and environmental factors can contribute to the observed phenotypic variability among individuals with CTEV, underscoring the importance of a multifactorial model in understanding the etiology of clubfoot [8, 17, 24]. Numerous studies have indeed demonstrated the substantial involvement of genetic factors in the development of CTEV. Although the exact mode of inheritance remains elusive, the recurrence patterns observed in families strongly suggest a multifactorial genetic basis characterized by both inherited and sporadic genetic mutations [25, 26]. However, it is important to acknowledge that research on the susceptibility genes of CTEV has made limited progress, and no specific gene has yet been identified as the CTEV gene [27, 28]. The heritability of CTEV is estimated to be approximately 65%, yet the precise patterns of heredity and penetrance remain unclear. Approximately 25% of all cases exhibit familial tendencies, and a sibling of an affected child has a 2–4% chance of developing CTEV [1, 29,30,31]. Furthermore, data from twin studies have revealed higher concordance rates among monozygotic twins (33%) compared to dizygotic twins (3%), indicating a more prominent role for heritability in the development of CTEV [13]. Presently, only one genome-wide association study (GWAS) has been conducted on CTEV, which identified genetic variants in the NCOR2 and ZNF664 genes as being associated with the pathogenesis of CTEV in non-Hispanic and European descent individuals with CTEV [32, 33].

The study of the genetics of CTEV is a captivating field of research. In humans, genetic variations in collagen IX genes are connected to the functional state of joint cartilages, while genetic abnormalities in collagen formation are associated with short stature, severe skeletal deformities, and chondrodysplasia [34]. Collagen IX is a heterodimer molecule composed of identical chains of α1, α2, and α3 (α1(IX)–α3(IX)), which are encoded by the genes collagen type 9 alpha 1 (COL9A1), alpha 2 (COL9A2), and alpha 3 (COL9A3) located on chromosome locus 6q13, respectively [35, 36]. Through the candidate cloning strategy, the COL9A1 gene in the candidate region was initially confirmed as an associated gene with the risk of CTEV. Furthermore, the transmission disequilibrium test analysis demonstrated a linkage between the COL9A1 gene and a predisposition to CTEV [37]. In abductor hallucis muscle samples of CTEV cases, an increased expression of COL9A1 was observed, indicating that the COL9A1 gene may be a significant predisposition gene for CTEV [38, 39]. Convincing evidence from case-control and family-based linkage studies has indicated an association between COL9A1 polymorphisms and the risk of CTEV. However, due to the sample sizes of individual studies, the results continue to be debatable. To the best of our knowledge, no relevant meta-analysis has examined the association of COL9A1 polymorphisms with the risk of CTEV. Therefore, the objective of this meta-analysis is to determine the association between genetic variants at the COL9A1 gene and the risk of CTEV.

Materials and methods

Search approaches

The ethical endorsement was not obligatory for the present investigation, as it constitutes a systematic review and meta-analysis. This study was conducted in accordance with the guidelines outlined in the Preferred Reporting Items for Systematic Reviews and Meta-Analyses. The bibliographic databases, including PubMed, Web of Science, Europe PMC, ResearchGate, Elsevier, Cochrane Library, EMBASE, SciELO, Google Scholar, Chinese National Knowledge Infrastructure (CNKI), Wanfang Data Company, Chaoxing, China/Asia On Demand (CAOD), Chinese Medical Citation Index (CMCI), VIP Information Consulting Company (VIP), Chinese Medical Current Contents (CMCC), Chinese Biomedical Database (CBD), and Weipu Periodical Database, were utilized to search for relevant studies that assessed the correlation between COL9A1 polymorphisms and the risk of Congenital Talipes Equinovarus (CTEV). The search period spanned from the inception of the databases until November 15, 2023. The search was conducted using a combination of specific keywords and Medical Subject Headings (MeSH) terms, including ‘Congenital Talipes Equinovarus’, ‘Talipes Equinovarus’, ‘CTEV’, ‘Clubfoot’, ‘Collagen alpha-1(IX)’, ‘COL9A1’, ‘Gene’, ‘Association’, ‘Correlation’, ‘Genetic’, ‘Predisposition’, ‘Polymorphism’, ‘DNA Sequence’, ‘Single-Nucleotide Polymorphism’, ‘SNPs’, ‘Genotype’, ‘Frequency’, ‘Mutation’, ‘Mutant’, ‘Allele’, ‘Variation’, and ‘Variant’. Additionally, a manual search of the references cited in the retrieved articles and reviews was conducted to identify any other relevant publications. The search was limited to human studies and did not impose any language restrictions.

Including and excluding criteria

All studies included in the analysis met the following criteria: (1) they were original studies with either a case-control or cohort design, published in various languages; (2) they investigated the correlation between COL9A1 polymorphisms and the risk of CTEV; (3) they provided sufficient and accessible data to calculate the odds ratios (ORs) and 95% confidence intervals (CIs). The primary exclusion criteria were: (1) studies without a control group; (2) case reports, case series, commentaries, letters, editorials, reviews, animal experiments, in vitro cell experiments, oral presentations, and previous meta-analyses; (3) studies with incomplete literature data or where the original text could not be obtained after contacting the original author; (4) inadequate or inaccessible data for analysis; (5) overlapping or duplicate data from the same population.

Data extraction

Two investigators independently examined the bibliographies, gathered data, and cross-referenced them according to the aforementioned inclusion and exclusion criteria. In case of discrepancies, resolution would be sought through debate or a meeting involving a third researcher. Additionally, contact with the original author via email would be made when deemed necessary. During the literature screening process, the title and abstract were initially reviewed to eliminate obviously irrelevant sources. Subsequently, the full text was carefully examined to determine eligibility for inclusion. The literature that met the required standards yielded the following key information: the first author’s name, ethnicity, publication date, country of origin, genotyping methods, the total number of cases and controls for CTEV, genotype frequencies of CTEV cases and healthy controls for each available COL9A1 polymorphism, Hardy-Weinberg equilibrium (HWE) test results, and minor allele frequencies (MAFs) in the healthy control group. In cases where multiple studies were conducted by the same investigator(s) and involved duplicated or overlapping records, only the most recent published data or the study with the largest sample size was considered for incorporation into the analysis.

Statistical analysis

The HWE in healthy subjects for each individual study was assessed using the Pearson’s chi-square Statistic, employing free online software. A p-value below 0.05 was deemed to be statistically significant. The association between COL9A1 polymorphisms and susceptibility to CTEV was quantified using odds ratios (ORs) with 95% confidence intervals (CIs). To determine the statistical significance of the combined data, the Z-test was applied to ascertain the difference between the population mean and sample mean [40]. This meta-analysis incorporated five genetic models, including allelic (B vs. A), heterozygote (BA vs. AA), homozygote (BB vs. AA), recessive (BB vs. BA + AA), and dominant (BB + BA vs. AA). The chi-square test is commonly utilized to evaluate the significance of heterogeneity, with a significance level of p < 0.05. Furthermore, based on Cochrane, the level of heterogeneity between studies was classified on a scale of 0 to 100% [41, 42]. The degree of heterogeneity was quantified using the I2 index, which represents the proportion of total variation between studies attributable to their distinct characteristics. Therefore, if I2 exceeded 50%, a random-effect model (DerSimonian-Laird method) was employed. Otherwise, fixed-effect models (Mantel-Haenszel method) were utilized for interpretation [42, 43]. Sensitivity analysis was performed by systematically excluding one study at a time to assess the robustness of the findings. To assess publication bias, Begg’s test was employed by plotting the standard error of each study against its odds ratio, and Egger’s test was conducted alongside visual examination of the funnel plot for asymmetry. If publication bias was detected, the trim-and-fill method was employed to adjust the conclusions accordingly [44, 45]. The primary studies were synthesized using the Comprehensive meta-analysis software (Version 4.0) developed by Biostat, USA. A p-value below 0.05 on either end was deemed to possess statistical importance.

Results

Study characteristics

Detailed flowcharts illustrating the complex selection process are presented in Fig. 1. Initially, a comprehensive search was conducted across multiple electronic databases, accompanied by a thorough manual review of the complete contents on a page-by-page basis, leading to the identification of 328 papers. Subsequently, by carefully evaluating the research titles and abstracts, 117 duplicate documents and 142 irrelevant articles were eliminated, resulting in 69 articles for detailed examination. Ultimately, a total of eight case-control studies, involving 833 cases with CTEV and 1280 healthy individuals, were selected, meeting our inclusion criteria, as detailed in four publications on the rs1135056 (c.1862 A > G, p.Gln621Arg), rs35470562 (c.2050G > A, p.Glu441Lys), and rs592121 (c.1015T > C, p.Ser339Pro) genetic variations. Among these investigations, four focused on rs1135056, with 432 CTEV cases and 603 controls, while two studies examined rs35470562, including 189 CTEV cases and 378 controls, and another two studies investigated rs592121, encompassing 212 CTEV cases and 299 controls. The selected studies were published between 2011 and 2022 in English and Chinese, and were conducted on populations of Asian descent in China and India. The genotyping technique employed for assessing the genetic variations in all studies included in this analysis utilized the polymerase chain reaction restriction fragment length polymorphism (PCR-RFLP) method. It is noteworthy that the genotype distribution for rs1135056 in the control cohort within a study deviated from Hardy-Weinberg Equilibrium principles (p = 0.005). Furthermore, Table 1 delineates the frequencies of the minor allele for each genetic variation examined in this investigation within the control cohort.

Fig. 1
figure 1

Flow diagram of the study selection process

Full size image
Table 1 Characteristics of the studies included in meta-analysis
Full size table

Quantitative synthesis

The key findings regarding the association between COL9A1 polymorphisms and susceptibility to CTEV are concisely presented in Table 2. Overall, a notable connection was observed between the COL9A1 rs1135056 polymorphism and the susceptibility to CTEV under the homozygote model (GG vs. AA: OR = 1.487, 95% CI 1.023–2.160, p = 0.037, Fig. 2). Importantly, subgroup analysis based on the country of origin revealed a significant correlation among Chinese individuals across four genetic models: allele (G vs. A: OR = 1.431, 95% CI 1.158–1.769, p = 0.001, Fig. 3A), homozygote (GG vs. AA: OR = 1.885, 95% CI 1.233–2.882, p = 0.003, Fig. 3B), heterozygote (GA vs. AA: OR = 1.565, 95% CI 1.129–2.170, p = 0.007, Fig. 3C), and dominant (GG + GA vs. AA: OR = 1.639, 95% CI 1.210–2.221, p = 0.001, Fig. 3D). Regarding rs35470562, a significant correlation was observed across three genetic models: allele (A vs. G: OR = 0.667, 95% CI 1.038–2.131, p = 0.031), homozygote (AA vs. GG: OR = 2.600, 95% CI 1.156–5.850, p = 0.021), and recessive (AA vs. AG + GG: OR = 2.225, 95% CI 1.054–4.699, p = 0.036). However, the combined data analysis showed no significant correlation between COL9A1 rs592121 and the susceptibility to CTEV across all five genetic models in the overall population.

Table 2 Summary risk estimates for association between COL9A1 polymorphisms and risk of CTEV
Full size table
Fig. 2
figure 2

Forest plots for the association of COL9A1 rs1135056 polymorphism and risk of CTEV under the homozygote model (GG vs. AA)

Full size image
Fig. 3
figure 3

Forest plots for the association of COL9A1 rs1135056 polymorphism and risk of CTEV among Chinese. A: allele (G vs. A); B: homozygote (GG vs. AA); C: heterozygote (GA vs. AA); and D: dominant (GG + GA vs. AA)

Full size image

Sensitivity analysis

A sensitivity analysis was conducted to assess whether the combined data were significantly affected by the presence of any individual study. This was achieved by systematically removing each study and recalculating the significance of the findings. The results showed that the fixed-effects and/or random-effects estimates remained largely consistent before and after the removal of each study, suggesting that the combined odds ratio estimates were stable. Additionally, a sensitivity analysis was performed by excluding studies in which the controls did not adhere to HWE. The analysis revealed that there was no heterogeneity observed before or after the exclusion of these studies.

Heterogeneity test

As presented in Table 2, the current meta-analysis findings indicate a moderate to high level of variation among the included studies for the rs1135056 polymorphism under three genetic models overall. However, this variation was not observed for the rs35470562 and rs592121 polymorphisms. Consequently, we conducted a subgroup analysis based on the country of origin, revealing that the country of origin is a significant source of the variation observed for this specific polymorphism (Table 2). This study provides an interesting perspective and serves as a compelling example of how ethnicity can influence the variation observed between studies.

Publication bias

The intuitive evaluation of publication bias was carried out by constructing a Begg’s funnel plot, while the quantitative analysis was performed using Egger’s test. This statistical method involved plotting the standard error of the natural logarithm of the odds ratios against the odds ratios for each study included in the meta-analysis. As per the commonly acknowledged interpretation, the identification of selection bias would be reflected in an asymmetric plot, indicating a significant influence on the overall meta-analysis. Nevertheless, the shapes of the Begg’s funnel plots concerning the rs1135056, rs35470562, and rs592121 polymorphisms, across all five genetic models, did not display any signs of publication bias. This lack of bias was further validated through the implementation of the Egger test, as illustrated in Table 2.

Discussion

There exists a multitude of theories regarding the etiology and pathogenesis of CTEV. These theories primarily revolve around abnormal bone development, neuromuscular disease, soft tissue contracture, vascular abnormalities, intrauterine growth retardation, maternal smoking, conception month, birth season, and other factors. However, it should be noted that there is a dearth of substantial evidence supporting these theories [1, 46, 47]. It has been reported that the development of clubfoot involves various genes, such as PITX1, HOXD9, and TBX4, thereby emphasizing the intricate nature of genetic contributions to this condition [48, 49]. These findings serve to underscore that CTEV is influenced by multiple genetic variants, highlighting the significance of considering a comprehensive range of genetic factors when attempting to comprehend its genetic architecture. In several prior studies, the ICTEV susceptibility gene was mapped to the 6q12-13 region. Within this candidate region, several known genes are present, including COL9A1, PGM3, BAI3, COL19A1, and COL12A1. Notably, the COL9A1 gene is situated among these candidate genes within the aforementioned region [50, 51].

The Type IX collagen exhibits the capability to form associations with various molecules such as tissue inhibitors of matrix metalloproteinases, surface receptors, and growth factors present on the chondrocyte membrane. The significance of its function lies in its ability to stabilize the internal environment of articular cartilage and provide protection to it [52]. The α1 chain, which is a major component of type IX collagen, is found in hyaline cartilage [53]. The COL9A1 gene was originally identified by Kimura et al. in 1989 [54] and its location was determined to be at 6q12-q13 by Warman et al. in 1993 using fluorescence in situ hybridization [55]. Two transcripts of different lengths exist, with full lengths of 3704 bp and 2985 bp, and containing 38 and 32 exons, respectively. Current research suggests that mutations in the type IX collagen gene can lead to mild achondroplasia. Additionally, mutations in exon 9 of the COL9A1 gene have been linked to Kashin-Beck disease, congenital epiphyseal dysplasia, osteoarthritis, and Stickler syndrome. A study conducted by Shi et al. in 2015 involving 274 cases of Kashin-Beck disease and 248 disease-free children revealed a significant correlation between COL9A1 rs6910140 and the risk of this pathology [56]. However, other studies have reported inconsistent findings. For instance, Czarny-Ratajczak et al. in 2001 failed to establish a substantial correlation between genetic variants in the COL9A1 gene and the risk of multiple epiphyseal dysplasia in the Finnish population [57]. Similarly, Jakkula and colleagues (2005) also reported no statistically significant association between various cartilage collagen genes, such as COL2A1, COL9A1, COL9A2, and COL9A3, and the occurrence of osteoarthritis in infants [58].

In the present investigation, a collective sum of eight case-control studies, comprising a total of 833 cases of CTEV and 1280 control individuals, were incorporated. Four of the studies examined rs1135056, with a total of 432 cases of CTEV and 603 controls. Another two studies specifically investigated rs35470562, with 189 CTEV cases and 378 controls. Additionally, two studies focused on rs592121, comprising 212 CTEV cases and 299 controls. Our amalgamated data implies that the COL9A1 rs1135056 and rs35470562 polymorphisms exhibit a statistically significant association with susceptibility to CTEV. Nevertheless, no significant correlation was found between the rs592121 polymorphism and the risk of CTEV. Previously, various studies have reported that polymorphisms in the COL9A1 gene are associated with the risk of CTEV. In particular, the COL9A1 rs1135056 has been found to have a significant relationship with CTEV in both children and their biological mothers [54]. Furthermore, it has been observed that individuals with the GG genotype for rs592121 have a higher likelihood of developing CTEV compared to those with other genotypes [59]. Another study has also indicated that the AA genotype of COL9A1 rs35470562 is significantly linked to an increased risk of CTEV when compared to the GG genotype [60]. However, no significant association was established between the COL9A1 rs1135056 and rs592121 polymorphisms and CTEV in any of the genetic models tested [38]. These findings suggest that COL9A1 polymorphisms, particularly rs1135056 and rs35470562, may play a role in the development of CTEV.

In 2011, Liu et al. conducted an examination on the correlation between COL9A1 polymorphism and susceptibility to CTEV. In their analysis, which included 118 CTEV cases and 100 normal cases, the authors proclaimed that the COL9A1 protein is substantially overexpressed in patients with CTEV, and the rs1135056 polymorphism is associated with its pathogenesis [38]. In a distinct investigation carried out by Liu et al. in the year 2007, two variations of genetic makeup, specifically rs592121 and rs1135056, were evaluated within the COL9A1 gene in a total of 84 nuclear pedigrees related to individuals with ICTEV. The findings of the study revealed that there was an observed transmission disequilibrium in the pedigrees for the rs592121 and rs1135056 polymorphisms within the COL9A1 gene. Furthermore, the expression levels of COL9A1 mRNA were significantly elevated in patients with ICTEV compared to individuals without the condition [61]. In the year 2022, Raikwar et al. conducted a study on a dyad consisting of a mother and child to investigate the impact of two significant polymorphisms, namely rs592121 and rs1135056, at the COL9A1 gene on the development of CTEV. The results of their research revealed a potential link between the rs592121 gene polymorphism and an increased susceptibility to CTEV in patients [54]. In 2016, the study conducted by Zhao et al. investigated the possible association between three genetic polymorphisms, namely COL9A1 rs1135056, rs35470562, and rs592121, and the risk of CTEV in a sample of 87 cases and 174 controls recruited from the Fourth People’s Hospital of Shaanxi and the First Hospital of Yulin. The authors performed a conditional regression analysis and found that individuals with the AA genotype of COL9A1 rs35470562 exhibited a significantly higher risk of developing CTEV compared to those with the GG genotype (OR = 2.60, 95% CI = 1.06–6.32). Furthermore, when using a recessive model, it was observed that carriers of the rs35470562 AA genotype are at a greater risk for CTEV in comparison to individuals with either the GG or GA genotypes (OR = 2.23, 95%CI = 1.03–5.04). However, no significant relationship was found between the COL9A1 rs1135056 and rs592121 polymorphisms and CTEV in Chinese children [59]. In another investigation involving 692 cases and 1513 healthy controls within the Han Chinese population, Zhao et al. conducted an assessment on the impact of COL9A1 gene variations on the vulnerability to CTEV. The findings of the study demonstrated that while there was no correlation observed between the genotyped SNPs and CTEV, an indication of a gene-environment interaction emerged between the SNP rs6455357 and maternal alcohol consumption [62]. In 2016, Zhao et al. reported a lack of association between the COL9A1 rs1135056 and rs592121 polymorphisms and the risk of CTEV, yet a significant correlation was identified between the rs35470562 polymorphism and this particular disease [59]. Additionally, several investigations have genotyped the COL9A1 rs1135056, rs35470562, and rs592121 using RFLP-PCR and have proposed a positive correlation between COL9A1 polymorphisms and the risk of CTEV within the Chinese population. For instance, two studies published in 2016 have demonstrated that the COL9A1 rs1135056 polymorphism may have an influence on the development of CTEV in the Chinese population [63, 64].

To the best of our knowledge, this meta-analysis represents the initial aggregation of data on the role of COL9A1 polymorphisms in the development of CTEV. The analysis was conducted through a comprehensive search of the available literature. However, there are some limitations to this meta-analysis. Firstly, the analysis was constrained by a relatively small number of studies, which limited the statistical power and generalizability of the findings. Additionally, the quality of the included studies varied, with some having small sample sizes and potential biases. Furthermore, the analysis was limited to published studies, which raises the possibility of publication bias. Secondly, most of the studies focused on Chinese individuals, limiting our ability to evaluate the potential impact of COL9A1 polymorphisms on other ethnic groups. This weakness weakens the true association between COL9A1 polymorphisms and the risk of CTEV. Therefore, it is crucial to conduct large, multi-center studies with sufficient power to thoroughly investigate these associations in various ethnic groups. Thirdly, this meta-analysis only included published papers available in the online network, and it is possible that some unpublished studies with null results were overlooked, potentially biasing the pooled data. Fourthly, we were unable to obtain data on various confounding factors known to influence the development of CTEV, and a more rigorous analysis could be performed if the selected studies provided adjusted evaluations. Lastly, despite significant advancements, there are still numerous challenges in elucidating the genetic basis of CTEV. The complex interplay between genetic and environmental factors, as well as the phenotypic variability, pose a formidable challenge. The meta-analysis focused solely on the association between COL9A1 polymorphisms and CTEV, and did not investigate the potential influence of other genetic or environmental factors on the development of this condition due to limited original data from the included studies.

Our aggregated data indicates that the COL9A1 rs1135056 (c.1862 A > G, p.Gln621Arg) and rs35470562 (c.2050G > A, p.Glu441Lys) variations are significantly associated with susceptibility to CTEV. However, the rs592121 (c.1015T > C, p.Ser339Pro) polymorphism does not show the same association. Nevertheless, caution is warranted when interpreting the findings of this meta-analysis, as our results only allow for the examination of the effects of COL9A1 polymorphisms in Asian children with CTEV. Understanding the genetic basis of CTEV has important implications for clinical practice. Genetic testing can facilitate early diagnosis and risk assessment for families with a history of clubfoot. Identifying the responsible genes opens up the possibility of gene-targeted therapies and personalized treatment options, leading to improved patient outcomes. Therefore, further studies with larger sample sizes and sufficient statistical power across diverse ethnic groups are necessary to gain a better understanding of the role of COL9A1 polymorphisms in the development of CTEV and to consider the combined impact of genetic and environmental factors on its etiology.

Data availability

No datasets were generated or analysed during the current study.

References

  1. Mustari MN, Faruk M, Bausat A, Fikry A. Congenital talipes equinovarus: a literature review. Annals Med Surg. 2022;81:104394.

    Article  Google Scholar 

  2. Panza R, Albano F, Casto A, Del Vecchio C, Laforgia N, Dibello D. Incidence and prevalence of congenital clubfoot in Apulia: a regional model for future prospective national studies. Italian J Pediatr 2023. 2023;49(1):49:1–6.

    Google Scholar 

  3. Cady R, Hennessey TA, Schwend RM. Clinical reports: diagnosis and treatment of idiopathic congenital clubfoot. Pediatrics. 2022;149:e2021055555.

    Article  PubMed  Google Scholar 

  4. Tileston K, Baskar D, Frick SL. What is new in pediatric orthopaedic foot and ankle. J Pediatr Orthop. 2022;42:E448–52.

    Article  PubMed  Google Scholar 

  5. Sun M, Tao B, Gao G, Shang A. Tethered cord syndrome associated with spina bifida and intradural angiolipoma coexisting with intramedullary arteriovenous malformation. Clin Neuroradiol. 2022;32:1127–30.

    Article  PubMed  Google Scholar 

  6. Trout SM, Whitaker AT. Management issues of congenital talipes equinovarus in the neonatal intensive care unit: a systematic review. Foot Ankle Surg. 2021;27:480–5.

    Article  PubMed  Google Scholar 

  7. Mathias RG, Lule JK, Waiswa G, Naddumba EK, Pirani S. Incidence of clubfoot in Uganda. Can J Public Health. 2010;101:341–4.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Smythe T, Rotenberg S, Lavy C. The global birth prevalence of clubfoot: a systematic review and meta-analysis. eClinicalMedicine. 2023;63:102178.

    Article  PubMed  PubMed Central  Google Scholar 

  9. Smythe T, Kuper H, Macleod D, Foster A, Lavy C. Birth prevalence of congenital talipes equinovarus in low- and middle-income countries: a systematic review and meta-analysis. Trop Med Int Health. 2017;22:269–85.

    Article  PubMed  Google Scholar 

  10. van Wijck SFM, Oomen AM, van der Heide HJL. Feasibility and barriers of treating clubfeet in four countries. Int Orthop. 2015;39:2415–22.

    Article  PubMed  Google Scholar 

  11. Todd MRLC, Kodovaru SS, Antoniou G, Cundy PJ. Clubfoot deformity in the Solomon Islands: melanesian versus Polynesian ethnicity, a retrospective cohort study. J Children’s Orthop. 2020;14:281–5.

    Article  Google Scholar 

  12. Vakili Ojarood M, Khanghah AS, Belalzadeh M. Gangrenous ischemic colitis due to acute promyelocytic leukaemia, and myelofibrosis in a 62-year-old man suffering from ESRD; Case report. Int J Surg Case Rep. 2021;89:106663.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Shabtai L, Segev E, Yavor A, Wientroub S, Hemo Y. Prolonged use of foot abduction brace reduces the rate of surgery in Ponseti-treated idiopathic club feet. J Children’s Orthop. 2015;9:177.

    Article  CAS  Google Scholar 

  14. Vakili M, Shahbazzadegan B. Esophageal rupture; a case report. J Isfahan Med School. 2016;33:2102–4.

    Google Scholar 

  15. Sadler B, Gurnett CA, Dobbs MB. The genetics of isolated and syndromic clubfoot. J Children’s Orthop. 2019;13:238.

    Article  CAS  Google Scholar 

  16. Haghighikian M, Khaleghian M, Saberi A, Nasiri M. A survey of knowledge level about pediatric oral/dental health among pediatricians. Indian J Forensic Med Toxicol. 2021;15:2686–90.

    Google Scholar 

  17. Pavone V, Chisari E, Vescio A, Lucenti L, Sessa G, Testa G. The etiology of idiopathic congenital talipes equinovarus: a systematic review. J Orthop Surg Res. 2018;13:1–11.

    Article  Google Scholar 

  18. Harsanyi S, Zamborsky R, Krajciova L, Kokavec M, Danisovic L. Developmental dysplasia of the hip: a review of etiopathogenesis, risk factors, and genetic aspects. Med (Kaunas Lithuania). 2020;56:153.

    Google Scholar 

  19. Chou DTS, Ramachandran M. Prevalence of developmental dysplasia of the hip in children with clubfoot. J Children’s Orthop. 2013;7:263.

    Article  Google Scholar 

  20. Cardy AH, Barker S, Chesney D, Sharp L, Maffulli N, Miedzybrodzka Z. Pedigree analysis and epidemiological features of idiopathic congenital talipes equinovarus in the United Kingdom: a case-control study. BMC Musculoskelet Disord. 2007;8:62.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Sharp L, Miedzybrodzka Z, Cardy AH, Inglis J, Madrigal L, Barker S, et al. The C677T polymorphism in the methylenetetrahydrofolate reductase gene (MTHFR), maternal use of folic acid supplements, and risk of isolated clubfoot: a case-parent-triad analysis. Am J Epidemiol. 2006;164:852–61.

    Article  PubMed  Google Scholar 

  22. Barker S, Chesney D, Miedzybrodzka Z, Maffulli N. Genetics and epidemiology of idiopathic congenital talipes equinovarus. J Pediatr Orthop. 2003;23:265–72.

    Article  PubMed  Google Scholar 

  23. Chesney D, Barker S, Miedzybrodzka Z, Haites N, Maffulli N. Epidemiology and genetic theories in the etiology of congenital talipes equinovarus. Bulletin: Hosp Joint Dis. 1999;58:59–64.

    CAS  Google Scholar 

  24. Hedley PL, Lausten-Thomsen U, Conway KM, Hinds K, Romitti PA, Christiansen M. Trends in congenital clubfoot prevalence and co-occurring anomalies during 1994–2021 in Denmark: a nationwide register-based study. BMC Musculoskelet Disord. 2023;24:784.

    Article  PubMed  PubMed Central  Google Scholar 

  25. Muhammad H, Haryana SM, Magetsari R, Kurniawan A, Baikuni B, Saraswati PA. Genetic role in recurrence of idiopathic CTEV: a systematic review. Orthop Res Reviews. 2023;15:19.

    Google Scholar 

  26. Vakili Ojarood M, Farzan R. Neural stem cell-conditioned medium and burn wound: a hopeful therapeutic approach to heal burn wounds. Burns. 2024;50:776–8.

    Article  PubMed  Google Scholar 

  27. Sadler B, Haller G, Antunes L, Nikolov M, Amarillo I, Coe B, et al. Rare and de novo duplications containing SHOX in clubfoot. J Med Genet. 2020;57:851.

    Article  CAS  PubMed  Google Scholar 

  28. Vakili Ojarood M, Yaghoubi T, Farzan R. Machine learning for prehospital care of patients with severe burns. Burns. 2024;24:00056–1.

    Google Scholar 

  29. Leyne E, Anselem O, Jordan P, Vivanti AJ, Benachi A, Salomon L, et al. Prenatal diagnosis of isolated bilateral clubfoot: is amniocentesis indicated? Acta Obstet Gynecol Scand. 2023;103:51–8.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Hasanpour Dargah M, Samadi N, Vakili J, Isazadefar K, Kebar SM, Zade ARM, et al. Comparative analysis of the effects of vasoperssin and norepinephrine on the renal function in patients undergoing CABG; a randomized clinical trial. Iran Red Crescent Med J. 2018;20:e67026.

    Google Scholar 

  31. Hoseini M, Negahi A, Vosough F, Farazmand B, Vosough F, Farazmand B, et al. Solitary fibrous tumor in pelvis extended to transverse mesocolon and peritoneum. Res J Pharm Technol. 2020;13:1941–8.

    Google Scholar 

  32. Zhang TX, Haller G, Lin P, Alvarado DM, Hecht JT, Blanton SH, et al. Genome-wide association study identifies new disease loci for isolated clubfoot. J Med Genet. 2014;51:334–9.

    Article  CAS  PubMed  Google Scholar 

  33. Kargar S, Dalimi A, Eslami G, Hajimohammadi B, Shirinzadeh A, Amouei A, et al. Cystic echinococcosis in Central Iran: G1 and G6 genotypes in patients. Surg Infect. 2022;23:451–7.

    Article  Google Scholar 

  34. Karimian F, Ghaderi H, Aminian A, Haghighikian SM, Mirjafari SA. Replantation of traumatic limb amputation above the elbow: a report of 4 cases. Tehran Univ Med Sci J. 2014;71:745–51.

    Google Scholar 

  35. Heilig J, Dietmar HF, Brachvogel B, Paulsson M, Zaucke F, Niehoff A. Collagen IX deficiency leads to premature vascularization and ossification of murine femoral heads through an imbalance of pro- and antiangiogenic factors. Osteoarthr Cartil. 2020;28:988–99.

    Article  CAS  Google Scholar 

  36. Asadollahi S, Neamatzadeh H, Namiranian N, Sobhan MR. Genetics of Legg-Calvé-Perthes disease: a review study. J Pediatr Rev. 2021;9:301–8.

    Article  Google Scholar 

  37. Ding J, Liang Z, Feng W, Cai Q, Zhang Z. Integrated bioinformatics analysis reveals potential pathway biomarkers and their interactions for clubfoot. Med Sci Monitor: Int Med J Experimental Clin Res. 2020;26:e925249–1.

    Article  CAS  Google Scholar 

  38. Liu L-Y, Jin C-L, Jiang L, Lin C-K. Expression of COL9A1 gene and its polymorphism in children with idiopathic congenital talipes equinovarus. Zhongguo Dang Dai Er Ke Za Zhi. 2011;13:478–81.

    CAS  PubMed  Google Scholar 

  39. Ahmad S, Reza E, Abdolhamid A, Jamal J, REZA SH, Shadi K, et al. Severe main bronchus obstruction due to pulmonary schwannoma: a case report. Iran J Ped Hematol Oncol 2022. 2022;12:140–4.

    Google Scholar 

  40. Gohari M, Neámatzadeh H, Jafari MA, Mazaheri M, Zare-Shehneh M, Abbasi-Shavazi E. Association between the p53 codon 72 polymorphism and primary open-angle glaucoma risk: meta-analysis based on 11 case-control studies. Indian J Ophthalmol. 2016;64:756–61.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Mazaheri M, Shahdadi V, Boron AN. Molecullar and biochemical effect of alcohlic extract of alpinia galanga on rat spermatogenesis process. Iran J Reprod Med. 2014;12:765–70.

    PubMed  PubMed Central  Google Scholar 

  42. Babakhanzadeh E, Khodadadian A, Nazari M, Tezerjani MD, Aghaei SM, Ghasemifar S, et al. Deficient expression of DGCR8 in human testis is related to spermatogenesis dysfunction, especially in meiosis I. Int J Gen Med. 2020;13:185–92.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Mirjalili SA, Moghimi M, Aghili K, Jafari M, Abolbaghaei SM, Neamatzadeh H, et al. Association of promoter region polymorphisms of interleukin-10 gene with susceptibility to colorectal cancer: a systematic review and meta-analysis. Arq Gastroenterol. 2018;55:306–13.

    Article  PubMed  Google Scholar 

  44. Motamedi S, Majidzadeh K, Mazaheri M, Anbiaie R, Mortazavizadeh SMR, Esmaeili R. Tamoxifen resistance and CYP2D6 copy numbers in breast cancer patients. Asian Pac J Cancer Prev. 2012;13:6101–4.

    Article  PubMed  Google Scholar 

  45. Karimi-Zarchi M, Sayad S, Vakili M, Shirinzadeh-Dastgiri A, Naseri A, Antikchi MH, et al. A collect of recommendations and guidelines for management and treatment of underlying malignancies during the COVID-19 pandemic. Acta Medica Iranica. 2023;61:443–8.

    Google Scholar 

  46. Cai G, Yang X, Chen T, Jin F, Ding J, Wu Z. Integrated bioinformatics analysis of potential pathway biomarkers using abnormal proteins in clubfoot. PeerJ. 2020;2020:e8422.

    Article  Google Scholar 

  47. Khosravi M, Jahanshahi F, Tafreshian A, Kian SMH. Medical evaluation of mortality and complications of thoracoscopic surgery with the help of primary video in patients with thoracic trauma with stable hemodynamics, 2019–2020. J Med Chem Sci. 2022;5:361–8.

    Google Scholar 

  48. Razmpoosh E, Safi S, Mazaheri M, Salehi-Abargouei A, Abdollahi N, Nazari M, et al. Effects of oral Nigella sativa oil on the expression levels and serum concentrations of adiponectin, PPAR-γ, and TNF-α in overweight and obese women: a study protocol for a crossover-designed, double-blind, placebo-controlled randomized clinical trial. Trials. 2019;20:512.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Bianco AM, Ragusa G, Di Carlo V, Faletra F, Di Stazio M, Racano C, et al. What is the exact contribution of PITX1 and TBX4 genes in clubfoot development? An Italian study. Genes. 2022;13:1958.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Li-ying L, Chun-lian J, Chang-kun L. Genetic polymorphism of 7 STR loci in idiopathic congenital talipes equinovarus predisposing genes. Chin J Public Health. 2010;26:1544–5.

    Google Scholar 

  51. Li-ying L, Chun-lian J, Dong-Hua C. Susceptibility gene location of idiopathic congenital talipes equinovarus. Chin J Public Health. 2010;26:433–5.

    Google Scholar 

  52. Michelacci YM, Baccarin RYA, Rodrigues NNP. Chondrocyte homeostasis and differentiation: transcriptional control and signaling in healthy and osteoarthritic conditions. Life. 2023;13:1460.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Alcaide-Ruggiero L, Molina-Hernández V, Granados MM, Domínguez JM. Main and minor types of collagens in the articular cartilage: the role of collagens in repair tissue evaluation in chondral defects. Int J Mol Sci. 2021;22:13329.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Raikwar A, Singh A, Mahdi AA, Yadav M. Assessment of the COL9A1 single nucleotide polymorphism with severity of clubfoot in a paediatric population along with their biological mothers. J Clin Orthop Trauma. 2022;35:102064.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Warman ML, Tiller GE, Polumbo PA, Seldin MF, Rochelle JU, Knoll JHM, et al. Physical and linkage mapping of the human and murine genes for the α1 chain of type ix collagen (col9a1). Genomics. 1993;17:694–8.

    Article  CAS  PubMed  Google Scholar 

  56. Wei SX, Feng Z, Yan LZ, Li LA, Xiong G, Wei SX, et al. Polymorphism in rs2229783 of the alpha 1(XI) collagen gene is associated with susceptibility to but not severity of kashin-beck disease in a northwest. Biomed Environ Sci. 2018;31:322–6.

    Google Scholar 

  57. Czarny-Ratajczak M, Lohiniva J, Rogala P, Kozlowski K, Perälä M, Carter L, et al. A mutation in COL9A1 causes multiple epiphyseal dysplasia: further evidence for locus heterogeneity. Am J Hum Genet. 2001;69:969–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Jakkula E, Melkoniemi M, Kiviranta I, Lohiniva J, Räinä SS, Perälä M, et al. The role of sequence variations within the genes encoding collagen II, IX and XI in non-syndromic, early-onset osteoarthritis. Osteoarthr Cartil. 2005;13:497–507.

    Article  CAS  Google Scholar 

  59. Zhao XL, Wang YJ, Wu YL, Han WH. Role of COL9A1 genetic polymorphisms in development of congenital talipes equinovarus in a Chinese population. Genet Mol Res. 2016;15:gmr15048773.

    Google Scholar 

  60. Dim E, Edagha I, Peter A, Umoh I, Ituen A, Oforjigha Dim C, et al. Congenital talipes equinovarus: a review. J Orthop Surg Techniques. 2022;5:474–86.

    Article  Google Scholar 

  61. Liu LY, Jin CL, Cao DH, Zhao N, Lin CK, Sun KL. Analysis of association between COL9A1 gene and idiopathic congenital talipes equinovarus. Yi Chuan = Hereditas. 2007;29:427–32.

    Article  CAS  PubMed  Google Scholar 

  62. Zhao J, Cai F, Liu P, Wei J, Chen Q. Gene environment interactions between the COL9A1 gene and maternal drinking of alcohol contribute to the risk of congenital talipes equinovarus. Genetic Test Mol Biomarkers. 2021;25:48–54.

    Article  CAS  Google Scholar 

  63. Li F, Fu J, Zhang X, Ping D, Song S, Zhao Z, et al. COL9A1 rs1135056 genetic polymorphism contributes to the susceptibility to congenital talipes equinovarus in a Chinese population. Int J Clin Exp Med. 2016;9:18740–5.

    CAS  Google Scholar 

  64. Li H, Sun Y, Chen F, Zhao F, Zhang J, Wu J, et al. Role of polymorphisms in COL9A1 gene in the development of congenital talipes equinovarus in a Chinese population. Int J Clin Exp Pathol. 2016;9:5712–6.

    CAS  Google Scholar 

Download references

Funding

There is no funding source.

Author information

Author notes

  1. Kamran Alijanpour and Reza Bahrami are equally attributed as corresponding authors.

Authors and Affiliations

  1. Department of Pediatrics, Islamic Azad University of Yazd, Yazd, Iran

    Mohammad Golshan-Tafti

  2. Department of Medical Genetics, School of Medicine, Shiraz University of Medical Sciences, Shiraz, Iran

    Seyed Alireza Dastgheib

  3. General Practitioner, Babol University of Medical Sciences, Babol, Iran

    Kamran Alijanpour

  4. Neonatal Research Center, Shiraz University of Medical Sciences, Shiraz, Iran

    Reza Bahrami

  5. Mother and Newborn Health Research Center, Shahid Sadoughi University of Medical Sciences, Yazd, Iran

    Mahta Mazaheri & Hossein Neamatzadeh

Authors
  1. Mohammad Golshan-Tafti
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Seyed Alireza Dastgheib
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Kamran Alijanpour
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Reza Bahrami
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Mahta Mazaheri
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Hossein Neamatzadeh
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Study concept and design: MGT, and RB. Data analysis and interpretation: SAD. Drafting of the manuscript: KA. Critical review of the manuscript: MM, HN. Statistical analysis: SAD, HN. Administrative and technical support: RB, MM. Study supervision: MGT, MM. All authors read and approved the final manuscript.

Corresponding authors

Correspondence to Kamran Alijanpour or Reza Bahrami.

Ethics declarations

Ethics approval

This article does not contain any studies with human participants or animals performed by any of the authors.

Consent for publication

Not applicable for this manuscript.

Consent to participate

Not applicable for this manuscript.

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.

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

Golshan-Tafti, M., Dastgheib, S.A., Alijanpour, K. et al. A thorough analysis of data on the correlation between COL9A1 polymorphisms and the susceptibility to congenital talipes equinovarus: a meta-analysis. J Orthop Surg Res 19, 345 (2024). https://doi.org/10.1186/s13018-024-04834-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-024-04834-5

Keywords

Journal of Orthopaedic Surgery and Research

ISSN: 1749-799X