Skip to main content

Evaluation of genetic susceptibility of common variants in SOX9 in patients with congenital talipes equinovarus in the Han Chinese population



Congenital talipes equinovarus (CTEV) is a common birth defect that causes severe deformities of one or both feet. Genetics have been proven to play a key role in the risk of CTEV. Our study aimed to evaluate the genetic susceptibility of common variants in the SOX9 gene to CTEV in a Han Chinese population.


In this study, we recruited 2,205 study participants, including 692 CTEV patients and 1513 healthy controls. A total of seven selected single-nucleotide polymorphisms (SNPs) within the SOX9 gene were genotyped, and environmental variables, including maternal smoking and alcoholic drinking habits, were assessed. In addition, bioinformatics analyses were performed to explore the potential biological functions of the associated SNPs.


The SNP rs73354570 was identified to be significantly associated with the risk of CTEV (OR = 1.53, P = 2.11 × 10−5), and the C allele was associated with an increased risk of CTEV. A dose-dependent pattern could be observed in genotypic analyses. The OR for individuals with AC genotypes was 1.37 (95% CI 1.09–1.71), and the OR for individuals with CC homozygotes was 1.47 (95% CI 1.18–1.82). Further analyses identified that rs73354570 is located within a region of multiple binding proteins, including CEBPB and POLR2A, which suggested that this SNP was also part of genetic motifs that are found within several cell types.


Our results provide evidence supporting the important role of the SOX9 gene in the contribution to the risk of CTEV.


Congenital talipes equinovarus (CTEV), or clubfoot, is a common birth defect that occurs in 1 in every 1000 live births [1]. As a severe birth defect, patients with this disorder have difficulty walking because of severe deformities of one or both feet [1]. The etiology of CTEV is still unclear, but genetics have been proven to play a key role in the risk of CTEV [2,3,4,5]. In a family-based study, Gurnett et al. showed that clubfoot segregates in an autosomal dominant fashion with incomplete penetrance [6]. Previous genetic studies have identified several genes that are related to clubfoot, including PITX1, TBX4, and MYH3. The only previous genome-wide association study (GWAS) of clubfoot identified a significant single-nucleotide polymorphism (SNP) between NCOR2 and ZNF664, and there was evidence suggesting that FOXN3, SORCS1, and MMP7/TMEM123 were also associated with clubfoot [7].

SOX9 encodes a transcription factor that can recognize and bind to the sequence of CCTTGAG. Early studies have shown that the products of SOX9 can bind to COL9A1 and regulate its transcription [8]. COL9A1 encodes an α chain of type IX collagen. Early studies have shown that mice with a knockout of the Col9a1 gene do not produce any collagen IX polypeptides, which indicates that COL9A1 is essential for the synthesis of type IX collagen. Recent studies have shown that variations in collagen IX-related genes can cause alterations in bone marrow hyperplasia and might therefore be related to articular cartilage-related diseases [8, 9]. Zhou et al. investigated the association between common SNPs of COL9A1 and clubfoot in a small sample of individuals with Chinese Han ancestry, and a potentially related SNP was identified [10]. Based on these multiple lines of evidence connecting SOX9 and COL9A1, we hypothesize that SOX9 might be involved in the regulation of articular cartilage development. In a recent study, Wang et al. detected significant differences in the mRNA and protein expression levels of SOX9 between idiopathic CTEV muscle samples and controls [11]. Nevertheless, no evidence has been identified to connect the genetic polymorphisms of SOX9 and the risk of CTEV in the human population.

In this study, we aimed to investigate the potential genetic association between genetic polymorphisms of SOX9 and the risk of CTEV based on thousands of study participants with Chinese Han ancestry. In addition to single marker-based association analyses, we also examined the effects of some environmental factors, including maternal smoking and alcoholic drinking habits. Gene-by-environment interactions were also explored.


Study participants

