ALL Metrics
-
Views
Get PDF
Get XML
Cite
Export
Track
Method Article

Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours

[version 1; peer review: 3 approved, 1 approved with reservations]
PUBLISHED 23 Apr 2018
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

Abstract

Transcriptomic profiling of skin disease using next generation sequencing allows for detailed information on aspects of RNA biology including gene expression, non-coding regulatory elements and gene splicing. The application of RNA sequencing to human skin disease and cancer is often hampered by degraded RNA. Here we describe a protocol that allows for consistently intact RNA to be extracted from snap frozen skin biopsy samples, which has been validated in a clinical trial setting.
Human skin tumour punch biopsies (n=28) ranging from 4-6mm in diameter were obtained from 14 patients with an inherited skin tumour syndrome (CYLD cutaneous syndrome) and frozen in liquid nitrogen prior to being stored at -80°C. These samples were then subject to cyrostat sectioning, allowing for histological assessment, and were homogenised using a bead-based lysis platform. RNA extraction was performed using a silica column-based system. RNA concentration was measured using fluorescent quantitation and RNA integrity assessed using microfluidic gel electrophoresis. We also processed normal skin biopsies using the same protocol (n=10). The mean RNA integrity score of the tumour and normal samples was 9.5, and the quantity of RNA obtained from the small amounts of tissue used exceeded requirements for RNA-seq library generation.
We propose that the method of RNA extraction suggested here allows for transcriptomic profiling from small pieces of human tissue without the need for PCR amplification during library preparation. This protocol could be utilised in healthy and diseased skin to improve mechanistic understanding in a range of human skin diseases.

Keywords

Skin biopsy, tumour, CYLD, RNA sequencing, cancer, clinical trial, transcriptomics.

Introduction

Skin has evolved as a protective barrier against environmental stresses such as physical damage, ultraviolet radiation exposure and pathogenic infections. It is a complex tissue, consisting of many cell types, including keratinocytes, immune cells, melanocytes, fibroblasts, adipocytes, nerve fibres, smooth muscle and endothelial cells (Cole et al., 2001). Skin diseases may involve any of these constituent cells, and transcriptomic profiling of diseased human skin can inform understanding of disease pathogenesis. The accessible nature of the skin, unlike tissues such as the brain or liver, and the low patient impact of skin biopsy, make the use of transcriptomics in human skin particularly attractive. Whole punch biopsies have been studied using transcriptome-wide approaches, giving insights into the inflammatory skin disease psoriasis (Zolotarenko et al., 2016).

Reliably obtaining high quality RNA of sufficient quantity from small samples typical of those derived from punch biopsies is, however, problematic. The reasons for this are diverse and include suboptimal tissue disruption and endogenous RNase activity, resulting in skin being a tissue that is amongst the most challenging to obtain intact RNA from (Walker et al., 2016)(Kaufmann et al., 1980). RNA is also a highly unstable molecule. Once extracted from cells or tissues it has a very short half-life and is easily degraded (Brooks, 1998; Tan & Yiap, 2009). RNA is chemically unstable due to the presence of a hydroxyl group at the 2’ and thus is highly susceptible to hydrolysis by ribonuclease enzymes. The hydrolysis or cleavage of RNA can also occur spontaneously, without the presence of a catalyst or enzyme (Elliott & Ladomery, 2011). Together, these challenges may discourage the inclusion of transcriptomics in clinical trials involving human skin, where such information may offer mechanistic insights.

Methods of RNA extraction from human skin have been evolving. Early approaches required large amounts of tissue, and were not amenable to medium/high throughput lysis, as each sample had to be disrupted in liquid nitrogen (Hipfel et al., 1998; Kaufmann et al., 1980). Whole punch biopsies of the skin have been successfully disrupted by low throughput devices, such as the Kinematica Polytron 1300D homogenizer and FastPrep120 bead beater (Berglund et al., 2007). 4 mm skin biopsies have also been subjected to ammonium thiocyanate-induced dermo-epidermal separation and subsequently homogenized using bead-based lysis (Clemmensen et al., 2009); however the RNA integrity may have been suboptimal due to tissue processing. Importantly none of these approaches allow for histological assessment of the skin sample, which is relevant in interpreting sequencing results of diseased skin. Laser capture microdissection is an alternative method with histological information, with the limitations being cost of laser capture microscopy equipment, small amount of RNA obtained, and RINs that may be lower (circa 7.0) (Butler et al., 2016). An additional amplification stage is also often required when starting from small amounts of RNA during library generation, which may introduce an element of bias (Jackson et al., 2014).

