Comprehensive genome methylation and whole genome expression analysis in penile carcinoma: Uncovering new molecular markers

Background: Penile carcinoma (PeCa) is frequently associated with high morbidity rates. Unlikely of the vast majority of tumors, there is no molecular markers described that are able to assist in diagnosis and prognosis or with potential to be therapeutic targets in PeCa. Patients and methods: DNA methylation status (244K Human DNA Methylation Microarray platform, Agilent Technologies) and large-scale expression analysis (4x44K Whole Human Genome Microarray, Agilent Technologies) were performed in 35 and 37 PeCa, respectively. Quantitative bisulfite pyrosequencing (qBP) and RT-qPCR were used to validate the findings in 93 samples. HPV status was assessed using the Linear Array HPV Genotyping kit (Roche Molecular Diagnostics, CA, USA). Results: Methylome analysis revealed 171 hypermethylated and 449 hypomethylated CpGs sites and the transcriptome profiling showed 2986 down- and 2817 over-expressed genes. HPV positivity was found in 32.7% of the cases, mainly the HPV16. The integrative analysis in 32 PeCa revealed a panel of 96 genes with inverse correlation between methylation and gene expression levels. The CpG hypermetlylation and gene downexpression, was confirmed for TWIST1, RSOP2, SOX3, SOX17, CD133, OTX2, HOXA3 and MEIS. In addition, BIRC5, DNMT1 and DNMT3B presented low levels of methylation and overexpression. The comparison of the results with clinical findings revealed that LIN28A, NKX2.2, NKX2.3, LHX5, BDNF, FOXA1 and CDX2 were associated with poor prognosis features. Conclusion: Putative prognostic markers were detected revealing that DNA methylation modulates the expression of several genes in PeCa. These data may prove instrumental for biomarker discovery in clinics and molecular epidemiology of PeCa.

Genome-wide methylation and transcriptome analysis in penile carcinoma: uncovering new molecular markers

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Whole genome methylation array analysis reveals new aspects in Balkan endemic nephropathy etiology

Background Balkan endemic nephropathy (BEN) represents a chronic progressive interstitial nephritis in striking correlation with uroepithelial tumours of the upper urinary tract. The disease has endemic distribution in the Danube river regions in several Balkan countries. DNA methylation is a primary epigenetic modification that is involved in major processes such as cancer, genomic imprinting, gene silencing, etc. The significance of CpG island methylation status in normal development, cell differentiation and gene expression is widely recognized, although still stays poorly understood. Methods We performed whole genome DNA methylation array analysis on DNA pool samples from peripheral blood from 159 affected individuals and 170 healthy individuals. This technique allowed us to determine the methylation status of 27 627 CpG islands throughout the whole genome in healthy controls and BEN patients. Thus we obtained the methylation profile of BEN patients from Bulgarian and Serbian endemic regions. Results Using specifically developed software we compared the methylation profiles of BEN patients and corresponding controls and revealed the differently methylated regions. We then compared the DMRs between all patient-control pairs to determine common changes in the epigenetic profiles. SEC61G, IL17RA, HDAC11 proved to be differently methylated throughout all patient-control pairs. The CpG islands of all 3 genes were hypomethylated compared to controls. This suggests that dysregulation of these genes involved in immunological response could be a common mechanism in BEN pathogenesis in both endemic regions and in both genders. Conclusion Our data propose a new hypothesis that immunologic dysregulation has a place in BEN etiopathogenesis. Keywords: Epigenetics; Whole genome array analysis; Balkan endemic nephropathy

Understanding the role of transforming growth factor beta signalling and epigenomics in osteoarthritis

