Characterization of In Vitro Differentiation of Human Primary Keratinocytes by RNA-Seq Analysis

Podsumowanie

Presented here is a stepwise procedure for in vitro differentiation of human primary keratinocytes by contact inhibition followed by characterization at the molecular level by RNA-seq analysis.

Streszczenie

Human primary keratinocytes are often used as in vitro models for studies on epidermal differentiation and related diseases. Methods have been reported for in vitro differentiation of keratinocytes cultured in two-dimensional (2D) submerged manners using various induction conditions. Described here is a procedure for 2D in vitro keratinocyte differentiation method by contact inhibition and subsequent molecular characterization by RNA-seq. In brief, keratinocytes are grown in defined keratinocyte medium supplemented with growth factors until they are fully confluent. Differentiation is induced by close contacts between the keratinocytes and further stimulated by excluding growth factors in the medium. Using RNA-seq analyses, it is shown that both 1) differentiated keratinocytes exhibit distinct molecular signatures during differentiation and 2) the dynamic gene expression pattern largely resembles cells during epidermal stratification. As for comparison to normal keratinocyte differentiation, keratinocytes carrying mutations of the transcription factor p63 exhibit altered morphology and molecular signatures, consistent with their differentiation defects. In conclusion, this protocol details the steps for 2D in vitro keratinocyte differentiation and its molecular characterization, with an emphasis on bioinformatic analysis of RNA-seq data. Because RNA extraction and RNA-seq procedures have been well-documented, it is not the focus of this protocol. The experimental procedure of in vitro keratinocyte differentiation and bioinformatic analysis pipeline can be used to study molecular events during epidermal differentiation in healthy and diseased keratinocytes.

Wprowadzenie

Human primary keratinocytes derived from the human skin are often used as a cellular model to study the biology of the epidermis^1,2,3,4. The stratification of the epidermis can be modeled by keratinocyte differentiation, either in a 2D submerged monolayer fashion or 3D air-lift organotypic model^2,3,5,6,7. Although 3D models have become increasingly important to assess the epidermal structure and function, 2D differentiation models are still widely used, due to their convenience and the possibility to generate large numbers of cells for analyses.

Various conditions have been applied for inducing keratinocyte differentiation in 2D, including addition of serum, high concentration of calcium, lower temperature and inhibition of epidermal growth factor receptors^2,3. Each of these methods has been validated by a number of keratinocyte differentiation marker genes and shown to be effective in assessing keratinocyte differentiation, including under pathological conditions. However, these induction conditions also show differences in their differentiation efficiency and kinetics when specific panels of marker genes are examined^2,3.

One of these methods involves keratinocyte contact inhibition and depletion of growth factors in the culture medium⁸. It has been shown that keratinocytes can differentiate spontaneously when cells reach full density.  Excluding growth factors in the culture medium can further enhance differentiation. The method combining contact inhibition and depleting growth factors has been shown to generate differentiated keratinocytes with gene expression patterns similar to the normal stratified epidermis when using several epidermal markers³, suggesting that this model is suitable for studying normal keratinocyte differentiation. Recently, two comprehensive gene expression analyses of keratinocyte differentiation using this model have been reported^9,10. Researchers validated this model at the molecular level and showed that it can be used to study normal and diseased keratinocyte differentiation.

This protocol describes the procedure for the in vitro differentiation method and molecular analysis of differentiated cells using RNA-seq. It also illustrates characterization of the transcriptome of cells on differentiation day 0 (proliferation stage), day 2, day 4, and day 7 (early, middle, and late differentiation, respectively). It is shown that differentiated keratinocytes display gene expression patterns that largely resemble cells during epidermal stratification. To examine whether this method can be used for studying skin pathology, we applied the same experimental and analysis pipeline to investigate keratinocytes carrying mutations of the transcription factor p63 that are derived from patients with ectrodactyly, ectodermal displasia and cleft lip/palate (EEC) syndrome^11,12. This protocol focuses on the in vitro differentiation of keratinocytes as well as subsequent bioinformatic analysis of RNA-seq. Other steps in the complete procedure such as RNA extraction, RNA-seq sample preparation and library construction, are well documented and can be easily followed, especially when using many commonly used commercial kits. Therefore, these steps are only briefly described in the protocol. The data show that this pipeline is suitable for studying molecular events during epidermal differentiation in healthy and diseased keratinocytes.