Here we report an efficient method of RNA extraction from skin tumours obtained using punch biopsies under local anaesthetic in a clinical trial involving patients with an inherited skin tumour syndrome. We also validated this method in normal human skin biopsies. Downstream analysis demonstrated high yield and RNA integrity, allowing for transcriptomic profiling using RNA sequencing.

Methods

Ethics

Regulatory approvals for the present study were sought and obtained from an ethics review committee (National Research Ethics Service Committee North East-Tyne and Wear Ref:14/NE//080 and 06/1059) and the Medicines Health Regulatory authority (EudraCT: 2014-001342-21). Patients provided written, informed consent for biopsy and use of their tissue samples for research. This study used samples collected as part of a clinical trial, registered at ISCRTN 75715723. The study was performed in accordance with the Declaration of Helsinki.

Sample collection

Tumours scheduled for biopsy as part of the trial protocol were identified using tumour maps and photographs. Tumour identifiers were labelled onto cryovials in advance of the procedure. Skin biopsies were performed under local anaesthetic, and one sample was taken at a time. The biopsy was transferred to a labelled cryovial, with identifiers checked together with a research nurse. The sample was then immediately immersed in liquid nitrogen. All samples were logged in accordance with standard operating procedures and stored in a -80°C freezer within a designated Human Tissue Authority freezer that is subject to temperature monitoring.

Sample processing

A full working protocol is detailed in Supplementary File 1, and the key steps are outlined here. Standard precautions to prevent contamination with RNAses were employed. The sample was removed from the freezer and transferred on dry ice to the cryostat. The sample identifier was used to label slides and preweighed lysis tubes, which were supplied prefilled with beads (Precellys). Skin punch biopsy samples of tumours (4–6 mm diameter) were mounted on a cork piece using cutting compound (OCT), and then serially sectioned (Figure 1). Each skin biopsy was cryosectioned at two levels and ten 30-micron curls were taken from each level for RNA extraction; two were taken from each biopsy. Sections were then taken for histology and applied to a standard slide, and subject to haematoxylin and eosin staining. Following cryosectioning, material was immediately placed into the bead lysis tube, weighed again and then cold lysis buffer (RLT) containing beta-mecaptoethanol was added (Qiagen RNeasy Micro Kit No 74004). The tube was frozen on dry ice and then returned to the -80°C freezer until RNA extraction was performed.

88c74a95-7ac5-4879-a54f-ea6cef3edbca_figure1.gif

Figure 1. A method for extracting high quality RNA from small skin biopsy samples which is suitable for RNA sequencing.

(a) Diagram indicating key steps in the work flow, and the need to keep samples cold throughout. (A) Human skin punch biopsy (B) Freezing of sample within 30 seconds. (C–D) Cryostat sectioning of punch biopsy and curls obtained. (E–F) Addition of lysis buffer to single use bead tube, followed by homogenisation in the bead-based lysis machine where up to 24 samples can be processed at once. (G) Silica spin column based RNA extraction performed at 4 degrees. (b) Microfluidic gel analysis of total RNA demonstrates distinct 18S and 28S ribosomal bands, consistent with the high RNA integrity scores (c) demonstrated across samples. (d) FASTQc assessment of reads indicate high quality reads in libraries developed from this dataset.

RNA extraction and quality control

The protocol supplied with the silica spin column kit (Qiagen RNeasy Micro Kit No 74004, Qiagen UK) was followed with modifications, as indicated in the working protocol (Supplementary File 1). Briefly, bead tubes were taken from the -80°C freezer on dry ice. Tissue was then homogenized in a Precellys Evolution homogeniser (Bertin, France) for 20 seconds at 5500 bpm. After homogenization, the lysate was applied to the column, and then wash steps were performed. On-column DNase digestion was performed for 7 minutes at room temperature. Further washes were performed before RNA was eluted from the column with RNase-free water. All protocol steps were performed in a cold centrifuge at 4°C at 15000 rpm apart from the DNase incubation step stated above. Eluted RNA was measured using the Qubit BR assay kit (Thermofisher, UK). RNA quality was measured using a microfluidic gel electrophoresis chip (Bioanalyer RNA 6000 Nano Chip, Agilent UK). RNA integrity numbers were obtained with the software provided (2100 Expert Software: Revision B.02.09 (SR1)) with the Agilent 2100 Bioanalyser (Agilent, UK) (Schroeder et al., 2006).

Transcriptomic analyses

To obtain the cytokeratin signatures of differentially expressed genes in cylindroma and spiradenoma compared to control skin, three tumour transcript files were compared to three control skin sample files. RNA from each sample was used to generate sequencing libraries using the Illumina Truseq stranded mRNA kit. Libraries were sequenced using an Illumina Hiseq 2500, giving 45 million paired end reads per sample which were 100 bp in length. FASTQ files were checked for quality using FASTQc and aligned using the splice aware aligner program STAR (v. 2.5.2b) to generate alignment files (Dobin et al., 2013). The read counts for each sample file were obtained using the R package Subread (v.1.28.1) (Liao et al., 2013). Differential gene expression analysis was carried out using DeSeq2 v.1.18.1 (Love et al., 2014). Cytokeratin genes that were differentially expressed with a false discovery rates of <0.05 after correction for multiple hypothesis testing are listed in Table 3.

Results

Consistent high-quality RNA is extracted from multiple skin tumours

28 skin biopsies were cryosectioned and material with an average weight of 10.4mg (range 5 –14 mg) was obtained at two levels in each biopsy Figure 1a. RNA was obtained from a total of 56 levels, with total yield exceeding the requirements for library preparation for RNA sequencing (yields typically >500ng) in the majority of samples. Sections taken for histology were stained using haematoxylin and eosin, and confirmed adjacent curls were taken from cylindroma and spiradenoma in 25 out of 28 of biopsies; three biopsies demonstrated trichoepithelioma. In 3 out of 56 levels, RNA was not obtained and this correlated with histology reflecting relatively acellular dermis beneath the level of the tumour. The range of concentrations and quality of RNA extracted from tumours are indicated in Table 1 and Figure 1b and c, indicating a mean RIN of 9.5 (range 8.5–10). Control skin sample yields and integrity (mean RIN 9.5; range 8.8–10) are indicated in Table 2.

Table 1. RNA concentrations and integrity in the 28 skin tumour biopsies studied, with 2 samples taken per biopsy.

Qualitative and quantitative measurements of the total RNA isolated from normal skin punch biopsies.

TumourSamplesConcentration (ng/ul)RNA integrity
11129.69.7
217.2-
2363.29.3
4too low/undetectable-
353269.3
62089.3
471189.3
8769
5920810
101429.8
611115.69.9
1229.28.5
71327210
145610
815182.89.6
16169.29.6
917147.49.5
1816.6-
101999.69.4
20105.69.3
1121989.5
2212.6-
122341.49.1
2482.89.7
132528210
262289.9
14273089.9
2830410
152925810
3037010
163180010
32112010
1733197.89.7
342449.6
18352869.6
361489.5
19372349.4
382509.6
20393629.7
402549.6
214125210
425510
224349.69.8
44269.5
234517.29.2
46219.8
2447258.6
48too low/undetectable-
2549360.89.7
50194.89.6
2651110.69.3
5250.69
275315.88.7
5416.69
285526.68.8
56too low/undetectable-
Average:180.89.5

Table 2. RNA concentrations and integrity of 10 normal skin samples studied.

Qualitative and quantitative measurements of the total RNA isolated from normal skin punch biopsies.

SamplesConcentration
(ng/ul)
RNA
Integrity
125.29.6
2178.8
318.69.2
42310
538.69.5
620.29.5
7329.4
836.69.4
9409.8
1021.89.7
Average:27.39.5

Table 3. Expression of known cytokeratin signatures of differentially expressed genes in cylindroma and spiradenoma compared to control skin.

Genelog2 Fold ChangeProb. FDRDescription
KRT133.941.37E-11Keratin 13
KRT17P22.440.011702576Keratin 17 pseudogene 2
KRTCAP31.050.012251062Keratinocyte associated protein 3
TCHP-0.830.022666568Trichoplein keratin filament binding
KRT15-0.890.031124672Keratin 15
KRT3-1.560.016500418Keratin 3
KRT9-1.710.002468716Keratin 9
KRT5-1.740.000610354Keratin 5
KRT8P26-1.960.012588263Keratin 8 pseudogene 26
KRT80-1.960.017781638Keratin 80
KRT19-2.070.035571416Keratin 19
KPRP-2.200.013678232Keratinocyte proline rich protein
KRTDAP-2.350.001196344Keratinocyte differentiation
associated protein
KRT78-2.390.00676354Keratin 78
KRT10-2.460.00126146Keratin 10
KRT14-2.480.00000538Keratin 14
KRT1-2.550.000588345Keratin 1
KRT77-2.670.00000744Keratin 77
KRT72-2.710.000336864Keratin 72
KRT27-2.760.005439881Keratin 27
KRT73-3.373.44E-08Keratin 73
KRT31-3.894.65E-11Keratin 31
KRT2-4.020.000000475Keratin 2

RNA sequencing validates tumour transcriptomic signature

Reads derived from next generation sequencing were subject to quality control using FASTQc (Figure 1d), and this demonstrated high quality scores, consistent with the high integrity RNA used. Differential gene expression studies focussed on known differentially expressed signature cytokeratin genes in cylindroma and spiradenoma, skin tumours seen in CYLD cutaneous syndrome (Brown et al., 2018). This confirmed that the recognised cytokeratin signature was expressed in these tumour transcriptomes (Table 3).

Discussion

Due to its high sensitivity and resolution, RNA sequencing is a powerful tool for investigation of skin disease. It can provide mechanistic insights behind disease pathogenesis, which may offer prognostic information, prediction of response to treatments and the potential for developing novel therapies (Gilmore, 2013). To overcome the inherent difficulties of obtaining high integrity RNA, there are in general, four crucial steps required during extraction: cells or tissues should be completely homogenised; nucleoprotein complexes should be disrupted; RNase should be inactivated; contamination including carbohydrate, lipids, protein and other nucleic acid should be avoided (Buckingham & Flaws, 2007; Doyle & Doyle, 1990).

We describe a novel method, incorporating cryostat sectioning, resulting in partial disruption of the tissue, whilst simultaneously allowing us to take 8 μm sections for histological assessment. We used bead-based tubes that allowed for up to 24 samples to be processed simultaneously. Importantly, we were able to work with small amounts of skin, which were robustly disrupted using the bead-based lysis system we employed. This prevented carryover of incompletely homogenised material and obstruction of spin columns used for RNA extraction. The bead-based lysis system we employed had a range of disruption settings and after optimisation with a range of bead sizes, we chose a medium bead mix (CK28 Mix), and lysis of 20 seconds at 5500 rpm in lysis buffer as this gave optimal RNA extraction and RINs.

After skin homogenisation, there are several established methods for RNA extraction. Phenol/chloroform extraction, one of the commonly used techniques, has potential for phenol contamination of the samples, which is often reported (Oñate-Sánchez & Vicente-Carbajosa, 2008). Therefore, we employed a silica column-based RNA extraction methodology (Sellin Jeffries et al., 2014). This approach suited our protocol, which required small volumes of lysis buffer, and elution of RNA performed in 14ul of water. This typically resulted in highly concentrated samples and thus satisfied the requirements of most RNA library preparation kits for next generation sequencing.

Given RNA’s sensitivity to temperature, we developed a protocol where the skin sample is kept as cold as possible throughout. The tumour samples were first snap frozen in liquid nitrogen, and then cryosectioned at a temperature of -20°C. During rapid homogenisation in our protocol, the tumour material was exposed to ambient temperatures for only a few minutes. This was immediately followed by RNA extraction in a precooled centrifuge at 4°C. Spin columns were incubated on wet ice between centrifugation steps and eluted RNA was immediately frozen after RNA extraction was complete.

In summary, we developed a method for efficient extraction of RNA from small cryosectioned pieces of human skin with histological data from adjacent tissue sections. We validated this method in human skin tumours samples in a clinical trial setting. This protocol could be utilised in healthy and diseased skin to improve mechanistic understanding in a range of human skin disease and cancer.

Data availability

FASTQ files (controls [n=3] and tumours [n=3] used to generate validation signatures), sample description files and uncropped gel images are available at Open Science Framework: http://doi.org/10.17605/OSF.IO/5YX96 (Rajan, 2018).

Data are available under the terms of the Creative Commons Zero "No rights reserved" data waiver (CC0 1.0 Public domain dedication).

Comments on this article Comments (0)

Version 1
VERSION 1 PUBLISHED 23 Apr 2018
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
Wellcome Open Research - -
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Danilenko M, Stones R and Rajan N. Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours [version 1; peer review: 3 approved, 1 approved with reservations] Wellcome Open Res 2018, 3:45 (https://doi.org/10.12688/wellcomeopenres.14360.1)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 23 Apr 2018
Views
0
Cite
Reviewer Report 20 Jun 2018
Tarl W. Prow, Future Industries Institute, University of South Australia, Adelaide, SA, Australia 
Approved
VIEWS 0
The article describes a practical and highly effective means to isolate RNA from fresh clinical skin samples for downstream transcriptome analysis.

The major issue that I see with the manuscript is the lack of text describing previous ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Prow TW. Reviewer Report For: Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours [version 1; peer review: 3 approved, 1 approved with reservations]. Wellcome Open Res 2018, 3:45 (https://doi.org/10.21956/wellcomeopenres.15623.r32973)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
0
Cite
Reviewer Report 01 Jun 2018
Verena Wally, EB House Austria, Research Program for Molecular Therapy of Genodermatoses, Department of Dermatology, University Hospital Salzburg, Paracelsus Medical University, Salzburg, Austria 
Mo Wimmer, EB House Austria, Research Program for Molecular Therapy of Genodermatoses, Department of Dermatology, University Hospital Salzburg, Paracelsus Medical University, Salzburg, Austria 
Approved
VIEWS 0
The paper "Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours" talks about a novel method to isolate high quality RNA from skin biopsies in a clinical trial ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Wally V and Wimmer M. Reviewer Report For: Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours [version 1; peer review: 3 approved, 1 approved with reservations]. Wellcome Open Res 2018, 3:45 (https://doi.org/10.21956/wellcomeopenres.15623.r33203)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
0
Cite
Reviewer Report 24 May 2018
Irene M. Leigh, Centre for Cutaneous Research, Barts and the London School of Medicine and Dentistry, Queen Mary University of London, London, UK 
Approved
VIEWS 0
The results show very good RIN figures by careful attention to detail to inhibit skin RNAases and prevent degradation of RNA etc. The transcriptomic profiling of skin disease for research purposes is well established and it would improve the article ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Leigh IM. Reviewer Report For: Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours [version 1; peer review: 3 approved, 1 approved with reservations]. Wellcome Open Res 2018, 3:45 (https://doi.org/10.21956/wellcomeopenres.15623.r33121)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
Views
0
Cite
Reviewer Report 24 Apr 2018
Veronica A. Kinsler, Pediatric Dermatology, Great Ormond Street Hospital for Children, London, UK;  Genetics and Genomic Medicine, UCL Institute of Child Health, London, UK 
Approved with Reservations
VIEWS 0
Overall this is an excellent contribution to technical knowledge of handling of skin samples, and a well-presented methods paper in so far as the optimised extraction method goes.

I have a few points to be addressed:
    ... Continue reading
    CITE
    CITE
    HOW TO CITE THIS REPORT
    Kinsler VA. Reviewer Report For: Transcriptomic profiling of human skin biopsies in the clinical trial setting: A protocol for high quality RNA extraction from skin tumours [version 1; peer review: 3 approved, 1 approved with reservations]. Wellcome Open Res 2018, 3:45 (https://doi.org/10.21956/wellcomeopenres.15623.r32972)
    NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.

    Comments on this article Comments (0)

    Version 1
    VERSION 1 PUBLISHED 23 Apr 2018
    Comment
    Alongside their report, reviewers assign a status to the article:
    Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
    Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
    Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions

    Are you a Wellcome-funded researcher?

    If you are a previous or current Wellcome grant holder, sign up for information about developments, publishing and publications from Wellcome Open Research.

    You must provide your first name
    You must provide your last name
    You must provide a valid email address
    You must provide an institution.

    Thank you!

    We'll keep you updated on any major new updates to Wellcome Open Research

    Sign In
    If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

    The email address should be the one you originally registered with F1000.

    Email address not valid, please try again

    You registered with F1000 via Google, so we cannot reset your password.

    To sign in, please click here.

    If you still need help with your Google account password, please click here.

    You registered with F1000 via Facebook, so we cannot reset your password.

    To sign in, please click here.

    If you still need help with your Facebook account password, please click here.

    Code not correct, please try again
    Email us for further assistance.
    Server error, please try again.