Variants in KIAA0825 underlie autosomal recessive postaxial polydactyly

Identification of truncated variants in GLI family zinc finger 3 (GLI3) associated with polydactyly

BACKGROUND

Polydactyly is a prevalent congenital anomaly with an incidence of 2.14 per 1000 live births in China. GLI family zinc finger 3 (GLI3) is a classical causative gene of polydactyly, and serves as a pivotal transcription factor in the hedgehog signaling pathway, regulating the development of the anterior-posterior axis in limbs.

METHODS

Three pedigrees of polydactyly patients were enrolled from Hunan Province, China. Pathogenic variants were identified by whole-exome sequencing (WES) and Sanger sequencing.

RESULTS

Three variants in GLI3 were identified in three unrelated families, including a novel deletion variant (c.1372del, p.Thr458GlnfsTer44), a novel insertion-deletion (indel) variant (c.1967_1968delinsAA, p.Ser656Ter), and a nonsense variant (c.2374 C > T, p.Arg792Ter). These variants were present exclusively in patients but not in healthy individuals.

CONCLUSIONS

We identified three pathogenic GLI3 variants in polydactyly patients, broadening the genetic spectrum of GLI3 and contributing significantly to genetic counseling and diagnosis for polydactyly.

Variants in EFCAB7 underlie nonsyndromic postaxial polydactyly

Polydactyly is the most common limb malformation that occurs in 1.6-10.6 per one thousand live births, with incidence varying with ancestry. The underlying gene has been identified for many of the ~100 syndromes that include polydactyly. While for the more common form, nonsydromic polydactyly, eleven candidate genes have been reported. We investigated the underlying genetic cause of autosomal recessive nonsyndromic postaxial polydactyly in four consanguineous Pakistani families. Some family members with postaxial polydactyly also present with syndactyly, camptodactyly, or clinodactyly. Analysis of the exome sequence data revealed two novel homozygous frameshift deletions in EFCAB7: [c.830delG;p.(Gly277Valfs*5)]; in three families and [c.1350_1351delGA;p.(Asn451Phefs*2)] in one family. Sanger sequencing confirmed that these variants segregated with postaxial polydactyly, i.e., family members with postaxial polydactyly were found to be homozygous while unaffected members were heterozygous or wild type. EFCAB7 displays expressions in the skeletal muscle and on the cellular level in cilia. IQCE-EFCAB7 and EVC-EVC2 are part of the heterotetramer EvC complex, which is a positive regulator of the Hedgehog (Hh) pathway, that plays a key role in limb formation. Depletion of either EFCAB7 or IQCE inhibits induction of Gli1, a direct Hh target gene. Variants in IQCE and GLI1 have been shown to cause nonsyndromic postaxial polydactyly, while variants in EVC and EVC2 underlie Ellis van Creveld and Weyers syndromes, which include postaxial polydactyly as a phenotype. This is the first report of the involvement of EFCAB7 in human disease etiology.

Identification of GLI1 and KIAA0825 Variants in Two Families with Postaxial Polydactyly

Polydactyly is a rare autosomal dominant or recessive appendicular patterning defect of the hands and feet, phenotypically characterized by the duplication of digits. Postaxial polydactyly (PAP) is the most common form and includes two main types: PAP type A (PAPA) and PAP type B (PAPB). Type A involves a well-established extra digit articulated with the fifth or sixth metacarpal, while type B presents a rudimentary or poorly developed superfluous digit. Pathogenic variants in several genes have been identified in isolated and syndromic forms of polydactyly. The current study presents two Pakistani families with autosomal recessive PAPA with intra- and inter-familial phenotype variability. Whole-exome sequencing and Sanger analysis revealed a novel missense variant in KIAA0825 (c.3572C>T: p.Pro1191Leu) in family A and a known nonsense variant in GLI1 (c.337C>T: p.Arg113*) in family B. In silico studies of mutant KIAA0825 and GLI1 proteins revealed considerable structural and interactional modifications that suggest an abnormal function of the proteins leading to the disease phenotype. The present study broadens the mutational spectrum of KIAA0825 and demonstrates the second case of a previously identified GLI1 variant with variable phenotypes. These findings facilitate genetic counseling in Pakistani families with a polydactyly-related phenotype.

Genetics of congenital anomalies of the hand