Protokół

Skin biopsies were taken from the trunk of healthy volunteers or patients with p63 mutations, to set up the primary keratinocyte culture. All procedures regarding establishing human primary keratinocytes were approved by the ethical committee of the Radboud University Nijmegen Medical Centre (“Commissie Mensgebonden Onderzoek Arnhem-Nijmegen”). Informed consent was obtained.

1. Human primary keratinocyte differentiation by contact inhibition

1. Prepare keratinocyte growth medium (KGM) from Keratinocyte Basal Medium (Table of Materials) when necessary.
2. Make 500 mL proliferation medium using pre-prepared stocks (KGM-pro; see Table 1).
3. Make 500 mL differentiation medium (KGM-dif; see Table 2).
   NOTE: KGM with supplements can be kept at 4 °C for two weeks. For longer storage, it can be aliquoted and stored at -20°C.
4. Seed primary keratinocytes in the density of 5.0–20 x 10³ cells/cm², depending on the cell proliferative capacity. Add sufficient KGM-pro medium to cover cells.
   NOTE: 1) Details of regular cell culture within KGM medium can be found in the website of Lonza¹³. 2) The seeding density should be tested for each cell line. For example, a normal primary keratinocyte line can be seeded at a density of 5.0 x 10³ cells/cm², whereas a line carrying mutations in the transcription factor p63 that is less proliferative should be seeded at a density of 20 x 10³ cells/cm². 3) Depending on the experimental set up, several dishes or wells of cells should be seeded to allow collecting samples on different differentiation days in replicas.
5. Refresh cells with the KGM-pro medium for two days after seeding (seeding on day 0, refreshing for the 1^st time on day 3). Afterwards, refresh with KGM-pro medium every other day. Check cells regularly, at least every other day.
6. Induce cell differentiation when cells are more than 90% confluent by changing the medium to KGM-dif. The day of changing medium to KGM-dif is defined as differentiation day 0. Collect cells for further RNA analyses by first washing 2x with DPBS, afterwards adding in the lysis buffer (from an RNA extraction kit).
   NOTE: With the abovementioned seeding density, cells should reach confluency within 7–10 days. Cells can be seeded with a higher initial density to reach confluency sooner. If cells are not proliferative, longer waiting may not give rise to fully confluent cells. A higher seeding density may be required.
7. Refresh cells with the KGM-dif medium every day and collect cells on differentiation day 2, 4. and 7 for further RNA extraction.

2. RNA extraction

1. Isolate total RNA. This can be performed using a commercial RNA isolation kit (such as Zymo Quick-RNA MicroPrep RNA), or using phenol¹⁴. However, an isolation kit containing a DNAse treatment is highly recommended for RNA-seq samples to prevent DNA contamination and phenol. An example of a commercial RNA-isolation kit is given in the Table of Materials.
2. Measure RNA concentration with a spectrometer. Both DNA and RNA have an absorption peak at 260 nm.
   NOTE: 1) The ratio between the absorbances at 260 nm and 280 nm is used to assess the purity of the RNA. A ratio of around 2.0 is generally accepted as ‘pure’ RNA. If the ratio is substantially lower than 2.0, the sample might be contaminated with contaminants that absorb at 280 nm such as proteins, phenols, or other substances. 2) The ratio between the absorbances at 260 nm and 230 nm measures nucleic acid purity. A ratio between 2.0 and 2.2 is expected. If the ratio is substantially lower, the sample might be contaminated with contaminants that absorb at 230 nm such as EDTA, ethanol, or other substances.

3. RNA quality check