Osteoarthritis (OA) is the most common form of arthritis with a high socioeconomic burden, with an incompletely understood etiology. Evidence suggests a role for the transforming growth factor beta (TGF-ß) signalling pathway and epigenomics in OA. The aim of this thesis was to understand the involvement of the TGF-ß pathway in OA and to determine the DNA methylation patterns of OA-affected cartilage as compared to the OA-free cartilage. First, I found that a common SNP in the BMP2 gene, a ligand in the Bone morphogenetic protein (BMP) subunit of TGF-ß pathway, was associated with OA in the Newfoundland population. I also showed a genetic association between SMAD3 - a signal transducer in the TGF-ß subunit of the TGF-ß signalling pathway - and the total radiographic burden of OA. I further demonstrated that SMAD3 is over-expressed in OA cartilage, suggesting an over activation of the TGF-ß signalling in OA. Next, I examined the connection of these genes in the regulation of matrix metallopeptidase 13 (MMP13) - an enzyme known to destroy extracellular matrix in OA cartilage - in the context of the TGF-ß signalling. The analyses showed that TGF-ß, MMP13, and SMAD3 were overexpressed in OA cartilage, whereas the expression of BMP2 was significantly reduced. The expression of TGF-ß was positively correlated with that of SMAD3 and MMP13, suggesting that TGF-ß signalling is involved in up-regulation of MMP13. This regulation, however, appears not to be controlled by SMAD3 signals, possibly due to the involvement of collateral signalling, and to be suppressed by BMP regulation in healthy cartilage, whose levels were reduced in end-stage OA. In a genome-wide DNA methylation analysis, I reported CpG sites differentially methylated in OA and showed that the cartilage methylome has a potential to distinguish between OA-affected and non-OA cartilage. Functional clustering analysis of the genes harbouring differentially methylated loci revealed that they are enriched in the skeletal system morphogenesis pathway, which could be a potential candidate for further OA studies. Overall, the findings from the present thesis provide evidence that the TGF-ß signalling pathway is associated with the development of OA, and epigenomics might be involved as a potential mechanism in OA.

State of genomics and epigenomics research in the perspective of HIV cure.

PURPOSE OF REVIEW: One of the seven key scientific priorities identified in the road map on HIV cure research is to 'determine the host mechanisms that control HIV replication in the absence of therapy'. This review summarizes the recent work in genomics and in epigenetic control of viral replication that is relevant for this mission. RECENT FINDINGS: New technologies allow the joint analysis of host and viral transcripts. They identify the patterns of antisense transcription of the viral genome and its role in gene regulation. High-throughput studies facilitate the assessment of integration at the genome scale. Integration site, orientation and host genomic context modulate the transcription and should also be assessed at the level of single cells. The various models of latency in primary cells can be followed using dynamic study designs to acquire transcriptome and proteome data of the process of entry, maintenance and reactivation of latency. Dynamic studies can be applied to the study of transcription factors and chromatin modifications in latency and upon reactivation. SUMMARY: The convergence of primary cell models of latency, new high-throughput quantitative technologies applied to the study of time series and the identification of compounds that reactivate viral transcription bring unprecedented precision to the study of viral latency.

Low relative humidity triggers RNA-directed de novo DNA methylation and suppression of genes controlling stomatal development

Environmental cues influence the development of stomata on the leaf epidermis, and allow plants to exert plasticity in leaf stomatal abundance in response to the prevailing growing conditions. It is reported that Arabidopsis thaliana ‘Landsberg erecta’ plants grown under low relative humidity have a reduced stomatal index and that two genes in the stomatal development pathway, SPEECHLESS and FAMA, become de novo cytosine methylated and transcriptionally repressed. These environmentally-induced epigenetic responses were abolished in mutants lacking the capacity for de novo DNA methylation, for the maintenance of CG methylation, and in mutants for the production of short-interfering non-coding RNAs (siRNAs) in the RNA-directed DNA methylation pathway. Induction of methylation was quantitatively related to the induction of local siRNAs under low relative humidity. Our results indicate the involvement of both transcriptional and post-transcriptional gene suppression at these loci in response to environmental stress. Thus, in a physiologically important pathway, a targeted epigenetic response to a specific environmental stress is reported and several of its molecular, mechanistic components are described, providing a tractable platform for future epigenetics experiments. Our findings suggest epigenetic regulation of stomatal development that allows for anatomical and phenotypic plasticity, and may help to explain at least some of the plant’s resilience to fluctuating relative humidity.

DNA repair gene excision repair cross complementing-group 1 (ERCC1) in head and neck squamous cell carcinoma: analysis of methylation and polymorphism (G19007A), protein expression and association with epidemiological and clinicopathological factors