Congenital anomalies of the hand are malformations occurring during the development of the human limb, and present as isolated disorders or as a part of a syndrome. During the last years, molecular analysis techniques have offered increasing knowledge about the molecular basis of hand malformations. Disturbances in the signaling pathways during the development of the upper limb result in malformations of the upper extremity. At present, several genes have been identified as responsible for hand anomalies and other have been recognized as suspect genes related to them. Different and new high throughput methods have been introduced for the identification of the gene mutations. In the current editorial, we summarize concisely the current molecular status of isolated hand genetic disorders and the recent progress in molecular genetics, including the genes related to the disorder. This progress improves the knowledge of these disorders and has implications on genetic counselling and prenatal diagnosis.

The molecular genetics of human appendicular skeleton

Disorders that result from de-arrangement of growth, development and/or differentiation of the appendages (limbs and digit) are collectively called as inherited abnormalities of human appendicular skeleton. The bones of appendicular skeleton have central role in locomotion and movement. The different types of appendicular skeletal abnormalities are well described in the report of "Nosology and Classification of Genetic skeletal disorders: 2019 Revision". In the current article, we intend to present the embryology, developmental pathways, disorders and the molecular genetics of the appendicular skeletal malformations. We mainly focused on the polydactyly, syndactyly, brachydactyly, split-hand-foot malformation and clubfoot disorders. To our knowledge, only nine genes of polydactyly, five genes of split-hand-foot malformation, nine genes for syndactyly, eight genes for brachydactyly and only single gene for clubfoot have been identified to be involved in disease pathophysiology. The current molecular genetic data will help life sciences researchers working on the rare skeletal disorders. Moreover, the aim of present systematic review is to gather the published knowledge on molecular genetics of appendicular skeleton, which would help in genetic counseling and molecular diagnosis.

Prenatal Detection of Novel Compound Heterozygous Splice Site Variants of the KIAA0825 Gene in a Fetus with Postaxial Polydactyly Type A

Postaxial polydactyly (PAP) is a common abnormality characterized by extra digits on hands and/or feet. To date, sequence variants in seven genes have been identified in non-syndromic PAP. In the present study, a fetus manifesting non-syndromic postaxial polydactyly type A (PAPA) was found by fetal ultrasonography. To better evaluate fetal prognosis, SNP array analysis and trio whole-exome sequencing (trio-WES) were performed to identify the underlying etiology. Although SNP array analysis revealed no abnormality, trio-WES identified compound heterozygous splice site variants in KIAA0825, c.-1-2A>T and c.2247-2A>G in intron 2 and intron 12, respectively. These two splice site variants were absent in control databases and were predicted to influence splicing by in silico analysis. To confirm the potential pathogenicity of the variants, in vitro splicing assays using minigene and RNA from peripheral leukocytes of the heterozygous parents were conducted. Minigene and RT-PCR assays demonstrated that the c.-1-2A>T variant led to the loss of the initiation codon, and the c.2247-2A>G variant mainly resulted in exon 13 skipping. Prenatal WES and subsequent functional studies are important approaches for defining the genetic etiology of fetuses with PAPA and are also essential for accurate genetic counseling and decision making. Taken together, this study expands the spectrum of KIAA0825 variations in PAPA patients and increases the knowledge of the molecular consequences of KIAA0825 splice site variants.

AutozygosityMapper: Identification of disease-mutations in consanguineous families

With the shift from SNP arrays to high-throughput sequencing, most researchers studying diseases in consanguineous families do not rely on linkage analysis any longer, but simply search for deleterious variants which are homozygous in all patients. AutozygosityMapper allows the fast and convenient identification of disease mutations in patients from consanguineous pedigrees by focussing on homozygous segments shared by all patients. Users can upload multi-sample VCF files, including WGS data, without any pre-processing. Genome-wide runs of homozygosity and the underlying genotypes are presented in graphical interfaces. AutozygosityMapper extends the functions of its predecessor, HomozygosityMapper, to the search for autozygous regions, in which all patients share the same homozygous genotype. We provide export of VCF files containing only the variants found in homozygous regions, this usually reduces the number of variants by two orders of magnitude. These regions can also directly be analysed with our disease mutation identification tool MutationDistiller. The application comes with simple and intuitive graphical interfaces for data upload, analysis, and results. We kept the structure of HomozygosityMapper so that previous users will find it easy to switch. With AutozygosityMapper, we provide a fast web-based way to identify disease mutations in consanguineous families. AutozygosityMapper is freely available at https://www.genecascade.org/AutozygosityMapper/.

