Genome-wide methylation and expression differences in HPV(+) and HPV(?) squamous cell carcinoma cell lines are consistent with divergent mechanisms of carcinogenesis

Oncogenic human papillomaviruses (HPV) are associated with nearly all cervical cancers and are increasingly important in the etiology of oropharyngeal tumors. HPV-associated head and neck squamous cell carcinomas (HNSCC) have distinct risk profiles and appreciate a prognostic advantage compared to HPV-negative HNSCC. Promoter hypermethylation is widely recognized as a mechanism in the progression of HNSCC, but the extent to which this mechanism is consistent between HPV(+) and HPV(−) tumors is unknown. To investigate the epigenetic regulation of gene expression in HPV-induced and carcinogen-induced cancers, we examined genome-wide DNA methylation and gene expression in HPV(+) and HPV(−) SCC cell lines. We used two platforms: the Illumina Infinium Methylation BeadArray and tiling arrays, and confirmed illustrative examples with pyrosequencing and quantitative PCR. These analyses indicate that HPV(+) cell lines have higher DNA methylation in genic and LINE-1 regions than HPV(−) cell lines. Differentially methylated loci between HPV(+) and HPV(−) cell lines significantly correlated with HPV-typed HNSCC primary tumor DNA methylation levels. Novel findings include higher promoter methylation of polycomb repressive complex 2 target genes in HPV(+) cells compared to HPV(−) cells and increased expression of DNMT3A in HPV(+) cells. Additionally, CDKN2A and KRT8 were identified as interaction hubs among genes with higher methylation and lower expression in HPV(−) cells. Conversely, RUNX2, IRS-1 and CCNA1 were major hubs with higher methylation and lower expression in HPV(+) cells. Distinct HPV(+) and HPV(−) epigenetic profiles should provide clues to novel targets for development of individualized therapeutic strategies.Key words: epigenetics, human papillomavirus, HNSCC, DNA methylation, squamous cell carcinoma, gene expression, microarrays, illumina infinium humanmethylation27 beadarray     

Methylated Host Cell Gene Promoters and Human Papillomavirus Type 16 and 18 Predicting Cervical Lesions and Cancer

Change in the host and/or human papillomavirus (HPV) DNA methylation profile is probably one of the main factors responsible for the malignant progression of cervical lesions to cancer. To investigate those changes we studied 173 cervical samples with different grades of cervical lesion, from normal to cervical cancer. The methylation status of nine cellular gene promoters, CCNA1, CDH1, C13ORF18, DAPK1, HIC1, RARβ2, hTERT1, hTERT2 and TWIST1, was investigated by Methylation Specific Polymerase Chain Reaction (MSP). The methylation of HPV18 L1-gene was also investigated by MSP, while the methylated cytosines within four regions, L1, 5’LCR, enhancer, and promoter of the HPV16 genome covering 19 CpG sites were evaluated by bisulfite sequencing. Statistically significant methylation biomarkers distinguishing between cervical precursor lesions from normal cervix were primarily C13ORF18 and secondly CCNA1, and those distinguishing cervical cancer from normal or cervical precursor lesions were CCNA1, C13ORF18, hTERT1, hTERT2 and TWIST1. In addition, the methylation analysis of individual CpG sites of the HPV16 genome in different sample groups, notably the 7455 and 7694 sites, proved to be more important than the overall methylation frequency. The majority of HPV18 positive samples contained both methylated and unmethylated L1 gene, and samples with L1-gene methylated forms alone had better prognosis when correlated with the host cell gene promoters’ methylation profiles. In conclusion, both cellular and viral methylation biomarkers should be used for monitoring cervical lesion progression to prevent invasive cervical cancer.     

Genome-wide methylation and expression differences in HPV(+) and HPV(-) squamous cell carcinoma cell lines are consistent with divergent mechanisms of carcinogenesis