Aims: To evaluate the associations of excision repair cross complementing-group 1 (ERCC1) (DNA repair protein) (G19007A) polymorphism, methylation and immunohistochemical expression with epidemiological and clinicopathological factors and with overall survival in head and neck squamous cell carcinoma (HNSCC) patients. Methods and results: The study group comprised 84 patients with HNSCC who underwent surgery and adjuvant radiotherapy without chemotherapy. Bivariate and multivariate analyses were used. The allele A genotype variant was observed in 79.8% of the samples, GG in 20.2%, GA in 28.6% and AA in 51.2%. Individuals aged more than 45 years had a higher prevalence of the allelic A variant and a high (83.3%) immunohistochemical expression of ERCC1 protein [odds ratio (OR) = 4.86, 95% confidence interval (CI): 1.2-19.7, P = 0.027], which was also high in patients with advanced stage (OR= 5.04, 95% CI: 1.07-23.7, P = 0.041). Methylated status was found in 51.2% of the samples, and was higher in patients who did not present distant metastasis (OR = 6.67, 95% CI: 1.40-33.33, P = 0.019) and in patients with advanced stage (OR = 5.04, 95% CI: 1.07-23.7, P = 0.041). At 2 and 5 years, overall survival was 55% and 36%, respectively (median = 30 months). Conclusion: Our findings may reflect a high rate of DNA repair due to frequent tissue injury during the lifetime of these individuals, and also more advanced disease presentation in this population with worse prognosis.

Epigenetic regulation of somatic angiotensin-converting enzyme by DNA methylation and histone acetylation

Somatic angiotensin-converting enzyme (sACE) is crucial in cardiovascular homeostasis and displays a tissue-specific profile. Epigenetic patterns modulate genes expression and their alterations were implied in pathologies including hypertension. However, the influence of DNA methylation and chromatin condensation state on the expression of sACE is unknown. We examined whether such epigenetic mechanisms could participate in the control of sACE expression in vitro and in vivo. We identified two CpG islands in the human ace-1 gene 3 kb proximal promoter region. Their methylation abolished the luciferase activity of ace-1 promoter/reporter constructs transfected into human liver (HepG2), colon (HT29), microvascular endothelial (HMEC-1) and lung (SUT) cell lines (p < 0.001). Bisulphite sequencing revealed a cell-type specific basal methylation pattern of the ace-1 gene -1,466/+25 region. As assessed by RT-qPCR, inhibition of DNA methylation by 5-aza-2'-deoxycytidine and/or of histone deacetylation by trichostatin A highly stimulated sACE mRNA expression cell-type specifically (p < 0.001 vs. vehicle treated cells). In the rat, in vivo 5-aza-cytidine injections demethylated the ace-1 promoter and increased sACE mRNA expression in the lungs and liver (p = 0.05), but not in the kidney. In conclusion, the expression level of somatic ACE is modulated by CpG-methylation and histone deacetylases inhibition. The basal methylation pattern of the promoter of the ace-1 gene is cell-type specific and correlates to sACE transcription. DNMT inhibition is associated with altered methylation of the ace-1 promoter and a cell-type and tissue-specific increase of sACE mRNA levels. This study indicates a strong influence of epigenetic mechanisms on sACE expression.

DNA methylation and adaptive response in forest tree species

Progressive increase of temperatures as well as longer seasonal drought periods revealed by climate studies correspond to fast environmental changes that forest species face with their actual genetic background. Natural selective processes cannot develop an adaptive response within this time frame. Thus the capability of forest tree species to adapt to the new environments will depend on their genetic background, but also rely on their phenotypic plasticity. Several reports have shown the involvement of epigenetic modifiers as the basis of the phenotypic plasticity, and in particular to the adaptation to abiotic stresses. DNA methylation (methylation of cytosine residues)is one the most important epigenetic modification in eukaryotes. Itis involved in specific biological processes such as gene transcription regulation, gene silencing, mobile element control or genome imprinting.Therefore, there is a great interest in analyzing cytosine methylation levels and distribution within the genome

A link between DNA methylation and epigenetic silencing in transgenic Volvox carteri

