Fig. 2

Immune gene network expression by HPV16⁺ stratified human keratinocytes from different anatomical sites vary significantly. Pairwise analysis, normalized to all medians (p-value cutoff < 0.05) was performed individually to comparing SHKs expression at a single mucosal site before and after HPV16 transfection. Differentially expressed genes were selected from a larger list of genes (immune responses/regulation) for the analysis. Data were analyzed using Ingenuity pathway analysis software (IPA®, QIAGEN Inc., https://www.qiagenbioinformatics.com/products/ingenuitypathway-analysis). Networks were generated using the Core analysis tool (a-e) (Kramer et al. 2014), and a Venn/Euler diagram was generated (f) to visualize common vs. unique gene expression by applying the data sets of SHKs (different sites, resting vs HPV16⁺ cells) to software at http://bioinformatics.psb.ugent.be/webtools/Venn/. The number of genes that were uniquely or commonly expressed by cells from each anatomical site are shown within each intersection (f)