We recruited 692 CTEV patients and 1513 unrelated healthy controls from Honghui Hospital of Xi’an Jiaotong University and Xi’an Children Hospital from January 2012 to April 2017. All participants included in the study were randomly chosen genetically unrelated Han Chinese individuals. All patient diagnoses were confirmed by X-ray examinations. All controls without foot deformities were matched with patients by age and gender. The demographic and clinical characteristics of all participants were obtained from questionnaires and medical records. This study was performed in accordance with the ethical guidelines of the Helsinki Declaration of 1975 (revised in 2008) and was approved through the Ethics Committee of Honghui Hospital of Xi’an Jiaotong University. Written informed consent was obtained from all participants.

SNP selection and genotyping

Tagged SNPs located within gene regions with minor allele frequency (MAF) > 0.01 in SOX9 in 1000 Chinese Han genomes were chosen for genotyping. Algorithm Tagger integrated in Haploview [12] was used for SNP tagging, and the r2 criterion used for tagging was 0.8 for both gene regions. The r2 criterion was used for tagging. A total of 7 candidate SNPs were selected for genotyping in this study.

Genomic DNA was isolated from peripheral blood leukocytes according to the manufacturer’s protocol (Genomic DNA kit, Axygen Scientific, Inc., CA, USA). SNP genotyping was performed using the high-throughput Sequenom MassARRAY platform with iPLEX GOLD chemistry (Sequenom, San Diego, CA, USA) based on the manufacturer’s protocols. The results were processed using Sequenom Typer 4.0 software, and genotype data were generated from the samples [13]. The case and control sample results were blinded for quality control during genotyping processes [14], and 5% of samples were randomly processed with a concordance of 100%.

Statistical methods

Genetic association analyses were performed at both genotypic and allelic levels using Plink v1.9 [15]. A linkage disequilibrium (LD) plot was made using Haploview [12]. In addition to ordinary genetic association analyses, we conducted gene-by-environment interaction analyses by logistic models. For gene-by-environment analyses, two environmental factors, i.e., maternal smoking and maternal alcoholic drinking, were included. The smoking and drinking status were recorded as their habit in general but not status during their pregnancy. Each of these seven selected SNPs was examined by pairing with each of the two environmental factors. Both maternal smoking and alcoholic drinking were coded as 0, 1, and 2, indicating never, occasionally, and often, respectively. Age and gender were included in the logistic models described above. Bonferroni corrections were applied to address multiple comparisons. In addition to statistical methods, we conducted bioinformatics analyses to explore the potential biological functions of our targeted SNPs. RegulomeDB, which annotates SNPs using ENCODE data, was used to investigate the potential regulatory roles of significant SNPs [16]. We also examined the potential eQTL patterns of those significant SNPs using GTEx [17]. Moreover, we investigated the gene-gene network of SOX9 using the STRING database, which is a database of known and predicted protein-protein interactions.


Demographic and clinical characteristics of our study participants

Of the patients included in this study, 8.1% had a family history of CTEV compared with 1.4% of the controls. There was a significant difference in this variable between the study groups (P < 0.001, Table 1). However, no significant differences were found in age, gender, maternal smoking, or maternal drinking between groups (Table 1).

Table 1 Demographic characteristics of CTEV patients and healthy controls

Genetic associations signal

In our study, all seven SNPs selected were in Hardy-Weinberg equilibrium. Basic information, including MAF and the results of Hardy-Weinberg equilibrium tests for these 7 SNPs, is included in Table 2. Weak LD could be identified among these 7 SNPs (Fig. 1). Rs73354570 of SOX9 was the only SNP identified to be significantly associated (odds ratio = 1.53, P = 2.11 × 10−5) with the disease status of CTEV (Table 3). The C allele was associated with an increased risk of CTEV. A dose-dependent pattern could be observed in genotypic analyses. The ORs for individuals with AC genotypes and CC homozygotes were 1.37 (with 95% confidence interval 1.09–1.71) and 1.47 (with 95% confidence interval 1.18–1.82), respectively. The only interaction signal that achieved nominal significance was identified between SNP rs73354570 and maternal alcoholic drinking habit (P = 0.02). However, this signal was no longer significant if we correct for multiple comparisons (threshold of P value is 0.007). The full results are summarized in Table 4.

