Genomic coevolution of papillomavirus and immune system in placental mammals indicates the role of IFN-γ in the emergence of new variants.

Papillomaviruses (PVs) are causative agents for warts and cancers in different parts of the body in the mammalian lineage. Therefore, these viruses are proposed as model organisms to study host immune responses to pathogens causing chronic infections. The virus-associated cancer progression depends on two integral processes namely angiogenesis and immune response (AIR). The angiogenesis process aids in tumour progression through vessel formation and maturation but the host immune response, in contrast, makes every attempt to eliminate pathogens and thereby maintain healthy tissues. However, the evolutionary contribution of individual viral genes and host AIR genes in carcinogenesis is yet to be explored. Here, we applied the evolutionary genomics approach to find correlated evolution between six PV genes and 23 host AIR-related genes. We estimated that IFN-γ is the only host gene evolving in a correlated manner with all six PV genes under study. Furthermore, three papillomavirus genes, L2, E6, and E7, are found to interact with two third of host AIR-related genes. Moreover, a combined differential gene expression analysis and network analysis showed that inflammatory cytokine IFN-γ is the key regulator of hub genes in the PPI network of the differentially expressed genes. Functional enrichment of these hub genes is consistent with their established role in different cancers and viral infections. Overall, we conclude that IFN-γ maintains selective pressure on mammalian PV genes and seems to be a potential biomarker for PV-related cancers. This study demonstrates the evolutionary importance of IFN-γ in deciding the fate of carcinogenic PV variants.