Novel C16orf57 mutations in patients with Poikiloderma with Neutropenia: bioinformatic analysis of the protein and predicted effects of all reported mutations
Orphanet Journal of Rare Diseases volume 7, Article number: 7 (2012)
Poikiloderma with Neutropenia (PN) is a rare autosomal recessive genodermatosis caused by C16orf57 mutations. To date 17 mutations have been identified in 31 PN patients.
We characterize six PN patients expanding the clinical phenotype of the syndrome and the mutational repertoire of the gene. We detect the two novel C16orf57 mutations, c.232C>T and c.265+2T>G, as well as the already reported c.179delC, c.531delA and c.693+1G>T mutations. cDNA analysis evidences the presence of aberrant transcripts, and bioinformatic prediction of C16orf57 protein structure gauges the mutations effects on the folded protein chain.
Computational analysis of the C16orf57 protein shows two conserved H-X-S/T-X tetrapeptide motifs marking the active site of a two-fold pseudosymmetric structure recalling the 2H phosphoesterase superfamily. Based on this model C16orf57 is likely a 2H-active site enzyme functioning in RNA processing, as a presumptive RNA ligase.
According to bioinformatic prediction, all known C16orf57 mutations, including the novel mutations herein described, impair the protein structure by either removing one or both tetrapeptide motifs or by destroying the symmetry of the native folding.
Finally, we analyse the geographical distribution of the recurrent mutations that depicts clusters featuring a founder effect.
In cohorts of patients clinically affected by genodermatoses with overlapping symptoms, the molecular screening of C16orf57 gene seems the proper way to address the correct diagnosis of PN, enabling the syndrome-specific oncosurveillance.
The bioinformatic prediction of the C16orf57 protein structure denotes a very basic enzymatic function consistent with a housekeeping function. Detection of aberrant transcripts, also in cells from PN patients carrying early truncated mutations, suggests they might be translatable. Tissue-specific sensitivity to the lack of functionally correct protein accounts for the main cutaneous and haematological clinical signs of PN patients.
Biallelic mutations of the C16orf57 (OMIM*613276) gene underlie Poikiloderma with Neutropenia (PN; OMIM#604173), an inherited genodermatosis characterized by early onset poikiloderma, pachyonychia, palmo-plantar hyperkeratosis, skeletal defects and non-cyclic neutropenia. This condition results in recurrent infections during infancy and childhood, primarily of a pulmonary nature, and contributes to the postnatal growth delay in weight and height of the patients.
The syndrome was first described by Clericuzio in Navajo Indians  and subsequently in patients of Caucasian ancestry from different geographic areas [2–5]. Following discovery of the causative gene  and molecular evidence for distinct genetic control between PN and Rothmund-Thomson syndrome (RTS; OMIM#268400) , 31 PN patients have been tested and found to bear 17 different mutations in the responsible C16orf57 gene, 84% of the patients were found in the homozygous state [6, 8–12] and only six were compound heterozygous [6, 11, 13]. Interestingly this cohort includes patients previously diagnosed as affected with Dyskeratosis Congenita (DC; OMIM#224230) and with Rothmund-Thomson syndrome illustrating significant phenotypic overlap among these entities [10, 13]. All the identified mutations, six nonsense, six frameshifts and five splicing mutations including the apparent missense change, c.502A>G , lead to loss-of-function. Despite the limited mutational repertoire a few recurrent mutations have been identified, including c.496delA, common among those of Athabaskan ancestry , c.531delA in Turkish families [10, 12] and c.179delC in patients of North African origin [8, 12], consistent with founder mutations restricted to different geographic areas.
Genotyping further patients could better define the spectrum, the type and the geographical distribution of the C16orf57 sequence changes.
The C16orf57 gene is phylogenetically conserved and ubiquitously expressed, suggesting a housekeeping function .
Whenever cDNA analysis has been performed, mutant alleles have been detected [6, 11] hinting they might be translatable. The function of the unidentified C16orf57 protein remains obscure, so we focused on bioinformatic tools to predict the possible structure and ancestry of C16orf57 and to gauge the functional consequences of the different mutations on the folded protein chain, a start point to address genotype-phenotype correlations in the patients.
In this study we characterize six PN patients and define at the DNA and cDNA level the underlying C16orf57 mutations. The mutations all present in the homozygous state include a novel early truncating mutation, a novel IVS2 splice-site mutation and an IVS6 splice-site mutation already described but in the heterozygous condition .
We also present the computational analysis of the C16orf57 protein chain, showing that two conserved H-X-S/T-X tetrapeptide motifs (where X is a hydrophobic residue) likely mark the active site of a two-fold pseudosymmetric structure related to the 2H phosphoesterase superfamily [14, 15]. Based on this analysis we predict that C16orf57 belongs to this enzyme class and the effect of all C16orf57 mutations reported so far cause disruptions of both protein fold and catalytic site that hold the critical function of the protein.
Six patients, three males and three females, referred to us by clinical geneticists and dermatologists, were enrolled in this study. Patients and their parents provided appropriate informed consent.
DNA extraction and mutation analysis
Genomic DNA was isolated from buccal swabs of patient #26 with NucleoSpin Tissue (Macherey-Nagel, Bethlehem, PA 18020, USA) and from whole peripheral blood of patients #11, #16, #17a, #21, #25 and their relatives using the Wizard Genomic DNA Purification Kit (Promega Corporation, Madison, WI, USA).
About 100 ng of DNA were amplified using GoTaq polymerase (Promega) with previously published primers and conditions . Amplicons were sequenced using Big Dye Terminator v.1.1 Cycle Sequencing Kit according to the manufacture's protocol on the ABI PRISM 3130 sequencer (Applied Biosystems, Foster City, CA, USA). Electropherograms were analyzed with ChromasPro software 1.42 (Technelysium Pty Ltd, Tewantin QLD, Australia) using the wild type sequence of C16orf57 gene [ENSG00000103005] as reference.
RNA extraction and RT-PCR analysis
EBV-transformed lymphoblastoid cell lines (LCLs) were established for patients #21, #17a and their parents. LCLs were cultured in RPMI 1640 medium with 2 mM L-glutamine (EuroClone, Milano, Italy) supplemented with 20% foetal bovine serum (Lonza, Walkersville, MD, USA) and 1% Penicillin, Streptomycin and Ampicillin in 37°C humidified incubator with 5% CO2.
Total RNA was extracted from LCLs using TRI-Reagent RNA Isolation reagent (Sigma-Aldrich, Saint-Louis, MI, USA) and from whole blood of patient #26 using PAXgene Blood RNA Kit (PreAnalitix, Hombrechtikon, Swiss).
500 ng of total RNA was reverse-transcribed into cDNA using High Capacity cDNA Reverse Transcription Kit (Applied Biosystems) with random hexamers. PCR amplification of the C16orf57 gene transcripts was performed for patients and positive controls using GoTaq polymerase (Promega). Specific primers and conditions for C16orf57 transcripts are listed in additional file 1, Table S1. Nucleotide sequences were compared to the major C16orf57 transcript reference sequence [ENST00000219281].
An evolutionary and structural profile of C16orf57 was constructed by iterative PsiBLAST searches of Genbank  followed by PsiPRED secondary structure prediction of the resulting multiple sequence alignment . This alignment was also used by HHrep  to look for internally repeated segments of sequence and predicted structure, and by HHpred to locate potentially related folds from the PDB by sensitive HMM-HMM comparison . Comparative modeling of the C16orf57 three-dimensional structure from fold recognition-derived templates was performed by MODELLER  and I-TASSER . Surface patterns of conservation and variability in the resulting model were derived by ConSurf . Structure manipulation and electrostatic potential surface visualization of C16orf57 and related structures were done with Pymol http://www.pymol.org.
The features of the six patients with a suspected PN clinical diagnosis are listed in Table 1.
The cohort comprises two previously described cases, one Italian girl (#21) with the association of Osteopetrosis with Poikiloderma  and one young adult male (#26) from Turkey diagnosed as Rothmund-Thomson , and four patients of novel description from North Africa (#11), US (#25) and Turkey (#16, #17a), respectively. The affected sister of case #17a is included as #17b in Table 1 although her clinical signs are partially recorded due to premature death.
Pictures of four cases are provided in Figure 1, from left to right, according to age ranging from early infancy (#25, panel A), to childhood (#16, panel B), to teen age (#17a, panel C) to young adulthood (#26, panel D).
The acute phase of the rash is visible in the 5 months old infant (#25) who displayed, firstly on arms and legs and subsequently on buttocks and face (Figure 1 A1, A2), flesh-colored papular changes superimposed on background scaly erythema. These early skin findings which have been described at times as ichthyosiform  later evolve into a more characteristic chronic poikiloderma visible in all the elder patients (#16, Figure 1 B1, B2; #17a, Figure 1 C1; #26, Figure 1 D1).
As outlined in Table 1 all our patients have nail dystrophy, most often pachyonychia, but anonychia was observed in one case (#26, Figure 1 D2, D3). The nail dystrophy progresses from multiple nails with distal onycholysis involving fingers and toes, as shown in the infant (#25, Figure 1 A3, A4) to clear-cut pachyonychia which is present in the two older male patients (#16, Figure 1 B3; #17a, Figure 1 C2, C3).
Supplemental skin findings include thin hair (hypotrichosis in #21) and actinic-hyperkeratosis of the helix (#17a).
Bone alterations (Table 1), detectable on radiograph [23, 24], include osteopenia (#26), severe toe and fingers dysplasia (#26, Figure 1 D2, D3), osteopetrosis (#21), calcification and multiple fractures in the long bones (#11). Delayed skeletal maturation (#21) and small stature (#16, #17a, #21, #26) are also observed.
Laboratory findings highlight severe, persistent and non-cyclic neutropenia, leading to recurrent respiratory infections, which are sometimes resistant to standard antibiotic treatments. Myelodysplasia and splenomegaly developed in patient #17a. With regard to sexual development, two Turkish males (#17a, #26) presented with hypogonadism. Cardiac abnormalities were observed in case #25 and demonstrated both structural (pulmonic stenosis) and functional abnormalities (ECG alterations). The family history of case #17a included a miscarriage and a sibling (#17b) born with a congenital heart defect.
C16orf57 mutations and transcript analysis
The splicing mutation affecting the donor splice site of IVS6, c.693+1G>T, is detected in homozygous state in the US infant (#25, Figure 2G).
The two Turkish males, #16 and #17a, carry the same homozygous deletion c.531delA within exon 5 (Figure 2H). cDNA amplification showed a PCR product of the expected size and its sequencing confirmed the frameshift mutation c.531delA predicting the aberrant protein p.H179MfsX86 (Figure 2 L, O).
The Italian girl #21 carries the novel homozygous transition c.232C>T within exon 2 (Figure 2I). Direct sequencing of the amplified cDNA fragment revealed the c.232C>T change leading to a premature stop codon p.R78X (Figure 2 M, P).
The third Turkish patient #26 shows the novel homozygous c.265+2T>G transversion affecting the donor splice site of IVS2 (Figure 2J). Direct sequencing of the cDNA fragment evidenced 18 intronic nucleotides added to exon 2 following the activation of a new cryptic splice site in IVS2 and predicting a premature stop codon p.Y89WfsX3 (Figure 2 N, Q).
Additional file 2 Table S2 sums up the location within the gene, the type of mutation and the predicted effect on the C16orf57 protein of the five homozygous mutations detected in the six investigated patients.
Recurrence and geographic distribution of C16orf57 mutations
A general overview of all 19 C16orf57 mutations detected so far in 37 molecularly tested PN patients depicts all recurrent mutations and their geographic distribution (Figure 3 A, B); 31 out of 37 patients carry mutations in homozygous state (84%) and only six are compound heterozygous.
The most frequent recurrent c.531delA mutation is recorded in the homozygous state in seven patients, including our patients #16 and #17a from two apparently unrelated Turkish families  which originate from the same rural area. This deletion has been also reported in five patients of three different families featuring a cluster of Caucasian ethnicity  (Figure 3B). Three more PN patients are of Turkish ancestry [3, 25].
The second most frequent mutation c.496delA is identified in five patients, four homozygous and one heterozygous, from the Athabaskan inbred ethnic group .
Our eldest patient of Algerian origin shares the homozygous mutation c.179delC with a described Moroccan kindred depicting a North African origin of this mutation , while the youngest patient of our cohort, from the US, shares the mutation c.693+1G>T with two fraternal twins of Scottish origin  (Figure 3B).
Three other mutations have been described in unrelated patients: the early truncating c.243G>A mutation in two Mediterranean families in homozygous state [9, 10, 26] and in a patient recently reported by a French group in the heterozygous state  (Figure 3B); the nonsense c.541C>T mutation reported in two unrelated homozygous patients [10, 11] and c.489_492del4 reported in the heterozygous state in one Navajo/Caucasian and two Caucasian sibs  (Figure 3B).
Bioinformatic prediction of the structure of C16orf57 protein
The unidentified C16orf57 protein sequence lacks revealing motifs or homologs through database searches and does not disclose any internal domains by SMART analysis . As protein structure is better conserved than sequence, we focused on more sensitive fold recognition and structure prediction tools  as a way to find a fold and related biochemical function for C16orf57. Iterative PsiBLAST searches  were first used to assemble a phylogenetically diverse set of C16orf57 paralogs from Genbank. Subsequently, a PsiPRED-derived secondary structure profile of the resulting alignment  suggested the presence of a α+β-rich domain in the C-terminal 185 residues of the C16orf57 chain. Within this conserved domain, the HHrep web server  detected a faint but significant duplication of sequence and secondary structure patterns, prominently anchored by conserved H-X-T/S-X tetramotifs that correspond to H 120 LS 122 L and H 208 LS 210 L sequences in human C16orf57 (Figure 4A). Accordingly, HMM-HMM comparisons by HHpred  located a string of significant matches between the C16orf57 domain (residues 80-265) and diverse members of the 2H phosphoesterase fold superfamily that are characterized by a pseudo-two-fold symmetric α+β fold with a central, dual histidines (2H) active site .
The implication from these fold recognition results is that C16orf57 is likely a 2H-active site enzyme (Additional file 3 Figure S1) functioning in RNA processing pathways, but the low degree of sequence identity (8-14%) in the HHpred alignments do not unambiguously cluster C16orf57 into one of the four major 2H phosphoesterase families: bacterial 2'-5' RNA ligases, fungal RNA ligases, plant and yeast 1',2'-cyclic nucleotide phosphodiesterases, or mammalian 2',3'-cyclic nucleotide phosphodiesterases (CNPase) [14, 29, 30]. Still, the highest scoring structural matches to C16orf57 are to: a Pyrococcus horikoshii 2'-5' RNA ligase (PDB file 1VJG) , a central domain from human AKAP18 with a 2H active site (PDB file 2VFK)  that binds AMP, lacks CNPase activity and resembles bacterial RNA ligase structures; and the Bacillus subtilis YjcG protein, a putative 2'-5' RNA ligase (PDB file 2D4G) .
To test the structure-based assertion that C16orf57 is a presumptive RNA ligase, the human C16orf57 2H phosphoesterase fold was modelled with the top HHpred match (1VJG) as the template for MODELLER , and also by I-TASSER drawing from multiple templates . The resulting comparative models are quite similar (1.42 Å RMSD over 177 aligned Cα positions), with differences parcelled out to loops and short helices. The characteristic geometry of the 2H active site is preserved in C16orf57, with catalytic H 120 and H 208 residues poised across an active site cleft between the two symmetric α+β lobes (Additional file 3 Figure S1). Drawing from the common catalytic mechanism of 2H phosphoesterase structures [29, 30], these enzymes are involved in the hydrolysis of 2',3'- or 1',2'-cyclic phosphates to 2'- or 1'-phosphates. The histidine residues signature functions either as a nucleophilic attacking group for the cyclic phosphate (C-terminal His), or protonate the leaving oxygen in the reaction (N-terminal His). In bacteria, the 2'-5' RNA ligase acts to ligate half-tRNA molecules (one half with a 2',3'-cyclic phosphate, the other half presenting a 5'-OH terminal). However, similar to the case of the RNA-ligase-like structure of the human AKAP18 2H domain , the true substrates of the C16orf57 enzyme remain unknown and await experimental elucidation.
Predicted effects of all described C16orf57 mutations according to the 3D model
The 2H phosphoesterase fold of human C16orf57 is particularly sensitive to the effects of mutations that cause deletions or truncations in the protein chain because of the nature of the internal sequence repeats that are interdigitated between the two structural repeats or lobes (Figure 4B). These lobes frame the active site groove and position the two catalytic His residues, but since the chain topology swaps β-strands and α-helices between lobes, the linear sequence repeats are individually incapable of forming a well-folded lobe-like domain. Since most of the C16orf57 mutations catalogued in the present study cause premature terminations of the protein chain (Figure 5), the resulting protein fragment is either comprised of (a) part of the N-terminal sequence repeat, (b) exactly the N-terminal repeat or (c) the complete N-terminal repeat plus some part of the C-terminal repeat. In all three cases, the mutant protein is unlikely to fold into a well-structured lobe; in addition, the frequent loss of the C-terminal nucleophilic His (Figure 5) further ensures the destruction of the 2H active site. The two instances of deletions of a protein segment that encompasses β-strand C (in the N-terminal sequence repeat) would also be severely disruptive to proper folding of an active protein (Figure 4B).
Since the discovery of the causative gene, Poikiloderma with Neutropenia syndrome can now be confirmed by molecular diagnostic testing for mutations in the C16orf57 gene and thereby differentiated from the phenotypically similar clinical entities Rothmund-Thomson syndrome and Dyskeratosis Congenita. Based on the molecular test, patients previously diagnosed on clinical grounds as having RTS or DC have been successfully reclassified .
We have identified biallelic mutations of the C16orf57 gene in six patients including #26, previously described as atypical Rothmund-Thomson  and #21, reported as affected with Poikiloderma associated with Osteopetrosis . All the patients herein molecularly confirmed as PN display the consensus PN clinical signs, i.e. early-onset poikiloderma involving the extremities and extending to the face, pachyonychia, palmo-plantar hyperkeratosis and non-cyclic neutropenia (Table 1). It is worth mentioning that our patient #17a shows poikiloderma also on the ear helix, as previously reported in other patients [3, 9].
Craniofacial dysmorphisms have been reported only in a few PN patients [5, 32]; facial features such as a saddle nose and midfacial hypoplasia are displayed by four out of our six patients, suggesting they are quite common and hence should be recorded during PN clinical evaluation. Dental defects are common in our cases, in agreement with the literature [4, 13, 26, 33]. Hypogonadism was identified in male patients #17a and #26 and previously reported [3, 34]. In contrast to RECQL4-positive RTS, where bony changes constitute a major diagnostic sign , skeletal involvement has only seldom been described in PN patients [3, 13, 26]. Indeed, apart from our patients #17a and the infant #25, all herein described PN patients display overt skeletal signs including zygodactyly between the second and third digit (#26), multiple bone fractures (#11), intermediate osteopetrosis (#21) or X-ray detectable skeletal findings (#16) (Table 1). A relationship may be envisaged between zygodactyly and the swan neck hand hyperflexibility noticed in a few described patients [3, 5, 33]. More generally delayed bone maturation has been recorded  in the girl subsequently confirmed to carry two distinct C16orf57 mutations , and diffuse osteosclerosis has been underlined by Porter  in the patient found to harbour C16orf57 mutations . This patient is one of three cases reported as RTS  with defects in the bone marrow, as reported for other clinically diagnosed RTS cases [32–34]. We believe these cases, previously classified as affected by RTS, more correctly represent cases of PN as confirmed by molecular testing of some of them [6, 10]. All our patients display non-cyclic neutropenia, the hallmark of PN, and this is a fundamental sign differentiating PN from both RTS or DC, which should be searched also in suspected cases where it remains silent . Neutropenia leads to recurrent infections, which are widely documented in the paediatric histories of PN cases and confers upon them a high risk to develop myelodysplasia from the second decade of life, as attested in at least ten C16orf57- positive patients [5, 10] and our patient #17a. When considering bone marrow hypocellularity, increased myeloid precursors and delayed neutrophils maturation, a higher percentage of C16orf57-positive cases displays this feature [4, 5, 10, 11] highlighting the sensitivity of the myeloid lineage to C16orf57 mutations. Evolution to acute myeloid leukemia has been reported in a few PN patients [10, 26] and in other untested cases [25, 33–35]. This evidence features PN as a cancer predisposing syndrome affecting the myeloid compartment and connects PN to RTS, a syndrome with an increased risk for osteosarcoma and skin cancer , and DC which predisposes to a wide variety of haematological and solid tumors . Molecular characterization is thus compulsory for assigning appropriate oncological surveillance.
As depicted in Figure 3 the recurrence of the c.531delA, c.496delA and c.179delC mutations delineates three clusters according to the geographic origin of the patients. A common ancestor is the likely hypothesis to explain the recurrence of these mutations restricted to specific ethnic groups, considering the very low frequency of the PN syndrome and the prevalence of patients with homozygous mutations (31 out of 37). Genetic analysis and reconstruction of ancient genetic links through haplotype segregation analysis could confirm this assumption.
Expansion of the C16orf57 mutational repertoire and validation of the observed geographical distribution may allow assessment of clinical variability of PN phenotype in distinct founder mutation cohorts, as it has been described for another autosomal recessive developmental disorder .
Lack of information of the function of C16orf57 protein makes it difficult to establish a link between mutations and the onset and evolution of syndromic presentation.
Bioinformatic analysis can provide a preliminary tool to predict the severity of specific mutations.
The structure-based inference that human C16orf57 is a member of the 2H phosphoesterase superfamily, despite very little sequence identity, rests on the recognition of a common protein fold by sensitive algorithms that weave together evolutionary information, in the form of sequence patterns that are conserved across C16orf57 orthologs, and accurate secondary structure predictions to comprehensively scan structural databases. The most prominent pattern in C16orf57 is a two-fold repeated sequence and structural segment with signature H-X-T/S-X tetramotifs; these conserved features are precisely mirrored in the fold-recognition-derived matches with bilobal 2H phosphoesterase folds that rely on the symmetrically poised His residues for catalytic activity. The 2H phosphoesterase superfamily is a diverse grouping of enzymes with a common core architecture, a basic hydrolytic focus for cyclic phosphates, and some functional variability in both substrates and reactions [14, 15]. This degree of functional diversification is observed in other structurally-assembled superfamilies of enzymes . Of the various biological tasks performed by 2H active site enzymes, we argue (by closer structural resemblance) that C16orf57 is perhaps an RNA ligase though the actual targets or substrates of this activity remain unknown.
All the mutations reported in our patients are in the homozygous state, which facilitates transcript analysis and understanding the effect of the mutations.
Indeed with the exception of the Algerian patient #11, transcripts have been tested, detected and sequenced in all other patients (Figure 2). With regards to the c.693+1G>T mutation in patient #25, transcript analysis has been reported  suggesting exon 6 skipping according to the size of the aberrant band; indeed a misspliced 84 base shorter transcript was found associated with the c.683_693+1del12 which affects the same IVS6 donor splice site . The general emerging feature is that the aberrant transcripts are detectable even for early truncating mutations, such as those of patients #21 and #26, pointing out they are relatively stable and translatable.
It is known that Dyskeratosis Congenita, Rothmund-Thomson and Poikiloderma with Neutropenia have many overlapping features. Starting from the clinical presentation it is very difficult to assess the correct diagnosis, and in fact most of the molecularly confirmed PN patients have a long history of wrong diagnosis, as DC or RTS. The availability of C16orf57 molecular testing allows the correct diagnosis which is compulsory for retargeting syndrome-specific oncosurveillance.
As the present study shows, all the 19 C16orf57 mutations linked to disease involve the predicted enzymatic domain of C16orf57 protein. By destroying the native fold, all mutations should cause a drastic loss of enzymatic activity. Future studies of protein presence and activity in PN patients could confirm whether and how the aberrant transcripts, which have been always detected whenever assayed, may be translated.
The onset of the poikiloderma, nail dystrophy and teeth malformations at early infancy reflects the perturbed morphogenesis of skin and cutaneous annexes, while neutropenia results from impaired homeostasis of the highly C16orf57 expressing myeloid cells [6, 10]. The life long risk of myelodysplastic syndrome features the increased tendency to apoptosis and leukemic transformation of C16orf57-defective myeloid progenitor cells .
However, whether different mutations impact differently on the clinical phenotype and on the risk of myelodysplastic syndrome awaits further clinical and molecular characterization of PN patients, along with delineation of patients' clinical expressivity in distinct geographical areas and dissection of the biological function of the C16orf57 protein.
Poikiloderma with Neutropenia
EBV-transformed lymphoblastoid cell lines
Clericuzio C, Hoyme HE, Asse JM: Immune deficient poikiloderma: a new genodermatosis. AJHG. 1991, 49 (Suppl): A661.
Wang LL, Gannavarapu A, Clericuzio CL, Erickson RP, Irvine AD, Plon SE: Absence of RECQL4 mutations in Poikiloderma with Neutropenia in Navajo and Non-Navajo patients. A J Med Genet A. 2003, 118A: 299-301. 10.1002/ajmg.a.10057.
Van Hove JLK, Jaeken J, Proesmans M, De Boeck K, Minner K, Matthijs G, Verbeken E, Demunter A, Boogaerts M: Clericuzio type Poikiloderma With Neutropenia is distinct from Rothmund-Thomson syndrome. A J Med Genet A. 2005, 132A: 152-158. 10.1002/ajmg.a.30430.
Mostefai R, Morice-Picard F, Boralevi F, Sautarel M, Lacombe D, Stasia MJ, McGrath J, Taïeb A: Poikiloderma with neutropenia, Clericuzio type, in a family from Morocco. Am J Med Genet. 2008, 146A: 2762-2769. 10.1002/ajmg.a.32524.
Concolino D, Roversi G, Muzzi GL, Sestito S, Colombo EA, Volpi L, Larizza L, Strisciuglio P: Clericuzio type Poikiloderma with Neutropenia Syndrome in three sibs with mutations in the C16orf57 gene: delineation of the phenotype. Am J Med Genet A. 2010, 152A: 2588-2594. 10.1002/ajmg.a.33600.
Volpi L, Roversi G, Colombo EA, Leijsten N, Concolino D, Calabria A, Mencarelli MA, Fimiani M, Macciardi F, Pfundt R, Schoenmakers EFPM, Larizza L: Targeted next-generation sequencing appoints C16orf57 as Clericuzio-type Poikiloderma with Neutropenia gene. AJHG. 2010, 86: 72-76. 10.1016/j.ajhg.2009.11.014.
Larizza L, Roversi G, Volpi L: Rothmund-Thomson syndrome. Orphanet J Rare Dis. 2010, 5: 2-10.1186/1750-1172-5-2.
Tanaka A, Morice-Picard F, Lacombe D, Nagy N, Hide M, Taïeb A, McGrath J: Identification of a homozygous deletion mutation in C16orf57 in a family with Clericuzio-type poikiloderma with neutropenia. Am J Med Genet A. 2010, 152A: 1347-1348.
Arnold AW, Itin PH, Pigors M, Kohlhase J, Bruckner-Tuderman L, Has C: Poikiloderma with neutropenia: a novel C16orf57 mutation and clinical diagnostic criteria. British J Dermatology. 2010, 163: 866-869. 10.1111/j.1365-2133.2010.09929.x.
Walne AJ, Vulliamy T, Beswick R, Kirwan M, Dokal I: Mutations in C16orf57 and normal-length telomeres unify a subset of patients with dyskeratosis congenita, poikiloderma with neutropenia and Rothmund-Thomson syndrome. Hum Mol Genet. 2010, 19: 4453-4461. 10.1093/hmg/ddq371.
Clericuzio C, Harutyunyan K, Jin W, Erickson RP, Irvine AD, McLean WH, Wen Y, Bagatell R, Griffin TA, Shwayder TA, Plon SE, Wang LL: Identification of a novel C16orf57 mutation in Athabaskan patients with Poikiloderma with Neutropenia. Am J Med Genet A. 2011, 155A: 337-342.
Colombo EA, Roversi G, Elcioglu N, Fremont G, Volpi L, Larizza L: Poikiloderma with neutropenia: Delineation of novel mutations and geographic clusters [abstract]. 2010, ASHG meeting 2380/T. http://www.ashg.org/cgi-bin/2010/showdetail.pl?absno=20843
Piard J, Holder-Espinasse M, Aral B, Gigot N, Rio M, Tardieu M, Puzenat E, Goldenberg A, Toutain A, Franques J, MacDermot K, Bessis D, Boute O, Callier P, Gueneau L, Huet F, Vabres P, Catteau B, Faivre L, Thauvin-Robinet C: Systematic search for neutropenia should be part of the first screening in patients with poikiloderma. Eur J Med Genet.
Mazumder R, Iyer LM, Vasudevan S, Aravind L: Detection of novel members, structure-function analysis and evolutionary classification of the 2H phosphoesterase superfamily. Nucleic Acids Res. 2002, 30: 5229-5243. 10.1093/nar/gkf645.
Sakamoto Y, Tanaka N, Ichimiya T, Kurihara T, Nakamura KT: Crystal structure of the catalytic fragment of human brain 2',3'-cyclic-nucleotide 3'-phosphodiesterase. J Mol Biol. 2005, 346: 789-800. 10.1016/j.jmb.2004.12.024.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Research. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics. 2000, 16: 404-405. 10.1093/bioinformatics/16.4.404.
Söding J, Remmert M, Biegert A: HHrep: de novo protein repeat detection and the origin of TIM barrels. Nucleic Acids Res. 2006, 34: W137-142. 10.1093/nar/gkl130.
Söding J: Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005, 21: 951-960. 10.1093/bioinformatics/bti125.
Sali A, Blundell TL: Comparative protein modelling by satisfaction of spatial restraints. J Mol Biol. 1993, 234: 779-815. 10.1006/jmbi.1993.1626.
Roy A, Kucukural A, Zhang Y: I-TASSER: a unified platform for automated protein structure and function prediction. Nat Protoc. 2010, 5: 725-738. 10.1038/nprot.2010.5.
Ashkenazy H, Erez E, Martz E, Pupko T, Ben-Tal N: ConSurf2010: calculating evolutionary conservation in sequence and structure of proteins and nucleic acids. Nucleic Acids Res. 2010, 38: W529-533. 10.1093/nar/gkq399.
Migliaccio S, Luciani M, Taranta A, De Rossi G, Minisola S, El Hachem M, Bosman C, De Felice L, Boldrini R, Corsi A, Bianco P, Teti A: Association of intermediate osteopetrosis with poikiloderma. J Bone Miner Res. 1999, 14: 834-836. 10.1359/jbmr.1918.104.22.1684.
Altunay I, Fisek N, Gokdemir G, Sakız D, Cetincelik U: Therapy-resistant leg ulcer in a patient with Rothmund-Thomson syndrome. Int Wound J. 2010, 7: 531-535. 10.1111/j.1742-481X.2010.00734.x.
Ilhan I, Arikan U, Buyukpamukçu M: Myelodysplastic syndrome and RTS. Pediatr Hematol Oncol. 1996, 13: 197-10.3109/08880019609030815.
Porter WM, Hardman CM, Abdalla SH, Powles AV: Haematological disease in siblings with Rothmund-Thomson syndrome. Clin Exp Dermatol. 1999, 24: 452-454. 10.1046/j.1365-2230.1999.00530.x.
Letunic I, Doerks T, Bork P: SMART 6: recent updates and new developments. Nucleic Acids Res. 2009, 37: D229-232. 10.1093/nar/gkn808.
Zhang Y: Protein structure prediction: when is it useful?. Current Opinion in Structural Biology. 2009, 19: 145-155. 10.1016/j.sbi.2009.02.005.
Gao YG, Yao M, Okada A, Tanaka I: The structure of Pyrococcus horikoshii 2'-5' RNA ligase at 1.94 Å resolution reveals a possible open form with a wider active-site cleft. Acta Cryst. 2006, F62: 1196-1200.
Gold MG, Smith FD, Scott JD, Barford D: AKAP18 Contains a phosphoesterase domain that binds AMP. Journal of Molecular Biology. 2008, 375: 1329-1343. 10.1016/j.jmb.2007.11.037.
Li D, Chan C, Liang YH, Zheng X, Li L, Su XD: Preparation, crystallization and preliminary X-ray analysis of YjcG protein from Bacillus subtilis. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2005, 61: 496-498. 10.1107/S1744309105011012.
Pianigiani E, De Aloe G, Andreassi A, Rubegni P, Fimiani M: Rothmund-Thomson syndrome (Thomson-type) and myelodysplasia. Pediatr Dermatol. 2001, 18: 422-425. 10.1046/j.1525-1470.2001.01971.x.
Narayan S, Fleming C, Trainer AH, Craig JA: Rothmund-Thomson syndrome with myelodysplasia. Pediatr Dermatol. 2001, 18: 210-212. 10.1046/j.1525-1470.2001.018003210.x.
Rizzari C, Bacchiocchi D, Rovelli A, Biondi A, Cantu'-Rajnoldi A, Uderzo C, Masera G: Myelodysplastic syndrome in a child with Rothmund-Thomson syndrome: a case report. J Pediatr Hematol Oncol. 1996, 18: 96-97. 10.1097/00043426-199602000-00020.
Knoell KA, Sidhu-Malik NK, Malik RK: Aplastic anemia in a patient with Rothmund-Thomson syndrome. J Pediatr Hematol Oncol. 1999, 21: 444-446. 10.1097/00043426-199909000-00021.
Walne AJ, Dokal I: Advances in the understanding of dyskeratosis congenita. Br J Haematol. 2009, 145: 164-172. 10.1111/j.1365-2141.2009.07598.x.
Douzgou S, Petersen MB: Clinical variability of genetic isolates of Cohen syndrome. Clin Genet. 2011, 79: 501-506. 10.1111/j.1399-0004.2011.01669.x.
Todd AE, Orengo CA, Thornton JM: Evolution of function in protein superfamilies, from a structural perspective. J Mol Biol. 2001, 307: 1113-1143. 10.1006/jmbi.2001.4513.
Bejar R, Levine R, Ebert BL: Unraveling the molecular pathophysiology of myelodysplastic syndromes. J Clin Oncol. 2011, 29: 504-515. 10.1200/JCO.2010.31.1175.
We thank the patients and their relatives for intensive cooperation in the study; Galliera Genetic Bank-Network of Telethon Genetic Biobank-project GTB07001- for providing us lymphoblastoid cell lines from the affected siblings of the PN families and DNA of patient #21; we thank Dr. Guillemette Frémont (Departement of Dermatology, Hôpital Saint-Louis, Paris, France) for referring to us patient #11.
As regards all patients, we followed the guidelines of the ethical committee of the Università degli Studi di Milano http://www.unimi.it/cataloghi/comitato_etico/CE_Rec_4_2006_HBMs.pdf.
This work was supported by Associazione Italiana per la Ricerca sul Cancro (grant 2008-2009/4217 to L.L.), CARIPLO N.O.B.E.L. (project 2007-2009 to L.L), Nando Peretti Foundation (grant 2007-2011 to L.V.) and the "Dote di Ricerca": FSE, Regione Lombardia granted to Elisa Adele Colombo, Dipartimento di Medicina Chirurgia e Odontoiatria, Università degli Studi di Milano for the research project "Functional characterization of the new C16orf57 gene".
L.V. and L.L.: equal contribution.
The authors declare that they have no competing interests.
EAC performed molecular analysis and interpreted the predicted effects of mutations, drafted the manuscript and approved the final version. JFB performed and interpreted computational analysis, drafted and revised the article, approved the final manuscript. GN performed and interpreted molecular analysis, mined the literature and approved the final manuscript. CG supported the molecular work and software access. NE and DY contributed with clinical data and biological samples of patients#16 and #17a, provided helpful feedback and approved the final manuscript. NA and UC contributed with clinical data and sample of patient #26 and approved the final manuscript. AT, ADF and ML contributed with clinical data and sample of patient #21 and approved the final manuscript. SKS and ACY contributed with clinical data and sample of patient #25, revised and approved the manuscript. LV contributed to recruit patients, designed the geographical analysis, drafted, revised and approved the manuscript. LL designed and coordinated the study; drafted, revised and approved the manuscript.
Electronic supplementary material
Additional file 1: Primer sequences, amplicons size, annealing temperatures (Ta) for C16orf57 cDNA analysis (ENST00000219281 Ensembl database) and analysed patients. Table S1 provides technical information on PCR conditions used to analyze cDNA of PN patients #17a, #21 and #26. (PDF 9 KB)
Additional file 2: C16orf57 homozygous mutations in the set of investigated PN patients. Table S2 provides a summary of identified mutations of PN patients including location within the gene, mutation type and predicted effect on the protein. All mutations are in the homozygous state. (PDF 9 KB)
Additional file 3: Electrostatic potential surface of human C16orf57 predicted protein. The solvent-accessible surface of the I-TASSER-derived  structural model of human C16orf57 was displayed in Pymol http://www.pymol.org and coloured according to the electrostatic potential (ESP) ranging from blue (positively charged or basic) to red (negatively charged or acidic). In a top view that looks directly down at the active site groove, the corresponding ESP surface is quite negatively charged, which is similar in nature to 2H phosphoesterase structures of the RNA ligase class that interact with positively charged substrates [29, 30]. (PDF 114 KB)
About this article
Cite this article
Colombo, E.A., Bazan, J.F., Negri, G. et al. Novel C16orf57 mutations in patients with Poikiloderma with Neutropenia: bioinformatic analysis of the protein and predicted effects of all reported mutations. Orphanet J Rare Dis 7, 7 (2012). https://doi.org/10.1186/1750-1172-7-7