Table 2 Genetic information of the selected SNPs
Fig. 1

Linkage disequilibrium plot made by the genotyped SNPs. Values of r2 were showed in each cell

Table 3 Genetic association between genetic polymorphisms of SOX9 and CTEV
Table 4 Results of gene by environmental interactions between selected SNPs and smoking and alcoholic drinking

Bioinformatics analyses

We explored the potential functional significance of SNP rs73354570 in RegulomeDB and GTEx. The RegulomeDB score for rs73354570 was 2b. The scoring system of RegulomeDB ranges from 1 to 6, and a lower score often indicates that the SNP has a more significant role in biological functions. A further examination of this SNP identified that rs73354570 is located within a region of multiple binding proteins, including CEBPB and POLR2A. This SNP was also part of genetic motifs found within several cell types. We did not identify significant eQTL signals from GTEx for SNP rs73354570 on SOX9 (Supplemental Table S1). The gene-gene network of SOX9 is shown in Fig. 2. According to the STRING database, the protein product of SOX9 interacts experimentally with several proteins encoded by the PRKG2, COL2A1, RUNX2, AMH, NR5A1, FOXL2, and CTNNB1 genes. In addition, SOX9 was also predicted to be connected with other proteins encoded by genes, including COL10A1, ACAN, and SHH.

Fig. 2

Gene-gene interaction network constructed based on protein-protein interaction (PPI) data. Blue and pink lines mean experimentally determined interactions, and the rest color lines mean the predicted interactions


In this study, we identified an SNP, i.e., rs73354570, that was significantly associated with the disease status of CTEV in a large Chinese Han population-based sample. To the best of our knowledge, our study is the first to report a significant association between SOX9 and CTEV. SOX9 was not reported to be significant in the first and only previous GWAS that focused on clubfoot, which was performed by Zhang et al. [7]. A GWAS with a stringent P value threshold (5 × 10−8) might miss some true positive hits. Our candidate gene-based study design could avoid this limitation by testing fewer markers of preselected genes. To date, the genes SOX9 and SOX9-AS1 (gene encoding anti-sense RNA1 of SOX9) have been identified to be associated with several traits and disorders, including angiotensin-converting enzyme inhibitor intolerance [18], height [19], liver enzyme levels (gamma-glutamyl transferase) [20], lung function [21], nose morphology [22], and thyroid hormone levels [23]. However, no GWAS has established the connection between SOX9 and CTEV-related disorders.

Our bioinformatics analyses have identified that SNP rs73354570 might have significant biological functions in the gene expression of SOX9. This SNP would likely affect the protein binding of its surrounding areas and thus might affect the gene expression of SOX9. However, this information extracted from RegulomeDB could not be replicated by GTEx eQTL data. No significant differences could be found in the gene expression of SOX9 in different genotype groups of SNP rs73354570 in multiple human tissues. Thus, the eQTL data extracted from publicly available databases should be interpreted carefully. Gene expression levels could be very different in target tissues of CTEV patients compared to samples collected in GTEx. Therefore, we believe that it is still too early to reach conclusions on the functional significance of the identified SNPs, although there is a good chance that rs73354570 is a surrogate of another underlying marker or a set of markers. No significant signals have been reported in candidate-based association studies. Wang et al. reported that SOX9 overexpression plays a potential role in idiopathic congenital talipes equinovarus [11]. In the present study, we identified the C allele of rs73354570 to be significantly associated with an increased risk of CTEV. In the future, an eQTL study for rs73354570 on the gene expression of SOX9 based on the target tissue of CTEV patients should be conducted to evaluate the allelic effects on gene expression.

