A homozygous HOXD13 missense mutation causes a severe form of synpolydactyly with metacarpal to carpal transformation

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

The pathogenic mechanism of syndactyly type V identified in a Hoxd13Q50R knock-in mice

Syndactyly type V (SDTY5) is an autosomal dominant extremity malformation characterized by fusion of the fourth and fifth metacarpals. In the previous publication, we first identified a heterozygous missense mutation Q50R in homeobox domain (HD) of HOXD13 in a large Chinese family with SDTY5. In order to substantiate the pathogenicity of the variant and elucidate the underlying pathogenic mechanism causing limb malformation, transcription-activator-like effector nucleases (TALEN) was employed to generate a Hoxd13Q50R mutant mouse. The mutant mice exhibited obvious limb malformations including slight brachydactyly and partial syndactyly between digits 2-4 in the heterozygotes, and severe syndactyly, brachydactyly and polydactyly in homozygotes. Focusing on BMP2 and SHH/GREM1/AER-FGF epithelial mesenchymal (e-m) feedback, a crucial signal pathway for limb development, we found the ectopically expressed Shh, Grem1 and Fgf8 and down-regulated Bmp2 in the embryonic limb bud at E10.5 to E12.5. A transcriptome sequencing analysis was conducted on limb buds (LBs) at E11.5, revealing 31 genes that exhibited notable disparities in mRNA level between the Hoxd13Q50R homozygotes and the wild-type. These genes are known to be involved in various processes such as limb development, cell proliferation, migration, and apoptosis. Our findings indicate that the ectopic expression of Shh and Fgf8, in conjunction with the down-regulation of Bmp2, results in a failure of patterning along both the anterior-posterior and proximal-distal axes, as well as a decrease in interdigital programmed cell death (PCD). This cascade ultimately leads to the development of syndactyly and brachydactyly in heterozygous mice, and severe limb malformations in homozygous mice. These findings suggest that abnormal expression of SHH, FGF8, and BMP2 induced by HOXD13Q50R may be responsible for the manifestation of human SDTY5.

Novel Heredity Basis of the Four-Horn Phenotype in Sheep Using Genome-Wide Sequence Data

Horns are an important breeding trait for sheep. However, no widely recognized viewpoint on the regulatory genes and mechanisms of horns is available, and the genetic basis of the four-horn phenotype (FHP) is unclear. This work conducted a genome-wide association study with 100 sheep genomes from multiple breeds to investigate the genetic basis of the FHP. The results revealed three significant associations (corrected as p < 1.64 × 10-8) of the InDels (CHR2: g.133,742,709delA, g.133,743,215insC, and g.133,743,940delT) for FHP in the intergenic sequence (IGS) between the MTX2 and the LOC105609047 of CHR2. Moreover, 14 significant associations (corrected as p < 1.42 × 10-9) of SNPs with the FHP phenotype were identified in CHR2 and CHR16, including five (e.g., CHR16: g.40,351,378G > A and g.40,352,577G > A) located in the intron of the ADAMTS12 gene, eight (e.g., CHR2: g.133,727,513C > T and g.133,732,145T > G) in the IGS between MTX2 and LOC105609047, and only one (CHR2: g.133,930,761A > G) in the IGS between HOXD1 and MTX2. Obvious divergence was also observed in genotype patterns between the FHP and others (two horns and hornless) in the HOXD1 and ADAMTS12 gene regions. An extremely significant linkage also occurred between Loci I and Loci II within 100 individuals (LD = -156.02186, p < 0.00001). In summary, our study indicated that the genomic sequences from CHR2 and CHR16 contributed to the FHP in sheep, specifically the key candidate genes HOXD1 and ADAMTS12. These results improved our understanding of the Mendelian genetic basis of the FHP in sheep.

A novel Hoxd13 mutation causes synpolydactyly and promotes osteoclast differentiation by regulating pSmad5/p65/c-Fos/Rank axis

