Skip to main content

Expression profiling of white sponge nevus by RNA sequencing revealed pathological pathways

An Erratum to this article was published on 23 September 2015

Abstract

Background

White sponge nevus (WSN) is a rare periodontal hereditary disease. To date, almost all WSN studies have focused on case reports or mutation reports. Thus, the mechanism behind WSN is still unclear. We investigated the pathogenesis of WSN using expression profiling.

Methods

Sequence analysis of samples from a WSN Chinese family revealed a mutation (332 T > C) in the KRT13 gene that resulted in the amino acid change Leu111Pro. The pathological pathway behind the WSN expression profile was investigated by RNA sequencing (RNA-seq).

Results

Construction of a heatmap revealed 24 activated genes and 57 reduced genes in the WSN patients. The ribosome structure was damaged in the WSN patients. Moreover, the translation rate was limited in the WSN patients, whereas ubiquitin-mediated proteolysis was enhanced.

Conclusions

Our results suggest that the abnormal degradation of the KRT13 protein in WSN patients may be associated with keratin 7 (KRT7) and an abnormal ubiquitination process.

Background

White sponge nevus (WSN) is a rare periodontal hereditary disease that was first described by Hyde [1] and coined by Cannon [2]. It is characterized by white, thickened, folded and spongy lesions of the oral mucosa, although the esophageal, laryngeal, nasal and anogenital mucosa might also be affected [3]. Recently, KRT4 [4] and KRT13 [5] gene mutations were shown be the underlying cause of WSN. WSN is an autosomal dominant genetic disease in the oral mucosa. Although WSN patients do not experience significant physical pain, they often complain of an altered texture of the mucosa or changes in their physical appearance induced by the lesions. To develop better therapeutic strategies for WSN, it is important to understand the consequences of the associated genetic mutations and molecular changes. To date, many WSN cases have been reported around the world, including China [68], Italy [911], Japan [12, 13], the UK [14, 15], Spain [16], Scotland [17], Iran [18], Brazil [19], and Turkey [20]. However, the exact pathogenic mechanism behind WSN remains unclear. With the advent of next-generation sequencing technologies, RNA-seq has become a useful tool for defining the transcriptomes of cells. Moreover, this technology may be useful for analyzing gene expression at the exon level and delineating novel splicing variants [2125]. Early applications of RNA-seq included expression profiling of yeast [26], mouse brains, liver tissues, skeletal muscle tissues [27], human embryonic kidneys, B cell lines [28] and early embryos [29]. RNA-seq offers several advantages over other expression profiling technologies, including higher sensitivity and the ability to detect splicing isoforms and somatic mutations [30]. To date, an RNA-seq analysis of WSN has not been published. Therefore, we applied RNA-seq technology to analyze WSN.

Recently, there have been many reports of new cases and mutations. However, the use of the modularity of transcriptional networks as a principle approach to understand this complex pathway has not been adequately explored. Elucidating the molecular mechanisms is a key factor for the development of successful treatments for WSN. In this study, we investigated the pathological pathway behind the WSN expression profile using RNA sequencing (RNA-seq).

This study provided a new direction for investigations into the mechanisms behind WSN, prenatal diagnosis and clinical treatment. Furthermore, our results provide a frame of reference and instructions for treating periodontal disease and other keratin-related diseases.

Methods

Ethics approval

Oral epidermis tissues were obtained from patients who were well informed of all of the purposes that the tissues might be used for in this study. This study was approved by Fudan University’s ethics committee.

Clinical report

The proband in the WSN family was a 44-year-old male Chinese patient from Hunan province who was affected by white asymptomatic oral plaques that were clinically diagnosed as WSN. In this six-generation-family from Hunan province, there were 120 members. Using pedigree analysis, we determined that the genetic modes of the disease for 28 WSN patients from this family were autosomal dominant disorders. The incidence of WSN was 23.3 % in this family. The major lesionsin these patients were white plaques of the tongue and the buccal mucosa on both sides. The diagnosis of WSN was supported by the family history and the clinical and histopathological findings.

Establishment of the cell lines