Oncogenic human papillomaviruses (HPV) are associated with nearly all cervical cancers and are increasingly important in the etiology of oropharyngeal tumors. HPV-associated head and neck squamous cell carcinomas (HNSCC) have distinct risk profiles and appreciate a prognostic advantage compared to HPV-negative HNSCC. Promoter hypermethylation is widely recognized as a mechanism in the progression of HNSCC, but the extent to which this mechanism is consistent between HPV(+) and HPV(-) tumors is unknown. To investigate the epigenetic regulation of gene expression in HPV-induced and carcinogen-induced cancers, we examined genome-wide DNA methylation and gene expression in HPV(+) and HPV(-) SCC cell lines. We used two platforms: the Illumina Infinium Methylation BeadArray and tiling arrays, and confirmed illustrative examples with pyrosequencing and quantitative PCR. These analyses indicate that HPV(+) cell lines have higher DNA methylation in genic and LINE-1 regions than HPV(-) cell lines. Differentially methylated loci between HPV(+) and HPV(-) cell lines significantly correlated with HPV-typed HNSCC primary tumor DNA methylation levels. Novel findings include higher promoter methylation of polycomb repressive complex 2 target genes in HPV(+) cells compared to HPV(-) cells and increased expression of DNMT3A in HPV(+) cells. Additionally, CDKN2A and KRT8 were identified as interaction hubs among genes with higher methylation and lower expression in HPV(-) cells. Conversely, RUNX2, IRS-1 and CCNA1 were major hubs with higher methylation and lower expression in HPV(+) cells. Distinct HPV(+) and HPV(-) epigenetic profiles should provide clues to novel targets for development of individualized therapeutic strategies.     

Leishmania donovani Infection Causes Distinct Epigenetic DNA Methylation Changes in Host Macrophages

Infection of macrophages by the intracellular protozoan Leishmania leads to down-regulation of a number of macrophage innate host defense mechanisms, thereby allowing parasite survival and replication. The underlying molecular mechanisms involved remain largely unknown. In this study, we assessed epigenetic changes in macrophage DNA methylation in response to infection with L. donovani as a possible mechanism for Leishmania driven deactivation of host defense. We quantified and detected genome-wide changes of cytosine methylation status in the macrophage genome resulting from L. donovani infection. A high confidence set of 443 CpG sites was identified with changes in methylation that correlated with live L. donovani infection. These epigenetic changes affected genes that play a critical role in host defense such as the JAK/STAT signaling pathway and the MAPK signaling pathway. These results provide strong support for a new paradigm in host-pathogen responses, where upon infection the pathogen induces epigenetic changes in the host cell genome resulting in downregulation of innate immunity thereby enabling pathogen survival and replication. We therefore propose a model whereby Leishmania induced epigenetic changes result in permanent down regulation of host defense mechanisms to protect intracellular replication and survival of parasitic cells.     

Epigenetic status of cell cycle regulation genes in the placenta of human embryos with chromosomal mosaicism

Epigenetic mechanisms of cell cycle regulation providing for maintenance of genome stability and precise transmission of hereditary information to the daughter cells attract considerable interest. Up-to-date molecular technologies allow the genome-wide methylation pattern to be determined in a single experiment. In this work, we pioneered in determining the epigenetic status of the placental tissues of the human embryos with mosaic karyotype using a genome-wide Illumina Infinium HumanMethylation27 BeadChip array (Illumina, United States), comprising 27578 CpG dinucleotides contained in 14475 genes. The groups of genes associated with the cell cycle and its regulation that contain differentially methylated CpG sites in their promoter regions have been determined. It was demonstrated that the methylation level most frequently changed in oncogenes (ARHGEF1 and FGF5), tumor suppressors (APC2, BRACA2, DCC, GRLF1, RB1, TP73, TSPYL2, and VHL), and the genes involved in the regulation of chromosome segregation (CNTROB, GMNN, PROCR, and TACC1).     