The seminal plasma proteins Peptidyl arginine deaminase 2, rRNA adenine N (6)-methyltransferase and KIAA0825 are linked to better motility post thaw in stallions

Seminal plasma plays an important role in sperm physiology. Seminal plasma proteins vehiculated in microvesicles, carry RNAs and proteins with a potential role in early embryo development. Additionally, proteins present in seminal plasma participate in redox regulation and energy metabolism. In view of these facts, we hypothesized that differences in protein composition of the seminal plasma among stallions may help to explain differences in freeze-ability seen among them. Three independent ejaculates from 10 different stallions of varying breeds were frozen using standard protocols in our laboratory. Aliquots of the ejaculate were separated and stored at -80 °C until further proteomic analysis. Semen analysis was performed using computer assisted sperm analysis and flow cytometry. Significant differences in proteome composition of seminal plasma were observed in the group of stallions showing better motility post thaw. 3116 proteins were identified, and of these, 34 were differentially expressed in stallions with better motility post thaw, 4 of them were also differentially expressed in stallions with different percentages of linearly motile sperm post thaw and 1 protein, Midasin, was expressed in stallions showing high circular velocity post thaw. Seminal plasma proteins may play a major role in sperm functionality; being vehiculated through extracellular vesicles and participating in sperm physiology. Bioinformatic analysis identifies discriminant proteins able to predict the outcome of cryopreservation, identifying potential new biomarkers to assess ejaculate quality.

A Novel Homozygous Missense Mutation in the Zinc Finger DNA Binding Domain of GLI1 Causes Recessive Post-Axial Polydactyly

Background: Polydactyly is a prevalent digit abnormality characterized by having extra digits/toes. Mutations in eleven known genes have been associated to cause nonsyndromic polydactyly: GLI3, GLI1, ZRS regulating LMBR1, IQCE, ZNF141, PITX1, MIPOL1, FAM92A, STKLD1, KIAA0825, and DACH1. Method: A single affected family member (IV-4) was subjected to whole-exome sequencing (WES) to identify the causal gene. Bi-directional Sanger sequencing was performed to segregate the identified variant within the family. In silico analysis was performed to investigate the effect of the variant on DNA binding properties. Results: whole-exome sequencing identified a bi-allelic missense variant (c.1010C > T; p. Ser337Leu) in exon nine of GLI1 gene located on chromosome 12q13.3. With the use of Sanger sequencing, the identified variant segregated perfectly with the disease phenotype. Furthermore, in silico analysis of this DNA binding protein revealed that the variant weakened the DNA binding interaction, resulting in indecorous GLI1 function. Conclusion: Herein, we report a novel variant in GLI1 gene, causing autosomal recessive post-axial polydactyly type A (PAPA) type 8. This confirms the critical role of GLI1 in digit development and might help in genotype-phenotype correlation in the future.

Identification of the genetic basis of sporadic polydactyly in China by targeted sequencing

Purpose

Polydactyly is a highly heterogeneous group of skeletal deformities in clinical and genetic background. The variation spectrum in Chinese sporadic polydactyly has not been comprehensively analyzed. To elucidate genetic variation spectrum and genotype-phenotype correlations in Chinese patients with polydactyly, we conducted comprehensive genetic analysis of patients nationwide using targeted sequencing.

Methods

A total of 181 patients diagnosed with polydactylies were recruited. We designed a targeted capture panel for sequencing 721 genes that are associated with the pathogenesis of skeletal dysplasia. We performed rigorous variant- and gene-level filtrations to identify potentially damaging variants, followed by enrichment analysis and gene prioritization.

Results

A total of 568 deleterious variants of 293 genes were identified in 173 of 181 patients with a positive rate of 95.6% by targeted sequencing. For each sample, an average of 3.17 deleterious variants were identified. Especially, 14 pathogenic or likely pathogenic variants were identified in 10 genes in 14 patients out of the 181 patients, providing a positive molecular diagnostic rate of 7.7%.

Conclusion

Targeted sequencing analysis provides a high efficiency approach for the genetic diagnosis of polydactyly. This is the largest next generation sequencing study performed to date in patients with polydactyly and represents the genetic basis of polydactyly typically encountered in genetics clinics.

Role of DZIP1-CBY-FAM92 transition zone complex in the basal body to membrane attachment and ciliary budding