Epigenetic silencing of foreign genes introduced into plants poses an unsolved problem for transgenic technology. Here we have used the simple multicellular green alga Volvox carteri as a model to analyse the relation of DNA methylation to transgenic silencing. Volvox DNA contains on average 1.1% 5-methylcytosine and 0.3% N6-methyladenine, as revealed by electrospray mass spectrometry and phosphoimaging of chromatographically separated 32P-labelled nucleotides. In two nuclear transformants of V.carteri, produced in 1993 by biolistic bombardment with a foreign arylsulphatase gene (C-ars), the transgene is still expressed in one (Hill 181), but not in the other (Hill 183), after an estimated 500–1000 generations. Each transformant clone contains multiple intact copies of C-ars, most of them integrated into the genome as tandem repeats. When the bisulphite genomic sequencing protocol was applied to examine two select regions of transgenic C-ars, we found that the inactivated copies (Hill 183) exhibited a high-level methylation (40%) of CpG dinucleotides, whereas the active copies (Hill 181) displayed low-level (7%) CpG methylation. These are average values from 40 PCR clones sequenced from each DNA strand in the two portions of C-ars. The observed correlation of CpG methylation and transgene inactivation in a green alga will be discussed in the light of transcriptional silencing.

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

A novel variant genotype C of hepatitis B virus identified in isolates from Australian Aborigines: complete genome sequence and phylogenetic relatedness

There have been no reports of DNA sequences of hepatitis B virus (HBV) strains from Australian Aborigines, although the hepatitis B surface antigen (HBsAg) was discovered among them. To investigate the characteristics of DNA sequences of HBV strains from Australian Aborigines, the complete nucleotide sequences of HBV strains were determined and subjected to molecular evolutionary analysis. Serum samples positive for HBsAg were collected from five Australian Aborigines. Phylogenetic analysis of the five complete nucleotide sequences compared with DNA sequences of 54 global HBV isolates from international databases revealed that three of the five were classified into genotype D and were most closely related in terms of evolutionary distance to a strain isolated from a healthy blood donor in Papua New Guinea. Two of the five were classified into a novel variant genotype C, which has not been reported previously, and were closely related to a strain isolated from Polynesians, particularly in the X and Core genes. These two strains of variant genotype C differed from known genotype C strains by 5.9-7.4% over the complete nucleotide sequence and 4.0-5.6 % in the small-S gene, and had residues Arg(122), Thr(127) and Lys(160) characteristic of serotype ayw3, which have not been reported previously in genotype C. In conclusion, this is the first report of the characteristics of complete nucleotide sequences of HBV from Australian Aborigines. These results contribute to the investigation of the worldwide spread of HBV, the relationship between serotype and genotype and the ancient common origin of Australian Aborigines.

Selection of patients with germline MLH1 mutated Lynch syndrome by determination of MLH1 methylation and BRAF mutation.

Lynch syndrome is one of the most common hereditary colorectal cancer (CRC) syndrome and is caused by germline mutations of MLH1, MSH2 and more rarely MSH6, PMS2, MLH3 genes. Whereas the absence of MSH2 protein is predictive of Lynch syndrome, it is not the case for the absence of MLH1 protein. The purpose of this study was to develop a sensitive and cost effective algorithm to select Lynch syndrome cases among patients with MLH1 immunohistochemical silencing. Eleven sporadic CRC and 16 Lynch syndrome cases with MLH1 protein abnormalities were selected. The BRAF c.1799T> A mutation (p.Val600Glu) was analyzed by direct sequencing after PCR amplification of exon 15. Methylation of MLH1 promoter was determined by Methylation-Sensitive Single-Strand Conformation Analysis. In patients with Lynch syndrome, there was no BRAF mutation and only one case showed MLH1 methylation (6%). In sporadic CRC, all cases were MLH1 methylated (100%) and 8 out of 11 cases carried the above BRAF mutation (73%) whereas only 3 cases were BRAF wild type (27%). We propose the following algorithm: (1) no further molecular analysis should be performed for CRC exhibiting MLH1 methylation and BRAF mutation, and these cases should be considered as sporadic CRC; (2) CRC with unmethylated MLH1 and negative for BRAF mutation should be considered as Lynch syndrome; and (3) only a small fraction of CRC with MLH1 promoter methylation but negative for BRAF mutation should be true Lynch syndrome patients. These potentially Lynch syndrome patients should be offered genetic counselling before searching for MLH1 gene mutations.