- Research article
- Open Access
- Open Peer Review
A novel mutation in homeobox DNA binding domain of HOXC13 gene underlies pure hair and nail ectodermal dysplasia (ECTD9) in a Pakistani family
BMC Medical Geneticsvolume 18, Article number: 42 (2017)
Pure hair and nail ectodermal dysplasia (PHNED) is a congenital disorder of hair abnormalities and nail dysplasia. Both autosomal recessive and dominant inheritance fashion of PHNED occurs. In literature, to date, five different forms of PHNED have been reported at molecular level, having three genes known and two loci with no gene yet.
In this study, a four generations consanguineous family of Pakistani origin with autosomal recessive PHNED was investigated. Affected members exhibited PHNED phenotypes with involvement of complete hair loss and nail dysplasia. To screen for mutation in the genes (HOXC13, KRT74, KRT85), its coding exons and exons-intron boundaries were sequenced. The 3D models of normal and mutated HOXC13 were predicted by using homology modeling.
Through investigating the family to known loci, the family was mapped to ectodermal dysplasia 9 (ECTD9) loci with genetic address of 12q13.13. Mutation screening revealed a novel missense mutation (c.929A > C; p.Asn310Thr) in homeobox DNA binding domain of HOXC13 gene in affected members of the family. Due to mutation, loss of hydrogen bonding and difference in potential energy occurs, which may resulting in alteration of protein function.
This is the first mutation reported in homeodomain, while 5th mutation reported in HOXC13 gene causing PHNED.
Pure hair and nail ectodermal dysplasia (PHNED; MIM 602032) is a rare genetic condition, characterized by sparse to complete absence of hairs and nail dystrophy without involvement of any other ectodermal and non-ectodermal manifestations. Hair abnormalities range its mild hair loss to complete absence on the scalp, or may be sparse/absence on rest of the body parts like eyebrows, eyelashes, axillary hair and pubic hair. Nail dystrophy affecting mostly all 20 nails varying in phenotypes including irregular, fragile and spoon shaped.
Inheritance mode of PHNED is both autosomal recessive and dominant [1, 2]. No gene has been reported for autosomal dominant PHNED. Based on OMIM classification, there are five autosomal recessive pure hair and nail ectodermal dysplasia types have been reported. These includes ectodermal dysplasia hair and nail type 4 (ECTD4; MIM 602032) lying on chromosome 12q13.13 with mutations in KRT85 gene , ectodermal dysplasia hair and nail type 5 (ECTD5; MIM 614927) located on chromosome 10q24.32-q25.1 , ectodermal dysplasia hair and nail type 6 (ECTD6; MIM 614928) with chromosomal address 17p12-q21.2 , ectodermal dysplasia hair and nail type 7 (ECTD7; MM 614929) lying on chromosome 12q13.13 with mutation in KRT74 gene  and ectodermal dysplasia hair and nail type 9 (ECTD9; MIM 614931) harboring HOXC13 gene on chromosome 12q13.13 [7–9].
KRT85 and KRT74 genes mutations have been reported only in Pakistani families, while HOXC13 gene mutations have been reported in few other ethnicities including Chinese, Afghani and Syrian along with Pakistani families. Both ECTD5 and ECTD6 linked families are also of Pakistani origin, however the causative genes are yet to be known.
In the present study, a consanguineous Pashto speaking origin family was recruited from Bannu district of Khyber Pakhtunkhwa province of Pakistan, segregating autosomal recessive form of PHNED. Genotyping data and sequence analysis identified a novel missense mutation (c.929A > C) in HOXC13 gene.
In the present investigation, a consanguineous Pakistani family of Pashtun origin with PHNED was recruited (Fig. 1a). To undertake the study, ethical approval was obtained from the Institutional Review Board (IRB) of Kohat University of Science & Technology (KUST), Khyber Pakhtunkhwa, Pakistan. The study was conducted after obtaining written informed consent from all participating subjects including affected individuals, parents and other normal members of the family.
Blood samples were collected from six individuals including three affected members including a male (IV-1, age: 37 years) and two females (IV-3, age: 31 years; IV-4, age: 25 years).
DNA was extracted from the blood samples collected from both affected and unaffected members of the family through phenol-chloroform method. For genotyping, highly polymorphic microsatellite markers were used mapped to the flanking region of the keratin cluster and HOXC cluster genes. After establishment of linkage to the region, based on phenotypes, three genes (HOXC13, KRT74 and KRT85) were prefer to sequence in both affected and unaffected individuals. Big Dye Terminator Cycle Sequencing Kit v3.1 (Applied Biosystems, CA USA) was used for sequencing and the product was analyzed on genetic analyzer ABI 3500 (Applied Biosystems, Inc., Foster City, CA, USA). Primer3 software [http://frodo.wi.mit.edu/primer3/] was used to design primer sequences. Sequence variants in the affected individuals were identified by BIOEDIT sequence alignment (editor version 6.0.7; lbis Biosciences Inc., Carlsbad, CA, USA). PolyPhen-2 [http://genetics.bwh.harvard.edu/pph/] was used to predict the possible effect of the mutation on function of protein. The protein sequences of the pathogenic sequence variant and the evolutionary conservation of amino acid was deduced through multiple sequence alignment tool ClustalW [www.ebi.ac.uk/clustalw/].
The crystal structure Human HOXC13 protein was not available in PDB. So, Homology modeling techniques were utilized to construct the three dimensional structure of Homeobox domain (241aa-322aa) of human HOXC13. The amino acid sequence of Homeobox domain was retrieved from the NCBI database [http://www.ncbi.nlm.nih.gov/], imported to pBlast search against PDB (Protein data bank). We first constructed the homology model of the wild type using the PDB structure 2L7Z as a template, and then built the mutant structure by mutating the selected residue . The model was built using the Molecular Operating Environment (MOE). A series of 10 independent models for protein was built using the Boltzmann weighted randomized procedure combined with specialized logic for the handling of sequence insertions and deletions . Out of 20 models, the model with best MOE packing score was selected for further mutation analysis. The protein structure was visually inspected using PYMOL viewer [http://www.pymol.org].
In the present investigation, three individuals in the family were affected including a brother and two sisters. They were born to healthy first cousin parents. All affected members of the family showed complete hair loss on scalp and rest of the body (Fig. 2a). Affected brother and sisters have dystrophic, irregularly shaped nails at the distal portion and distal onycholysis of the digits on hands while toe nails are dystrophic in all affected members (Fig. 2b-g). Symptoms were present since birth. Other ectodermal structures (teeth and sweat glands) were normal. No other associated abnormalities like mental retardation and skeletal involvement were observed. The parents and other normal siblings show normal hair texture of scalp and rest of the body as well as normal nails.
Based on phenotypes, linkage to the known locus was tested by genotyping highly polymorphic microsatellite markers in six members of the family. All the three affected members showed homozygosity with several markers (D12S1604, D12S325, D12S1724, D12S1632, D12S1298) linked to Keratin and HOXC genes cluster located on chromosome 12p11.1–q21.1. Three different genes (KRT74, KRT85, HOXC13) have been reported in the linked region to be responsible for PHNED [3, 6, 9].
In order to identify the causative gene responsible for this family harboring PHNED, we sequenced all the reported three genes. Sequence analysis revealed a novel homozygous A to C transversion mutation in HOXC13 gene at nucleotide position 929 (c.929A > C). This sequence change results in substitution of Aspargine residue with threonine at amino position 310 (p.Asn310Thr) of HOXC13 protein. Pathogenic sequence variant identified here is present in heterozygous conditions in both healthy siblings and a mother. This mutation is a pathogenic, and not present outside the family because we screened 102 ethnically match control individuals. The variant is absent both in “Exome variant server and Exome Aggregation consortium”. No pathogenic sequence variant was identified in KRT85 and KRT74.
The amide Asn residue in the mutant is replaced by Thr residue. The Asn310 in wild type is involved in the hydrogen bonding with nearby Ile306 while the Thr310 loss the contact (Fig. 3a-c). Hence, the mutant is less stable which is supported by increase in potential energy (−1504.414 kcal/mol) as compared to wild type HOXC13 (−1505.945 kcal/mol). As a consequence of this loss of hydrogen bonding, difference in potential energy and nature of amino acids and different surface area may alter the function of the protein (Fig. 3d-f).
The present study describes PHNED in a consanguineous Pakistani family. Clinical features of affected individuals observed in the present family is mostly similar to those reported earlier in ECTD9 families [7–9]. All affected members have complete hair loss and nail dysplasias. To date, five mutations (Tyr130*, Leu119Trpfs∗20, His68Glnfs*84, Ser135*, 27.6 kb deletion) have been reported in this gene causing PHNED [7–9], of these two mutations have been reported in Pakistani families [http://pakmutation.com/nonsyndromicsearch.aspx]. Based on genotyping and subsequently Sanger sequencing data, we identified a novel mutation (c.929A > C; p.Asn310Thr) in HOXC13 gene.
HOXC13 gene is composed of two exons spanning 2.423 kbs of genomic DNA with chromosomal address of 12q13.3, while encoding a protein of 330 amino acids (GenBank accession number, NM 017410). It contains a cluster of 61 amino acids (amino acids 258–318) forming the DNA binding homeodomain [http://www.uniprot.org/uniprot/Q86TI1]. This homeodomain is involved in the regulation of transcriptional activity of several hair keratin genes and FOXN1 in hair follicles and nails. The mutated aspargine residue at position 310 (identified in the present family), lying in the homeodomain region is highly conserved in other HOXC13 orthologs (Fig. 1e). Analysis of the protein sequence by protein prediction tool PolyPhen2 [http://genetics.bwh.harvard.edu/pph2/] revealed that the substitution of aspargine by threonine (p.Asn310Thr) could potentially have a damaging effect on HOXC13 structure. The mutation was also tested on mutation taster, predicting disease causing.
HOXC gene cluster is located near keratin type II gene cluster. These genes encode highly conserved homeobox family members, which are involved in morphogenesis of several organs. HOXC13 is involved in the development of hair and nails . It also regulates the expression of several other genes (FOXN1, DSG4, CRISP1, FOXQ1) involved in hair follicles and/or nail units.
We have reported a novel mutation in the HOXC13 gene, results in pure hair and nail dysplasia. A homozygous missense mutation reported here lying in the most conserve homeobox DNA binding domain, which confirms the significant role assigned to HOXC13 in ectodermal development. This study further supports the previously reported findings that homozygous mutations in the HOXC13 gene cause ECTD9.
Ectodermal dysplasia type 9
Pure hair nail ectodermal dysplasia
Barbareschi M, Cambiaghi S, Crupi AC, Tadini G. Family with “pure” hair-nail ectodermal dysplasia. Am J Med Genet. 1997;72:91–3.
Calzavara-Pinton P, Carlino A, Benetti A, De Panfilis G. Pili torti and onychodysplasia. Report of a previously undescribed hidrotic ectodermal dysplasia. Dermatologica. 1991;182:184–7.
Naeem M, Wajid M, Lee K, Leal SM, Ahmad W. A mutation in the hair matrix and cuticle keratin KRTHB5 gene causes ectodermal dysplasia of hair and nail type. J Med Genet. 2006;43:274–9.
Rafiq MA, Faiyaz-Ul-Haque M, Ud Din MA, Malik S, Sohail M, Anwar M, Haque S, Paterson AD, Tsui LC, Ahmad W. A novel locus of ectodermal dysplasia maps to chromosome 10q24.32-q25.1. J Invest Dermatol. 2005;124(2):338–42.
Naeem M, Jelani M, Lee K, Ali G, Chishti MS, Wali A, Gul A, John P, Hassan MJ, Leal SM, Ahmad W. Ectodermal dysplasia of hair and nail type: mapping of a novel locus to chromosome 17p12–q21.2. Br J Dermatol. 2006;155:1184–90.
Raykova D, Klar J, Azhar A, Khan TN, Malik NA, Iqbal M, Tariq M, Baig SM, Dahl N. Autosomal recessive transmission of a rare KRT74 variant causes hair and nail ectodermal dysplasia: allelism with dominant woolly hair/hypotrichosis. PLoS One. 2014;9(4):e93607.
Lin Z, Chen Q, Shi L, Giehl KA, Tang Z, Wang H, Zhang J, Yin J, Wu L, Xiao R, Liu X, Dai L, Zhu X, Li R, Betz RC, Zhang X, Yang Y. Loss-of-function mutations in HOXC13 cause pure hair and nail ectodermal dysplasia. Am J Hum Genet. 2012;91:906–11.
Farooq M, Kurban M, Fujimoto A, Fujikawa H, Abbas O, Nemer G, Saliba J, Sleiman R, Tofaili M, Kibbi AG, Ito M, Shimomura Y. A homozygous frameshift mutation in the HOXC13 gene underlies pure hair and nail ectodermal dysplasia in a Syrian family. Hum Mutat. 2013;34:578–81.
Ali RH, Habib R, Ud-Din N, Khan MN, Ansar M, Ahmad W. Novel mutations in the gene HOXC13 underlying pure hair and nail ectodermal dysplasia in consanguineous families. Br J Dermatol. 2013;169(2):478–80.
Zhang Y, Larsen CA, Stadler HS, Ames JB. Structural basis for sequence specific DNA binding and protein dimerization of HOXA13. PLoS One. 2011;6:e23069–e23069.
Fechteler T, Dengler U, Schomburg D. Prediction of protein three-dimensional structures in insertion and deletion regions: a procedure for searching data bases of representative protein fragments using geometric scoring criteria. J Mol Biol. 1995;253(1):114–31.
Shang L, Pruett ND, Awgulewitsch A. Hoxc12 expression pattern in developing and cycling murine hair follicles. Mech Dev. 2002;113(2):207–10.
We are grateful to the research volunteers for their participation and cooperation in this study.
Financially the work presented here was supported through NRPU project (No.4857/NRPU/R&D/HEC 2014) by the Higher Education Commission (HEC), Islamabad, Pakistan.
Availability of data and materials
All data generated or analysed during this study are included in this article. No any other data available for supplementary materials.
AKK design the study, performed molecular analysis and participated in manuscript writing. NM, SAK, MAK and KS located and studied the family, collected blood samples and extracted DNA. AA and AN performed bioinformatics analysis. SK analyzed the data, participated in manuscript preparation and collected funds for the study. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Patients and other family members provide written consent for publishing this study.
Ethical approval & consent to participate
Ethical approval was obtained from the Institutional Review Board (IRB) of Kohat University of Science & Technology (KUST), Khyber Pakhtunkhwa, Pakistan. Patients and other family members provide written consent for participation in the study.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.