Molecular profiling of high-level athlete skeletal muscle after acute endurance or resistance exercise - A systems biology approach.

OBJECTIVE: Long-term high-level exercise training leads to improvements in physical performance and multi-tissue adaptation following changes in molecular pathways. While skeletal muscle baseline differences between exercise-trained and untrained individuals have been previously investigated, it remains unclear how training history influences human multi-omics responses to acute exercise. METHODS: We recruited and extensively characterized 24 individuals categorized as endurance athletes with >15 years of training history, strength athletes or control subjects. Timeseries skeletal muscle biopsies were taken from M. vastus lateralis at three time-points after endurance or resistance exercise was performed and multi-omics molecular analysis performed. RESULTS: Our analyses revealed distinct activation differences of molecular processes such as fatty- and amino acid metabolism and transcription factors such as HIF1A and the MYF-family. We show that endurance athletes have an increased abundance of carnitine-derivates while strength athletes increase specific phospholipid metabolites compared to control subjects. Additionally, for the first time, we show the metabolite sorbitol to be substantially increased with acute exercise. On transcriptional level, we show that acute resistance exercise stimulates more gene expression than acute endurance exercise. This follows a specific pattern, with endurance athletes uniquely down-regulating pathways related to mitochondria, translation and ribosomes. Finally, both forms of exercise training specialize in diverging transcriptional directions, differentiating themselves from the transcriptome of the untrained control group. CONCLUSIONS: We identify a "transcriptional specialization effect" by transcriptional narrowing and intensification, and molecular specialization effects on metabolomic level Additionally, we performed multi-omics network and cluster analysis, providing a novel resource of skeletal muscle transcriptomic and metabolomic profiling in highly trained and untrained individuals.

Defining super-enhancers by highly ranked histone H4 multi-acetylation levels identifies transcription factors associated with glioblastoma stem-like properties.

BACKGROUND: Super-enhancers (SEs), which activate genes involved in cell-type specificity, have mainly been defined as genomic regions with top-ranked enrichment(s) of histone H3 with acetylated K27 (H3K27ac) and/or transcription coactivator(s) including a bromodomain and extra-terminal domain (BET) family protein, BRD4. However, BRD4 preferentially binds to multi-acetylated histone H4, typically with acetylated K5 and K8 (H4K5acK8ac), leading us to hypothesize that SEs should be defined by high H4K5acK8ac enrichment at least as well as by that of H3K27ac. RESULTS: Here, we conducted genome-wide profiling of H4K5acK8ac and H3K27ac, BRD4 binding, and the transcriptome by using a BET inhibitor, JQ1, in three human glial cell lines. When SEs were defined as having the top ranks for H4K5acK8ac or H3K27ac signal, 43% of H4K5acK8ac-ranked SEs were distinct from H3K27ac-ranked SEs in a glioblastoma stem-like cell (GSC) line. CRISPR-Cas9-mediated deletion of the H4K5acK8ac-preferred SEs associated with MYCN and NFIC decreased the stem-like properties in GSCs. CONCLUSIONS: Collectively, our data highlights H4K5acK8ac's utility for identifying genes regulating cell-type specificity.

xcore: an R package for inference of gene expression regulators.

BACKGROUND: Elucidating the Transcription Factors (TFs) that drive the gene expression changes in a given experiment is a common question asked by researchers. The existing methods rely on the predicted Transcription Factor Binding Site (TFBS) to model the changes in the motif activity. Such methods only work for TFs that have a motif and assume the TF binding profile is the same in all cell types. RESULTS: Given the wealth of the ChIP-seq data available for a wide range of the TFs in various cell types, we propose that gene expression modeling can be done using ChIP-seq "signatures" directly, effectively skipping the motif finding and TFBS prediction steps. We present xcore, an R package that allows TF activity modeling based on ChIP-seq signatures and the user's gene expression data. We also provide xcoredata a companion data package that provides a collection of preprocessed ChIP-seq signatures. We demonstrate that xcore leads to biologically relevant predictions using transforming growth factor beta induced epithelial-mesenchymal transition time-courses, rinderpest infection time-courses, and embryonic stem cells differentiated to cardiomyocytes time-course profiled with Cap Analysis Gene Expression. CONCLUSIONS: xcore provides a simple analytical framework for gene expression modeling using linear models that can be easily incorporated into differential expression analysis pipelines. Taking advantage of public ChIP-seq databases, xcore can identify meaningful molecular signatures and relevant ChIP-seq experiments.