1. Run RNA either on a bioanalyzer RNA Pico chip to validate RNA Integrity Number (RIN) quantitatively, or on a gel for a qualitative measure. In general, a RIN of at least eight is highly advisable. However, when extracting RNA from tissues the RIN may be lower.
   NOTE: Optionally, quality control can be performed by qPCR on the RNA material. 

4. RNA-seq library preparation

NOTE: The RNA-seq library preparation is often performed with a commercial kit or under commercial settings. The described protocol is adapted from a commercial kit, KAPA RNA HyperPrep Kit with RiboErase (Illumina), with a brief description of all required steps: rRNA depletion with oligo hybridization to human ribosomal RNAs, RNA fragmentation, first-strand synthesis, second-strand synthesis and A-tailing, and cleanup after each step¹⁵. Other library preparation kits can also be used for this purpose. It is recommended to perform this step using a commercially available kit, as the quality of generated cDNA library is often more consistent. 2) The following steps are described for 1x library preparation. If preparing several samples, make master mixes with 10% extra volume.

1. Oligo hybridization and rRNA depletion
   1. Prepare oligo hybridization master mix (total = 11 µL, with 4.4 µL of hybridization buffer, 4.4 µL of hybridization oligos, and 2.2 µL of RNase-free water) and depletion master mix (total = 5.5 µL, with 3.3 µL of depletion buffer and 2.2 µL of RNase H).
   2. Set up the PCR reaction program on a thermocycler: 95 °C for 2 min; ramp down to 45 °C at -0.1 °C/s; 45 °C pause; 45 °C for 30 min; 4 °C forever.
      NOTE: Consider starting with double the amount of RNA than necessary for amplification. In the first attempt, use one-half of the amount to continue with the following steps. This is to make sure that there is still material to repeat these steps with a different number of cycles (see step 4.7.2), if the cycles of amplification turn out to be insufficient or overamplified.
   3. Use between 25 ng and 1 µg of total RNA in 10 µL of RNAse-free water, add 10 μL of oligo hybridization master mix. Place samples in the pre-programmed thermocycler and start the program.
   4. When the program reaches the pause step at 45 °C, add 5 µL of depletion master mix to the 20 µL hybridization reaction without removing it from the thermocycler. Mix thoroughly by pipetting up and down several times.
   5. Resume the thermocycler program to continue with the depletion step (45 °C for 30 min).
      NOTE: 1) Put the beads for rRNA depletion (e.g., KAPA pure beads) at room temperature (RT) during the depletion step for the next part of the protocol. 2) Consider preparing the DNase digestion master mix (for section 4.2), since the reagents are needed directly after the rRNA depletion cleanup.
   6. Perform a 2.2x bead-based KAPA Pure Beads cleanup: combine the RNA mix (25 μL) and KAPA Pure Beads (55 μL), by thoroughly resuspend the beads in the RNA mix by pipetting up and down several times.
   7. Incubate at RT for 5 min to allow RNA binding to the beads, and place the tube(s) on a magnet rack to capture the beads until the liquid is clear, and carefully remove and discard 60 µL of supernatant.
   8. Keeping the tube(s) on the magnet, wash 2x with 200 µL of 80% ethanol by incubating the beads with the ethanol for ≥30 s and discard the ethanol. Try to remove all residual ethanol without disturbing the beads after the second wash.
   9. Dry the beads at RT for 3–5 min or until all ethanol has evaporated.
      NOTE: Over-drying the beads may result in reduced yield.
2. DNase digestion
   1. Prepare DNase digestion master mix (total = 22 µL, with 2.2 µL of DNase buffer, 2 µL of DNase, and 17.8 µL of RNase-free water) .
   2. Resuspend the beads in DNAse digestion master mix (20 μL) by pipetting up and down several times. Incubate the tube(s) at RT for 3 min to elute the RNA off the beads.
   3. Place the tube(s) on a magnet rack to capture the beads, until the liquid is clear, and carefully transfer 20 µL of supernatant into a clean tube.
   4. Incubate the tube(s) with supernatant at 37 °C for 30 min.
      NOTE: Consider preparing the RNA elution, fragmentation, and priming master mixes (for section 4.3), since the reagents are directly needed after the DNase digestion cleanup.
   5. Perform a 2.2x bead-based cleanup by following 4.1.6–4.1.9.