Cilia play important signaling or motile functions in various organisms. In Human, cilia dysfunctions are responsible for a wide range of diseases, called ciliopathies. Cilia assembly is a tightly controlled process, which starts with the conversion of the centriole into a basal body, leading to the formation of the ciliary bud that protrudes inside a ciliary vesicle and/or ultimately at the cell surface. Ciliary bud formation is associated with the assembly of the transition zone (TZ), a complex architecture of proteins of the ciliary base which plays critical functions in gating proteins in and out of the ciliary compartment. Many proteins are involved in the assembly of the TZ, which shows structural and functional variations in different cell types or organisms. In this review, we discuss how a particular complex, composed of members of the DZIP1, CBY and FAM92 families of proteins, is required for the initial stages of cilia assembly leading to ciliary bud formation and how their functional hierarchy contributes to TZ assembly. Moreover, we summarize how evidences in Drosophila reveal functional differences of the DZIP1-CBY-FAM92 complex in the different ciliated tissues of this organism. Whereas it is essential for proper TZ assembly in the two types of ciliated tissues, it is involved in stable anchoring of basal bodies to the plasma membrane in male germ cells. Overall, the DZIP1-CBY-FAM92 complex reveals a molecular assembly pathway required for the initial stages of ciliary bud formation and that is conserved from Drosophila to Human.

A Frameshift Variant in KIAA0825 Causes Postaxial Polydactyly

Postaxial polydactyly (PAP) is characterized by counterproductive 5th digit (pinky finger) duplication on hands and/or feet which often leads to functional complications. To date, at least 11 genes involved in causing various types of nonsyndromic polydactylies have been reported. In the present study, a consanguineous family of Sindhi origin with a segregating nonsyndromic form of PAP in an autosomal recessive manner was clinically and genetically evaluated. Genotyping, using polymorphic microsatellite markers, established linkage in the family on chromosome 5q15 harboring the KIAA0825 gene (MIM 617266). Sequence analysis of the gene revealed a novel frameshift variant leading to a premature stop codon [c.143delG, p.(Cys48Serfs*28)]. This is only the 4th novel variant in the KIAA0825 gene that leads to PAP type A10 (PAPA10) (MIM 618498). Identification of variants in the PAP causative gene will support the diagnosis of patients with limb malformations in the Pakistani population.

Novel heterozygous sequence variant in the GLI1 underlies postaxial polydactyly