Germline-derived DNA methylation and early embryo epigenetic reprogramming: The selected survival of imprints

DNA methylation is an essential epigenetic mechanism involved in many essential cellular processes. During development epigenetic reprograming takes place during gametogenesis and then again in the pre-implantation embryo. These two reprograming windows ensure genome-wide removal of methylation in the primordial germ cells so that sex-specific signatures can be acquired in the sperm and oocyte. Following fertilization the majority of this epigenetic information is erased to give the developing embryo an epigenetic profile coherent with pluripotency. It is estimated that ∼65% of the genome is differentially methylated between the gametes, however following embryonic reprogramming only parent-of-origin methylation at known imprinted loci remains. This suggests that trans-acting factors such as Zfp57 can discriminate imprinted differentially methylated regions (DMRs) from the thousands of CpG rich regions that are differentially marked in the gametes. Recently transient imprinted DMRs have been identified suggesting that these loci are also protected from pre-implantation reprograming but succumb to de novo remethylation at the implantation stage. This highlights that “ubiquitous” imprinted loci are also resilient to gaining methylation by protecting their unmethylated alleles. In this review I examine the processes involved in epigenetic reprograming and the mechanisms that ensure allelic methylation at imprinted loci is retained throughout the life of the organism, discussing the critical differences between mouse and humans.This article is part of a Directed Issue entitled: Epigenetics Dynamics in development and disease.     

Mycoplasma CG- and GATC-specific DNA methyltransferases selectively and efficiently methylate the host genome and alter the epigenetic landscape in human cells

Aberrant DNA methylation is frequently observed in disease, including many cancer types, yet the underlying mechanisms remain unclear. Because germline and somatic mutations in the genes that are responsible for DNA methylation are infrequent in malignancies, additional mechanisms must be considered. Mycoplasmas spp., including Mycoplasma hyorhinis, efficiently colonize human cells and may serve as a vehicle for delivery of enzymatically active microbial proteins into the intracellular milieu. Here, we performed, for the first time, genome-wide and individual gene mapping of methylation marks generated by the M. hyorhinis CG- and GATC-specific DNA cytosine methyltransferases (MTases) in human cells. Our results demonstrated that, upon expression in human cells, MTases readily translocated to the cell nucleus. In the nucleus, MTases selectively and efficiently methylated the host genome at the DNA sequence sites free from pre-existing endogenous methylation, including those in a variety of cancer-associated genes. We also established that mycoplasma is widespread in colorectal cancers, suggesting that either the infection contributed to malignancy onset or, alternatively, that tumors provide a favorable environment for mycoplasma growth. In the human genome, ∼11% of GATC sites overlap with CGs (e.g., CGAT^mCG); therefore, the methylated status of these sites can be perpetuated by human DNMT1. Based on these results, we now suggest that the GATC-specific methylation represents a novel type of infection-specific epigenetic mark that originates in human cells with a previous exposure to infection. Overall, our findings unveil an entirely new panorama of interactions between the human microbiome and epigenome with a potential impact in disease etiology.     

Genome‐wide DNA methylation signatures of infection status in Trinidadian guppies (Poecilia reticulata)

Epigenetic modification, especially DNA methylation, can play an important role in mediating gene regulatory response to environmental stressors and may be a key process affecting phenotypic plasticity and adaptation. Parasites are potent stressors with profound physiological and ecological effects on their hosts, yet it remains unclear how parasites influence host methylation patterns. Here, we used a well‐studied host–parasite system, the guppy Poecilia reticulata and its ectoparasitic monogenean Gyrodactylus turnbulli to gain mechanistic insight into the dynamics of DNA methylation in host–parasite interactions. To explore this, we quantitatively measured genome‐wide DNA methylation in guppy skin tissue using reduced representation bisulphite sequencing and characterized differential methylation patterns in guppies during distinct phases of infection. We identified 365, 313, and 741 differentially methylated regions (DMRs) between infected and control fish in early infection, peak infection and recovery phases, respectively. The magnitude of the methylation difference was moderate in DMRs, with an average of 29% (early infection), 27% (peak infection) and 30% (recovery) differential methylation per DMR. Approximately 50% of DMRs overlapped with CpG islands, and over half of the DMRs overlapped with gene bodies, several of which encode proteins relevant to immune response. These findings provide the first evidence of an epigenetic signature of infection by ectoparasites and demonstrate the changing relationship between epigenetic variation and immune response in distinct phases of infection.     