With the rapid development of sequencing technology, numerous susceptibility loci contributing to complex diseases have been reported, such as schizophrenia [24,25,26]. Considering that analyses of only some SNPs are not sufficient to draw conclusions [27,28,29,30,31], we conducted gene-by-environment interaction analyses to investigate potential interactions between our selected SNPs and two environmental factors along with genetic association analyses. However, no significant interaction signals were identified. Nevertheless, we believe that it is not necessary to overinterpret this negative result. Since smoking and alcoholic drinking have been proven to be significantly associated with multiple birth defects [32, 33], it is probable that multiple environmental factors could combine with genetic factors to play a role in the process of CTEV onset. Notably, the most significant interaction signal was obtained from the significant hit of association analysis (rs73354570). Although this signal failed to achieve genome-wide significance, it is still worth investigating this interaction signal further in the future because it might partly explain its single marker-based association signal. A potential limitation of the present study is that smoking and drinking status were recorded as habits in general but not status during pregnancy. Additional studies with appropriately measured environmental factors should be conducted to further investigate the combined effects of both environmental factors and genetic factors on CTEV.

Our study has several limitations. The most important factor is the lack of replication. Future studies are needed to replicate our findings regarding SOX9 in Chinese Han and other populations. Moreover, we did not perform any procedures to adjust for population stratification in the study. However, we tried our best to at least partially control this potential confounding factor because geographic location is a good indicator for genetic matching in the Han Chinese population [34, 35]. Another limitation is that we only included SNPs located within gene regions. However, several important regulatory regions are located at up- or downstream regions that are not within gene regions, and a significant portion of the GWAS panel markers cannot be mapped to any gene regions (but are several kb away from the targeted genes). It is very difficult to claim that seven preselected SNPs could represent most of the genetic information of SOX9. In addition, as a candidate gene-based study focusing on the effects of common SNPs, we did not examine any rare or low-frequency variants. However, a recent study showed that low-frequency and rare variants might play an important role in the onset and development of CTEV [36]. Based on Chinese family samples with CTEV, Zhang et al. identified two pathogenic variations from mediator complex subunit 13L (MED13L) and transforming growth factor-β receptor 2 (TGFBR2). Therefore, a sequencing-based study might provide more information about the genetic etiology of CTEV.

In this study, we identified potential links between genetic polymorphisms of SOX9 and the risk of CTEV. Our study could improve our understanding of the genetic architecture of CTEV and provide a basis for novel intervention plans. Replication studies involving Chinese Han and other populations are still needed in the future.

Availability of data and materials

Please contact the authors for reasonable requests.



Congenital talipes equinovarus


Genome-wide association study


Minor allele frequency


Hardy-Weinberg equilibrium


Linkage disequilibrium


Expression quantitative trait loci


Odds ratios