Inhibition of BCL2A1 by STAT5 inactivation overcomes resistance to targeted therapies of FLT3-ITD/D835 mutant AML.

Tyrosine kinase inhibitors (TKIs) are established drugs in the therapy of FLT3-ITD mutated acute myeloid leukemia (AML). However, acquired mutations, such as D835 in the tyrosine kinase domain (FLT3-ITD/D835), can induce resistance to TKIs. A cap analysis gene expression (CAGE) technology revealed that the gene expression of BCL2A1 transcription start sites was increased in primary AML cells bearing FLT3-ITD/D835 compared to FLT3-ITD. Overexpression of BCL2A1 attenuated the sensitivity to quizartinib, a type II TKI, and venetoclax, a selective BCL2 inhibitor, in AML cell lines. However, a type I TKI, gilteritinib, inhibited the expression of BCL2A1 through inactivation of STAT5 and alleviated TKI resistance of FLT3-ITD/D835. The combination of gilteritinib and venetoclax showed synergistic effects in the FLT3-ITD/D835 positive AML cells. The promoter region of BCL2A1 contains a BRD4 binding site. Thus, the blockade of BRD4 with a BET inhibitor (CPI-0610) downregulated BCL2A1 in FLT3-mutated AML cells and extended profound suppression of FLT3-ITD/D835 mutant cells. Therefore, we propose that BCL2A1 has the potential to be a novel therapeutic target in treating FLT3-ITD/D835 mutated AML.

Inducing human retinal pigment epithelium-like cells from somatic tissue.

Regenerative medicine relies on basic research outcomes that are only practical when cost effective. The human eyeball requires the retinal pigment epithelium (RPE) to interface the neural retina and the choroid at large. Millions of people suffer from age-related macular degeneration (AMD), a blinding multifactor genetic disease among RPE degradation pathologies. Recently, autologous pluripotent stem-cell-derived RPE cells were prohibitively expensive due to time; therefore, we developed a faster reprogramming system. We stably induced RPE-like cells (iRPE) from human fibroblasts (Fibs) by conditional overexpression of both broad plasticity and lineage-specific transcription factors (TFs). iRPE cells displayed critical RPE benchmarks and significant in vivo integration in transplanted retinas. Herein, we detail the iRPE system with comprehensive single-cell RNA sequencing (scRNA-seq) profiling to interpret and characterize its best cells. We anticipate that our system may enable robust retinal cell induction for basic research and affordable autologous human RPE tissue for regenerative cell therapy.

Development of p53 knockout U87MG cell line for unbiased drug delivery testing system using CRISPR-Cas9 and transcriptomic analysis.

Antibody-drug conjugates offers many advantages as a drug delivery platform that allows for highly specific targeting of cell types and genes. Ideally, testing the efficacy of these systems requires two cell types to be different only in the gene targeted by the drug, with the rest of the cellular machinery unchanged, in order to minimize other potential differences from obscuring the effects of the drug. In this study, we created multiple variants of U87MG cells with targeted mutation in the TP53 gene using the CRISPR-Cas9 system, and determined that their major transcriptional differences stem from the loss of p53 function. Using the transcriptome data, we predicted which mutant clones would have less divergent phenotypes from the wild type and thereby serve as the best candidates to be used as drug delivery testing platforms. Further in vitro and in vivo assays of cell morphology, proliferation rate and target antigen-mediated uptake supported our predictions. Based on the combined analysis results, we successfully selected the best qualifying mutant clone. This study serves as proof-of-principle of the approach and paves the way for extending to additional cell types and target genes.

Functional annotation of human long noncoding RNAs via molecular phenotyping.

Long noncoding RNAs (lncRNAs) constitute the majority of transcripts in the mammalian genomes, and yet, their functions remain largely unknown. As part of the FANTOM6 project, we systematically knocked down the expression of 285 lncRNAs in human dermal fibroblasts and quantified cellular growth, morphological changes, and transcriptomic responses using Capped Analysis of Gene Expression (CAGE). Antisense oligonucleotides targeting the same lncRNAs exhibited global concordance, and the molecular phenotype, measured by CAGE, recapitulated the observed cellular phenotypes while providing additional insights on the affected genes and pathways. Here, we disseminate the largest-to-date lncRNA knockdown data set with molecular phenotyping (over 1000 CAGE deep-sequencing libraries) for further exploration and highlight functional roles for ZNF213-AS1 and lnc-KHDC3L-2.

