An N-terminal G11A mutation in HOXD13 causes synpolydactyly and interferes with Gli3R function during limb pre-patterning

Overlapping binding sites underlie TF genomic occupancy

Sequence-specific DNA binding by transcription factors (TFs) is a crucial step in gene regulation. However, current high-throughput in vitro approaches cannot reliably detect lower affinity TF-DNA interactions, which play key roles in gene regulation. Here, we developed PADIT-seq ( p rotein a ffinity to D NA by in vitro transcription and RNA seq uencing) to assay TF binding preferences to all 10-bp DNA sequences at far greater sensitivity than prior approaches. The expanded catalogs of low affinity DNA binding sites for the human TFs HOXD13 and EGR1 revealed that nucleotides flanking high affinity DNA binding sites create overlapping lower affinity sites that together modulate TF genomic occupancy in vivo . Formation of such extended recognition sequences stems from an inherent property of TF binding sites to interweave each other and expands the genomic sequence space for identifying noncoding variants that directly alter TF binding.

One-Sentence Summary

Overlapping DNA binding sites underlie TF genomic occupancy through their inherent propensity to interweave each other.

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.

Identification of a HOXD13 variant in a Mongolian family with incomplete penetrance syndactyly by exon sequencing

Background

Syndactyly (SD) refers to a deformity caused by the fusion and limb differentiation disorder of soft tissues and/or skeletons to varying extents between adjacent fingers (toes). The main features of this disease are phenotypic heterogeneity and genetic heterogeneity. In this study, we examined four generations of a Chinese Mongolian with different phenotypes of syndactylia and analysed and identified the pathogenic genetic variants of SD by exon sequencing.

Methods

The clinical phenotypes of patients were analysed, and the hands and feet were examined by X-ray. The pedigree was drawn, and the family data were analysed. Peripheral blood was collected from the family members, and genomic DNA was extracted. The candidate genes of SD were identified by exon sequencing, and the mutation sites of the captured candidate genes were amplified by PCR and verified by Sanger sequencing.

Results

The family has congenital syndactyly, which is an autosomal dominant disease. At present, this condition has been passed down for 4 generations and was identified in 9 patients, including 4 males and 5 females. Five patients, I₂, II₄, III₅, III_,7 and III₁₀, had unilateral syndactyly, and four patients, III₁₆, IV₃, IV₆ and IV₇, had bilateral finger syndactyly. All of their toes were unaffected. The proband and the other patients in this family had a c.917G > A (p.R306Q) mutation, which is located at position 917 of the second exon of the HOXD13 gene. This mutation results in a change in the amino acid at position 306, in which arginine is changed to glutamine. This mutation cosegregates in unaffected individuals and affected patients in this family. Moreover, 201 Mongolian genome databases and a thousand human genome databases were referenced to further confirm that the pathogenic genetic variant that causes syndactyly in this family is found in HOXD13.

Conclusion

This study found that the mutation site of the pathogenic gene in this family was HOXD13, c.917G > A (p.R306Q). The phenotype of the family member III₁₂ was normal, but this member was also a carrier of the pathogenic genetic variant. This indicates that the disease of this family has incomplete penetrance characteristics. Our results further enrich the expression profile of the HOXD13 gene.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01360-3.

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 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.

A heterozygous duplication variant of the HOXD13 gene caused synpolydactyly type 1 with variable expressivity in a Chinese family

Background

Synpolydactyly type 1 (SPD1), also known as syndactyly type II, is an autosomal dominant limb deformity generally results in webbing of 3rd and 4th fingers, duplication of 4th or 5th toes. It is most commonly caused by mutation in HOXD13 gene. In this study, a five-generation Chinese family affected with SPD1 disease were collected. We tried to identify the pathogenic variations associated with SPD1 involved in the family.

Methods

We used the whole genome sequencing (WGS) to identify the pathogenic variant in this family which was later confirmed by PCR-Sanger sequencing. The genetic variation were evaluated with the frequencies in the 1000 Genome Project and Exome Aggregation Consortium (ExAC) dataset. The significance of variants were assessed using different mutation predictor softwares like Mutation Taster, PROVEAN and SIFT. The classification of variants was assessed according to American College of Medical Genetics and Genomics (ACMG) guidelines.

Results