Oral epithelial cells from the normal subjects and the WSN patients were cultured in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10 % fetal bovine serum (FBS), 100 units/ml penicillin, and 100 μg/ml streptomycin at 37 °C in 5 % CO2.

RNA isolation and library construction

Total cellular RNA from the normal subjects and the WSN patients was extracted using the Trizol Reagent (Invitrogen, USA). Library construction was performed following the Illumina manufacturer’s suggestions. The libraries were sequenced on the Illumina Hiseq 2000 platform. Sequencing reads that contained polyA, low quality, and adapters were pre-filtered prior to mapping. The filtered reads were mapped on to the hg19 genome and the mm9 genome using default parameters with BWA aligner29; reads that failed to map to the genome were re-mapped to their respective mRNA sequences to capture reads that spanned exons.

RNA-Seq analysis and pathogenic pathway analysis

We generated signal networks to identify and visualize the hub genes. We expected that module genes would have significant positive module membership values. Functional annotation was performed with the Database for Annotation, Visualization and Integrated Discovery (DAVID) Bioinformatics Resource. We focused our bioinformatics analysis on Ingenuity Pathways Analysis (IPA), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Set Enrichment Analysis (GSEA). RNA-seq represented an advanced method to investigate disease pathogenesis.

Results

Clinical report and establishment of the cell line

During the oral clinical examination, the WSN patient had white lesions located bilaterally on the lips, the lateral margin of the tongue and the bilateral buccal mucosa (Fig. 1a). The arrow in Fig. 1b indicates the position of the 332 T > C mutation identifiedin this family, this mutation predicts the amino acid change Leu111Pro. To investigate the pathogenic mechanism behind WSN, we collected tissues from the oral mucosa of WSN patients and normal controls from the same family for cell culture. The cells of the two cell lines attached and showed typical morphology for oral mucosa cells. Then, we established a stable mesenchymal stem cell line (MSC) from the subjects’ gum tissues and analyzed the cells by FACS (Fig. 1c).

Fig. 1
figure 1

Mutation analysis and establishment of cell lines. a: White spongy oral plaques in the buccal mucosa of the normal subjects and the WSN patients; b: Partial DNA sequences of exon 1A of the KRT13 gene from a Chinese family. The arrow indicates the position of the mutation 332 T > C. This mutation predicts the amino acid change L111P in the KRT13 polypeptide from the WSN patient; c: FACS analysis of human mesenchymal stem cells (MSC) isolated from gum tissues

RNA-Seq analysis and pathogenic pathway analysis

The scatterplot in Fig. 2a depicts the number of activated (red) and reduced (green) genes in the normal subjects compared to the WSN subjects. The heatmap demonstrated the presence of 24 activated genes and 57 reduced genes in the WSN patient relative to the control (Fig. 2b). Using the IPA software package, we identified 19 significant bio-function terms and 10 significant canonical pathways through KEGG analysis (Fig. 3a and b) (Additional files 1 and 2).The strongest enriched gene ontology terms found in the GSEA/GO analysis indicated that genes involved in the structural constituents of ribosomes are expressed at reduced levels in WSN patients; this finding is in agreement with the observation that the structure of the ribosome is damaged in WSN patients (Fig. 4a). The snapshot of the KEGG enrichment pathways showed that protein degradation levels and ubiquitin-mediated proteolysis were enhanced in the WSN patients (Fig. 4b). Next, we chose the top 20 upregulated classes between the normal subjects and the WSN patients. There were significant differences between the size, enrichment score (ES) and NES (Tables 1 and 2) between the 2 groups. The WSN enrichment area was mainly focused on monooxygenase activity, the JAK Stat cascade, and extracellular regions. In contrast, enrichment areas from the normal patient included the detection of stimuli and G-protein signaling.

Fig. 2
figure 2

RNA-seq analysis. a: The distribution of differentially expressed genes between the normal subjects and the WSN patients; b: Scatterplot showing the number of activated (red) and reduced (green) genes in the normal subjects compared to the WSN patients; c: Heatmap showing the relative expression of activated genes in the WSN patients (n = 24)

Fig. 3
figure 3

