GSM9462239: RNA-seq for count.Ctrl3; Homo sapiens; RNA-Seq
Total RNA was isolated and purified using TRIzol reagent (Invitrogen, Carlsbad, CA, USA) following the manufacturer's procedure. The RNA amount and purity of each sample was quantified using NanoDrop ND-1000 (NanoDrop, Wilmington, DE, USA). fastp software (https://github.com/OpenGene/fastp) were used to remove the reads that contained adaptor contamination, low quality bases and undetermined bases with default parameter. Then sequence quality was also verified using fastp. We used HISAT2 (https://ccb.jhu.edu/software/hisat2) to map reads to the reference genome of Homo sapiens GRCh38. The mapped reads of each sample were assembled using StringTie (https://ccb.jhu.edu/software/stringtie) with default parameters. Then, all transcriptomes from all samples were merged to reconstruct a comprehensive transcriptome using gffcompare (https://github.com/gpertea/gffcompare/). After the final transcriptome was generated, StringTie and was used to estimate the expression levels of all transcripts. StringTie was used to perform expression level for mRNAs by calculating FPKM (FPKM = [total_exon_fragments / mapped_reads(millions) × exon_length(kB)]). The differentially expressed mRNAs were selected with fold change > 2 or fold change < 0.5 and with parametric F-test comparing nested linear models (p value < 0.05) by R package edgeR (https://bioconductor.org/packages/release/bioc/html/edgeR.html).