Porokeratosis (PK) is a heterogeneous group of keratinization disorders. No causal genes except MVK have been identified, even though the disease was linked to several genomic loci. Here, we performed massively parallel sequencing and exonic CNV screening of 12 isoprenoid genes in 134 index PK patients (61 familial and 73 sporadic) and identified causal mutations in three novel genes (PMVK, MVD, and FDPS) in addition to MVK in the mevalonate pathway. Allelic expression imbalance (AEI) assays were performed in 13 lesional tissues. At least one mutation in one of the four genes in the mevalonate pathway was found in 60 (98%) familial and 53 (73%) sporadic patients, which suggests that isoprenoid biosynthesis via the mevalonate pathway may play a role in the pathogenesis of PK. Significantly reduced expression of the wild allele was common in lesional tissues due to gene conversion or some other unknown mechanism. A G-to-A RNA editing was observed in one lesional tissue without AEI. In addition, we observed correlations between the mutations in the four mevalonate pathway genes and clinical manifestations in the PK patients, which might support a new and simplified classification of PK under the guidance of genetic testing.https://doi.org/10.7554/eLife.06322.001
Porokeratosis refers to a group of around twenty skin conditions that involve a build-up of a protein called keratin in skin cells. Keratin forms the tough fibres that give strength to hair and nails, and people suffering from porokeratosis develop hardened skin lesions. Porokeratosis is an uncommon condition; most cases are inherited and have been linked to exposure to ultraviolet light and having a weakened immune system.
Mutations in one gene called MVK are known to cause two forms of the disorder, but it is suspected that other genetic causes of porokeratosis will also be identified. The MVK gene encodes an enzyme that is involved in making chemicals called isoprenoids. This large and diverse class of chemicals provides the building blocks for making many other important molecules in all living species. Zhang, Li et al. have now analysed genetic material from 134 different porokeratosis patients to search for mutations in other genes involved in the production of isoprenoids. The patients examined include 61 people with a family history of the disorder, and 73 cases in which the condition seems to be a one-off occurrence.
This search identified mutations in three additional genes (called PMVK, MVD and FDPS) that are all linked to porokeratosis. Further analysis of these three genes and MVK revealed that about half of the patients with mutations in the MVK gene developed large lesions (that were over 5 centimetres in diameter). However, those with mutations in the other three genes did not develop such large lesions. Mutations in some of the newly identified genes were instead linked to porokeratosis affecting specific areas of the body; for example, PMVK and MVD mutations are linked to porokeratosis localized to the genitals and around the eyes, respectively. This means that, in the future, doctors might be able to simplify the diagnosis of the different varieties of porokeratosis based on information gained via genetic tests.https://doi.org/10.7554/eLife.06322.002
Porokeratosis (PK, MIM 175800) is a heterogeneous group of keratinization disorders that exhibit an autosomal dominant mode of inheritance. PK is also a skin-specific autoinflammatory disease which was often inherited and linked to ultraviolet light exposure and immunosuppression (Schamroth et al., 1997; Abramovits and Oquendo, 2013). For example, eruptive pruritic papular porokeratosis exemplifies the inflammatory manifestation, and complications to inflammatory conditions such as localized cutaneous amyloidosis are seen in PK patients (Biswas, 2015). Indeed, PK and psoriasis share some features at both clinical and molecular levels and sometimes coexist in the same patients (Zhang et al., 2008).
As a histological hallmark that unifies all variants of PK, cornoid lamella (CL) is a vertical ‘column’ of parakeratosis. The pattern of CL can be slender, broad, or confluent, which is related to epidermal hyperplasia and dermal inflammation. However, CL is not a unique feature of PK because it can be seen in some inflammatory and inherited cutaneous disorders and also as an incidental finding (Biswas, 2015).
PK is currently classified according to the clinical manifestations, such as number, size, morphology, and distribution of the histological lesions. A better system of classification is expected because some variants of PK are fraught with confusing terminology (Schamroth et al., 1997; Sertznig et al., 2012; Biswas, 2015). For example, it is sometimes hard to completely differentiate disseminated superficial actinic porokeratosis (DSAP) from disseminated superficial porokeratosis (DSP) by age of onset and sun-exposed areas.
In addition to the heterogeneity in clinical manifestations, genetic heterogeneity is also observed in PK. At least five linkage loci (i.e., 12q23.2-24.1, 15q25.1-26.1, 18p11.3, 1p31.3-p31.1, 16q24.1-24.3) have been reported for the disseminated forms of PK which include DSAP, DSP, porokeratosis palmaris et plantaris disseminata (PPPD), and immunosuppression-induced porokeratosis (ISIP) (Schamroth et al., 1997; Luan et al., 2011). However, only one causal gene, the mevalonate kinase gene (MVK) at 12q24, has been identified in DSAP and porokeratosis of Mibille (PM) (Zhang et al., 2012; Li and Zhang, 2014; Zeng et al., 2014). Here we performed genetic analysis in 134 index patients with PK to identify additional causal genes and to establish new parameters for better differential diagnosis of PK. Allelic expression imbalance (AEI) and cDNA mutation were analyzed in 13 pairwise lesional tissues (LTs) and neighboring normal-appearing skin (NNS) to investigate the underlying pathogenicity of PK.
We previously identified a linkage locus for DSAP on chromosome 16q24.1-24.3 in a four-generation Chinese DSAP family (Luan et al., 2011). We applied the candidate gene approach and hypothesized that the MVD gene in the 16q24.1-24.3 region was the causal gene of DSAP in the PK family. We performed Sanger sequencing of all MVD exons and identified a c.746T>C (p.Phe249Ser) mutation in MVD. The c.746T>C mutation displayed 100% co-segregation with the PK phenotype in this family (Figure 1).
Since both MVD and MVK are involved in isoprenoid biosynthesis via the mevalonate pathway (Thurnher et al., 2013), we hypothesized that mutations in other members of the mevalonate pathway exist in PK and screened for mutations in 12 genes in the mevalonate pathway including MVD and MVK in 134 PK patients (61 familial and 73 sporadic). The 12 genes are: acetyl-CoA acetyltransferase (ACAT) 1 and 2, 3-hydroxy-3-methylglutaryl-CoA synthase (HMGCS) 1 and 2, 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR), MVK, phosphomevalonate kinase (PMVK), MVD, isopentenyl-diphosphate delta isomerase (IDI) 1 and 2, farnesyl diphosphate synthase (FDPS), and geranylgeranyl diphosphate synthase 1 (GGPS1) (Kanehisa et al., 2012) (Figure 2). A cycled primer extension and ligation-dependent amplification (CPELA) reaction for enrichment of target genomic fragments and next generation sequencing (NGS) were performed in each sample. The NGS data were analyzed for rare single nucleotide variants (SNVs) with no record in dbSNP or a minor allele frequency of <1% in Chinese from 1000 genome database. After Sanger sequencing validation, 60 out of 61 SNVs called from NGS data were demonstrated to be true SNVs, which were predicted to impact protein function. Of these, 12 non-pathogenic SNVs were excluded from subsequent analysis (Supplementary file 1). In total, 48 mutation sites were identified in MVK, PMVK, MVD, and FDPS (Figure 3, Figure 3—source data 1).
For 27 samples with no confirmed point mutation in MVK, PMVK, MVD, and FDPS genes, large deletion or duplication mutations were interrogated by analyzing the copy number of promoters and all exons in these genes. As a result, three deletion mutations were identified by the CNVplex assay. These three deletions were identified in four (F-26, F-40, S-69 and S-21), one (S-56) and one (F-6) cases, respectively. The breakpoints of three large deletions were also determined (Figure 3, Figure 3—figure supplement 1).
We determined the genomic variations by Sanger sequencing in other family members apart from one with an unknown mutation. In 20 pedigrees with more than two blood samples collected, all patients carried the same mutation with the proband. On the other hand, no mutation was detected in any of the healthy family members. Except for one common MVD mutation of c.746T>C, all deleterious candidate SNVs and large deletions were undetected in healthy family members and 270 unrelated ethnically matched controls. In total, we identified 13 mutations in MVD, six mutations in PMVK, four mutations in FDPS, and 28 mutations in MVK in 113 of the 134 PK patients (Supplementary file 2). Notably, 21 of the 28 MVK mutations are novel (Li and Zhang, 2014). These mutations include 30 missense (58%), two start codon (4%), five nonsense (10%), seven indels (14%), four splicing defects (8%), and three large deletions (6%), and predictably affect functions of the respective enzymes. The two MVD mutations (c.746T>C and c.875A>G) were identified in 50 unrelated patients, accounting for 81% of all patients with MVD mutations. Among the 60 familial patients in whom at least one mutation was found, 30 had MVD mutations, three had PMVK mutations, two had FDPS mutations, 24 had MVK mutations, and one had both MVK and MVD mutations. Among the 53 sporadic patients in whom at least one mutation was found, 31 had MVD mutations, six had PMVK mutations, two had FDPS mutations, and 14 had MVK mutations. No mutation was found in the remaining 21 PK patients.
For the 21 PK patients in whom no pathogenic mutation in MVK, PMVK, MVD, or FDPS was found, we sequenced all exons of the SLC17A9 gene by Sanger sequencing. Mutations in the SLC17A9 gene were recently reported to be another genetic cause for DSAP (Cui et al., 2014). However, no mutation in the SLC17A9 gene was found in the 21 PK patients studied.
To explore the genotype–phenotype correlations, we studied the clinical records of 134 index patients (92 males and 42 females) with PK. In all patients, the clinical diagnosis was supported by histopathological examinations. The onset of the disease ranged from birth to 78 years (mean age 31 years). The distribution of lesions could be categorized as localized and disseminated forms (Schamroth et al., 1997). Several classical variants were found to coexist in the same patient. However, no malignant degeneration was found in the lesions of these patients. Although all affected members in the same family carried the same mutation, each affected member showed different clinical manifestations and severity.
We observed some interesting correlations between gene mutations and clinical phenotypes of PK (Table 1 and Figure 4). First, giant plaque-type porokeratosis ptychotropica (PPt) with lesion diameters at least 5 cm appeared to be a unique phenotype associated with MVK mutations. This feature was observed in 50% (19/38) of index patients with MVK mutations, but not in any patients with MVD, PMVK, or FDPS mutations. In addition, patients with MVK mutations generally showed the widest range of phenotypes in terms of both the number and size of lesions. Second, localized genital PK (Chen et al., 2006) and porokeratoma (Walsh et al., 2007) appeared to be unique phenotypes associated with PMVK mutations. Third, in PK patients with MVD mutations, the age of onset spanned from 5 to 70 years, and the diameter of the lesions was generally less than 2 cm. Five females and one male with MVD mutations manifested mild solar facial PK (Sharquie and Al-Baghdady, 2003), which was not found in any patients with mutations in other genes. Fourth, in patients with FDPS mutations, the number of lesions was generally more than 500 and the diameter of the lesions was less than 1 cm.
As a whole, the lesions harboring mutations in MVD and FDPS genes tended to be more homogeneous and superficial than those carrying mutations in MVK and PMVK genes. The phenotypes of MVD and FDPS mutations on the natal cleft were mild and unobvious. However, some patients harboring MVK and PMVK mutations firstly presented with untypical lesions of genitogluteal porokeratosis (GGP).
The expression of the wild allele in 10 out of 13 LTs was significantly reduced at the 1% significance level (Figures 5–7, Supplementary file 3). In contrast, the expression of the mutant allele was significantly reduced in almost all NNS (except one from F-42) carrying a mutation resulting in nonsense, frameshift, or splicing defect. As for LT from F-31, the mutant to wild allelic ratio in genomic DNA was >3 compared with the ratio of about 1 in the corresponding NNS genomic DNA. However, no deletion or duplication of all MVK exons was detected in F-31, which might indicate that gene conversion of the wild to the mutant allele occurred in the LT (Figure 5). Gene conversion might also occur in the thigh LT of F-43 (Figure 7). No copy number change or CpG island methylation was detected at the corresponding mutated genes in all LTs. A mutation of c.1003G>A (p.Gly335Ser) on the wild allele (T) of MVK: c.1093T>A was identified in the cDNA of the LT without AEI but not the cDNA of NNS collected from the left forearm of F-38 (Figure 6). Surprisingly, this mutation of c.1003G>A was not identified in the genomic DNA of this LT either, which indicated that G-to-A RNA editing occurred on the wild allele of MVK: c.1093T>A. No mutation was observed in cDNA samples from other tissues.
In this present study we were able to identify mutations in four genes (MVK, PMVK, MVD, and FDPS) underlying PK by massively parallel sequencing and copy number variation analysis in 134 index patients with PK. In the 12 PK families where the affected members span three to six generations, we found mutations only in MVK or MVD but not in PMVK or FDPS (Figure 1, Figure 1—figure supplement 1). We also observed the correlation of genotypes and several clinical manifestations in the PK patients (Table 1), which might be helpful to simplify the classification of PK under the guidance of genetic testing.
Interestingly, two patients were identified to carry two mutations: c.302C>G (p.Pro101Arg) and c.683G>A (p.Arg228Gln) in MVD for S-62, and c.746T>C (p.Phe249Ser) in MVD and c.235G>A (p.Asp79Asn) in MVK for F-28. According to the prediction programs of bioinformatics, the c.235G>A in MVK and c.302C>G in MVD were less conservative than the other mutation carried in F-28 and S-62, respectively. The two mutations in MVD for S-62 were located in the trans position, for his daughter carried only one of them (Figure 1—figure supplement 2). S-62 had late-onset PK at the age of 50 and his 30-year-old daughter had no lesions to date. According to our observation, index patients F-28 and S-62 did not present with more severe phenotypes or earlier age of onset or any other unique clinical features.
In this study we did not detect pathogenic mutation in 21 patients: five cases with linear PK (LP), six cases with solitary plaque-type PK, and 10 cases with disseminated PK (only one had family history). There were three possible explanations: (1) these patients might have other non-classical or mosaic/somatic mutations that could not be detected by current mutation screening methods; (2) these patients might have mutations in unknown genes; (3) these patients might be acquired ISIP without mutation. Recently, mutations in the SLC17A9 gene were reported to be another genetic cause for DSAP (Cui et al., 2014). We sequenced all exons of SLC17A9 in these 21 patients by Sanger sequencing but found no mutations in SLC17A9.
RNA editing is a process of some discrete changes to specific nucleotide sequences of a RNA molecule which may be critical during normal development and diseases (Witkin, et al., 2015). G-to-A RNA editing in PK was first reported in our study. Interestingly, the RNA editing seemed to occur selectively on the one copy of MVK and this selective editing pattern seemed to be capable of being transmitted during mitosis. The mechanism underlying this observation should be investigated further.
Using AEI assays, we observed that the expression of wild allele was significantly reduced to various extents in most LTs (>77%). The variation in the expression reduction level might be due to the LT containing various amounts of normal-appearing keratinocytes and other different types of cells. Somatic loss of function of the wild allele in the epidermal stem cells might play a key role in the formation of lesions, and it might occur spontaneously and be induced by environmental factors such as ultraviolet light.
There was only one RNA editing identified in the LT, therefore the somatic loss of function of the wild allele in LTs of PK could be mainly caused by suppression of the wild allele through gene conversion in a minor proportion of cases or an unknown DNA methylation-independent epigenetic mechanism in most cases. Allelic imbalance of histone methylation modification unlinked to DNA methylation, whose pattern is able to be transmitted during cell mitosis, could be the DNA methylation-independent epigenetic mechanism (Xu et al., 2011). Further investigation would be necessary to test this hypothesis. We did not find reduced expression or somatic mutation of the wild allele in only two LTs. The reason for this could be that the LT might contain many normal cells affecting the sensitivity of AEI and mutation detection.
As a starter button of PK, gene conversion of the wild allele to the mutant allele was first identified in this study. This observation is generally consistent with the phenomenon known as ‘revertant mosaicism’, in which the mutant allele is corrected in some stem cells. Revertant mosaicism is mainly identified in blood diseases such as Wiskott–Aldrich syndrome and Fanconi’s anemia, and skin diseases such as ichthyosis and epidermolysis bullosa (May, 2011). Revertant mosaicism shows the bright side of gene conversion, but our observation revealed the opposite side. Therefore, gene conversion could be a natural ‘rescue angel’ in some recessive diseases or a ‘damaging devil’ in some dominant diseases, which could mostly affect tissues with extensive ability to self-renew for life.
Our findings support the notion that Mendelian diseases can be caused by multiple genes involved in the same metabolic or signal transduction pathway, or producing proteins on interaction. The four genes MVK, PMVK, MVD, and FDPS are members of the mevalonate pathway, suggesting that this pathway is involved in the pathogenesis of PK. The mevalonate pathway of isoprenoid biosynthesis provides precursors of isoprenoids, which are ubiquitous in living species and diverse in biological function (Smit and Mushegian, 2000) and serve as precursors of cholesterol, heme A, ubiquinones, dolichol, and isoprenylated proteins (e.g., RAS and Lamin B) (Figure 2). Moreover, it is also the subject of many pharmacological regulating drugs such as statins. It is well known that cholesterol is an important constituent of the cell membrane of most eukaryotic cells, and several inherited disorders have been linked to defects in cholesterol biosynthesis. The isoprenylated proteins play a role in the regulation of cell growth, division, and differentiation (Moir and Spann, 2001; Zwerger et al., 2011; Davidson and Lammerding, 2014) and are probably associated with the retained nuclei in the stratum corneum (i.e., parakeratosis). We propose that the accumulation of abnormal metabolites or shortage of isoprenoids might predispose patients to idiopathic inflammation of the skin, and that correction of the abnormal isoprenoid biosynthesis might be a novel therapeutic direction for PK. Further functional studies will be necessary to confirm this hypothesis.
All procedures followed the guidelines of the Helsinki Declaration and were approved by the ethics committee and by the Scientific Ethical Committee of Fudan University. Participating centers provided local Institutional Review Board approval for genetic analysis. Study participants provided informed consent for genetic testing.
From 2001 to 2013, peripheral blood samples were collected from 134 index patients with PK (92 males and 42 females) and 180 family members (103 males and 77 females). Thirteen pairwise LTs and NNS were dissected from 12 index patients with a MVK, MVD, or FDPS mutation (Supplementary file 3). In addition, blood samples were collected from 270 healthy adult individuals undergoing routine medical examination in the hospital. The 134 index patients were diagnosed by at least two experienced dermatologists, based on both clinical features and histological examinations. In view of the family history, these index patients were divided into family (F) and sporadic (S) cases. Among the 61 familial cases, 40 index patients were probands from each family. We examined members of 21 families from different provinces of China, whose generations ranged from 2 to 6.
Genomic DNA was extracted from the peripheral blood using the QIAamp DNA blood mini kit (Qiagen, Germany). Total RNA was isolated from the skin using the RNeasy Protect Mini Kit (Qiagen, Germany) according to the manufacturer's instructions.
PCR products were sequenced by BigDye3.1 (Applied Biosystems, Foster City, CA, USA). Raw data were collected on an ABI 3130XL sequencer and mutations were called by the PolyPhred program (Nickerson et al., 1997). The 60 SNVs in 270 control samples were genotyped using the ABI PRISM SNaPshot Multiplex Kit (Applied Biosystems, Foster City, CA, USA) and ABI3730xl sequencer.
A total of 162 fragments covering the promoter region, 5′UTR, coding sequence and splicing site of 12 genes in the mevalonate pathway (ACAT1, ACAT2, HMGCS1, HMGCS2, HMGCR, MVK, PMVK, MVD, IDI1, IDI2, FDPS, and GGPS1) were amplified using the EasyTarget® amplification kit (Genesky Biotechnologies Inc, Shanghai, China) which was developed according to the CPELA method. The CPELA was a new method developed by Genesky Biotechnologies for fast and simple enrichment of multiple gene regions for massively parallel sequencing. The principle of CPELA is described in Figure 3—figure supplement 2.
The amplification products by CEPLA were mixed and size-separated by 2% agarose gel electrophoresis. Products of length 200–400bp were recovered. The final concentration of the library mixture was determined by real-time quantitative PCR and the average concentration for each library was estimated. The libraries tagged with different index sequences from several projects were well mixed in an appropriate concentration for each library corresponding to its aimed sequencing depth. The final library mixture was sequenced on MiSeq sequencer (Illumina) using Miseq reagent kit v2. The sequencing reads were separated for each sample by running CASAVA (Illumina Inc, San Diego, CA, USA) and target reads were determined by comparing them with fragment reference sequences (hg19) using the Blat program (Kent, 2002).
In order to reduce the reading error from the sequencing reaction step, we first compared the two reads from the same cluster and integrated them into one correction read. The correction reads were then aligned to hg19 using the Burrows–Wheeler Aligner (BWA) (Li and Durbin, 2010). SNV calling was performed using both GATK (McKenna et al., 2010) and Varscan programs (Koboldt et al., 2012), and the called SNV data were then combined. The Annovar program was used for SNV annotation (Wang et al., 2010). The functional effect of non-synonymous SNVs was assessed by the PolyPhen-2, SIFT, and MutationTaster (Ng and Henikoff, 2003; Adzhubei et al., 2010; Schwarz et al., 2010). Non-synonymous SNVs with SIFT score of <0.05, Polyphen-2 score of >0.85 or MutationTaster score of >0.85 were considered as significant of not being benign. To sort potentially deleterious variants from benign polymorphisms, perl scripts were used to filter the SNVs against those of dbSNP135. Any SNV recorded in dbSNP135 and with a minor allele frequency of ≥1% in Chinese from 1000 genome database was considered as benign polymorphisms and therefore removed for subsequent analysis.
The copy number of the target regions was measured by a CNVplex assay, a high-throughput multiplex CNV analysis method recently developed by Genesky Biotechnologies. The principle of CNVplex technology is described in Figure 3—figure supplement 3. We utilized this technology for quantitative analysis of copy numbers of all 37 exons and upstream promoter regions in the MVK, MVD, PMVK, and FDPS genes for the blood DNA samples with no point mutation identified and the tissue DNA samples. Based on the copy number measurements for all target sequences, the breakpoints were estimated to be located between two neighboring probe target sites showing different copy numbers. Several primer sets flanking the two probe target sites were tested to amplify the target region from case and control DNA samples using a long PCR protocol. Specific PCR products from case samples were sequenced using the ABI BigDye3.1 and the breakpoints were determined by blasting the sequences with human reference genome assembly.
Both DNA and RNA were extracted from pairwise LT and NNS. Each RNA sample was reversely transcripted into cDNA twice using polyA and N9 primer mix and Reverse Transcriptase M-MLV (RNase H-) (New England Biolabs, England). For each pairwise tissue set, single nucleotide extension was used to quantitate the ratio of the mutant to the wild allele both in the tissue cDNA and DNA using the ABI PRISM SNaPshot Multiplex Kit (Applied Biosystems), followed by normalization to the ratio value in the corresponding patient’s blood DNA as a reference of 1:1. As for cDNA from F-42 and S-36 patients with splicing defect mutations of MVK: c.371+2T>A and FDPS: c.684+1G>A respectively, fluorescent PCR followed by capillary electrophoresis was used to quantitate the ratios of mutant to wild transcripts.
Total RNA was reversely transcribed into cDNA. Two or three sets of primers were designed to amplify the cDNA of MVK, MVD, or FDPS, covering all coding sequences for each gene. The sequences of PCR products were determined by Sanger sequencing. Two allele-specific PCRs for c.1093T>A followed by Sanger sequencing were performed in order to determine whether the mutant allele (A) of c.1003G>A was in the cis position with the wild allele (T) of c.1093T>A. DNA from each pairwise LT and NNS were subject to bisulfite conversion using the EZ DNA Methylation Kit (Zymo Research Corporation, Irvine, CA, USA). PCR products amplifying the bisulfite converted DNA at CpG islands of MVK, MVD, or FDPS were sequenced on MiSeq sequencer (Illumina) using Miseq reagent kit v2 or on the ABI 3130xl genetic analyzer using the ABI BigDye3.1.
Introduction to autoinflammatory syndromes and diseasesDermatologic Clinics 31:363–385.https://doi.org/10.1016/j.det.2013.04.010
Cornoid lamellation revisited: apropos of porokeratosis with emphasis on unusual clinicopathological variantsThe American Journal of Dermatopathology 37:145–155.https://doi.org/10.1097/DAD.0000000000000039
Genital porokeratosis: a series of 10 patients and review of the literatureThe British Journal of Dermatology 155:325–329.https://doi.org/10.1111/j.1365-2133.2006.07295.x
Broken nuclei—lamins, nuclear mechanics, and diseaseTrends in Cell Biology 24:247–256.https://doi.org/10.1016/j.tcb.2013.11.004
KEGG for integration and interpretation of large-scale molecular data setsNucleic Acids Research 40:D109–D114.https://doi.org/10.1093/nar/gkr988
The structure and function of nuclear lamins: implications for diseaseCellular and Molecular Life Sciences 58:1748–1757.https://doi.org/10.1007/PL00000814
SIFT: predicting amino acid changes that affect protein functionNucleic Acids Research 31:3812–3814.https://doi.org/10.1093/nar/gkg509
Porokeratosis of Mibelli. Overview and review of the literatureActa Dermato-Venereologica 77:207–213.
Solar facial porokeratosisThe Journal of Dermatology 30:216–221.https://doi.org/10.1111/j.1346-8138.2003.tb00374.x
Biosynthesis of isoprenoids via mevalonate in Archaea: the lost pathwayGenome Research 10:1468–1484.https://doi.org/10.1101/gr.145600
Regulation of mevalonate metabolism in cancer and immune cellsBiochimica et Biophysica Acta 1831:1009–1015.https://doi.org/10.1016/j.bbalip.2013.03.003
ANNOVAR: functional annotation of genetic variants from high-throughput sequencing dataNucleic Acids Research 38:e164.https://doi.org/10.1093/nar/gkq603
RNA editing, epitranscriptomics, and processing in cancer progressionCancer Biology & Therapy 16:21–27.https://doi.org/10.4161/15384047.2014.987555
Splicing mutation in MVK is a cause of porokeratosis of MibelliArchives of Dermatological Research 306:749–755.https://doi.org/10.1007/s00403-014-1465-7
Nuclear mechanics in diseaseAnnual Review of Biomedical Engineering 13:397–428.https://doi.org/10.1146/annurev-bioeng-071910-124736
Helen H HobbsReviewing Editor; University of Texas Southwestern Medical Center, United States
eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.
Thank you for sending your work entitled “Genomic Variations of the Mevalonate Pathway in Porokeratosis” for consideration at eLife. Your article has been favorably evaluated by Stylianos Antonarakis (Senior editor) and three reviewers, one of whom, Helen Hobbs, is a member of our Board of Reviewing Editors.
The Reviewing editor and the other reviewers discussed their comments before we reached this decision, and the Reviewing editor has assembled the following comments to help you prepare a revised submission.
This paper explores the genetic basis of porokeratosis, a genetically heterogeneous autosomal dominant disorder of keratinization that includes both inherited and sporadic forms. Mutations in mevalonate kinase (MVK), a key enzyme in sterol synthesis, were shown previously by these authors to cause the disease in a subset of individuals, but the etiology in most affected individuals has remained unknown. Here the authors test the hypothesis that the disease is due to defective isoprenoid generation and function. First, they used linkage analysis followed by candidate gene sequencing to identify a second defective gene in the mevalonate pathway that causes this disorder, mevalonate decarboxylase (MVD). They then sequenced 12 genes in the mevalonate and isoprenoid pathways in 134 index cases and provide evidence that mutations inactivating 3 more genes in the pathway (PMVK, MVD and FDPS) also cause this disease. Collectively, implicated variants in the 4 genes identified explain disease in 98% of familial cases and 73% of apparently sporadic cases. They also provide evidence for correlation between the underlying gene and selected aspects of clinical presentation.
The three reviewers agree that the paper is clearly written and the major conclusions of the paper are supported by the data presented. Although this represents an important advance in the molecular characterization of this disorder, and further implicates and refines the disrupted pathway, the paper does not provide sufficient mechanistic insights to warrant publication in eLife and without such data is better suited for a subspecialty journal. How do defects in mevalonate metabolism cause the keratinization defects – toxic intermediary metabolites vs. deficiency of metabolic end-products? Most individuals were heterozygous for a mutation. Since the genes all encode enzymes, is there a mutation in the normal allele, or possibly in another allele in the pathway, in the skin lesions? Providing data to address any one of these questions would suffice.
Other major concerns:
1) While the authors describe that complete co-segregation between each putative mutation and disease was observed (and representative pedigrees are shown), these data should be expanded. While it might not be plausible to show pedigrees for all 60 familial cases with a defined mutation, it would be worth noting the size distribution of families for each gene and the number of families that achieved an independent LOD score >2.0 or >3.0, for example.
2) Curiously, it is stated that the vast majority of cases associated with variation in PMVK or FDPS were sporadic. If familial segregation was ever seen for these 2 genes, it should be shown. Genetic evidence that at least some of the sporadic mutations are de novo should be provided. It would also be informative to perform a formal burden test for these 2 genes to document enrichment for novel/rare and predicted deleterious variation in association with phenotype.
3) There is no evidence of pathogenicity for a number of the missense variants listed in Table 1 aside from the fact that they are present in cases and not controls (or 1 control, in the case of variant 44, MVD c.746T>C in Table 1). Providing information regarding the presence of these mutations in the public databases would be informative. A formal burden test for each new gene should be performed.
4) While the authors acknowledge that MVK mutations have previously been observed in other disorders, including mevalonic aciduria (MA) and hyper IgD syndrome (HIGDS), the same appears not to be true for the other 3 genes. Curiously, patients with biallelic mutations causing MA or HIGDS do not show porokeratosis, nor do their heterozygous parents. This is particularly puzzling because this group has identified porokeratosis in individuals heterozygous for some of the same mutations. Variation in environmental exposure (e.g. UV light) has been invoked as a potential mechanism, however this seems less tenable given the complete penetrance seen for porokeratosis in the described families. Somatic loss of the other MVK allele in skin lesions was previously excluded in a handful of patients, but this should be revisited in this manuscript.https://doi.org/10.7554/eLife.06322.020
- Zhenghua Zhang
- Zhenghua Zhang
- Zhenghua Zhang
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
We thank all patients and family members who participated in this study. This work is sponsored by grants from the Shanghai Pujiang Program, the Research Special Fund for Public Welfare Industry of Health, Guideline-oriented Research in the Management of Some Common and Severe Skin Diseases (No. 201002016), Construction of National Key Clinical Specialty in Shanghai, Ministry of Health. Special thanks should be given to Dr Lihua Zheng for helpful discussion.
Human subjects: All procedures followed the guidelines of the Helsinki Declaration and were approved by the ethics committee and by the Scientific Ethical Committee of Fudan University. Participating centers provided local Institutional Review Board approval for genetic analysis. Study participants provided informed consent for genetic testing.
- Helen H Hobbs, Reviewing Editor, University of Texas Southwestern Medical Center, United States
- Received: January 4, 2015
- Accepted: June 26, 2015
- Version of Record published: July 23, 2015 (version 1)
© 2015, Zhang et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.