Functional analysis and pathway analysis. a: Functional annotation of differentially expressed genes; b: Pathway analysis

Fig. 4
figure 4

Gene Ontology and KEGG Enrichment analysis. a: The strongest enriched Gene Ontology terms in the GSEA/GO analysis; b: The snapshot of KEGG enrichment

Table 1 GSEA analysis of a normal subject
Table 2 GSEA analysis of a WSN patient

We used signaling networks to identify and visualize the hub genes. We expected the module genes to have significant positive module membership values. Then, we constructed a network using published data sets and generated an independent list of hub genes.

We compared the expression levels of different KRT proteins from the normal subjects and the WSN patients. The KRT7 expression level was lower in the WSN patients; however, there was almost no change in the KRT13 levels (Fig. 5a).The constructed pathway was in agreement with the module visualization of network connections and associated functions. The function of KRT13 was dependent on KRT7 and UBC (Fig. 5b). There was almost no changein the UBC expression levels between the normal subjects and the WSN samples (Fig. 5c). Thus, the abnormal degradation of the KRT13 protein in the WSN patients may be associated with an abnormal ubiquitination process.

Fig. 5
figure 5

Pathogenic mechanism analysis. a: The KRT expression levels between the normal subjects and the WSN patients; b: Module visualization of network connections and associated functions. Bioinformatics analysis of target genes and network analysis of these genes using the String 8.3 software indicated the central involvement of KRT 13 signaling; c: The UBC expression levels between the normal subjects and the WSN patients

Discussion

With the development of sequencing technology, there are two important bioinformatic gene expression profiling methods: Microarray and RNA-Seq. In oral research, microarrays have diverse applications in oral cancers, including early diagnosis of the transformation of premalignant lesions, identification of malignancy in tissue biopsies, drug discovery, identification of biomarkers, and subclassification of histologically identified tumors [31]. In premalignant lesions, such as leukoplakias and erythroplakias, microarrays have been used to identify genes that could serve as biomarkers for dysplastic lesions with the potential to progress to cancer [32]. For clinical applications, microarrays have been used for a longer period of time and will probably have regulatory approvals for diagnostic use prior to RNA-Seq obtaining approvals. RNA-Seq will eventually be used more routinely than microarray, but right now the techniques can be complementary to each other. Microarrays will not become obsolete but might be relegated to only a few uses. RNA-Seq clearly has a bright future in bioinformatic data collection [33].

WSN is a rare oral hereditary disease. Gene expression profiling of WSN patientshas not elucidated the mechanisms behind this disease. In this study, we analyzed the expression differences at the gene level between normal subjects and WSN patients and found 81 differentially expressed genes. These genes was divided into 8 categories according to their gene function and 10 canonical pathways, which provided important clues for understanding the molecular mechanisms behind WSN pathogenesis.

In summary, we demonstrated that the use of RNA-seq markedly improved the transcriptome quantification associated with WSN. We expect that RNA-seq will also be useful for quantitatively delineating the structures, isoforms, and specific expression patterns of both coding genes and non-coding regulatory RNAs. Furthermore, RNA-seq analysis has the advantage of providing quantitative results, which is in contrastto exome or genomic sequencing.

To the best of our knowledge, WSN is an autosomal dominant genetic disease. Therefore, there is no appropriate animal model or no effective therapeutic drugs have been developed. Although WSN patients do not exhibit significant physical pain, they often complain of an altered texture of the mucosa or changes to their physical appearance induced by the lesions. To develop better therapeutic strategies for WSN, it is important to understand the consequences of the associated genetic mutations and molecular changes. Therefore, we investigated the pathogenesis and signaling pathways involved in WSN using RNA-seq. Our results suggested that the KRT13 mutation may be associate with KRT7 and UBC. UBC is a protein-coding gene that encodes a ubiquitin protein that exists either covalently attached to another protein or free (unanchored). Our previous study found that the abnormal degradation of the KRT13 protein in WSN patients contributed to an abnormal ubiquitination process [6]. Therefore, a key target for the treatment of WSN is to prevent the degradation of the KRT13 protein.

