Abstract
BACKGROUND: Clubfoot, presenting as a rigid inward and downward turning of the foot, is one of the most common congenital musculoskeletal anomalies. The aetiology of clubfoot is poorly understood and variants in known clubfoot disease genes account for only a small portion of the heritability. METHODS: Exome sequence data were generated from 1190 non-syndromic clubfoot cases and their family members from multiple ethnicities. Ultra-rare variant burden analysis was performed comparing 857 unrelated clubfoot cases with European ancestry with two independent ethnicity-matched control groups (1043 in-house and 56 885 gnomAD controls). Additional variants in prioritised genes were identified in a larger cohort, including probands with non-European ancestry. Segregation analysis was performed in multiplex families when available. RESULTS: Rare variants in 29 genes were enriched in clubfoot cases, including PITX1 (a known clubfoot disease gene), HOXD12, COL12A1, COL9A3 and LMX1B. In addition, rare variants in posterior HOX genes (HOX9-13) were enriched overall in clubfoot cases. In total, variants in these genes were present in 8.4% (100/1190) of clubfoot cases with both European and non-European ancestry. Among these, 3 are de novo and 22 show variable penetrance, including 4 HOXD12 variants that segregate with clubfoot. CONCLUSION: We report HOXD12 as a novel clubfoot disease gene and demonstrate a phenotypic expansion of known disease genes (myopathy gene COL12A1, Ehlers-Danlos syndrome gene COL9A3 and nail-patella syndrome gene LMX1B) to include isolated clubfoot.