Polydactyly is one of the most common congenital abnormal phenotype of autopod, which is characterized by extra supernumerary digit in hands/feet with or without well-developed bony structure within the digits. Preaxial polydactyly (PPD), postaxial polydactyly (PAP), and meso-axial (central) polydactyly are three different isoforms of polydactyly. Genetically, at least 10 genes have been identified causing nonsyndromic polydactyly. In the present study, we have investigated a large family segregating autosomal dominant form of nonsyndromic polydactyly. Whole exome sequencing followed by Sanger sequencing revealed a novel heterozygous missense variant (NM_005269.3; c.1064C>A; p.(Thr355Asn) in the gene GLI1 segregating with the disease phenotype within the family. This study presents first familial case of autosomal dominant form of polydactyly caused by the GLI1 variant.

Identification of a novel biallelic missense variant in the KIAA0825 underlies postaxial polydactyly type A

Postaxial polydactyly (PAP) is characterized by development of extra digits, which mostly segregates in autosomal recessive pattern. The underlying genetic cause of recessive non-syndromic PAP type A has been associated with sequence variants in five different genes (ZNF141, IQCE, GLI1, FAM92A, KIAA0825). The present study was aimed to investigate clinical and genetic causes of PAPA in a consanguineous family of Pakistani origin. Microsatellite-based linkage analysis was used to search for the disease-causing gene. Linkage in the family was established at chromosome 5q15 harbouring a candidate gene KIAA0825. Subsequently, Sanger sequencing revealed a novel homozygous missense variant [c.50T>C; p. (Leu17Ser)] in the gene, which co-segregated with the disease within the family. Protein structural analysis predicted a substantial change in the secondary structure of the mutant protein affecting its function. This is the third disease causing variant identified in the KIAA0825. This has not only expanded spectrum of the mutations in the gene but also further substantiated its role in the limb development in human.

Deleterious Variants in WNT10A, EDAR, and EDA Causing Isolated and Syndromic Tooth Agenesis: A Structural Perspective from Molecular Dynamics Simulations

The dental abnormalities are the typical features of many ectodermal dysplasias along with congenital malformations of nails, skin, hair, and sweat glands. However, several reports of non-syndromic/isolated tooth agenesis have also been found in the literature. The characteristic features of hypohidrotic ectodermal dysplasia (HED) comprise of hypodontia/oligodontia, along with hypohidrosis/anhidrosis, and hypotrichosis. Pathogenic variants in EDA, EDAR, EDARADD, and TRAF6, cause the phenotypic expression of HED. Genetic alterations in EDA and WNT10A cause particularly non-syndromic/isolated oligodontia. In the current project, we recruited 57 patients of 17 genetic pedigrees (A-Q) from different geographic regions of the world, including Pakistan, Egypt, Saudi Arabia, and Syria. The molecular investigation of different syndromic and non-syndromic dental conditions, including hypodontia, oligodontia, generalized odontodysplasia, and dental crowding was carried out by using exome and Sanger sequencing. We have identified a novel missense variant (c.311G>A; p.Arg104His) in WNT10A in three oligodontia patients of family A, two novel sequence variants (c.207delinsTT, p.Gly70Trpfs*25 and c.1300T>G; p.Try434Gly) in EDAR in three patients of family B and four patients of family C, respectively. To better understand the structural and functional consequences of missense variants in WNT10A and EDAR on the stability of the proteins, we have performed extensive molecular dynamic (MD) simulations. We have also identified three previously reported pathogenic variants (c.1076T>C; p.Met359Thr), (c.1133C>T; p.Thr378Met) and (c.594_595insC; Gly201Argfs*39) in EDA in family D (four patients), E (two patients) and F (one patient), correspondingly. Presently, our data explain the genetic cause of 18 syndromic and non-syndromic tooth agenesis patients in six autosomal recessive and X-linked pedigrees (A-F), which expand the mutational spectrum of these unique clinical manifestations.

Novel IQCE variations confirm its role in postaxial polydactyly and cause ciliary defect phenotype in zebrafish

Polydactyly is one of the most frequent inherited defects of the limbs characterized by supernumerary digits and high-genetic heterogeneity. Among the many genes involved, either in isolated or syndromic forms, eight have been implicated in postaxial polydactyly (PAP). Among those, IQCE has been recently identified in a single consanguineous family. Using whole-exome sequencing in patients with uncharacterized ciliopathies, including PAP, we identified three families with biallelic pathogenic variations in IQCE. Interestingly, the c.895_904del (p.Val301Serfs*8) was found in all families without sharing a common haplotype, suggesting a recurrent mechanism. Moreover, in two families, the systemic phenotype could be explained by additional pathogenic variants in known genes (TULP1, ATP6V1B1). RNA expression analysis on patients' fibroblasts confirms that the dysfunction of IQCE leads to the dysregulation of genes associated with the hedgehog-signaling pathway, and zebrafish experiments demonstrate a full spectrum of phenotypes linked to defective cilia: Body curvature, kidney cysts, left-right asymmetry, misdirected cilia in the pronephric duct, and retinal defects. In conclusion, we identified three additional families confirming IQCE as a nonsyndromic PAP gene. Our data emphasize the importance of taking into account the complete set of variations of each individual, as each clinical presentation could finally be explained by multiple genes.

Heterozygous pathogenic variants in GLI1 are a common finding in isolated postaxial polydactyly A/B

Postaxial polydactyly (PAP) is a frequent limb malformation consisting in the duplication of the fifth digit of the hand or foot. Morphologically, this condition is divided into type A and B, with PAP-B corresponding to a more rudimentary extra-digit. Recently, biallelic truncating variants in the transcription factor GLI1 were reported to be associated with a recessive disorder, which in addition to PAP-A, may include syndromic features. Moreover, two heterozygous subjects carrying only one inactive copy of GLI1 were also identified with PAP. Herein, we aimed to determine the level of involvement of GLI1 in isolated PAP, a condition previously established to be autosomal dominantly inherited with incomplete penetrance. We analyzed the coding region of GLI1 in 95 independent probands with nonsyndromic PAP and found 11.57% of these subjects with single heterozygous pathogenic variants in this gene. The detected variants lead to premature termination codons or result in amino acid changes in the DNA-binding domain of GLI1 that diminish its transactivation activity. Family segregation analysis of these variants was consistent with dominant inheritance with incomplete penetrance. We conclude that heterozygous changes in GLI1 underlie a significant proportion of sporadic or familial cases of isolated PAP-A/B.