Our study utilized a next generation sequencing platform to comprehensively characterize the KRT13-related WSN transcriptome for the first time. It provided the basis for an understanding of the molecular mechanisms behind WSN pathogenesis at a system-wide level. Future research based on our findings may speed up the discovery of novel biomarkers and drug targets that can be used to improve the diagnosis and therapy of WSN. The results of our RNA-Seq analysis suggested that the abnormal degradation of the KRT13 protein in WSN patients may be associated with keratin 7 (KRT7) and an abnormal ubiquitination process. The structure of the ribosome was found to be damaged in WSN patients. Moreover, the translation rate was reduced in WSN patients, whereas ubiquitin-mediated proteolysis was enhanced. Therefore, the development of a valuable drug to reduce the degradation of KRT13is crucial for WSN patients. We are optimistic that the problem will be solved with future studies on WSN. We also expect the development of RNA-Seq to enable applications involved in determining the structural dynamics of the transcriptome and the pathogenic mechanism of disease.

Although WSN patients had no physical pain significantly, but they often complained of an altered texture of the mucosa or the bad looking created by the lesions. Many WSN patients made therapy treatments medication with nystatin, antihistamines, vitamins and mouth rinses. Azithromycin, tetracycline [5], chlorhexidine [34], Victoria A acid [35] and penicillin had succeed in the clinical progress. However, there is no standard treatment protocol for WSN till now. We applied RNA-seq technology to explor the WSN mechanism. The Human induced pluripotent stem cells (iPSCs) represent an excellent tool for many clinical trials [36]. The RNA-seq technology and induced pluripotent stem cells (iPSCs) technology may be a point way to the treatment of the rare disease. Meanwhile the WSN patient should perform a careful oral hygiene to reduce infection in the oral cavity. To lead to the proper diagnosis and treatment of this rare disease, it is great importance of collaboration between anamnesis, clinical examination and pathologic findings. With the further research of WSN, there is optimism that the problem will be solved in the next decade years.

Conclusions

The genetic disease WSN occurs infrequently. Experts have struggled to grasp factors contributing to its clinical symptoms, gene mutations and treatment; however, the mechanism behind the disease is still unclear. Here, we report the pathogenic mechanism of WSN. Our results suggest that the abnormal degradation of the KRT13 protein in WSN patients may be associated with KRT7 and an abnormal ubiquitination process. This finding may contribute to the development of a molecular therapy for WSN. Gene-based diagnosis and therapy for WSN patients may become available in the near future and may provide references and instructions for treating other keratin-associated diseases. This finding will hopefully improve the levels of prenatal diagnoses and treatment of rare diseases.

Abbreviations

KRT:

Keratin

KRTs:

Keratins

KRT4:

Keratin 4

KRT13:

Keratin 13

KRT7:

Keratin 7

RNA-seq:

RNA sequencing

WSN:

White sponge nevus

DAVID:

Database for Annotation Visualization and Integrated Discovery

IPA:

Ingenuity Pathways Analysis

KEGG:

Kyoto Encyclopedia of Genes and Genomes

GSEA:

Gene Set Enrichment Analysis

iPSCs:

induced Pluripotent Stem Cells.