Our results showed the mutation of 24-base pair duplication (c.183_206dupAGCGGCGGCTGCGGCGGCGGCGGC) in exon one of HOXD13 in heterozygous form which was predicted to result in eight extra alanine (A) residues in N-terminal domain of HOXD13 protein. The mutation was detected in all affected members of the family.

Conclusion

Based on our mutation analysis of variant c.183_206dupAGCGGCGGCTGCGGCGGCGGCGGC in HOXD13 and its cosegregation in all affected family members, we found this variant as likely pathogenic to this SPD1 family. Our study highlights variable expressivity of HOXD13 mutation. Our results also widen the spectrum of HOXD13 mutation responsible for SPD1.

Identification of gene variants in 130 Han Chinese patients with hypospadias by targeted next-generation sequencing

Background

Hypospadias is a common congenital malformation of male external genitalia, which mainly manifests as an abnormal urethral opening on the ventral side of the penis. The etiology and clinical phenotype of hypospadias is highly heterogeneous, and its clinical diagnosis is challenging. Currently, over 70% of patients have an unknown etiology. Here, we performed a targeted analysis of gene mutations in 130 patients with hypospadias of unknown etiology to find the precise genetic cause.

Methods

We developed a targeted next‐generation sequencing (NGS) panel, encompassing the exon coding regions of 105 genes involved in external genitalia and urogenital tract development and performed sequencing analysis on 130 children with hypospadias of unknown etiology.

Results

In total, 25 patients with hypospadias (19.2%) were found to have 20 mutations among the nine genes involved in external genitalia and urogenital tract development, including 16 reported and four novel mutation sites. Twenty‐two patients (16.9%) had diagnostic variants. Multiple genetic mutations were identified in three of the 25 patients. Hypospadias combined with micropenis was the most common phenotype (68%) in 25 patients.

Conclusions

Higher frequency mutations were identified in SRD5A2 (52%) and AR (24%) in our patient cohort. Middle or posterior hypospadias with micropenis may be significant indicators of genetic variations. Polygenic inheritance may be a rare genetic cause of hypospadias.

Pathogenic gene screening in 91 Chinese patients with short stature of unknown etiology with a targeted next-generation sequencing panel

Background

Dwarfism is a common severe growth disorder, but the etiology is unclear in the majority of cases. Recombinant human growth hormone may be a treatment option, but it has limited efficacy. The currently known laboratory assays do not meet the precision requirements for clinical diagnosis. Here, we have constructed a targeted next-generation sequencing (NGS) panel of selected genes that are suspected to be associated with dwarfism for genetic screening.

Methods

Genetic screening of 91 children with short stature of unknown etiology was performed with the help of the NGS panel. All the coding regions and exon-intron boundaries of 166 genes were included in the panel. To clarify the pathogenicity of these mutations, their clinical data were reviewed and analyzed.

Results

The assay identified p.A72G, p.I282V, and p.P491S variants of the PTPN11 gene and a p.I437T variant of the SOS1 gene in 4 cases with Noonan syndrome. A frameshift mutation (p.D2407fs) of the ACAN gene was identified in a case of idiopathic short stature with moderately advanced bone age. A p.R904C variant of the COL2A1 gene was found in a patient, who was accordingly diagnosed with Stickler syndrome. Severe short stature without limb deformity was associated with a p.G11A variant of HOXD13. In addition, we evaluated evidence that a p.D401N variant of the COMP gene may cause multiple epiphyseal dysplasia.

Conclusions

Our findings suggest that syndromes, particularly Noonan syndrome, may be overlooked due to atypical clinical features. This gene panel has been verified to be effective for the rapid screening of genetic etiologies associated with short stature and for guiding precision medicine-based clinical management.

Exome sequencing identifies a novel nonsense mutation of HOXD13 in a Chinese family with synpolydactyly

Synpolydactyly (SPD) is an autosomal dominant limb malformation with a distinctive combination of syndactyly and polydactyly. SPD is clinically heterogeneous and could be genetically classified into three types. The clinical phenotype of SPD is complicated by its variable expressivity. In the present study, whole exome sequencing (WES) was used to identify the affected gene(s) in a Chinese family with atypical SPD phenotype. Our results showed that a novel heterogenous nonsense mutation (c.556C > T, p.R186X) in HOXD13 was associated with this SPD case. Due to variable expressivity, the diagnosis of a clinical heterogenous disease such as SPD is usually difficult. Our results also suggested that WES is an efficient tool to assist with these diagnoses.