Influence of in vitro manipulation on the stability of methylation patterns in the Snurf/Snrpn-imprinting region in mouse embryonic stem cells

Recent work on embryonic stem (ES) cells showed that stem cell-derived tissues and embryos, cloned from ES cell nuclei, often fail to maintain epigenetic states of imprinted genes. This deregulation is frequently associated with in vitro manipulations and culture conditions which might affect the cells potential to develop into normal fetuses. Usually, epigenetic instability is reported in differentially methylated regions of mostly growth-related imprinted genes. However, little is known about the epigenetic stability of genes that function late in organogenesis. Hence, we set out to investigate the epigenetic stability of neuronal genes and analyzed DNA methylation patterns in the Snurf/Snrpn imprinted cluster in several cultured mouse ES cell lines. We also determined the effects of in vitro stress factors such as consecutive passaging, trypsination, mechanical handling, single cell cloning, centrifugation, staurosporine-induced neurogenesis and the insertion of viral (foreign) DNA into the host genome. Intriguingly, none of these in vitro manipulations interfered with the stability of the methylation patterns in the analyzed neuronal genes. These data imply that, in contrast to growth-related genes like Igf2, H19, Igf2r or Grb10, the methylation imprints of the analyzed neuronal genes in the Snurf/Snrpn cluster may be particularly stable in manipulated ES cells.     

Differential DNA methylation identified in the blood and retina of AMD patients

Age-related macular degeneration (AMD) is a major cause of blindness in the western world. While genetic studies have linked both common and rare variants in genes involved in regulation of the complement system to increased risk of development of AMD, environmental factors, such as smoking and nutrition, can also significantly affect the risk of developing the disease and the rate of disease progression. Since epigenetics has been implicated in mediating, in part, the disease risk associated with some environmental factors, we investigated a possible epigenetic contribution to AMD. We performed genome-wide DNA methylation profiling of blood from AMD patients and controls. No differential methylation site reached genome-wide significance; however, when epigenetic changes in and around known GWAS-defined AMD risk loci were explored, we found small but significant DNA methylation differences in the blood of neovascular AMD patients near age-related maculopathy susceptibility 2 (ARMS2), a top-ranked GWAS locus preferentially associated with neovascular AMD. The methylation level of one of the CpG sites significantly correlated with the genotype of the risk SNP rs10490924, suggesting a possible epigenetic mechanism of risk. Integrating genome-wide DNA methylation analysis of retina samples with and without AMD together with blood samples, we further identified a consistent, replicable change in DNA methylation in the promoter region of protease serine 50 (PRSS50). These methylation changes may identify sites in novel genes that are susceptible to non-genetic factors known to contribute to AMD development and progression.     

Characterization of Two Novel Gammapapillomaviruses,HPV179 and HPV184, Isolated from Common Warts of a Renal-Transplant Recipient