References

  1. Hyde J. An unusual naevus of the tongue in a five-year-old boy. J Cutan Dis. 1909;27:256.

    Google Scholar 

  2. Cannon A. White sponge nevus of the mucosa (naevusspongiosusalbus mucosae). Acta DermVenereol. 1935;31:365.

    Google Scholar 

  3. Jorgenson RJ, Levin LS. White sponge nevus. Arc Dermatol. 1981;117:73.

    Article  CAS  Google Scholar 

  4. Rugg E, McLean W, Allison W, Lunny D, Macleod R, Felix D, et al. A mutation in the mucosal keratin K4 is associated with oral white sponge nevus. Nat Genet. 1995;11:450–2.

    Article  CAS  PubMed  Google Scholar 

  5. Richard G, De Laurenzi V, Didona B, Bale SJ, Compton JG. Keratin 13 point mutation underlies the hereditary mucosal epithelia disorder white sponge nevus. Nat Genet. 1995;11:453–5.

    Article  CAS  PubMed  Google Scholar 

  6. Cai W, Chen Z, Jiang B, Yu F, Xu P, Wang M, et al. Keratin 13 mutations associated with oral white sponge nevus in two Chinese families. Meta Gene. 2014;2:374–83.

    Article  PubMed Central  PubMed  Google Scholar 

  7. Zhang J, Yang Z, Chen R, Gao P, Zhang Y, Zhang L. Two new mutations in the keratin 4 gene causing oral white sponge nevus in Chinese family. Oral Dis. 2009;15:100–5.

    Article  CAS  PubMed  Google Scholar 

  8. Chao S, Tsai Y, Yang M, Lee J. A novel mutation in the keratin 4 gene causing white sponge naevus. Br J Dermatol. 2003;148:1125–8.

    Article  CAS  PubMed  Google Scholar 

  9. Marrelli M, Tatullo M, Dipalma G, Inchingolo F. Oral infection by Staphylococcus aureus in patients affected by White Sponge Nevus: a description of two cases occurred in the same family. Int J Med Sci. 2012;9:47.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Terrinoni A, Candi E, Oddi S, Gobello T, Camaione DB, Mazzanti C, et al. A glutamine insertion in the 1A alpha helical domain of the keratin 4 gene in a familial case of white sponge nevus. J Invest Dermatol. 2000;114:388–91.

    Article  CAS  PubMed  Google Scholar 

  11. Lucchese A, Favia G. White sponge naevus with minimal clinical and histological changes: report of three cases. J Oral Pathol Med. 2006;35:317–9.

    Article  PubMed  Google Scholar 

  12. Kimura M, Nagao T, Machida J, Warnakulasuriya S. Mutation of keratin 4 gene causing white sponge nevus in a Japanese family. Int J Oral Maxillofac Surg. 2013;42:615–8.

    Article  CAS  PubMed  Google Scholar 

  13. Shimizu A, Yokoyama Y, Shimomura Y, Ishikawa O. White sponge nevus caused by a missense mutation in the keratin 4 gene. Eur J Dermatol. 2012;22:571–2.

    PubMed  Google Scholar 

  14. McDonagh A, Gawkrodger D, Walker A. White sponge nevus successfully treated with topical tetracycline. Clin Exp Dermatol. 1990;15:152–3.

    Article  CAS  PubMed  Google Scholar 

  15. Cutlan J, Saunders N, Olsen S, Fullen D. White sponge nevus presenting as genital lesions in a 28-year-old female. J Cutan Pathol. 2010;37:386–9.

    Article  CAS  PubMed  Google Scholar 

  16. Jornet PL. White sponge nevus: presentation of a new family. Pediatr Dermatol. 2008;25:116–7.

    Article  Google Scholar 

  17. Terrinoni A, Rugg E, Lane E, Melino G, Felix D, Munro C, et al. A novel mutation in the keratin 13 gene causing oral white sponge nevus. J Dent Res. 2001;80:919–23.

    Article  CAS  PubMed  Google Scholar 

  18. Aghbali A, Pouralibaba F, Eslami H, Pakdel F, Jamali Z. White sponge nevus: a case report. J Dent Res Dent Clin Dent Prospects. 2009;3:70.

    PubMed Central  PubMed  Google Scholar 

  19. Martins FP, Brasileiro B, Piva M, Trento C, Santos TS. Familial case of oral white sponge nevus: a rare hereditary condition. An Bras Dermatol. 2011;86:39–41.

    Article  Google Scholar 

  20. Songu M, Adibelli H, Diniz G. White Sponge Nevus: Clinical Suspicion and Diagnosis. Pediatric Dermatol. 2012;29:495–7.

    Article  Google Scholar 

  21. Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Mardis ER. The impact of next-generation sequencing technology on genetics. Trends Genet. 2008;24:133–41.

    Article  CAS  PubMed  Google Scholar 

  23. Marioni JC, Mason CE, Mane SM, Stephens M, Gilad Y. RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. Genome Res. 2008;18:1509–17.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Pawitan Y, Michiels S, Koscielny S, Gusnanto A, Ploner A. False discovery rate, sensitivity and sample size for microarray studies. Bioinformatics. 2005;21:3017–24.

    Article  CAS  PubMed  Google Scholar 

  25. Iizuka N, Oka M, Yamada-Okabe H, Mori N, Tamesa T, Okada T, et al. Comparison of gene expression profiles between hepatitis B virus-and hepatitis C virus-infected hepatocellular carcinoma by oligonucleotide microarray data on the basis of a supervised learning method. Cancer Res. 2002;62:3939–44.

    CAS  PubMed  Google Scholar 

  26. Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, et al. The transcriptional landscape of the yeast genome defined by RNA sequencing. Science. 2008;320:1344–9.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5:621–8.

    Article  CAS  PubMed  Google Scholar 

  28. Sultan M, Schulz MH, Richard H, Magen A, Klingenhoff A, Scherf M, et al. A global view of gene activity and alternative splicing by deep sequencing of the human transcriptome. Science. 2008;321:956–60.

    Article  CAS  PubMed  Google Scholar 

  29. Xue Z, Huang K, Cai C, Cai L, Jiang C-y, Feng Y, et al. Genetic programs in human and mouse early embryos revealed by single-cell RNA sequencing. Nature. 2013;500:593–7.

    Article  CAS  PubMed  Google Scholar 

  30. Berger MF, Levin JZ, Vijayendran K, Sivachenko A, Adiconis X, Maguire J, et al. Integrative analysis of the melanoma transcriptome. Genome Res. 2010;20:413–27.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Todd R, Wong DT. DNA hybridization arrays for gene expression analysis of human oral cancer. J Dent Res. 2002;81:89–97.

    Article  CAS  PubMed  Google Scholar 

  32. Carinci F, Lo Muzio L, Piattelli A, Rubini C, Palmieri A, Stabellini G, et al. Genetic portrait of mild and severe lingual dysplasia. Oral Oncol. 2005;41:365–74.

    Article  CAS  PubMed  Google Scholar 

  33. Mantione KJ, Kream RM, Kuzelova H, Ptacek R, Raboch J, Samuel JM, et al. Comparing Bioinformatic Gene Expression Profiling Methods: Microarray and RNA-Seq. Medical Science Monitor Basic Research. 2014;20:138–42.

    Article  PubMed Central  PubMed  Google Scholar 

  34. Moll R, Divo M, Langbein L. The human keratins: biology and pathology. Histochem Cell Biol. 2008;129:705–33.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. McGinnis JP, Turner JE. Ultrastructure of the white sponge nevus. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 1975;40:644–51.

    Article  Google Scholar 

  36. Wang S, Qu XB, Zhao CH. Clinical applications of mesenchymal stem cells. J Hematol Oncol. 2012;5:19.

    Article  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgements

