Library preparation and sequencing
cDNA libraries were prepared from the RNA and libraries were pooled equimolarly and sequenced on a Novaseq platform, generating 60 million paired-end reads of 150 base pairs per sample. Raw data underwent quality control, followed by read alignment to a reference genome. Gene expression quantification and subsequent differential expression analysis were performed to gain insights into PBMCs’ gene expression profiles. The fastq files were uploaded into Galaxy Online, and the samples were mapped to the hg38 RefSeq mRNA transcriptome (UCSC) using the Sailfish algorithm [