RNAseq in WT and Dnmt3a KO HSCs naive and after 1-month of infection with M. avium

5000-10000 HSCs (KL CD150+ CD48- CD34-) were sorted into lysis buffer from the pools of naive or 1-month infected WT or Dnmt3a-/- mice (n=10-12 per group). RNA was isolated with the NucleoSpin ® RNA Plus XS kit (Macherey Nagel). RNA-seq libraries were prepared by using SMARTer® Stranded Total RNA-Seq Kit v2 - Pico Input Mammalian (Takara Bio Usa). Illumina NovaSeq SP was used for sequencing with a paired-end sequencing length of 10bp. FASTQ files were preprocessed using HTS stream (https://github.com/ibest/HTStream) and the clean FASTQ file were aligned using STAR. Differential expression (DE) analysis of gene expression was performed using Limma-Voom. False discovery rate (FDR)<0.05 was considered statistically significant. We performed gene ontology analysis for differentially expressed genes with q value <0.05.