The mutations of HOXD13 gene have been involved in synpolydactyly (SPD), and the polyalanine extension mutation of Hoxd13 gene could lead to SPD in mice. In this study, a novel missense mutation of Hoxd13 (NM_000523: exon2: c.G917T: p.R306L) was identified in a Chinese family with SPD. The mice carrying the corresponding Hoxd13mutation were generated. The results showed that the homozygous mutation of Hoxd13 also caused SPD, but heterozygous mutation did not affect limbs development, which was different from that of SPD patients. With the increasing generation, the mice with homozygous Hoxd13 mutation presented more severe syndactyly. Western blotting showed that this mutation did not affect the protein expression of Hoxd13, suggesting that this mutation did not result in haploinsufficiency. Further analysis demonstrated that this homozygous Hoxd13mutation promoted osteoclast differentiation and bone loss, and enhanced the mRNA and protein expression of osteoclast-related genes Rank, c-Fos, and p65. Meanwhile, this homozygous Hoxd13 mutation elevated the level of phosphorylated Smad5 (pSmad5). Co-immunoprecipitation verified that this mutation attenuated the interaction between pSmad5 and HOXD13, suggesting that this mutation released more pSmad5. Inhibition of pSmad5 reduced the expression of Rank, c-Fos, and p65 despite in the mutation group. In addition, inhibition of pSmad5 repressed the osteoclast differentiation. ChIP assay confirmed that p65 and c-Fos could bind to the promoter of Rank. These results suggested that this novel Hoxd13 mutation promoted osteoclast differentiation by regulating Smad5/p65/c-Fos/Rank axis, which might provide a new insight into SPD development.

Diagnostic yield of exome sequencing in congenital vertical talus

BACKGROUND

Congenital vertical talus (CVT), also known as "rocker-bottom foot", is a rare foot deformity associated with a dislocation of the talonavicular joint. Although genetic causes of CVT have been described in single isolated and syndromic families, whole-exome sequencing (WES) of large cohorts have not yet been reported.

METHODS

In this study, 62 probands with CVT were evaluated for likely causative single nucleotide variants (SNVs) and copy number variants (CNVs) using WES. Segregation of variants within families was determined by Sanger sequencing.

RESULTS

In this cohort, CVT occurred as an isolated anomaly in 75.8% (47/62) and was familial in 19.3% (12/62) of cases. Analysis of WES data led to the identification of likely causative variants in known disease genes in 30.6% (19/62) of all CVT probands. More than one proband had likely causative SNVs in TSHZ1, GDF5, and LMX1B. Only two probands had likely causative CNVs: a chromosome 12q13.13 deletion of the 5' HOXC gene cluster, and a chromosome 18q22.3q23 deletion involving TSHZ1. Familial CVT was strongly predictive of identifying a molecular diagnosis [75% (9/12) of familial cases compared to 20% (10/50) of non-familial cases (Chi-square test, P-value = 0.0002)]. There was no difference in the solved rate based on isolated or syndromic presentation, unilateral or bilateral affectation, or sex.

CONCLUSIONS

CVT is genetically heterogeneous and more often caused by SNVs than CNVs. There is a high yield of WES in familial CVT cases (∼75%). Additional research is needed to identify the causes of sporadic CVT, which had much lower solved rates.

Recent Advances in Syndactyly: Basis, Current Status and Future Perspectives

A comprehensive summary of recent knowledge in syndactyly (SD) is important for understanding the genetic etiology of SD and disease management. Thus, this review article provides background information on SD, as well as insights into phenotypic and genetic heterogeneity, newly identified gene mutations in various SD types, the role of HOXD13 in limb deformities, and recently introduced modern surgical techniques for SD. This article also proposes a procedure for genetic analysis to obtain a clearer genotype-phenotype correlation for SD in the future. We briefly describe the classification of non-syndromic SD based on variable phenotypes to explain different phenotypic features and mutations in the various genes responsible for the pathogenesis of different types of SD. We describe how different types of mutation in HOXD13 cause various types of SD, and how a mutation in HOXD13 could affect its interaction with other genes, which may be one of the reasons behind the differential phenotypes and incomplete penetrance. Furthermore, we also discuss some recently introduced modern surgical techniques, such as free skin grafting, improved flap techniques, and dermal fat grafting in combination with the Z-method incision, which have been successfully practiced clinically with no post-operative complications.

Incomplete Penetrance and Variable Expressivity: From Clinical Studies to Population Cohorts

