scATAC-seq reads were aligned to the GRCh38 (hg38) reference genome and quantified using CellRanger-ATAC count (10x Genomics, v.1.0.0).