JQ1 affects BRD2-dependent and independent transcription regulation without disrupting H4-hyperacetylated chromatin states.

The bromodomain and extra-terminal domain (BET) proteins are promising drug targets for cancer and immune diseases. However, BET inhibition effects have been studied more in the context of bromodomain-containing protein 4 (BRD4) than BRD2, and the BET protein association to histone H4-hyperacetylated chromatin is not understood at the genome-wide level. Here, we report transcription start site (TSS)-resolution integrative analyses of ChIP-seq and transcriptome profiles in human non-small cell lung cancer (NSCLC) cell line H23. We show that di-acetylation at K5 and K8 of histone H4 (H4K5acK8ac) co-localizes with H3K27ac and BRD2 in the majority of active enhancers and promoters, where BRD2 has a stronger association with H4K5acK8ac than H3K27ac. Although BET inhibition by JQ1 led to complete reduction of BRD2 binding to chromatin, only local changes of H4K5acK8ac levels were observed, suggesting that recruitment of BRD2 does not influence global histone H4 hyperacetylation levels. This finding supports a model in which recruitment of BET proteins via histone H4 hyperacetylation is predominant over hyperacetylation of histone H4 by BET protein-associated acetyltransferases. In addition, we found that a remarkable number of BRD2-bound genes, including MYC and its downstream target genes, were transcriptionally upregulated upon JQ1 treatment. Using BRD2-enriched sites and transcriptional activity analysis, we identified candidate transcription factors potentially involved in the JQ1 response in BRD2-dependent and -independent manner.

Transcriptional landscape of Mycobacterium tuberculosis infection in macrophages.

Mycobacterium tuberculosis (Mtb) infection reveals complex and dynamic host-pathogen interactions, leading to host protection or pathogenesis. Using a unique transcriptome technology (CAGE), we investigated the promoter-based transcriptional landscape of IFNγ (M1) or IL-4/IL-13 (M2) stimulated macrophages during Mtb infection in a time-kinetic manner. Mtb infection widely and drastically altered macrophage-specific gene expression, which is far larger than that of M1 or M2 activations. Gene Ontology enrichment analysis for Mtb-induced differentially expressed genes revealed various terms, related to host-protection and inflammation, enriched in up-regulated genes. On the other hand, terms related to dis-regulation of cellular functions were enriched in down-regulated genes. Differential expression analysis revealed known as well as novel transcription factor genes in Mtb infection, many of them significantly down-regulated. IFNγ or IL-4/IL-13 pre-stimulation induce additional differentially expressed genes in Mtb-infected macrophages. Cluster analysis uncovered significant numbers, prolonging their expressional changes. Furthermore, Mtb infection augmented cytokine-mediated M1 and M2 pre-activations. In addition, we identified unique transcriptional features of Mtb-mediated differentially expressed lncRNAs. In summary we provide a comprehensive in depth gene expression/regulation profile in Mtb-infected macrophages, an important step forward for a better understanding of host-pathogen interaction dynamics in Mtb infection.

Integrative CAGE and DNA Methylation Profiling Identify Epigenetically Regulated Genes in NSCLC.

Lung cancer is the leading cause of cancer-related deaths worldwide. The majority of cancer driver mutations have been identified; however, relevant epigenetic regulation involved in tumorigenesis has only been fragmentarily analyzed. Epigenetically regulated genes have a great theranostic potential, especially in tumors with no apparent driver mutations. Here, epigenetically regulated genes were identified in lung cancer by an integrative analysis of promoter-level expression profiles from Cap Analysis of Gene Expression (CAGE) of 16 non-small cell lung cancer (NSCLC) cell lines and 16 normal lung primary cell specimens with DNA methylation data of 69 NSCLC cell lines and 6 normal lung epithelial cells. A core set of 49 coding genes and 10 long noncoding RNAs (lncRNA), which are upregulated in NSCLC cell lines due to promoter hypomethylation, was uncovered. Twenty-two epigenetically regulated genes were validated (upregulated genes with hypomethylated promoters) in the adenocarcinoma and squamous cell cancer subtypes of lung cancer using The Cancer Genome Atlas data. Furthermore, it was demonstrated that multiple copies of the REP522 DNA repeat family are prominently upregulated due to hypomethylation in NSCLC cell lines, which leads to cancer-specific expression of lncRNAs, such as RP1-90G24.10, AL022344.4, and PCAT7. Finally, Myeloma Overexpressed (MYEOV) was identified as the most promising candidate. Functional studies demonstrated that MYEOV promotes cell proliferation, survival, and invasion. Moreover, high MYEOV expression levels were associated with poor prognosis.Implications: This report identifies a robust list of 22 candidate driver genes that are epigenetically regulated in lung cancer; such genes may complement the known mutational drivers.Visual Overview: http://mcr.aacrjournals.org/content/molcanres/15/10/1354/F1.large.jpg Mol Cancer Res; 15(10); 1354-65. ©2017 AACR.