3. RNA elution, fragmentation, and priming
   1. Prepare the Fragment, Prime and Elute Buffer (1x) master mix with 11 µL of Fragment, Prime and Elute Buffer (2x) and 11 µL of RNase-free water.
   2. Thoroughly resuspend the beads with purified, DNase-treated RNA in 22 µL of Fragmentation, Prime and Elute Buffer (1x) by pipetting up and down several times.
   3. Incubate at room temperature for 3 min to elute RNA off the beads, then place the tube(s) on a magnet to capture the beads until the liquid is clear.
   4. Carefully transfer 20 µL of supernatant into a tube(s). Discard the tube(s) with beads.
      NOTE: This is a safe stopping point, as samples can be stored at -20 °C for ≤24 h.
   5. Place the tube(s) in a thermocycler and carry out the fragmentation and priming at 94 °C for 6 min. Resulting in approximate 200–300 bp long fragments.
      NOTE: Incubate for 8 min at 94 °C for 100–200 bp fragments. When using partially degraded RNA, fragment between 1–6 min at 85 °C depending on the severity of degradation.
   6. Place the tube(s) on ice and proceed immediately to first strand synthesis.
4. First strand synthesis
   1. Prepare  first strand synthesis master mix (total = 12 µL, with 11 µL of first strand synthesis buffer and 1 µL of KAPA script).
   2. Set up the PCR reaction program on a thermocycler: 25 °C for 10 min; 42 °C for 15 min; 70 °C for 15 min; 4 °C forever.
   3. On ice, combine 20 µL of fragmented primed RNA with 10 µL of first strand synthesis master mix. Keep the tube(s) on ice, mix thoroughly by gently pipetting the reaction up and down several times.
   4. Place the tube(s) in the pre-programmed thermocycler and start the program.
5. Second strand synthesis and A-tailing
   1. Prepare second strand synthesis and A-tailing master mix on ice (total = 33 µL, with 31 µL of second strand synthesis buffer and 2 µL of second strand synthesis & A-tailing enzyme mix).
   2. Set up the PCR reaction program on a thermocycler: 16 °C for 30 min; 62 °C for 10 min; 4 °C forever.
   3. Combine 30 µL of first strand synthesis product with 30 µL of second strand synthesis and A-tailing master mix, and mix thoroughly by pipetting up and down several times on ice.
   4. Place the tube(s) in the pre-programmed thermocycler and start the program.
6. Adapter ligation and post-ligation cleanup
   1. Dilute stock adapters (NEXTflex DNA barcodes, 25 µM) 3.57x in RNase-free water to 7 µM.
   2. Prepare adapter ligation master mix (total = 50 µL, with 40 µL of ligation buffer and 10 µL of DNA ligase).
   3. Combine 60 µL of second strand synthesis product, 45 µL of adapter ligation master mix, and 5 µL of diluted adapters. Mix thoroughly by pipetting up and down several times on ice.
   4. Incubate the tubes at 20 °C for 15 min and continue immediately with the first post ligation cleanup.
   5. Perform a 0.63x bead-based cleanup by combining adapter-ligated DNA (110 µL) and KAPA pure beads (70 µL), then thoroughly resuspend the beads in the DNA mix by pipetting up and down several times.
   6. Repeat steps 4.1.7–4.1.9.
   7. Remove the tubes from the magnet. Thoroughly resuspend the beads in 50 µL of 10 mM Tris HCL (pH = 8.0–8.5) and incubate the beads at RT for 2 min.
   8. Place the plate on the magnet and wait until the liquid is clear. Transfer 50 μL of the clear supernatant to a new tube.
      NOTE: Safe stopping point for less than 24 h at 4 °C.
   9. Perform a 0.7x bead-based cleanup by combining 50 µL of beads with purified adapter-ligated DNA with 35 µL of PEG/NaCL solution. Mix thoroughly by vortexing.
   10. Perform cleanup steps 4.1.7–4.1.9. Thoroughly resuspend the beads in 20 µL of 10 mM Tris HCL (pH = 8.0–8.5), and incubate the beads at RT for 2 min.
   11. Place the plate on the magnet; wait until the liquid is clear. Transfer 20 μL of the clear supernatant to a new tube and proceed to section 4.7.
       NOTE: Safe stopping point for less than 1 week at 4 °C or less than 1 month at -20 °C.
