Articles

  • the Korean Society for Brain and Neural Sciences

Article

Original Article

Exp Neurobiol 2014; 23(2): 169-172

Published online June 30, 2014

https://doi.org/10.5607/en.2014.23.2.169

© The Korean Society for Brain and Neural Sciences

Dynamic Transcriptional Events in Distal Sural Nerve Revealed by Transcriptome Analysis

Young Bin Hong1#, Sung-Chul Jung2#, Jinho Lee1, Heui-Soo Moon3, Ki Wha Chung4* and Byung-Ok Choi1*

1Department of Neurology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Seoul 135-710, 2Department of Biochemistry, Ewha Womans University School of Medicine, Seoul 158-710, 3Department of Neurology, Kangbuk Samsung Hospital, Sungkyunkwan University School of Medicine, Seoul 110-746, 4Department of Biological Science, Kongju National University, Gongju 314-701, Korea

Correspondence to: *To whom correspondence should be addressed.
Byung-Ok Choi
TEL: 82-2-3410-1296, FAX: 82-2-3410-0052
e-mail: bochoi77@hanmail.net
Ki Wha Chung
TEL: 82-41-850-8506, FAX: 82-41-850-0957
e-mail: kwchung@kongju.ac.kr
#Those authors contributed equally to this study

Received: March 1, 2014; Revised: April 6, 2014; Accepted: May 1, 2014

Compared with biochemical information available about the diseases in the central nervous system, that for peripheral neuropathy is quite limited primarily due to the difficulties in obtaining samples. Characterization of the core pathology is a prerequisite to the development of personalized medicine for genetically heterogeneous diseases, such as hereditary motor and sensory neuropathy (HMSN). Here, we first documented the transcriptome profile of distal sural nerve obtained from HMSN patients. RNA-seq analysis revealed that over 12,000 genes are expressed in distal sural nerve. Among them 4,000 transcripts are novel and 10 fusion genes per sample were observed. Comparing dataset from whole exome sequencing revealed that over 1,500 transcriptional base modifications occur during transcription. These data implicate that dynamic alterations are generated when genetic information are transitioned in distal sural nerve. Although, we could not find significant alterations associated with HMSN, these data might provide crucial information about the pathophysiology of HMSN. Therefore, next step in the development of therapeutic strategy for HMSN might be unveiling biochemical and biophysical abnormalities derived from those potent variation.

Keywords: hereditary motor and sensory neuropathy (HMSN), peripheral nervous system (PNS), distal sural nerve, transcriptome

Recent explosive advances in multi-omics have provided us with powerful tools to investigate the pathophysiological complexity of disease processes. In Mendelian diseases, the application of a single-omic approach, such as whole exome sequencing (WES), has dramatically reduced the time and effort required to identify causative mutations [1]. Thus, the next step, in which multi-omics are effectively applied to Mendelian diseases, might reveal biochemical and biophysical abnormalities that present therapeutic targets for personalized medicine.

Compared with the biochemical information available about the central nervous system (CNS), that available for the peripheral nervous system (PNS) is quite limited, thereby therapeutic approaches to the PNS diseases more difficult. Characterization of the core pathology, in conjunction with the identifications of genetic causes, is a prerequisite to the development of personalized medicine for genetically and clinically heterogeneous disorder, such as hereditary motor and sensory neuropathy (HMSN) [2].

HMSN is frequently classified as type 1, demyelinating form (HMSN1), or type 2, axonal form (HMSN2). HMSN1 exhibits markedly reduced nerve conduction velocities (NCV), whereas HMSN2 shows slightly reduced or even normal NCVs. HMSN2 is occasionally difficult to distinguish phenotypically from distal hereditary motor neuropathy (dHMN) which is closely resembles HMSN but without sensory abnormality, because sensory signs are often lacking in HMSN2 patients [3].

Recent monumental attempts to decipher the transcriptome in the human brain have undoubtedly contributed to our anatomical understanding of the transcriptional network and disease-related changes in inaccessible tissues [4, 5]. Subsequent proteomic or metabolomic analyses will provide further crucial information in the near future. Therefore, documentation of the transcriptome in the human distal sural nerve might provide primary knowledge as a basis for understanding the pathophysiology of PNS diseases. Here, we performed the transcriptome analysis using distal sural nerve of HMSN patients.

Patients