Single-nucleotide polymorphism


  1. 1.

    Miedzybrodzka Z. Congenital talipes equinovarus (clubfoot): a disorder of the foot but not the hand. J Anat. 2003;202:37–42.

    Article  Google Scholar 

  2. 2.

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

    CAS  PubMed  Google Scholar 

  3. 3.

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

    PubMed  Google Scholar 

  4. 4.

    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  Google Scholar 

  5. 5.

    Dobbs MB, Gurnett CA. Genetics of clubfoot. J Pediatr Orthop B. 2012;21:7–9.

    Article  Google Scholar 

  6. 6.

    Gurnett CA, Alaee F, Kruse LM, Desruisseau DM, Hecht JT, Wise CA, et al. Asymmetric lower-limb malformations in individuals with homeobox PITX1 gene mutation. Am J Hum Genet. 2008;83:616–22.

    CAS  Article  Google Scholar 

  7. 7.

    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.

    CAS  Article  Google Scholar 

  8. 8.

    Liu LY, Jin CL, Jiang L, Lin CK. 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 

  9. 9.

    Brachvogel B, Zaucke F, Dave K, Norris EL, Stermann J, Dayakli M, et al. Comparative proteomic analysis of normal and collagen IX null mouse cartilage reveals altered extracellular matrix composition and novel components of the collagen IX interactome. J Biol Chem. 2013;288:13481–92.

    CAS  Article  Google Scholar 

  10. 10.

    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.

  11. 11.

    Wang Z, Yan N, Liu L, Cao D, Gao M, Lin C, et al. SOX9 overexpression plays a potential role in idiopathic congenital talipes equinovarus. Mol Med Rep. 2013;7:821–5.

    CAS  Article  Google Scholar 

  12. 12.

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

    CAS  Article  Google Scholar 

  13. 13.

    Guan F, Zhang C, Wei S, Zhang H, Gong X, Feng J, et al. Association of PDE4B polymorphisms and schizophrenia in Northwestern Han Chinese. Hum Genet. 2012;131:1047–56.

    CAS  Article  Google Scholar 

  14. 14.

    Guan F, Zhang B, Yan T, Li L, Liu F, Li T, et al. MIR137 gene and target gene CACNA1C of miR-137 contribute to schizophrenia susceptibility in Han Chinese. Schizophr Res. 2014;152:97–104.

    Article  Google Scholar 

  15. 15.

    Chang CC, Chow CC, Tellier LC, Vattikuti S, Purcell SM, Lee JJ. Second-generation PLINK: rising to the challenge of larger and richer datasets. Gigascience. 2015;4:7.

    Article  Google Scholar 

  16. 16.

    Boyle AP, Hong EL, Hariharan M, Cheng Y, Schaub MA, Kasowski M, et al. Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 2012;22:1790–7.

    CAS  Article  Google Scholar 

  17. 17.

    Consortium GT. The Genotype-Tissue Expression (GTEx) project. Nat Genet. 2013;45:580–5.

    Article  Google Scholar 

  18. 18.

    Mahmoudpour SH, Veluchamy A, Siddiqui MK, Asselbergs FW, Souverein PC, de Keyser CE, et al. Meta-analysis of genome-wide association studies on the intolerance of angiotensin-converting enzyme inhibitors. Pharmacogenet Genomics. 2017;27:112–9.

    CAS  Article  Google Scholar 

  19. 19.

    Tachmazidou I, Suveges D, Min JL, Ritchie GRS, Steinberg J, Walter K, et al. Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits. Am J Hum Genet. 2017;100:865–84.

    CAS  Article  Google Scholar 

  20. 20.

    Chambers JC, Zhang W, Sehmi J, Li X, Wass MN, Van der Harst P, et al. Genome-wide association study identifies loci influencing concentrations of liver enzymes in plasma. Nat Genet. 2011;43:1131–8.

    CAS  Article  Google Scholar 

  21. 21.

    Kichaev G, Bhatia G, Loh PR, Gazal S, Burch K, Freund MK, et al. Leveraging Polygenic Functional Enrichment to Improve GWAS Power. Am J Hum Genet. 2019;104:65–75.

    CAS  Article  Google Scholar 

  22. 22.

    Cha S, Lim JE, Park AY, Do JH, 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:481.

    Article  Google Scholar 

  23. 23.

    Porcu E, Medici M, Pistis G, Volpato CB, Wilson SG, Cappola AR, et al. A meta-analysis of thyroid-related traits reveals novel loci and gender-specific differences in the regulation of thyroid function. PLoS Genet. 2013;9:e1003266.

    CAS  Article  Google Scholar 

  24. 24.

    Zhang T, Zhu L, Ni T, Liu D, Chen G, Yan Z, et al. Voltage-gated calcium channel activity and complex related genes and schizophrenia: A systematic investigation based on Han Chinese population. J Psychiatr Res. 2018;106:99–105.

    Article  Google Scholar 

  25. 25.

    Han W, Zhang T, Ni T, Zhu L, Liu D, Chen G, et al. Relationship of common variants in CHRNA5 with early-onset schizophrenia and executive function. Schizophr Res. 2019;206:407–12.

    Article  Google Scholar 

  26. 26.

    Guan F, Ni T, Han W, Lin H, Zhang B, Chen G, et al. Evaluation of the relationships of the WBP1L gene with schizophrenia and the general psychopathology scale based on a case–control study. Am J Med Genet B Neuropsychiatr Genet. 2020;183:164–71.

    CAS  Article  Google Scholar 

  27. 27.

    Zhu L, Li J, Dong N, Guan F, Liu Y, Ma D, et al. mRNA Changes in Nucleus Accumbens Related to Methamphetamine Addiction in Mice. Sci Rep. 2016;6:36993.

    CAS  Article  Google Scholar 

  28. 28.

    Sun H, Luo C, Chen X, Tao L. Assessment of Cognitive Dysfunction in Traumatic Brain Injury Patients: A Review. Forensic Sci Res. 2017;2:174–9.

    Article  Google Scholar 

  29. 29.

    Zhang Z, Gong Q, Feng X, Zhang D, Quan L. Astrocytic clasmatodendrosis in the cerebral cortex of methamphetamine abusers. Forensic Sci Res. 2017;2:139–44.

    Article  Google Scholar 

  30. 30.

    Li J, Zhu L, Guan F, Yan Z, Liu D, Han W, et al. Relationship between schizophrenia and changes in the expression of the long non-coding RNAs Meg3, Miat, Neat1 and Neat2. J Psychiatr Res. 2018;106:22–30.

    Article  Google Scholar 

  31. 31.

    Guan F, Zhang T, Han W, Zhu L, Ni T, Lin H, et al. Relationship of SNAP25 variants with schizophrenia and antipsychotic-induced weight change in large-scale schizophrenia patients. Schizophr Res. 2020;215:250–5.

    Article  Google Scholar 

  32. 32.

    McDonough M. Update on medicines for smoking cessation. Aust Prescr. 2015;38:106–11.

    Article  Google Scholar 

  33. 33.

    Yang J, Qiu H, Qu P, Zhang R, Zeng L, Yan H. Prenatal alcohol exposure and congenital heart defects: a meta-analysis. PLoS One. 2015;10:e0130681.

    Article  Google Scholar 

  34. 34.

    Jia X, Zhang T, Li L, Fu D, Lin H, Chen G, et al. Two-stage additional evidence support association of common variants in the HDAC3 with the increasing risk of schizophrenia susceptibility. Am J Med Genet B Neuropsychiatr Genet. 2016;171:1105–11.

    CAS  Article  Google Scholar 

  35. 35.

    Guan F, Zhang T, Li L, Fu D, Lin H, Chen G, et al. Two-stage replication of previous genome-wide association studies of AS3MT-CNNM2-NT5C2 gene cluster region in a large schizophrenia case-control sample from Han Chinese population. Schizophr Res. 2016;176:125–30.

    Article  Google Scholar 

  36. 36.

    Zhang J, Li S, Ma S, Liu Y, Wang X, Li Y. Whole-exome sequencing study identifies two novel rare variations associated with congenital talipes equinovarus. Mol Med Rep. 2020;21:2597–602.

    CAS  PubMed  Google Scholar 