Gammapapillomavirus (Gamma-PV) is a diverse and rapidly expanding PV-genus, currently consisting of 76 fully characterized human papillomavirus (HPV) types. In this study, DNA genomes of two novel HPV types, HPV179 and HPV184, obtained from two distinct facial verrucae vulgares specimens of a 64 year-old renal-transplant recipient, were fully cloned, sequenced and characterized. HPV179 and HPV184 genomes comprise 7,228-bp and 7,324-bp, respectively, and contain four early (E1, E2, E6 and E7) and two late genes (L1 and L2); the non-coding region is typically positioned between L1 and E6 genes. Phylogenetic analysis of the L1 nucleotide sequence placed both novel types within the Gamma-PV genus: HPV179 was classified as a novel member of species Gamma-15, additionally containing HPV135 and HPV146, while HPV184 was classified as a single member of a novel species Gamma-25. HPV179 and HPV184 type-specific quantitative real-time PCRs were further developed and used in combination with human beta-globin gene quantitative real-time PCR to determine the prevalence and viral load of the novel types in the patient’s facial warts and several follow-up skin specimens, and in a representative collection, a total of 569 samples, of HPV-associated benign and malignant neoplasms, hair follicles and anal and oral mucosa specimens obtained from immunocompetent individuals. HPV179 and HPV184 viral loads in patients’ facial warts were estimated to be 2,463 and 3,200 genome copies per single cell, respectively, suggesting their active role in the development of common warts in organ-transplant recipients. In addition, in this particular patient, both novel types had established a persistent infection of the skin for more than four years. Among immunocompetent individuals, HPV179 was further detected in low-copy numbers in a few skin specimens, indicating its cutaneous tissue tropism, while HPV184 was further detected in low-copy numbers in one mucosal and a few skin specimens, suggesting its dual tissue tropism.     

Sequential changes in genome-wide DNA methylation status during adipocyte differentiation

DNA methylation is an epigenetic mark on the mammalian genome. There are numerous tissue-dependent and differentially methylated regions (T-DMRs) in the unique sequences distributed throughout the genome. To determine the epigenetic changes during adipocyte differentiation, we investigated the sequential changes in DNA methylation status of 3T3-L1 cells at the growing, confluent, postconfluent and mature adipocyte cell stages. Treatment of 3T3-L1 cells with 5-aza-2′-deoxycytidine inhibited differentiation in a stage-dependent manner, supporting the idea that formation of accurate DNA methylation profile, consisting of methylated and unmethylated T-DMRs, may be involved in differentiation. Analysis by methylation-sensitive quantitative real-time PCR of the 65 known T-DMRs which contain NotI sites detected 8 methylations that changed during differentiation, and the changes in the patterns of these methylations were diverse, confirming that the differentiation process involves epigenetic alteration at the T-DMRs. Intriguingly, the dynamics of the methylation change vary depending on the T-DMRs and differentiation stages. Restriction landmark genomic scanning detected 32 novel T-DMRs, demonstrating that differentiation of 3T3-L1 cells involves genome-wide epigenetic changes by temporal methylation/demethylation, in addition to maintenance of a static methylated/demethylated state, and both depend on differentiation stage.     

UHRF1 Induces Methylation of the TXNIP Promoter and Down-Regulates Gene Expression in Cervical Cancer

DNA methylation, and consequent down-regulation, of tumour suppressor genes occurs in response to epigenetic stimuli during cancer development. Similarly, human oncoviruses, including human papillomavirus (HPV), up-regulate and augment DNA methyltransferase (DNMT) and histone deacetylase (HDAC) activities, thereby decreasing tumour suppressor genes (TSGs) expression. Ubiquitin-like containing PHD and RING finger domain 1 (UHRF1), an epigenetic regulator of DNA methylation, is overexpressed in HPV-induced cervical cancers. Here, we investigated the role of UHRF1 in cervical cancer by knocking down its expression in HeLa cells using lentiviral-encoded short hairpin (sh)RNA and performing cDNA microarrays. We detected significantly elevated expression of thioredoxin-interacting protein (TXNIP), a known TSG, in UHRF1-knockdown cells, and this gene is hypermethylated in cervical cancer tissue and cell lines, as indicated by whole-genome methylation analysis. Up-regulation of UHRF1 and decreased TXNIP were further detected in cervical cancer by western blot and immunohistochemistry and confirmed by Oncomine database analysis. Using chromatin immunoprecipitation, we identified the inverted CCAAT domain-containing UHRF1-binding site in the TXNIP promoter and demonstrated UHRF1 knockdown decreases UHRF1 promoter binding and enhances TXNIP expression through demethylation of this region. TXNIP promoter CpG methylation was further confirmed in cervical cancer tissue by pyrosequencing and methylation-specific polymerase chain reaction. Critically, down-regulation of UHRF1 by siRNA or UHRF1 antagonist (thymoquinone) induces cell cycle arrest and apoptosis, and ubiquitin-specific protease 7 (USP7), which stabilises and promotes UHRF1 function, is increased by HPV viral protein E6/E7 overexpression. These results indicate HPV might induce carcinogenesis through UHRF1-mediated TXNIP promoter methylation, thus suggesting a possible link between CpG methylation and cervical cancer.     