Transcriptome Analysis of Recurrently Deregulated Genes across Multiple Cancers Identifies New Pan-Cancer Biomarkers.

Genes that are commonly deregulated in cancer are clinically attractive as candidate pan-diagnostic markers and therapeutic targets. To globally identify such targets, we compared Cap Analysis of Gene Expression profiles from 225 different cancer cell lines and 339 corresponding primary cell samples to identify transcripts that are deregulated recurrently in a broad range of cancer types. Comparing RNA-seq data from 4,055 tumors and 563 normal tissues profiled in the The Cancer Genome Atlas and FANTOM5 datasets, we identified a core transcript set with theranostic potential. Our analyses also revealed enhancer RNAs, which are upregulated in cancer, defining promoters that overlap with repetitive elements (especially SINE/Alu and LTR/ERV1 elements) that are often upregulated in cancer. Lastly, we documented for the first time upregulation of multiple copies of the REP522 interspersed repeat in cancer. Overall, our genome-wide expression profiling approach identified a comprehensive set of candidate biomarkers with pan-cancer potential, and extended the perspective and pathogenic significance of repetitive elements that are frequently activated during cancer progression.

CAGE profiling of ncRNAs in hepatocellular carcinoma reveals widespread activation of retroviral LTR promoters in virus-induced tumors.

An increasing number of noncoding RNAs (ncRNAs) have been implicated in various human diseases including cancer; however, the ncRNA transcriptome of hepatocellular carcinoma (HCC) is largely unexplored. We used CAGE to map transcription start sites across various types of human and mouse HCCs with emphasis on ncRNAs distant from protein-coding genes. Here, we report that retroviral LTR promoters, expressed in healthy tissues such as testis and placenta but not liver, are widely activated in liver tumors. Despite HCC heterogeneity, a subset of LTR-derived ncRNAs were more than 10-fold up-regulated in the vast majority of samples. HCCs with a high LTR activity mostly had a viral etiology, were less differentiated, and showed higher risk of recurrence. ChIP-seq data show that MYC and MAX are associated with ncRNA deregulation. Globally, CAGE enabled us to build a mammalian promoter map for HCC, which uncovers a new layer of complexity in HCC genomics.

Batf2/Irf1 induces inflammatory responses in classically activated macrophages, lipopolysaccharides, and mycobacterial infection.

Basic leucine zipper transcription factor Batf2 is poorly described, whereas Batf and Batf3 have been shown to play essential roles in dendritic cell, T cell, and B cell development and regulation. Batf2 was drastically induced in IFN-γ-activated classical macrophages (M1) compared with unstimulated or IL-4-activated alternative macrophages (M2). Batf2 knockdown experiments from IFN-γ-activated macrophages and subsequent expression profiling demonstrated important roles for regulation of immune responses, inducing inflammatory and host-protective genes Tnf, Ccl5, and Nos2. Mycobacterium tuberculosis (Beijing strain HN878)-infected macrophages further induced Batf2 and augmented host-protective Batf2-dependent genes, particularly in M1, whose mechanism was suggested to be mediated through both TLR2 and TLR4 by LPS and heat-killed HN878 (HKTB) stimulation experiments. Irf1 binding motif was enriched in the promoters of Batf2-regulated genes. Coimmunoprecipitation study demonstrated Batf2 association with Irf1. Furthermore, Irf1 knockdown showed downregulation of IFN-γ- or LPS/HKTB-activated host-protective genes Tnf, Ccl5, Il12b, and Nos2. Conclusively, Batf2 is an activation marker gene for M1 involved in gene regulation of IFN-γ-activated classical macrophages, as well as LPS/HKTB-induced macrophage stimulation, possibly by Batf2/Irf1 gene induction. Taken together, these results underline the role of Batf2/Irf1 in inducing inflammatory responses in M. tuberculosis infection.