This study included a total of 14 individuals of Korean origin with HMSN or normal phenotype. Written informed consent was obtained from all the participants and from the parents of participants younger than 18 years, according to the protocol approved by the Institutional Review Board of Ewha Womans University, Mokdong Hospital (ECT 11-58-37).

Clinical assessments

The patients were examined for motor and sensory impairment, deep tendon reflexes, and muscle atrophy by two independent neurologists. The muscle strengths of the flexor and extensor were assessed manually using the Medical Research Council scale. Physical disabilities were expressed on two different scales: the Functional Disability Scale (FDS) [6] and Charcot-Marie-Tooth disease neuropathy score (CMTNS) [7] (see Table 1). Disease severity was determined using a nine-point FDS, which was based on the following criteria: 0, normal; 1, normal but with cramps and fatigability; 2, an inability to run; 3, walking difficulty but still possible unaided; 4, walking with a cane; 5, walking with crutches; 6, walking with a walker; 7, wheelchair bound; and 8, bedridden. In addition, we determined the CMTNS, and the patient was divided into mild (CMTNS≤10), moderate (CMTNS, 11-20), and severe (CMTNS≥21) categories.

Sural nerve biopsy

Sural nerves from 14 individuals were obtained by incision under sterile conditions and local anesthesia [8]. After obtaining about 10 mm length of nerve, sharp dissection followed to obtain pure sural nerve.

Transcriptome analysis