7. Library amplification and cleanup
   1. Prepare library amplification master mix (total = 33 µL, with 27.5 µL of 2x KAPA HiFi HotStart ReadyMix and 5.5 µL of 10x library amplification primer mix).
   2. Set up the PCR reaction program on a thermocycler using the following parameters. 98 °C for 45 s; N cycles (see the note below) of: 98 °C for 15 s; 60 °C for 30 s; 72 °C for 30 s; 72 °C for 1 min; and 4 °C forever.
      NOTE: The cycles for the library amplification depend on the input material. It can roughly be estimated as such: starting with RNA = 25–100 ng, N = 11–15 cycles; 100–250 ng, N = 9–12 cycles; 250–500 ng, N = 7–10 cycles.
   3. Perform a 0.8x bead-based cleanup by combining amplified library DNA (50 µL) and KAPA pure beads (40 µL), then thoroughly resuspend the beads in the DNA mix by pipetting up and down several times.
   4. Perform cleanup steps 4.1.7–4.1.9. Thoroughly resuspend the beads in 50 µL of 10 mM Tris HCL (pH = 8.0–8.5), and incubate the beads at RT for 2 min.
   5. Transfer 50 μL of the clear supernatant to a new tube and proceed to the second library amplification cleanup.
   6. Perform a 1x bead-based cleanup by combining amplified library DNA (50 µL) and KAPA pure beads (50 µL), then thoroughly resuspend the beads in the DNA mix by pipetting up and down several times.
   7. Perform cleanup steps 4.1.7–4.1.9. Thoroughly resuspend the beads in 22 µL of 10 mM Tris HCL (pH = 8.0–8.5), and incubate the beads at RT for 2 min.
   8. Transfer 20 μL of the clear supernatant to a new tube proceed to Qbit and bioanalyzer measurements.
8. Concentration and fragmentation size
   1. Measure DNA concentrations of the samples. Using a highly sensitive fluorescent dye-based kit (e.g., Denovix Qbit, see Table of Materials), or if the concentration seems to be below 0.5 ng/μL, requantify using a more sensitive qPCR approach (e.g., Kappa Quantification, see Table of Materials).
   2. Determine the fragment size of the library, using high sensitive electrophoresis (e.g., a bioanalyzer).
      NOTE: Optionally, quality control by qPCR can be performed before sequencing (Appendix) using a diluted prepared library instead of cDNA.
9. Sequencing
   1. Send the library for sequencing. For RNA-seq differential gene expression analysis, a sequencing depth of reads between 10–25 million per sample is generally sufficient.

5. Data pre-processing

1. Quality check Fastq files
   1. Download and install Trimgalore¹⁶ to remove low quality base pairs and empty reads within the Fastq files.
      NOTE: 1) Most software mentioned here only works in Linux. If limited to a Windows computer, it is possible to run analyses on a Linux server using the software MobaXterm or Putty. Keep in mind that for genome indexing a lot of random-access memory (RAM) is required (around 64 GB). 2) Installing all required software into a conda environment is highly recommended. This will allow easy software installation and package management. For more information on conda, visit <https://docs.conda.io>.
   2. Create a directory (folder) for the data e.g. the folder ‘RNA_seq_KC_diff’. Set this as the working directory by typing: “cd /home/user/RNA_seq_KC_diff/”.
      NOTE: For an overview of the folder structure, see ‘folder_structure.txt’ in the supplementary coding files.
   3. Create several folders in the working directory with the specific names: ‘Fastq’, ‘CRCh38_fasta’, ‘CRCh38’, ‘scripts’, and ‘Mapping’.
   4. Move the fasta files of the sequenced data into the fastq folder. Alternatively, generate soft links using the command: “ln –s home/user/old_location_fastq/FastqFile home/user/RNAseq_KC_diff/Fastq/Fastqfile” to the Fastq files to save disk space.
   5. Run Trim galore on the Fastq files, see the TrimGalore.txt file for the code to copy paste into bash. Change folder names and settings within the command where necessary.