Cancers of unknown primary origin (CUP) are characterized by chromosomal instability (CIN) compared to metastasis of know origin.

BACKGROUND: Cancers of unknown primary (CUPs) constitute ~5% of all cancers. The tumors have an aggressive biological and clinical behavior. The aim of the present study has been to uncover whether CUPs exhibit distinct molecular features compared to metastases of known origin. METHODS: Employing genome wide transcriptome analysis, Linear Discriminant Analysis (LDA) and Quadratic Discriminant Analysis (QDA), we defined the putative origins of a large series of CUP and how closely related a particular CUP was to corresponding metastases of known origin. LDA predictions were subsequently used to define a universal CUP core set of differentially expressed genes, that by means of gene set enrichment analysis was exploited to depict molecular pathways characterizing CUP. RESULTS: The analyses show that CUPs are distinct from metastases of known origin. CUPs exhibit inconsistent expression of conventional cancer biomarkers and QDA derived outlier scores show that CUPs are more distantly related to their primary tumor class than corresponding metastases of known origin. Gene set enrichment analysis showed that CUPs display increased expression of genes involved in DNA damage repair and mRNA signatures of chromosome instability (CIN), indicating that CUPs are chromosome unstable compared to metastases of known origin. CONCLUSIONS: CIN may account for the uncommon clinical presentation, chemoresistance and poor outcome in patients with CUP and warrant selective diagnostic strategies and treatment.

Down-Regulation of miR-129-5p and the let-7 Family in Neuroendocrine Tumors and Metastases Leads to Up-Regulation of Their Targets Egr1, G3bp1, Hmga2 and Bach1.

Expression of miRNAs in Neuroendocrine Neoplasms (NEN) is poorly characterized. We therefore wanted to examine the miRNA expression in Neuroendocrine Tumors (NETs), and identify their targets and importance in NET carcinogenesis. miRNA expression in six NEN primary tumors, six NEN metastases and four normal intestinal tissues was characterized using miRNA arrays, and validated by in-situ hybridization and qPCR. Among the down-regulated miRNAs miR-129-5p and the let-7f/let-7 family, were selected for further characterization. Transfection of miR-129-5p inhibited growth of a pulmonary and an intestinal carcinoid cell line. Analysis of mRNA expression changes identified EGR1 and G3BP1 as miR-129-5p targets. They were validated by luciferase assay and western blotting, and found robustly expressed in NETs by immunohistochemistry. Knockdown of EGR1 and G3BP1 mimicked the growth inhibition induced by miR-129-5p. let-7 overexpression inhibited growth of carcinoid cell lines, and let-7 inhibition increased protein content of the transcription factor BACH1 and its targets MMP1 and HMGA2, all known to promote bone metastases. Immunohistochemistry analysis revealed that let-7 targets are highly expressed in NETs and metastases. We found down-regulation of miR-129-5p and the let-7 family, and identified new neuroendocrine specific targets for these miRNAs, which contributes to the growth and metastatic potential of these tumors.

MicroRNA expression analysis and Multiplex ligation-dependent probe amplification in metastatic and non-metastatic uveal melanoma.

PURPOSE: To determine the association of microRNA expression and chromosomal changes with metastasis and survival in uveal melanoma (UM). METHODS: Thirty-six patients with UM were selected based on the metastatic status, and clinicopathological data were collected. Multiplex ligation-dependent probe amplification (MLPA) was used to identify chromosomal changes. Chromosomal changes and clinicopathological data were correlated with survival and metastasis. The microRNA expression was analysed in 26 of the 36 archived UM samples. Unsupervised analysis, differential expression analysis and Cox regression analysis were performed to determine the association with metastasis and survival. RESULTS: Metastasis and metastatic death occurred in 20 patients, two patients died of other causes and one patient of unknown causes. A significant association between increasing size category (p = 0.002, log-rank), extraocular extension (p = 0.001), chromosome 3 loss (p = 0.033) and 1p loss (p = 0.030) and development of metastases was observed. Tumour, node, metastasis (TNM) staging showed a significant association with survival (p < 0.0001, log-rank). Adjusting for gender and age TNM size category T4 (p = 0.016, Cox regression analysis), mixed (p = 0.029) and epithelioid (p = 0.0058) cell types, chromosome 3 loss (p = 0.014) and 8q gain (p = 0.010) were significant prognosticators for a poor survival. Hierarchical clustering divided the UM into three groups based on microRNA expression. The clusters showed no association with clinical or histopathological features, TNM staging, metastasis or survival. Differential expression analysis did not reveal microRNAs related to metastasis or survival. CONCLUSIONS: The prognostic significance of chromosome 3 loss and 8q gain identified by MLPA analysis was confirmed in archived UM samples. The value of microRNA expression as a predictor of metastasis and survival in UM could not be confirmed.