Download references


We would like to thank all the study participants for their cooperation.


This research was totally supported by Shaanxi Province Natural Science Foundation (No. 2018JQ3040). The funding sources had no role in the design of this study, the collection, analysis and interpretation of data, the writing of the report, or the decision to submit the paper for publication.

Author information




WF conceived and designed the study. JL was involved in the data search and selection of data, analyzed the data, and wrote the manuscript. ZW, DF, and WW contributed to the subject screening and the collection and preparation of control DNA samples. All authors recruited, diagnosed, and gathered patients. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Weilou Feng.

Ethics declarations

Ethics approval and consent to participate

Written informed consent was obtained from all participants prior to their participation. The research protocol was approved by the Ethics Committee of Honghui Hospital of Xi’an Jiaotong University. The ethical approval was consistent with the standards of the Declaration of Helsinki.

Consent for publication

Not applicable.

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: Table S1.

eQTL signals for SNP rs73354570 on gene SOX9.

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 The Creative Commons Public Domain Dedication waiver ( 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

Verify currency and authenticity via CrossMark

Cite this article

Li, J., Wang, Z., Feng, D. et al. Evaluation of genetic susceptibility of common variants in SOX9 in patients with congenital talipes equinovarus in the Han Chinese population. J Orthop Surg Res 15, 276 (2020).

Download citation


  • Congenital talipes equinovarus
  • SOX9 gene
  • Single-nucleotide polymorphisms
  • Genetic association
  • Case-control study