The same genetic variant found in different individuals can cause a range of diverse phenotypes, from no discernible clinical phenotype to severe disease, even among related individuals. Such variants can be said to display incomplete penetrance, a binary phenomenon where the genotype either causes the expected clinical phenotype or it does not, or they can be said to display variable expressivity, in which the same genotype can cause a wide range of clinical symptoms across a spectrum. Both incomplete penetrance and variable expressivity are thought to be caused by a range of factors, including common variants, variants in regulatory regions, epigenetics, environmental factors, and lifestyle. Many thousands of genetic variants have been identified as the cause of monogenic disorders, mostly determined through small clinical studies, and thus, the penetrance and expressivity of these variants may be overestimated when compared to their effect on the general population. With the wealth of population cohort data currently available, the penetrance and expressivity of such genetic variants can be investigated across a much wider contingent, potentially helping to reclassify variants that were previously thought to be completely penetrant. Research into the penetrance and expressivity of such genetic variants is important for clinical classification, both for determining causative mechanisms of disease in the affected population and for providing accurate risk information through genetic counseling. A genotype-based definition of the causes of rare diseases incorporating information from population cohorts and clinical studies is critical for our understanding of incomplete penetrance and variable expressivity. This review examines our current knowledge of the penetrance and expressivity of genetic variants in rare disease and across populations, as well as looking into the potential causes of the variation seen, including genetic modifiers, mosaicism, and polygenic factors, among others. We also considered the challenges that come with investigating penetrance and expressivity.

A Novel Missense Variant of HOXD13 Caused Atypical Synpolydactyly by Impairing the Downstream Gene Expression and Literature Review for Genotype-Phenotype Correlations

Synpolydactyly (SPD) is a hereditary congenital limb malformation with distinct syndactyly designated as SPD1, SPD2, and SPD3. SPD1 is caused by mutations of HOXD13, which is a homeobox transcription factor crucial for limb development. More than 143 SPD patients have been reported to carry HOXD13 mutations, but there is a lack of genotype-phenotype correlation. We report a novel missense mutation of c. 925A > T (p.I309F) in an individual with atypical synpolydactyly inherited from her father with mild clinodactyly and three other different alanine insertion mutations in HOXD13 identified by whole exome sequencing (WES) in 12 Chinese SPD families. Unlike polyalanine extension, which tends to form α-helix and causes protein aggregation in the cytoplasm as shown by molecular simulation and immunofluorescence, the c. 925A > T mutation impairs downstream transcription of EPHA7. We compiled literature findings and analyzed genotype-phenotype features in 173 SPD individuals of 53 families, including 12 newly identified families. Among the HOXD13-related individuals, mutations were distributed in three regions: polyalanine, homeobox, and non-homeobox. Polyalanine extension was the most common variant (45%), followed by missense mutations (32%) mostly in the homeobox compared with the loss-of-function (LOF) variants more likely in non-homeobox. Furthermore, a more severe degree and classic SPD were associated with polyalanine mutations although missense variants were associated with brachydactyly and syndactyly in hands and feet and LOF variants with clinodactyly in hands. Our study broadens the HOXD13 mutation spectrum and reveals the profile of three different variants and their severity of SPD, the genotype-phenotype correlation related to the HOXD13 mutation site provides clinical insight, including for genetic counseling.

A novel microdeletion upstream of HOXD13 in a Chinese family with synpolydactyly

Synpolydactyly (SPD) is a digital malformation with the typical clinical phenotype of the webbing of 3/4 fingers and/or 4/5 toes, and combined with polydactyly. In this study, we investigated a Chinese family with SPD and genetic analysis found that all of the affected individuals in the family carry a heterozygous 11,451 bp microdeletion at chr2:176933872-176945322 (GRCh37), which is located upstream of HOXD13 gene, the known disease gene for SPD1. All the affected individuals in the family carry the heterozygous deletion variant, and the variant co-segregated with SPD in the family. Thus, we speculate that the 11,451 bp microdeletion is the disease-causing variant in the family. To date, the microdeletion associating with SPD1 which we identified is the smallest deletion upstream of the HOXD13 gene and not altering the sequence of the HOXD13 gene.

Novel HOXD13 variants in syndactyly type 1b and type 1c, and a new spectrum of TP63-related disorders

Syndactyly is the most common limb defect depicting the bony and/or cutaneous fusion of digits. Syndactyly can be of various types depending on the digits involved in the fusion. To date, eight syndactyly-associated genes have been reported, of which HOXD13 and GJA1 have been explored in a few syndactyly but most of them have unknown underlying genetics. In the present study HOXD13, GJA1 and TP63 genes have been screened by resequencing in 24 unrelated sporadic cases with various syndactyly. The screening revealed two pathogenic HOXD13 variants, NM_000523:c.500 A > G [p.(Y167C)], and NM_000523:c.961 A > C [p.(T321P)] in syndactyly type 1b and type 1c, respectively. This is the first report to identify HOXD13 pathogenic variant in syndactyly type 1b and third report in syndactyly type 1c pathogenesis. Furthermore, this study also reports a TP63 pathogenic variant, NM_003722:c.953 G > A [p.(R318H)] in Ectrodactyly and Cleft lip and palate (ECLP). In conclusion, the current study expands the clinical spectrum of HOXD13 and TP63-related disorders.