RNA-Seq data were generated from all 14 samples. The cDNA library was prepared using the TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA). The sequencing protocol consisted of total RNA extraction from a biopsy of the distal sural nerve and polyA RNA extraction, RNA fragmentation, random-hexamer-primed reverse transcription, and 100 nt paired-end sequencing with HiSeq 2000 (Illumina). To estimate the expression levels and to identify alternatively spliced transcripts, the RNA-Seq reads were mapped to the human genome using TopHat (version 1.3.3, http://tophat.cbcb.umd.edu/index.html). Then all variants were annotated by the ANNOVAR (version 2011Feb20) program. To determine whether the variants were previously reported, dbSNP135 (http://www.ncbi.nlm.nih.gov) was used. The transcript levels were calculated and the relative transcript abundances were measured in fragments per kilobase of exon per million fragments mapped (FPKM) using Cufflinks (http://cufflinks.cbcb.umd.edu/index.html). The deFuse program (version 0.4.3, http://sourceforge.net/apps/mediawiki/defuse) was used to identify gene fusions from the RNA-Seq data because this program considered all alignments and all possible locations for fusion boundaries [9].

Whole-exome sequencing

Total DNA was purified from peripheral blood with a QIAamp Blood DNA Purification Kit (Qiagen, Hilden, Germany). The whole exome was captured using SeqCap EZ V3.0 (Roche-NimbleGen, Madison, WI), and sequencing was performed using the HiSeq 2000 Genome Analyzer (Illumina). The UCSC assembly hg19 (NCBI build 37.1) was used as the reference sequence. Variant calling was achieved for those with ≥20 SNP quality. The variants were annotated with the ANNOVAR program.

To investigate the transcriptional profile in the human distal sural nerve which is the most common tissue for observation of pathophysiology of PNS, we enrolled 14 Korean individuals. RNA-Seq and subsequent analyses were used to generate mRNA expression profiles in which fragments per kilobase of transcript per million fragments mapped (FPKM) was used to quantify each transcript. Using a threshold of FPKM>1, which corresponds to approximately one transcript per cell [10], we found that over 12,000 transcripts are expressed in the human distal sural nerve (Table 1 and Supplementary Material) is almost equivalent to the number expressed in human lymphoblastoid cells (11,000 active transcripts), which was calculated by a similar analytical parameter [11].

RNA-Seq also identified over 4,000 unreported transcripts (FPKM > 1) per sample. As genome tiling arrays revealed over 10,000 novel transcripts in the human liver [12], the transcriptional events in the human distal sural nerve are also dynamic. We also found that around 10 transcripts per sample were fused in both the controls and patients (Table 1). However, we have no clue to the function or biochemical significance of these unreported and fused transcripts in this tissue.

Our transcriptome analysis also revealed over 5,000 functionally significant variants (nonsynonymous splicing site, stop gain/loss, and coding indels) in each sample. Because these variants might primarily derive from changes in the genomic DNA, we analyzed the genomic variants in nine samples from peripheral blood DNA using WES. WES revealed that over 9,000 functionally significant variants are present per sample. Therefore, the variants of transcriptome are fewer than those of genome. Next, we calculated transcriptional base modifications (TBMs) by comparing the two data sets. TBM analysis revealed that there are over 1,500 functionally significant variants. Among them 77% (over 1,200 variants) were unreported to dbSNP135. Although these might include many false positives, arising from methodological or analytical limitations and the discordance of data sources, these data imply that dramatic alterations occur when genetic information is actively transitioned in specific tissues. The significance of TBM is still unclear, until its correlation with disease onset, progression, susceptibility, or penetrance is reported.

Based on this expression profile, we analyzed the changes in the expression patterns of HMSN-causative genes and tried to identify the genes that were differentially expressed in the patients and the controls. Several changes in putative mechanistic/metabolic pathways were identified in the HMSN patient groups (data not shown), but the heterogeneity and size of the samples limited the statistical validity of the results. Nonetheless, future well-designed in vitro experimental approaches might convert parts of these data into valid information, which will allow us to understand the pathophysiology of the disease and may suggest therapeutic strategies.

In summary, we first generated a transcriptional profile of the human distal sural nerve, providing basic information about the gene expression in this tissue based on pathophysiological studies of HMSN. Our data also suggest that considerable numbers of events occur beyond the information encoded in the DNA, which has been the main source of biomedical information in the pre-omics era. The contribution of multi-omics might accelerate the biomedical approaches in the PNS by presenting systematic information. Subsequent studies that characterize these novel transcripts and TBMs and validate the HMSN-related changes in the expression profile will provide crucial information about the pathological mechanism of HMSN and provide therapeutic options for each genetic causative.

Supplementary Material

Supplementary questionnaire can be found via http://en-journal.org/src/sm/en-23-169-s001.xlsx

Questionnaires for Peyronie disease.

en-23-169-s001.xlsx
Table. 1.

Table 1. Patients enrolled in this study and transcriptome analysis results

aAge at sural nerve biopsy; bFunctional disability scale [6]; cCharcot-Marie-Tooth disease neuropathy score [7]; dExpressed transcripts were determined as FPKM>1. Original data are available in Supplementary Material; eNovel transcripts determined with the Cufflinks program and FPKM>1; fGene fusion was analyzed with the deFuse software; gNonsynonymous splicing site, stop gain/loss, and coding indels; hVariants unreported in dbSNP135; iTranscriptional base modifications; jNumber of transcriptome variants that were not identified in the whole-exome analysis; kNonsynonymous splicing site, stop gain/loss, and coding indels that are not reported in dbSNP135; lcausative gene was not defined.


  1. Ng SB, Buckingham KJ, Lee C, Bigham AW, Tabor HK, Dent KM, Huff CD, Shannon PT, Jabs EW, Nickerson DA, Shendure J, Bamshad MJ. Exome sequencing identifies the cause of a mendelian disorder. Nat Genet 2010;42:30-35.
    Pubmed
  2. Park J, Hyun YS, Kim YJ, Nam SH, Kim SH, Hong YB, Park JM, Chung KW, Choi BO. Exome sequencing reveals a novel PRPS1 mutation in a family with CMTX5 without optic atrophy. J Clin Neurol 2013;9:283-288.
    Pubmed
  3. Nakhro K, Park JM, Kim YJ, Yoon BR, Yoo JH, Koo H, Choi BO, Chung KW. A novel Lys141Thr mutation in small heat shock protein 22 (HSPB8) gene in Charcot-Marie-Tooth disease type 2L. Neuromuscul Disord 2013;23:656-663.
    Pubmed
  4. Hawrylycz MJ, Lein ES, Guillozet-Bongaarts AL, Shen EH, Ng L, Miller JA, van de Lagemaat LN, Smith KA, Ebbert A, Riley ZL, Abajian C, Beckmann CF, Bernard A, Bertagnolli D, Boe AF, Cartagena PM, Chakravarty MM, Chapin M, Chong J, Dalley RA, Daly BD, Dang C, Datta S, Dee N, Dolbeare TA, Faber V, Feng D, Fowler DR, Goldy J, Gregor BW, Haradon Z, Haynor DR, Hohmann JG, Horvath S, Howard RE, Jeromin A, Jochim JM, Kinnunen M, Lau C, Lazarz ET, Lee C, Lemon TA, Li L, Li Y, Morris JA, Overly CC, Parker PD, Parry SE, Reding M, Royall JJ, Schulkin J, Sequeira PA, Slaughterbeck CR, Smith SC, Sodt AJ, Sunkin SM, Swanson BE, Vawter MP, Williams D, Wohnoutka P, Zielke HR, Geschwind DH, Hof PR, Smith SM, Koch C, Grant SG, Jones AR. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature 2012;489:391-399.
    Pubmed
  5. Konopka G, Friedrich T, Davis-Turak J, Winden K, Oldham MC, Gao F, Chen L, Wang GZ, Luo R, Preuss TM, Geschwind DH. Human-specific transcriptional networks in the brain. Neuron 2012;75:601-617.
    Pubmed
  6. Birouk N, LeGuern E, Maisonobe T, Rouger H, Gouider R, Tardieu S, Gugenheim M, Routon MC, Léger JM, Agid Y, Brice A, Bouche P. X-linked Charcot-Marie-Tooth disease with connexin 32 mutations: clinical and electrophysiologic study. Neurology 1998;50:1074-1082.
    Pubmed
  7. Shy ME, Blake J, Krajewski K, Fuerst DR, Laura M, Hahn AF, Li J, Lewis RA, Reilly M. Reliability and validity of the CMT neuropathy score as a measure of disability. Neurology 2005;64:1209-1214.
    Pubmed
  8. McPherson A, Hormozdiari F, Zayed A, Giuliany R, Ha G, Sun MG, Griffith M, Heravi Moussavi A, Senz J, Melnyk N, Pacheco M, Marra MA, Hirst M, Nielsen TO, Sahinalp SC, Huntsman D, Shah SP. deFuse: an algorithm for gene fusion discovery in tumor RNA-Seq data. PLoS Comput Biol 2011;7:e1001138.
    Pubmed
  9. Sommer CL, Brandner S, Dyck PJ, Harati Y, LaCroix C, Lammens M, Magy L, Mellgren SI, Morbin M, Navarro C, Powell HC, Schenone AE, Tan E, Urtizberea A, Weis J, Peripheral Nerve Society. Peripheral Nerve Society Guideline on processing and evaluation of nerve biopsies. J Peripher Nerv Syst 2010;15:164-175.
    Pubmed
  10. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods 2008;5:621-628.
    Pubmed
  11. Ju YS, Kim JI, Kim S, Hong D, Park H, Shin JY, Lee S, Lee WC, Kim S, Yu SB, Park SS, Seo SH, Yun JY, Kim HJ, Lee DS, Yavartanoo M, Kang HP, Gokcumen O, Govindaraju DR, Jung JH, Chong H, Yang KS, Kim H, Lee C, Seo JS. Extensive genomic and transcriptional diversity identified through massively parallel DNA and RNA sequencing of eighteen Korean individuals. Nat Genet 2011;43:745-752.
    Pubmed
  12. Bertone P, Stolc V, Royce TE, Rozowsky JS, Urban AE, Zhu X, Rinn JL, Tongprasit W, Samanta M, Weissman S, Gerstein M, Snyder M. Global identification of human transcribed sequences with genome tiling arrays. Science 2004;306:2242-2246.
    Pubmed
  13. Chung KW, Suh BC, Cho SY, Choi SK, Kang SH, Yoo JH, Hwang JY, Choi BO. Early-onset Charcot-Marie-Tooth patients with mitofusin 2 mutations and brain involvement. J Neurol Neurosurg Psychiatry 2010;81:1203-1206.
    Pubmed
  14. Lee SS, Lee HJ, Park JM, Hong YB, Park KD, Yoo JH, Koo H, Jung SC, Park HS, Lee JH, Lee MG, Hyun YS, Nakhro K, Chung KW, Choi BO. Proximal dominant hereditary motor and sensory neuropathy with proximal dominance association with mutation in the TRK-fused gene. JAMA Neurol 2013;70:607-615.
    Pubmed
  15. Choi BO, Park MH, Chung KW, Woo HM, Koo H, Chung HK, Choi KG, Park KD, Lee HJ, Hyun YS, Koo SK. Clinical and histopathological study of Charcot-Marie-Tooth neuropathy with a novel S90W mutation in BSCL2. Neurogenetics 2013;14:35-42.
    Pubmed