The interaction of high and low-risk human papillomavirus genotypes increases the risk of developing genital warts: A population-based cohort study

Cervical cancer is among the most common type of cancers in women and is associated with human papillomavirus (HPV) infection. Genital warts are also reported to be linked with HPV infection types 11 and 6. In turn, clinical characteristics and morphological features of warts may be useful in the prediction of prognosis and in making treatment decisions. Thus, we have investigated the association of high and low-risk HPVs genotype with genital wart risk, as well as pathological and cytological information in cases recruited from a population-based cohort study of 1380 patients. Patients infected with HPV genotype 6 or 11 had an increased risk of having warts, with OR of 2.34 (95% CI: 0.955-5.737, P = 0.06). Also, this association was enhanced in the presence of high plus low-risk HPV for having genital wart (OR: 2.814; 95%: 1.208-6.55, P = 0.017) and cases having high-risk HPV (OR: 2.329; 95% CI: 1.029-5.269, P = 0.042). Moreover, we observed patients with genital warts having CIN2/3, indicating the importance of informing the physician to the patient to prevent more severe lesions. Our data demonstrated that patients with both low/high-risk HPV types had an increased risk of developing genital warts and persistent infection with HPV was a necessary precursor for the increase in cervical lesions.     

Allele-Specific Transcriptome and Methylome Analysis Reveals Stable Inheritance and Cis-Regulation of DNA Methylation in Nasonia

Gene expression divergence between closely related species could be attributed to both cis- and trans- DNA sequence changes during evolution, but it is unclear how the evolutionary dynamics of epigenetic marks are regulated. In eutherian mammals, biparental DNA methylation marks are erased and reset during gametogenesis, resulting in paternal or maternal imprints, which lead to genomic imprinting. Whether DNA methylation reprogramming exists in insects is not known. Wasps of the genus Nasonia are non-social parasitoids that are emerging as a model for studies of epigenetic processes in insects. In this study, we quantified allele-specific expression and methylation genome-wide in Nasonia vitripennis and Nasonia giraulti and their reciprocal F1 hybrids. No parent-of-origin effect in allelic expression was found for >8,000 covered genes, suggesting a lack of genomic imprinting in adult Nasonia. As we expected, both significant cis- and trans- effects are responsible for the expression divergence between N. vitripennis and N. giraulti. Surprisingly, all 178 differentially methylated genes are also differentially methylated between the two alleles in F1 hybrid offspring, recapitulating the parental methylation status with nearly 100% fidelity, indicating the presence of strong cis-elements driving the target of gene body methylation. In addition, we discovered that total and allele-specific expression are positively correlated with allele-specific methylation in a subset of the differentially methylated genes. The 100% cis-regulation in F1 hybrids suggests the methylation machinery is conserved and DNA methylation is targeted by cis features in Nasonia. The lack of genomic imprinting and parent-of-origin differentially methylated regions in Nasonia, together with the stable inheritance of methylation status between generations, suggests either a cis-regulatory motif for methylation at the DNA level or highly stable inheritance of an epigenetic signal in Nasonia.