Context-dependent HOX transcription factor function in health and disease

During animal development, HOX transcription factors determine the fate of developing tissues to generate diverse organs and appendages. The power of these proteins is striking: mis-expressing a HOX protein causes homeotic transformation of one body part into another. During development, HOX proteins interpret their cellular context through protein interactions, alternative splicing, and post-translational modifications to regulate cell proliferation, cell death, cell migration, cell differentiation, and angiogenesis. Although mutation and/or mis-expression of HOX proteins during development can be lethal, changes in HOX proteins that do not pattern vital organs can result in survivable malformations. In adults, mutation and/or mis-expression of HOX proteins disrupts their gene regulatory networks, deregulating cell behaviors and leading to arthritis and cancer. On the molecular level, HOX proteins are composed of DNA binding homeodomain, and large regions of unstructured, or intrinsically disordered, protein sequence. The primary roles of HOX proteins in arthritis and cancer suggest that mutations associated with these diseases in both the structured and disordered regions of HOX proteins can have substantial functional effects. These insights lead to new questions critical for understanding and manipulating HOX function in physiological and pathological conditions.

A Review of the Phenotype of Synpolydactyly Type 1 in Homozygous Patients: Defining the Relatively Long and Medially Deviated Big Toe with/without Cupping of the Forefoot as a Pathognomonic Feature in the Phenotype

Synpolydactyly type 1 (SPD1, OMIM 186000) is inherited as autosomal dominant and is caused by HOXD13 mutations. The condition is rare and is known for its phenotypic heterogeneity. In the homozygous state, the phenotype is generally more severe and is characterized by three main features: a more severe degree of syndactyly, a more severe degree of brachydactyly, and the frequent loss of the normal tubular shape of the metacarpals/metatarsals. Due to the phenotypic heterogeneity and the phenotypic overlap with other types of syndactyly, no pathognomonic feature has been described for the homozygous phenotype of SPD1. In the current communication, the author reviews the literature on the phenotypes of SPD1 in homozygous patients. The review documents that not all homozygous patients show a severe hand phenotype. The review also defines the “relatively long and medially deviated big toe with/without cupping of the forefoot” as a pathognomonic feature in the phenotype. Illustration of this feature is done through a demonstrative clinical report in a multigeneration family with SPD1 and HOXD13 polyalanine repeat expansion. Finally, the pathogenesis of the clinical features is reviewed.

Identification of novel missense mutations associated with non-syndromic syndactyly in two vietnamese trios by whole exome sequencing

BACKGROUND AND METHODS

Syndactyly is a congenital disorder caused by an irregularity in limb formation during the embryonic development. Many studies have demonstrated the critical effect of genetic factor in controlling the outcome of non-syndromic syndactyly. However the signaling pathway causing this disease has not been fully understood. The aim of this study was to identify the genetic mutations that related to syndactyly type I-c and I-d by exome sequencing.

RESULTS

The exome sequence from two patients revealed two novel heterozygous missense mutations: GLI3: cG1622A pT541M and GJA1: cT274C p.Y92H. Sanger sequencing result confirmed that these mutations were present under heterozygous form in the affected mothers, but not in the unaffected fathers. In-silico analyses by SIFT, Polyphen-2, PredictSNP, PhD-SNP, and PROVEAN did confirm the damaging effect of these mutations in the structure and function of the proteins.

CONCLUSIONS

The result suggested that the two novel mutations may be pathogenic for the disease in these families under the dominant model, provided the initial data for further functional studies to investigate whether those mutations play a disturbing role in the molecular network of syndactyly.

A Review of the Genetics and Pathogenesis of Syndactyly in Humans and Experimental Animals: A 3-Step Pathway of Pathogenesis

Embryology of normal web space creation and the genetics of syndactyly in humans and experimental animals are well described in the literature. In this review, the author offers a 3-step pathway of pathogenesis for syndactyly. The first step is initiated either by the overactivation of the WNT canonical pathway or the suppression of the Bone Morphogenetic Protein (BMP) canonical pathway. This leads to an overexpression of Fibroblast Growth Factor 8 (FGF8). The final step is the suppression of retinoic acid in the interdigital mesenchyme leading to suppression of both apoptosis and extracellular matrix (ECM) degradation, resulting in syndactyly.