This work was supported, in whole or in part by the National Basic Research Program (973 Program, No. 2013CB967501), the Natural Science Foundation of Shanghai (No. 12ZR1434200), the Fundamental Research Funds for the Central Universities (No. 20120072110016) and the Key Program of Shanghai (No. 074119614).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Shangfeng Liu, Shouliang Zhao or Xiaoping Wang.

Additional information

Competing interests

The authors declare that they have no competing interests and sources of funding.

Authors’ contributions

WPC contributed to the content and drafted the manuscript. TNF, BZJ, JFX and ZHC analyzed and interpreted the data. WPC, RBW, and JJL designed and prepared the figures. SFL, XPW, SLZ and JHY were responsible for critical revision of the content. All authors read and approved the final manuscript.

Wenping Cai, Beizhan Jiang and Tienan Feng contributed equally to this work.

An erratum to this article is available at http://dx.doi.org/10.1186/s13023-015-0308-8.

Additional files

Additional file 1: Table S1.

Functional annotation of different expression genes.

Additional file 2: Table S2.

Pathway analysis.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Cai, W., Jiang, B., Feng, T. et al. Expression profiling of white sponge nevus by RNA sequencing revealed pathological pathways. Orphanet J Rare Dis 10, 72 (2015). https://doi.org/10.1186/s13023-015-0285-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13023-015-0285-y

Keywords