2. Mapping
   1. Download a genome to map the reads to, e.g., hg38 ensembl release 97 (unmasked version): ftp://ftp.ensembl.org/pub/release-97/fasta/homo_sapiens/dna/Homo_sapiens.GRCH38.primary_assembly.fa.gz ; and the corresponding gene annotation file: hg 38 gene annotation ensembl release 97, ftp://ftp.ensembl.org/pub/release-97/gtf/homo_sapiens/Homo_sapiens.GRCh38.97.gtf.gz. Transfer both the files to the CRCh38_fasta folder.
      NOTE: Alternatively use a different version of the genome, for example, the UCSC version of hg38, if mapping to transcript IDs rather than gene IDs is preferred.
   2. Install STAR 2.7.1¹⁷.
   3. Generate an in-house reference genome using STAR 2.7.1 by the generate genome script (see Supplemental Coding Files). Change the amount of threads to what the processer has (--runThreadN X). Change folder names and settings in this script where necessary. Run it by copy/pasting into bash.
   4. Install samtools 1.9¹⁸ for indexing the bam files and gzip to compress the wiggle files.
   5. Align the Trim Galore validated sequencing reads to the human genome assembly hg38 (ensembl release 97) using STAR 2.7.1 followed by indexing using samtools and gzip. This should be performed using the script map_fastq.txt (see Supplemental Coding Files). Change the folder names and settings in this script where necessary. Run it by copying/pasting into bash.
      NOTE: STAR mapping generates several output files, including a bam file, a wiggle file, and read counts table named *_ReadsPerGene.out.tab, which can directly be used concatenated by R to use as input for Deseq2.
3. Genome browser visualization
   1. Install wigToBigWig from the UCSC genome browser tools to generate bigwig files with the big2bw script¹⁹.
   2. Transfer the files ‘convertBigWigChroms.py’ and ‘CRCH38_ensembl2UCSC.txt’ from the supplementary coding files to a folder ‘scripts’ in the working directory.
   3. Make the convertBigWigChroms.py executable (using the command: ‘chmod +x scripts/convertBigWigChroms.py’ on a Linux machine).
   4. Generate bigWig files from the Wiggle files, using the script wig2bw.txt (see Supplemental Coding Files). Run it by copy/pasting into bash.
   5. Input the BigWig files on the UCSC genome browser or Integrative Genomics Viewer (IGV) for visualization.

6. RNA-seq data analysis

1. Input sample data into Deseq2
   1. Download and install Rstudio (version 1.1.456) and R (version 3.6).
   2. Install all the required R packages (see Supplemental coding files).
      NOTE: For a detailed R-markdown file showing all programming code and output, see the attached files: “RNA_seq_kc_differentiation_wt.html” and “RNA_seq_kc_differentiation_patient.html”. A general description of all steps is included below.
   3. Generate a count table from the ReadsPerGene.out.tab files.
   4. Write a sample data file, containing all filenames, the day of differentiation and other relevant sample data. For an example sample file, see “sample_data_example.csv” in the Supplemental Coding Files.
   5. Use the count table and sample data to generate a Deseq2²⁰ object containing both the countable and sample data.