HemaExplorer: a database of mRNA expression profiles in normal and malignant haematopoiesis.

The HemaExplorer (http://servers.binf.ku.dk/hemaexplorer) is a curated database of processed mRNA Gene expression profiles (GEPs) that provides an easy display of gene expression in haematopoietic cells. HemaExplorer contains GEPs derived from mouse/human haematopoietic stem and progenitor cells as well as from more differentiated cell types. Moreover, data from distinct subtypes of human acute myeloid leukemia is included in the database allowing researchers to directly compare gene expression of leukemic cells with those of their closest normal counterpart. Normalization and batch correction lead to full integrity of the data in the database. The HemaExplorer has comprehensive visualization interface that can make it useful as a daily tool for biologists and cancer researchers to assess the expression patterns of genes encountered in research or literature. HemaExplorer is relevant for all research within the fields of leukemia, immunology, cell differentiation and the biology of the haematopoietic system.

The rectal cancer microRNAome--microRNA expression in rectal cancer and matched normal mucosa.

PURPOSE: miRNAs play a prominent role in a variety of physiologic and pathologic biologic processes, including cancer. For rectal cancers, only limited data are available on miRNA expression profiles, whereas the underlying genomic and transcriptomic aberrations have been firmly established. We therefore, aimed to comprehensively map the miRNA expression patterns of this disease. EXPERIMENTAL DESIGN: Tumor biopsies and corresponding matched mucosa samples were prospectively collected from 57 patients with locally advanced rectal cancers. Total RNA was extracted, and tumor and mucosa miRNA expression profiles were subsequently established for all patients. The expression of selected miRNAs was validated using semi-quantitative real-time PCR. RESULTS: Forty-nine miRNAs were significantly differentially expressed (log(2)-fold difference >0.5 and P < 0.001) between rectal cancer and normal rectal mucosa. The predicted targets for these miRNAs were enriched for the following pathways: Wnt, TGF-beta, mTOR, insulin, mitogen-activated protein kinase, and ErbB signaling. Thirteen of these 49 miRNAs seem to be rectal cancer-specific, and have not been previously reported for colon cancers: miR-492, miR-542-5p, miR-584, miR-483-5p, miR-144, miR-2110, miR-652, miR-375, miR-147b, miR-148a, miR-190, miR-26a/b, and miR-338-3p. Of clinical impact, miR-135b expression correlated significantly with disease-free and cancer-specific survival in an independent multicenter cohort of 116 patients. CONCLUSION: This comprehensive analysis of the rectal cancer miRNAome uncovered novel miRNAs and pathways associated with rectal cancer. This information contributes to a detailed view of this disease. Moreover, the identification and validation of miR-135b may help to identify novel molecular targets and pathways for therapeutic exploitation.

Integrative analyses reveal novel strategies in HPV11,-16 and -45 early infection.

The interaction between human papillomavirus (HPV) and host cells is not well understood. We investigate the early stage of HPV infections by global expression profiling in a cell model, in which HaCaT cells were transfected with HPV11, HPV16 or HPV45 genomes. We report the differential expression of genes not previously implicated in HPV biology, such as the PSG family and ANKRD1, and of genes implicated in the biology of other viruses, e.g. MX1, IFI44 and DDX60. Carcinogenesis-related genes, e.g. ABL2, MGLL and CYR61, were upregulated by high-risk HPV16 and -45. The integrative analysis revealed the suppression of DNA repair by HPV11 and -16, and downregulation of cytoskeleton genes by all HPV types. Various signalling pathways were affected by the HPVs: IL-2 by HPV11; JAK-STAT by HPV16; and TGF-ß, NOTCH and tyrosine kinase signalling by HPV45. This study uncovered novel strategies employed by HPV to establish infection and promote uncontrolled growth.