FACS-sorted tTreg from WT-Foxp3RFP and Kcnk18G339R/Foxp3RFP thymus were used for bulk RNA sequencing. RNA was isolated with Quick-RNA Microprep Kit (Zymo Research) as described above. Quality and amount of RNA were verified by NanoDrop and Bioanalyzer RNA 6000 Nano Kit (Agilent). Samples with RIN values > 6.5 were used for RNA sequencing. NEBNext® rRNA depletion was performed followed by NEBNext directional Ultra RNA II Library preparation and sequencing on NextSeq500 (Illumina) platform (75 cycles, high output v2 kit). Raw sequencing data were analyzed by Linux bash tools following the analysis pipeline: (1) quality control (fastqc), (2) trimming (Trimmomatic 0.36),61 (link) (3) alignment to mouse genome (Hisat 2.1.0; build: mm10),62 (link) (4) aligned read sorting (Samtools)63 (link) and (5) read counting (HTseq 0.10.0.64 (link)) Expression analysis was performed with R/Bioconducter DESeq2.65 (link) Treg signature genes were assigned according to a previous study.66 (link) Significantly regulated genes (FDR < 0.05) were used for further analysis and visualization using pHeatmap (https://cran.r-project.org/package=pheatmap) package in R. GO term gene enrichment analysis was performed using the PANTHER classification system online tool.67 (link),68 (link)
Free full text: Click here