2. Gene expression normalization, sample distance and PCA
   1. Normalize the count table in the Deseq2 object, using either Deseq2 rld, or vst normalization. Rld normalization is preferred, but with many samples, vst normalization is much faster.
   2. Plot sample distance based on the normalized read count intensities using the “dist” function in R and subsequently performing “hclust” clustering based on sample distance. Plot the heatmap itself using the pheatmap package.
   3. Generate a PCA plot of the normalized read count intensities using the plotPCA function of Deseq2.
      NOTE: PCA serves as a tool in exploratory data analysis and can be used to visualize distance and relatedness between different samples.
3. Differential/highly variable gene expression analysis
   1. Calculate either the differential genes using the Deseq2 result function. However, when assessing changes over several timepoints in which pairwise comparisons are not applicable, use highly variable genes. In this case, extract the top 500 highly variable genes via ordering on the variance between samples from different time points with the rowVars function.
      NOTE: In our analysis, for the control samples, the top 500 highly variable genes (the genes with the highest standard deviation of their normalized intensity across days of differentiation) were used for further analysis. However, for the diseased vs control analysis differentially expressed genes between disease and healthy controls over time were calculated by Deseq2 using a multiple testing corrected p-value of 1 x 10^-4 as a cutoff. Another possible filtering step is to exclude differential genes changing less than a certain fold change cutoff.
   2. Perform kmeans clustering on the differential or highly variable genes to cluster them by different expression patterns.
   3. Visualize differential or highly variable genes in a heatmap using the pheatmap package. The intensity plotted in the heatmap is the Deseq2 normalized intensity with the median value subtracted.
4. Gene Ontology (GO) annotation enrichment analysis
   1. Generate a list of expressed background genes by taking all genes with more than 10 counts in a single sample.
   2. Perform GO analysis using online tools such as GOrilla²¹. Use the lists of the differential/highly variable genes in clusters as “gene list”, and the background genes as the background for comparison.
      NOTE: More advanced R-users can use the package clusterProfiler²² to automate Go-term enrichment analysis.

Wyniki

Normal keratinocyte differentiation and RNA-seq analysis
In this experiment, keratinocyte lines derived from five individuals were used for differentiation and RNA-seq analyses. Figure 1 summarizes the experimental procedure of differentiation and RNA-seq analysis results. An overview of in vitro differentiation procedures of normal keratinocytes and cell morphology changes during differentiation are illustrated in Figure 1A. Principle com...

Dyskusje

This work describes a method for inducing human keratinocyte differentiation and subsequent characterization using RNA-seq analyses. In the current literature, many studies on human keratinocyte differentiation use two other methods, with a high calcium concentration or with serum as methods to induce differentiation^2,3,23. A previous report carefully compared these three different methods³ and showe...

Materiały

Name	Company	Catalog Number	Comments
Bioanalyzer 2100	Agilent	G2929BA
Bovine pituitary extract (BPE)	Lonza		Part of the bulletKit
CFX96 Real-Time system	Bio-Rad		qPCR machine
Dulbecco's Phosphate-Buffered Saline (DPBS)	Sigma-Aldrich	D8537
Epidermal Growth Factor (EGF)	Lonza		Part of the bulletKit
Ethanolamine >= 98%	Sigma-Aldrich	E9508
High Sensitivity DNA chips	Agilent	5067-4626
Hydrocortison	Lonza		Part of the bulletKit
Insulin	Lonza		Part of the bulletKit
iQ SYBR Green Kit	BioRad	170-8886
iScript cDNA synthesis	Bio rad	1708890
KAPA Library Quant Kit	Roche	07960255001	Low concentration measure kit
KAPA RNA HyperPrep Kit with RiboErase	Roche	KK8540	RNAseq kit
KGM Gold Keratinocyte Growth Medium BulletKit	Lonza	192060
Nanodrop	deNovix	DS-11 FX (model)	Nanodrop and Qbit for DNA and RNA measurements
NEXTflex DNA barcodes -24	Illumnia	NOVA-514103	6 bp long primers
Penicillin-Streptomycin	Gibco	15140122
RNA Pico Chip	Agilent	5067-1513