PCSK5 mutation in a patient with the VACTERL association
BMC Research Notes volume 8, Article number: 228 (2015)
The VACTERL association is a typically sporadic, non-random collection of congenital anomalies that includes vertebral defects, anal atresia, cardiac defects, tracheoesophageal fistula with esophageal atresia, renal anomalies, and limb abnormalities. Although several chromosomal aberrations and gene muta tions have been reported as disease-causative, these findings have been sparsely replicated to date.
In the present study, whole exome sequencing of a case with the VACTERL association uncovered a novel frameshift mutation in the PCSK5 gene, which has been reported as one of the causative genes for the VACTERL association. Although this mutation appears potentially pathogenic in its functional aspects, it was also carried by the healthy father. Furthermore, a database survey revealed several other deleterious variants in the PCSK5 gene in the general population.
Further studies are necessary to clarify the etiological role of the PCSK5 mutation in the VACTERL association.
The VACTERL association (OMIM #192350) is a congenital disease characterized by the presence of at least three of the conditions of vertebral defects, anal atresia, cardiac defects, tracheoesophageal fistula, renal anomalies, and limb abnormalities. The VACTERL association occurs sporadically with a heterogeneous background, and its incidence is estimated at approximately 1 in 10,000 to 1 in 40,000 live-born infants . The etiology of the VACTERL association is not still fully understood, although HOXD13, ZIC3, PTEN, FANCB, FOXF1, and TRAP1 mutations have been reported as possible contributors to the VACTERL association or the VACTERL-like association as described previously [2–9]. Several microduplications and microdeletions are also reportedly associated with the VACTERL association [10, 11]. However, these findings have had sparsely be replicated.
Whole exome sequencing (WES) is a powerful tool in the investigation of genomic bases of human diseases. WES covers more than 95% of the exons of human genes, and 85% of the disease-causing mutations in Mendelian disorders are identified therein. In addition, many disease-associated single nucleotide polymorphisms (SNPs) are also located in this portion of the genome. Accordingly, WES is expected to provide information on nearly all functional, protein-coding regions in individuals with agnogenic congenital disorders [12, 13].
In this study, we analyzed a Japanese patient with the VACTERL association and his parents by WES and identified a novel frameshift mutation in the PCSK5 gene. Since several non-synonymous variants in the cysteine-rich motif (CRM) of the PCSK5 gene have been described in patients with the VACTERL association as a potentially responsible gene for this condition [14, 15], we further studied the causal role of the PCSK5 mutation in the VACTERL association.
A 10-year-old Japanese boy with the Fallot tetralogy (ventricular septal defect, branch pulmonary artery stenosis, right ventricular hypertrophy, and aortic overriding), congenital scoliosis affecting the eighth thoracic hemivertebra, unilateral renal anomaly of the right kidney, limb anomalies and dislocation of the right knee, and choroidal coloboma was referred to our institution. On the basis of four out of seven component features, the patient was diagnosed as having the VACTERL association by a board-certified pediatrician. The patient was born by vaginal delivery following a gestation period of 37 weeks and weighed 2,740 g at birth. He underwent heart surgery using a Rastelli procedure when he was 2 years of age. Since then, he has led a normal life apart from the abstinence of vigorous exercise. The parent’s family histories were unremarkable. There was no maternal or paternal history of congenital malformations, hypertension, diabetes, miscarriage, or consanguineous marriage.
We analyzed the case of a Japanese boy with the VACTERL association based on the presence of four out of seven component features as determined by a board-certified pediatrician at our institution. After informed consent was obtained, peripheral blood samples were taken from the patient and his parents. Our study was approved by the Ethical Review Board for Human Genome Studies at Shinshu University School of Medicine.
Exome capture and sequencing
A total of 3 μg of genomic DNA obtained from peripheral leukocytes of the patient and his parents was subjected to an exome capture procedure using the Agilent Sure Select Human All Exon Kit V5 (Agilent Technology, Tokyo) (target size: 50 megabases [Mb]) according to the manufacturer’s protocols. Briefly, genomic DNA was fragmented with Covaris (Covaris Inc., MA, USA) and purified using Agencourt AMPure XP beads (Beckman Coulter Inc., Tokyo). The quality of the fragmentation and purification was assessed with an Agilent 2100 Bioanalyzer. The fragment ends were repaired, an adenosine residue was added to the 3′ end of the fragments, and then SureSelect adaptors were ligated to the fragments using an Agilent SureSelect Library Prep Kit, ILM. At each step, the fragments were purified using Agencourt AMPure XP beads. The DNA libraries were amplified by PCR following quality confirmation. Exon-containing fragments were captured by biotinylated RNA library “baits” using an Agilent SureSelect Target Enrichment Kit. The captured DNA was purified with streptavidin-coated magnetic beads and re-amplified. The DNA libraries of the patient and his parents were sequenced using the illumina HiSeq system in 101-base-pair (bp) paired-end reads.
Alignment and variant calling
Sequence reads were aligned to the human reference genome (GRCh37/hg19 + decoy sequences) obtained from the 1000 Genomes FTP site (ftp://1000genomes.ebi.ac.uk/) using the Burrows-Wheeler Aligner (BWA) version 0.6.2 . Multiple identical reads from the exact same fragment were marked as duplicates and removed using Picard Tools version 1.8.3 (http://picard.sourceforge.net.). Exome target regions of the SureSelect V5 were downloaded from the Agilent SureDesign website (https://earray.chem.agilent.com/suredesign/) for calculation of read depth to assess coverage. Sequence data and coverage are shown in Additional file 1: Table S1. Local realignment around known indels, base quality score recalibration (BQSR), and variant calling using UnifiedGenotyper variant caller were performed with the Genome Analysis Toolkit (GATK) 2 version 2013.2 obtained from the Appistry website (http://www.appistry.com/gatk). Variant quality score recalibration (VQSR) using the GATK was performed to filter out false positive variants. This method builds a probabilistic model from a training set of known true mutations and assigns an accurate confidence score to each putative mutation call. The workflow of the GATK and its parameters were carried out as recommended by the GATK best practice guide (http://www.appistry.com/gatk).
Additional WES Data for variant filtering and VQSR
We previously sequenced the genomes of 30 Japanese people, including 19 apparently healthy individuals. We included these data in our variant analysis since: (1) mutations found in the genomes of healthy individuals are not likely disease-causing and therefore are usable for further variant filtering; and (2) the VQSR of GATK requires at least 30 WES samples to achieve optimal results according to the GATK best practice guide. With these extra subjects, a total of 33 samples were available for variant analysis.
Variant annotation, filtering, and family-based analysis
Functional annotations of the Ensembl database GRCh37.70  and the possible effects of variants were added using SnpEff version 3.2 . Using these annotations, VQSR-passing variants were filtered first for those that were predicted to alter amino acid sequences (missense, nonsense, and splice-site mutations and indels in coding regions), and then for those that were rare (<1.0% Minor Allele Frequencies [MAF] in the HapMap-JPT [Japanese in Tokyo, Japan] or the 1000 Genomes ASN [the East Asian population, composed mostly of Japanese and Chinese] databases). The filtering procedure of variants and their numbers are presented in Additional file 1: Table S1. Genetic mutations classified as de novo, dominant, or recessive (homozygotes or compound heterozygotes) were identified by trio-family-based analysis using an in-house script. The variants were also compared with the recently-released Human Genetic Variation Database (HGVD) (http://www.genome.med.kyoto-u.ac.jp/SnpDB/), which contained genetic variations determined by WES in 1,208 Japanese individuals.
The variants identified via WES were confirmed using the Sanger method. Primer sequences were designed using Primer-BLAST . PCR-amplified fragments were purified and then sequenced using an ABI3730xl DNA Analyzer (Life Technologies, Tokyo). Sequence data were analyzed with Sequence Scanner software version 1.0 (Life Technologies, Tokyo).
Variant filtering, family-based analysis, and validation by Sanger sequencing
After the removal of variants that did not affect protein amino acid sequences or were common (≥1.0% MAF in HapMap-JPT or 1000 Genomes ASN databases) [20, 21], we analyzed the patient with the VACTERL association and his parents to identify genetically functional variants using the following classifications: (1) de novo mutations that were non-inherited or unique; (2) homozygous recessive mutations that were derived from both parents; and (3) compound heterozygous mutations for which the patient had at least two mutations derived separately from parents in the same gene. Two de novo mutations in two genes, one homozygous recessive mutation in one gene, and four compound heterozygote mutations in one gene in the patient were identified and validated by following Sanger sequencing (Additional file 1: Table S1). These mutations were not found in our in-house database of previously seqenced WES data obtained from 19 healthy Japanese individuals.
The identified two de novo mutations included a nonsense mutation in the NT5C3L (5′-nucleotidase, cytosolic III-like) gene and a missense mutation in the TTLL9 (tubulin tyrosine ligase-like family, member 9) gene. However, since both genes are functionally unknown, we could not conclusively determine their contribution to the VACTERL association.
WES identified a novel frameshift mutation in the PCSK5 gene
Among the candidates for recessive mutations, we focused on the PCSK5 gene, which has been reported as one of the causative genes for the VACTERL association [9, 10]. We identified a novel frameshift mutation and three missense mutations as a compound heterozygote in the PCSK5 gene in the patient that were validated by Sanger sequencing (Additional file 2: Table S2). The PCSK5 gene has five transcript variants (splice variants/isoforms) according to the Ensembl database , and we called the mutations according to the longest transcript, PCSK5-201 (ENST00000545128). The frameshift mutation was c.4870_4871delTG (p.C1624fs) and a missense mutation was c.4876C>A (p.Q1626K) (Figure 1). The other two missense mutations only appeared in transcript PCSK5-003 (ENST00000376767) and we called them c.2027G>C (p.W676S) and c.2032G>C (p.G678R), respectively.
Two missense mutations in exon 14 of the PCSK5-003 transcript (p.W676S and p.G678R) were also present in the patient’s mother. The genomic region of these elements was located within a short simple repeat of (GAATG)n (Chromosome 9: 78790113-78790255, 143 bp) according to the RepeatMasker of the UCSC genome browser (http://www.repeatmasker.org/). These missense mutations were found to occur due to repeat number variations of 5-nucleotide repeats since there were nucleotide variations in each repeat unit (Figure 2a). The 5-nucleotide repeats were transcribed only in the PCSK5-003 isoform. The evolutional conservation of the elements in which the mutations were detected was analyzed by multiple alignments using the UCSC genome browser , and this region demonstrated no conservation across species (Figure 3a). Taken together, these mutations are unlikely to affect the function of the gene product. Whereas the dbSNP database contains a single submission for these mutations (rs62556590, rs77249767) , those of HapMap-JPT, 1000 Genomes (including ASN) , and HGVD have no such entries.
A frameshift mutation and the other missense mutation were also present in the boy’s father (Figure 2b). Although the dbSNP and 1000 Genomes databases do not contain records of the frameshift or the missense mutations in exons 34 (PCSK5-201) or 28 (PCSK5-002), the HGVD has entries for 4 and 2 identical mutant alleles, respectively, among a total of 858 alleles (allele frequency: 0.5 and 0.2%, respectively). Both mutations were located in the PCSK5-201 and PCSK5-002 isoforms. The frameshift mutation results in a truncated protein that disrupts a region highly conserved across species (Figure 3b), suggesting that it might cause a functional alteration of PCSK5. The position of the missense mutation (p.Q1626K) may not be significant because it is located downstream of the frameshift mutation (Figure 2b).
PCSK5 proteins are composed of 3 or 4 domains according to predictions by the InterPro database version 46.0  (Figure 1). The protein domain that included the frameshift mutation was that of an insulin-like growth factor binding protein (IPR009030) containing a cysteine-rich motif (CRM). CRMs are found in a variety of eukaryotic proteins that are involved in receptor tyrosine kinase signaling. CRMs are also considered to be responsible for interactions with cell surface heparan sulfate proteoglycans (HSPGs) via tissue inhibitors of metalloproteases (TIMPs) . The frameshift mutation converts a cysteine residue to a stop codon to produce a truncated protein with the loss of the C-terminal half of the CRM domain or a loss of protein production due to nonsense-mediated mRNA decay.
Since we identified the p.C1624fs (PCSK5-201) frameshift mutation in the database for the general population (HGVD), we next searched for other possible deleterious mutations in the PCSK5 gene (Additional file 2: Table S2). We identified two nonsense mutations, c.3052C>T (p.Q1018X, rs373172614) and c.4615C>T (p.R1539X, rs140709872), in the PCSK5 gene in the variant database from the general population. The frequency of these variants is 0.02%, which is several-fold higher than that of the VACTERL association. We also identified five frameshift mutations, c.1490_1491insACAC (p.N497fs rs138544337), c.1496_1497insC (p.R500fs, rs372197834), c.3915_3916insG (p.V1307fs, rs34898066), c.4665_4666insG (p.Y1557fs, rs34889598), and c.5095_5096insCC (p.D1699fs, rs138280866), in the general population. Altogether, it appears that the mutation in the PCSK5 gene might be innocuous and not associated with the VACTERL association.
Discussion and conclusion
PCSK5 is a proprotein convertase enzyme that cleaves prohormones at consensus sequences [23–25]. PCSK5 is expressed in somites, the dorsal surface ectoderm, and primordial vertebral cartilage , as well as in the skeletal regions of the developing vertebrae, limbs, and craniofacium . Recessive PCSK5 mutant mice induced by ethylnitrosourea exhibited hypoplastic hindlimbs, absent tail, cardiac malformations, palatal agenesis, tracheoesophageal malformation, pulmonary hypoplasia, exomphalos, and renal agenesis , and epiblast-specific PCSK5 conditional knockout mice using Meox2Cre showed similar phenotypes . These findings demonstrate that PCSK5 plays a pivotal role in skeletogenesis and organogenesis in mice. Meanwhile, non-synonymous PCSK5 mutations were implicated in 4 of the 36 cases of the VACTERL association in humans reported by Szumska et al. and in 3 of the 39 cases described by Winberg et al. [14, 15].
It has also been reported that GDF11 is involved in murine skeletogenesis . PCSK5 can cleave the GDF11 propeptide into its mature form , which suggests that PCSK5 activates GDF11 in mice. Mutant PCSK5 proteins were seen to not cleave the pro-protein of GDF11 . In addition, GDF11 knockout mice showed absent tail, palatal agenesis, renal agenesis, and increased numbers of thoracic vertebrae and ribs, all of which resembled the phenotypes of PCSK5-deficient mice . Therefore, it is likely that the mechanism by which such phenotypes are observed in mice is via an inability to produce active GDF11 .
We witnessed a novel frameshift mutation in the PCSK5 gene in our patient with the VACTERL association. Since three heterozygous missense mutations , four heterozygous missense mutations, and one homozygous missense mutation in the PCSK5 gene  were previously reported in patients with the VACTERL association, we first thought this might be a strong candidate gene that was responsible for the patient’s phenotype. However, this is improbable for several reasons: (1) the healthy father also carries this mutation, (2) databases for the general population include this frameshift mutation, and (3) databases for the general population include not only this mutation, but also other deleterious mutations in the PCSK5 gene, such as nonsense and frameshift mutations. Furthermore, the frequency of these variants is several-fold higher than that of the VACTERL association. A recessive trait was unlikely for the VACTERL association since there have been no previous reports on the simultaneous occurrence of this condition in siblings. Taken together, the PCSK5 mutation appears to be benign and unrelated to the etiology of the VACTERL association.
It remains possible that the effect of the PCSK5 frameshift mutation observed in this study might be dominant with incomplete penetrance. To date, Winberg and other groups have suggested that detected PCSK5 variants could represent pathogenic entities with reduced penetrance [14, 15]. It is also possible that the maternally transmitted missense mutation, although innocuous alone, can affect phenotype when combined with the paternal frameshift mutation. Further mutations or copy number variations in other genes might have affected phenotype as well. Moreover, since the etiology of the VACTERL association is believed to be multifactorial, both modifier genes and environmental factors, such as hyperglycemia, may have precluded the development of this malformation .
In addition to the PCSK5 mutations, we identified two de novo mutations in our patient. One was identified in the NT5C3L gene encoding a cytosolic nucleotidase-like protein, which might play a role in RNA metabolism. The other was harbored in the TTLL9 gene encoding a ciliary protein that may act in tubulin glutamylation [29, 30]. Since ciliary genes are involved in some human malformation syndromes, they represent another candidate gene responsible for the patient’s phenotype. Further studies are necessary to clarify the elusive etiology of the VACTERL association.
Solomon BD (2011) VACTERL/VATER Association. Orphanet J Rare Dis 6:56
Reardon W, Zhou XP, Eng C (2001) A novel germline mutation of the PTEN gene in a patient with macrocephaly, ventricular dilatation, and features of VATER association. J Med Genet 38:820e3
Garcia-Barceló MM, Wong KK, Lui VC, Yuan ZW, So MT, Ngan ES et al (2008) Identification of a HOXD13 mutation in a VACTERL patient. Am J Med Genet 146A:3181e5
Stankiewicz P, Sen P, Bhatt SS, Storer M, Xia Z, Bejjani BA et al (2009) Genomic and genic deletions of the FOX gene cluster on 16q24.1 and inactivating mutations of FOXF1 cause alveolar capillary dysplasia and other malformations. Am J Hum Genet 84:780–791
Wessels MW, Kuchinka B, Heydanus R, Smit BJ, Dooijes D, de Krijger RR et al (2010) Polyalanine expansion in the ZIC3 gene leading to X- linked heterotaxy with VACTERL association: a new polyalanine dis- order? J Med Genet 47:351–355
McCauley J, Masand N, McGowan R, Rajagopalan S, Hunter A, Michaud JL et al (2011) X-linked VACTERL with hydrocephalus syndrome: further delineation of the phenotype caused by FANCB mutations. Am J Med Genet 155A:2370–2380
Saisawat P, Kohl S, Hilger AC, Hwang DY, Yung Gee H, Dworschak GC et al (2014) Whole-exome resequencing reveals recessive mutations in TRAP1 in individuals with CAKUT and VACTERL association. Kidney Int 85:1310–1317
Chung B, Shaffer LG, Keating S, Johnson J, Casey B, Chitayat D (2011) From VACTERL-H to heterotaxy: variable expressivity of ZIC3-related disorders. Am J Med Genet A 155A:1123–1128
Schramm C, Draaken M, Bartels E, Boemers TM, Aretz S, Brockschmidt FF et al (2011) De novo microduplication at 22q11.21 in a patient with VACTERL association. Eur J Med Genet 54:9–13
Hilger A, Schramm C, Pennimpede T, Wittler L, Dworschak GC, Bartels E et al (2013) De novo microduplications at 1q41, 2q37.3, and 8q24.3 in patients with VATER/VACTERL association. Eur J Hum Genet 21:1377–1382
Dworschak GC, Draaken M, Marcelis C, de Blaauw I, Pfundt R, van Rooij IA et al (2013) De novo 13q deletions in two patients with mild anorectal malformations as part of VATER/VACTERL and VATER/VACTERL-like association and analysis of EFNB2 in patients with anorectal malformations. Am J Med Genet A 161A:3035–3041
Ng SB, Buckingham KJ, Lee C, Bigham AW, Tabor HK, Dent KM et al (2010) Exome sequencing identifies the cause of a mendelian disorder. Nat Genet 42:30–35
Rabbani B, Tekin M, Mahdieh N (2014) The promise of whole-exome sequencing in medical genetics. J Hum Genet 59:5–15
Winberg J, Gustavsson P, Papadogiannakis N, Sahlin E, Bradley F, Nordenskjöld E et al (2014) Mutation Screening and Array Comparative Genomic Hybridization Using a 180 K Oligonucleotide Array in VACTERL Association. PLoS One 9:e85313
Szumska D, Pieles G, Essalmani R, Bilski M, Mesnard D, Kaur K et al (2008) VACTERL/caudal regression/Currarino syndrome-like malformations in mice with mutation in the proprotein convertase Pcsk5. Genes Dev 22:1465–1477
Li H, Durbin R (2009) Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25:1754–1760
Flicek P, Amode MR, Barrell D, Beal K, Billis K, Brent S et al (2014) Ensembl 2014. Nucleic Acids Res 42 (Database issue): D749–D755
Cingolani P, Platts A, le Wang L, Coon M, Nguyen T, Wang L et al (2012) A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3”. Fly (Austin) 6:80–92
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL (2012) Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinform 13:134
Sherry ST, Ward MH, Kholodov M, Baker J, Phan L, Smigielski EM et al (2001) dbSNP: the NCBI database of genetic variation. Nucleic Acids Res 29:308–311
1000 Genomes Project Consortium, Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM et al (2012) An integrated map of genetic variation from 1,092 human genomes. Nature 491:56–65
Zdobnov EM, Apweiler R (2001) InterProScan—an integration platform for the signature-recognition methods in InterPro. Bioinformatics 17:847–848
Constam DB, Calfon M, Robertson EJ (1996) SPC4, SPC6, and the novel protease SPC7 are coexpressed with bone morphogenetic proteins at distinct sites during embryogenesis. J Cell Biol 134:181–191
Rancourt SL, Rancourt DE (1997) Murine subtilisin-like proteinase SPC6 is expressed during embryonic implantation, somitogenesis, and skeletal formation. Dev Genet 21:75–81
Essalmani R, Zaid A, Marcinkiewicz J, Chamberland A, Pasquato A, Seidah NG et al (2008) In vivo functions of the proprotein convertase PC5/6 during mouse development: Gdf11 is a likely substrate. Proc Natl Acad Sci USA 105:5750–5755
Essalmani R, Hamelin J, Marcinkiewicz J, Chamberland A, Mbikay M, Chrétien M et al (2006) Deletion of the gene encoding proprotein convertase 5/6 causes early embryonic lethality in the mouse. Mol Cell Biol 26:354–361
McPherron AC, Lawler AM, Lee SJ (1999) Regulation of anterior/posterior patterning of the axial skeleton by growth/differentiation factor 11. Nat Genet 22:260–264
Solomon BD (2011) VACTERL/VATER Association. Orphanet J Rare Dis. 16(6):56. doi:10.1186/1750-1172-6-56 (Review)
Wloga D, Rogowski K, Sharma N, Van Dijk J, Janke C, Eddé B et al (2008) Glutamylation on alpha-tubulin is not essential but affects the assembly and functions of a subset of microtubules in Tetrahymena thermophila. Eukaryot Cell 7:1362–1372
Kubo T, Yanagisawa HA, Yagi T, Hirono M, Kamiya R (2010) Tubulin polyglutamylation regulates axonemal motility by modulating activities of inner-arm dyneins. Curr Biol 20:441–445
SK, SS, MT, NI, HT, MM, FF, HI, TK, TK, HK, and HK participated in the design of the study and helped to draft the manuscript. YN conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
This study was supported by grants from the Takeda Science Foundation, the Japan Orthopaedics and Traumatology Foundation, Inc. (No.227), the Kanae Foundation for the Promotion of Medical Science, and the Naito Foundation to Y.N.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Consent for publication Written informed consent was obtained from the patient for publication of this Case report and any accompanying images. A copy of the written consent is available for review by the Editor of this journal.
About this article
Cite this article
Nakamura, Y., Kikugawa, S., Seki, S. et al. PCSK5 mutation in a patient with the VACTERL association. BMC Res Notes 8, 228 (2015). https://doi.org/10.1186/s13104-015-1166-0