Fig. 5

Integration of ATAC-seq and RNA-seq data. a) Log2 Fold change correlation analysis (Spearman R test) in the 313 genes with significant and concordant changes in both ATAC-seq (DARs) and RNA-seq (DEGs) (e.g., higher chromatin accessibility corresponding to higher expression levels in HBV infected cells, and vice versa) at 72 h p.i. (co-regulated genes). The 242 genes (77.3%) showing a reduced chromatin accessibility and reduced expression after HBV infection (negative Log2 fold change values) are shadowed in light blue. The 71 genes (22.7%) with increased chromatin accessibility and increased expression after HBV infection (positive Log2 fold change values) are shadowed in light red. b) MSigDB (blue bars), Hallmark MSigDB (green bars), Computational (Comp) MSigDB (orange bars) and KEGG (dark yellow bars) enrichment analysis (ShinyGO 0.76 tool) for21 the 313 concordant ATAC-seq (DARs) and RNA-seq (DEGs) genes at 72 h p.i. (see Table S8). c) Ingenuity Pathways Analysis (IPA) of the 313 concordant ATAC-seq (DARs) and RNA-seq (DEGs) genes at 72 h p.i. d) MSigDB (blue bars), Hallmark MSigDB (green bars), Computational (Comp) MSigDB (orange bars) and KEGG (dark yellow bars) enrichment analysis (ShinyGO 0.76 tool) for 643 genes impacted by HBV infection both in ATAC-sec (DARs) and RNA-seq RNA-seq (DEGs) (comodulated genes) at 72 h p.i. (see Table S9). e) Ingenuity Pathways Analysis (IPA) of the 643 co-modulated genes with significant changes in both RNA-seq (DEGs) and ATAC-seq (DARs) profiles at 72 h p.i