Variable expressivity of the phenotype in two families with brachydactyly type E, craniofacial dysmorphism, short stature and delayed bone age caused by novel heterozygous mutations in the PTHLH gene

Brachydactyly refers to shortening of digits due to hypoplasia or aplasia of bones forming the hands and/or feet. Isolated brachydactyly type E (BDE), which is characterized by shortened metacarpals and/or metatarsals, results in a small proportion of patients from HOXD13 or PTHLH mutations, although in the majority of cases molecular lesion remains unknown. BDE, like other brachydactylies, shows clinical heterogeneity with highly variable intrafamilial and interindividual expressivity. In this study, we investigated two Polish cases (one familial and one sporadic) presenting with BDE and additional symptoms due to novel PTHLH mutations. Apart from BDE, the affected family showed short stature, mild craniofacial dysmorphism and delayed bone age. Sanger sequencing of PTHLH revealed a novel heterozygous frameshift mutation c.258delC(p.N87Tfs*18) in two affected individuals and one relative manifesting mild brachydactyly. The sporadic patient, in addition to BDE, presented with craniofacial dysmorphism, normal stature and bone age, and was demonstrated to carry a de novo heterozygous c.166C>T(p.R56*) mutation. Our paper reports on the two novel truncating PTHLH variants, resulting in variable combination of BDE and other symptoms. Data shown here expand the knowledge on the phenotypic presentation of PTHLH mutations, highlighting significant clinical variability and incomplete penetrance of the PTHLH-related symptoms.

Functional classification and mutation analysis of a synpolydactyly kindred

The aim of the present study was to analyze a congenital syndactyly/polydactyly kindred and propose a new functional classification method of clinical significance. The modes of inheritance and mutational mechanisms were also determined using genetic analyses. Hand and foot anatomy and functions were measured using photographic images, X-ray imaging and grip ability tests. Genetic analysis comprised the genotyping of polymorphic microsatellite markers at known polydactyly-associated loci and the sequencing of the candidate gene. A functional classification system was devised to divide the clinical features into three types, which included mild, moderate or severe deformity. The family was concluded to have syndactyly type II with autosomal dominant inheritance. The microsatellites, D2S2310 and D2S2314, at the 2q31–32 chromosome, which have previously been associated with synpolydactyly type I, were found to be associated with the disorder in the current family. A 27-bp insertion mutation was identified in the affected individuals in the HOXD13 gene at this locus. The insertion added a further nine alanine residues to the polyalanine stretch within the encoded protein. In conclusion, the functional classification method described in the present study may be used to guide surgical approaches to treatment. A family was identified in whom expansion of the polyalanine tract in the HOXD13 gene causes autosomal dominant hereditary synpolydactyly.

Mutations in the homeodomain of HOXD13 cause syndactyly type 1-c in two Chinese families

Background

Syndactyly type 1 (SD1) is an autosomal dominant limb malformation characterized in its classical form by complete or partial webbing between the third and fourth fingers and/or the second and third toes. Its four subtypes (a, b, c, and d) are defined based on variable phenotypes, but the responsible gene is yet to be identified. SD1-a has been mapped to chromosome 3p21.31 and SD1-b to 2q34–q36. SD1-c and SD1-d are very rare and, to our knowledge, no gene loci have been identified.

Methods and Results

In two Chinese families with SD1-c, linkage and haplotype analyses mapped the disease locus to 2q31-2q32. Copy number variation (CNV) analysis, using array-based comparative genomic hybridization (array CGH), excluded the possibility of microdeletion or microduplication. Sequence analyses of related syndactyly genes in this region identified c.917G>A (p.R306Q) in the homeodomain of HOXD13 in family A. Analysis on family B identified the mutation c.916C>G (p.R306G) and therefore confirmed the genetic homogeneity. Luciferase assays indicated that these two mutations affected the transcriptional activation ability of HOXD13. The spectrum of HOXD13 mutations suggested a close genotype-phenotype correlation between the different types of HOXD13-Syndactyly. Overlaps of the various phenotypes were found both among and within families carrying the HOXD13 mutation.

Conclusions

Mutations (p.R306Q and p.R306G) in the homeodomain of HOXD13 cause SD1-c. There are affinities between SD1-c and synpolydactyly. Different limb malformations due to distinct classes of HOXD13 mutations should be considered as a continuum of phenotypes and further classification of syndactyly should be done based on phenotype and genotype.