Methylation of NR3C1 is related to maternal PTSD, parenting stress and maternal medial prefrontal cortical activity in response to child separation among mothers with histories of violence exposure.

Prior research has shown that mothers with Interpersonal violence-related posttraumatic stress disorder (IPV-PTSD) report greater difficulty in parenting their toddlers. Relative to their frequent early exposure to violence and maltreatment, these mothers display dysregulation of their hypothalamic pituitary adrenal axis (HPA-axis), characterized by hypocortisolism. Considering methylation of the promoter region of the glucocorticoid receptor gene NR3C1 as a marker for HPA-axis functioning, with less methylation likely being associated with less circulating cortisol, the present study tested the hypothesis that the degree of methylation of this gene would be negatively correlated with maternal IPV-PTSD severity and parenting stress, and positively correlated with medial prefrontal cortical (mPFC) activity in response to video-stimuli of stressful versus non-stressful mother-child interactions. Following a mental health assessment, 45 mothers and their children (ages 12-42 months) participated in a behavioral protocol involving free-play and laboratory stressors such as mother-child separation. Maternal DNA was extracted from saliva. Interactive behavior was rated on the CARE-Index. During subsequent fMRI scanning, mothers were shown films of free-play and separation drawn from this protocol. Maternal PTSD severity and parenting stress were negatively correlated with the mean percentage of methylation of NR3C1. Maternal mPFC activity in response to video-stimuli of mother-child separation versus play correlated positively to NR3C1 methylation, and negatively to maternal IPV-PTSD and parenting stress. Among interactive behavior variables, child cooperativeness in play was positively correlated with NR3C1 methylation. Thus, the present study is the first published report to our knowledge, suggesting convergence of behavioral, epigenetic, and neuroimaging data that form a psychobiological signature of parenting-risk in the context of early life stress and PTSD.

Expression of the hyaluronan-mediated motility receptor RHAMM in tumor budding cells identifies aggressive colorectal cancers.

Expression of the hyaluronan-mediated motility receptor (RHAMM, CD168) predicts adverse clinicopathological features and decreased survival for colorectal cancer (CRC) patients. Using full tissue sections, we investigated the expression of RHAMM in tumor budding cells of 103 primary CRCs to characterize the biological processes driving single-cell invasion and early metastatic dissemination. RHAMM expression in tumor buds was analyzed with clinicopathological data, molecular features and survival. Tumor budding cells at the invasive front of CRC expressed RHAMM in 68% of cases. Detection of RHAMM-positive tumor budding cells was significantly associated with poor survival outcome (P = .0312), independent of TNM stage and adjuvant therapy in multivariate analysis (P = .0201). RHAMM-positive tumor buds were associated with frequent lymphatic invasion (P = .0007), higher tumor grade (P = .0296), and nodal metastasis (P = .0364). Importantly, the prognostic impact of RHAMM expression in tumor buds was maintained independently of the number of tumor buds found in an individual case (P = .0246). No impact of KRAS/BRAF mutation, mismatch repair deficiency and CpG island methylation was observed. RHAMM expression identifies an aggressive subpopulation of tumor budding cells and is an independent adverse prognostic factor for CRC patients. These data support ongoing efforts to develop RHAMM as a target for precision therapy.

H19 lncRNA alters DNA methylation genome wide by regulating S-adenosylhomocysteine hydrolase

DNA methylation is essential for mammalian development and physiology. Here we report that the developmentally regulated H19 lncRNA binds to and inhibits S-adenosylhomocysteine hydrolase (SAHH), the only mammalian enzyme capable of hydrolysing S-adenosylhomocysteine (SAH). SAH is a potent feedback inhibitor of S-adenosylmethionine (SAM)-dependent methyltransferases that methylate diverse cellular components, including DNA, RNA, proteins, lipids and neurotransmitters. We show that H19 knockdown activates SAHH, leading to increased DNMT3B-mediated methylation of an lncRNA-encoding gene Nctc1 within the Igf2-H19-Nctc1 locus. Genome-wide methylation profiling reveals methylation changes at numerous gene loci consistent with SAHH modulation by H19. Our results uncover an unanticipated regulatory circuit involving broad epigenetic alterations by a single abundantly expressed lncRNA that may underlie gene methylation dynamics of development and diseases and suggest that this mode of regulation may extend to other cellular components.

An Sp1/Sp3 site polymorphism associated with hypermethylation of the candidate tumor suppressor gene RIL in cancer

Gene silencing due to promoter methylation is an alternative to mutations and deletions, which inactivate tumor suppressor genes (TSG) in cancer. We identified RIL by Methylated CpG Island Amplification technique as a novel aberrantly methylated gene. RIL is expressed in normal tissues and maps to the 5q31 region, frequently deleted in leukemias. We found methylation of RIL in 55/80 (69%) cancer cell lines, with highest methylation in leukemia and colon. We also observed methylation in 46/80 (58%) primary tumors, whereas normal tissues showed substantially lower degrees of methylation. RIL expression was lost in 13/16 cancer cell lines and was restored by demethylating agent. Screening of 38 cell lines and 13 primary cancers by SSCP revealed no mutations in RIL, suggesting that methylation and LOH are the primary inactivation mechanisms. Stable transfection of RIL into colorectal cancer cells resulted in reduction in cell growth, clonogenicity, and increased apoptosis upon UVC treatment, suggesting that RIL is a good candidate TSG. ^ In searching for a cause of RIL hypermethylation, we identified a 12-bp polymorphic sequence around the transcription start site of the gene that creates a long allele containing 3CTC repeat. Evolutionary studies suggested that the long allele appeared late in evolution due to insertion. Using bisulfite sequencing, in cancers heterozygous for RIL, we found that the short allele is 4.4-fold more methylated than the long allele (P = 0.003). EMSA results suggested binding of factor(s) to the inserted region of the long allele, but not to the short. EMSA mutagenesis and competition studies, as well as supershifts using nuclear extracts or recombinant Sp1 strongly indicated that those DNA binding proteins are Sp1 and Sp3. Transient transfections of RIL allele-specific expression constructs showed less than 2-fold differences in luciferase activity, suggesting no major effects of the additional Sp1 site on transcription. However, stable transfection resulted in 3-fold lower levels of transcription from the short allele 60 days post-transfection, consistent with the concept that the polymorphic Sp1 site protects against time-dependent silencing. Thus, an insertional polymorphism in the RIL promoter creates an additional Sp1/Sp3 site, which appears to protect it from silencing and methylation in cancer. ^

Involvement of histone H3 methylation in Tup1-Ssn6 repression

The Tup1-Ssn6 complex regulates the expression of diverse classes of genes in Saccharomyces cerevisiae including those regulated by mating type, DNA damage, glucose, and anaerobic stress. The complex is recruited to target genes by sequence-specific repressor proteins. Once recruited to particular promoters, it is not completely clear how it functions to block transcription. Repression probably occurs through interactions with both the basal transcriptional machinery and components of chromatin. Tup1 interactions with chromatin are strongly influenced by acetylation of histories H3 and H4. Tup1 binds to underacetylated histone tails and requires multiple histone deacetylases (HDACs) for its repressive functions. Like acetylation, histone methylation is involved in regulation of gene expression. The possible role of histone methylation in Tup1 repression is not known. Here we examined possible roles of histone methyltransferases in Tup1-Ssn6 functions. We found that like other genes, Tup1-Ssn6 target genes exhibit increases in the levels of histone H3 lysine 4 methylation upon activation. However, deletion of individual or multiple histone methyltransferases (HMTs) and other SET-domain containing genes has no apparent effect on Tup1-Ssn6 mediated repression of a number of well-defined targets. Interestingly, we discovered that Ssn6 interacts with Set2. Since deletion of SET2 does not affect Tup1-Ssn6 repression, Ssn6 may utilize Set2 in other contexts to regulate gene repression. In order examine if the two components of the Tup1-Ssn6 complex have independent functions in the cell, we identified genes differentially expressed in tup1Δ and ssn6Δ mutants using DNA microarrays. Our data indicate that ∼4% of genes in the cell are regulated by Ssn6 independently of Tup1. In addition, expression of genes regulated by Tup1-Ssn6 seems to be differently affected by deletion of Ssn6 and deletion of Tup1, which indicates that these components might have separate functions. Our data shed new light on the classical view of Tup1-Ssn6 functions, and indicate that Ssn6 might have repressive functions as well. ^

Mechanism of methylation gene reactivation in human cancer cells

Epigenetic silencing of tumor suppressor genes by DNA hypermethylation at promoter regions is a common event in carcinogenesis and tumor progression. Abrogation of methylation and reversal of epigenetic silencing is a very potent way in cancer treatment. However, the reactivation mechanisms are poorly understood. In this study, we first developed a cell line model system named YB5, derived from SW48 cancer cell line, which bears one copy of stably integrated EGFP gene on Chromosome 1p31.1 region. The GFP gene expression is transcriptionally silenced due to the hypermethylated promoter CMV. However, the GFP expression can be restored using demethylating agent 5-aza-2' deoxycytidine (DAC), and detected by FACS and fluorescent microscopy. Using this system, we observed the heterogeneous reactivation induced by DAC treatment. After flow sorting, GFP negative cells exhibited similar level of incomplete demethylation compared to GFP positive cells on repetitive LINE1 element, tumor suppressor genes such as P16, CDH13, and RASSF1a, and CMV promoter as well. However, the local chromatin of CMV-GFP locus altered to an open structure marked by high H3 lysine 9 acetylation and low H3 lysine 27 tri-methylation in GFP positive cells, while the GFP negative cells retained mostly the original repressive marks. Thus, we concluded that DAC induced DNA hypomethylation alone does not directly determine the level of re-expression, and the resetting of the local chromatin structure under hypomethylation environment is required for gene reactivation. Besides, a lentivirus vector-based shRNA screening was performed using the YB5 system. Although it is the rare chance that vector lands in the neighboring region of GFP, we found that the exogenous vector DNA inserted into the upstream region of GFP gene locus led to the promoter demethylation and reactivated the silenced GFP gene. Thus, epigenetic state can be affected by changing of the adjacent nucleic acid sequences. Further, this hypermethylation silenced system was utilized for epigenetic drug screening. We have found that DAC combined with carboplatin would enhance the GFP% yield and increase expression of other tumor suppressor genes than DAC alone, and this synergistic effect may be related to DNA repair process. In summary, these studies reveal that reversing of methylation silencing requires coordinated alterations of DNA methylation, chromatin structure, and local microenvironment. ^

TAZ AS A REGULATOR OF MESENCHYMAL TRANSFORMATION AND CLINICAL AGGRESSIVENESS IN GLIOMAS

Glioblastoma multiforme (GBM) is an aggressive, high grade brain tumor. Microarray studies have shown a subset of GBMs with a mesenchymal gene signature. This subset is associated with poor clinical outcome and resistance to treatment. To establish the molecular drivers of this mesenchymal transition, we correlated transcription factor expression to the mesenchymal signature and identified transcriptional co-activator with PDZ-binding motif (TAZ) to be highly associated with the mesenchymal shift. High TAZ expression correlated with worse clinical outcome and higher grade. These data led to the hypothesis that TAZ is critical to the mesenchymal transition and aggressive clinical behavior seen in GBM. We investigated the expression of TAZ, its binding partner TEAD, and the mesenchymal marker FN1 in human gliomas. Western analyses demonstrated increased expression of TAZ, TEAD4, and FN1 in GBM relative to lower grade gliomas. We also identified CpG islands in the TAZ promoter that are methylated in most lower grade gliomas, but not in GBMs. TAZ-methylated glioma stem cell (GSC) lines treated with a demethylation agent showed an increase in mRNA and protein TAZ expression; therefore, methylation may be another novel way TAZ is regulated since TAZ is epigenetically silenced in tumors with a better clinical outcome. To further characterize the role of TAZ in gliomagenesis, we stably silenced or over-expressed TAZ in GSCs. Silencing of TAZ decreased invasion, self-renewal, mesenchymal protein expression, and tumor-initiating capacity. Over-expression of TAZ led to an increase in invasion, mesenchymal protein expression, mesenchymal differentiation, and tumor-initiating ability. These actions are dependent on TAZ interacting with TEAD since all these effects were abrogated with TAZ could not bind to TEAD. We also show that TAZ and TEAD directly bind to mesenchymal gene promoters. Thus, TAZ-TEAD interaction is critically important in the mesenchymal shift and in the aggressive clinical behavior of GBM. We identified TAZ as a regulator of the mesenchymal transition in gliomas. TAZ could be used as a biomarker to both estimate prognosis and stratify patients into clinically relevant subgroups. Since mesenchymal transition is correlated to tumor aggressiveness, strategies to target and inhibit TAZ-TEAD and the downstream gene targets may be warranted in alternative treatment.

RNA-SEQUENCING APPLICATIONS: GENE EXPRESSION QUANTIFICATION AND METHYLATOR PHENOTYPE IDENTIFICATION

My dissertation focuses on two aspects of RNA sequencing technology. The first is the methodology for modeling the overdispersion inherent in RNA-seq data for differential expression analysis. This aspect is addressed in three sections. The second aspect is the application of RNA-seq data to identify the CpG island methylator phenotype (CIMP) by integrating datasets of mRNA expression level and DNA methylation status. Section 1: The cost of DNA sequencing has reduced dramatically in the past decade. Consequently, genomic research increasingly depends on sequencing technology. However it remains elusive how the sequencing capacity influences the accuracy of mRNA expression measurement. We observe that accuracy improves along with the increasing sequencing depth. To model the overdispersion, we use the beta-binomial distribution with a new parameter indicating the dependency between overdispersion and sequencing depth. Our modified beta-binomial model performs better than the binomial or the pure beta-binomial model with a lower false discovery rate. Section 2: Although a number of methods have been proposed in order to accurately analyze differential RNA expression on the gene level, modeling on the base pair level is required. Here, we find that the overdispersion rate decreases as the sequencing depth increases on the base pair level. Also, we propose four models and compare them with each other. As expected, our beta binomial model with a dynamic overdispersion rate is shown to be superior. Section 3: We investigate biases in RNA-seq by exploring the measurement of the external control, spike-in RNA. This study is based on two datasets with spike-in controls obtained from a recent study. We observe an undiscovered bias in the measurement of the spike-in transcripts that arises from the influence of the sample transcripts in RNA-seq. Also, we find that this influence is related to the local sequence of the random hexamer that is used in priming. We suggest a model of the inequality between samples and to correct this type of bias. Section 4: The expression of a gene can be turned off when its promoter is highly methylated. Several studies have reported that a clear threshold effect exists in gene silencing that is mediated by DNA methylation. It is reasonable to assume the thresholds are specific for each gene. It is also intriguing to investigate genes that are largely controlled by DNA methylation. These genes are called “L-shaped” genes. We develop a method to determine the DNA methylation threshold and identify a new CIMP of BRCA. In conclusion, we provide a detailed understanding of the relationship between the overdispersion rate and sequencing depth. And we reveal a new bias in RNA-seq and provide a detailed understanding of the relationship between this new bias and the local sequence. Also we develop a powerful method to dichotomize methylation status and consequently we identify a new CIMP of breast cancer with a distinct classification of molecular characteristics and clinical features.

Methylation index of Branched Tetraethers (MBT) and Cyclization ratio of Branched Tetraethers (CBT), annual mean air temperature and soil pH of sediment core GeoB6518-1 in the Congo River basin

We analyzed the distribution of branched tetraether membrane lipids derived from soil bacteria in a marine sediment record that was recovered close to the Congo River outflow, and the results enabled us to reconstruct large-scale continental temperature changes in tropical Africa that span the past 25,000 years. Tropical African temperatures gradually increased from ~21° to 25°C over the last deglaciation, which is a larger warming than estimated for the tropical Atlantic Ocean. A direct comparison with sea-surface temperature estimates from the same core revealed that the land-sea temperature difference was, through the thermal pressure gradient, an important control on central African precipitation patterns.

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

Incidence and functional consequences of hMLH1 promoter hypermethylation in colorectal carcinoma

Inactivation of the genes involved in DNA mismatch repair is associated with microsatellite instability (MSI) in colorectal cancer. We report that hypermethylation of the 5′ CpG island of hMLH1 is found in the majority of sporadic primary colorectal cancers with MSI, and that this methylation was often, but not invariably, associated with loss of hMLH1 protein expression. Such methylation also occurred, but was less common, in MSI− tumors, as well as in MSI+ tumors with known mutations of a mismatch repair gene (MMR). No hypermethylation of hMSH2 was found. Hypermethylation of colorectal cancer cell lines with MSI also was frequently observed, and in such cases, reversal of the methylation with 5-aza-2′-deoxycytidine not only resulted in reexpression of hMLH1 protein, but also in restoration of the MMR capacity in MMR-deficient cell lines. Our results suggest that microsatellite instability in sporadic colorectal cancer often results from epigenetic inactivation of hMLH1 in association with DNA methylation.

Sequence motifs in adenoviral DNA block immune activation by stimulatory CpG motifs

Unmethylated CpG dinucleotides in particular base contexts (CpG-S motifs) are relatively common in bacterial DNA but are rare in vertebrate DNA. B cells and monocytes have the ability to detect such CpG-S motifs that trigger innate immune defenses with production of Th1-like cytokines. Despite comparable levels of unmethylated CpG dinucleotides, DNA from serotype 12 adenovirus is immune-stimulatory, but serotype 2 is nonstimulatory and can even inhibit activation by bacterial DNA. In type 12 genomes, the distribution of CpG-flanking bases is similar to that predicted by chance. However, in type 2 adenoviral DNA the immune stimulatory CpG-S motifs are outnumbered by a 15- to 30-fold excess of CpG dinucleotides in clusters of direct repeats or with a C on the 5′ side or a G on the 3′ side. Synthetic oligodeoxynucleotides containing these putative neutralizing (CpG-N) motifs block immune activation by CpG-S motifs in vitro and in vivo. Eliminating 52 of the 134 CpG-N motifs present in a DNA vaccine markedly enhanced its Th1-like function in vivo, which was increased further by the addition of CpG-S motifs. Thus, depending on the CpG motif, prokaryotic DNA can be either immune-stimulatory or neutralizing. These results have important implications for understanding microbial pathogenesis and molecular evolution and for the clinical development of DNA vaccines and gene therapy vectors.

In utero exposure to cigarette chemicals induces sex-specific disruption of one-carbon metabolism and DNA methylation in the human fetal liver

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated. Acknowledgements: We thank Ms Margaret Fraser, Ms Samantha Flannigan, and Dr Wing Yee Kwong for their expert assistance. The staff at Grampian NHS Pregnancy Counselling Service were essential for collecting fetuses. We thank Professor Geoffrey Hammond and Dr Marc Simard, University of British Colombia for helpful comments on the manuscript. Supported by grants as follows: Scottish Senior Clinical Fellowship (AJD); Chief Scientist Office (Scottish Executive, CZG/1/109 to PAF, & CZG/4/742 (PAF & PJOS); NHS Grampian Endowments 08/02 (PAF, SB & PJOS); the European Community’s Seventh Framework Programme (FP7/2007-2013) under grant agreement no 212885 (PAF & SMR); the Medical Research Council grants MR/L010011/1 (PAF & PJOS) and MR/K018310/1 (AJD). None of the funding bodies played any role in the design, collection, analysis, and interpretation of data, in the writing of the manuscript, nor in the decision to submit the manuscript for publication

Hepatic IGF1 DNA methylation is influenced by gender but not by intrauterine growth restriction in the young lamb

Acknowledgements: The authors would like to thank Chaz Mein and Eva Wozniak at the Genome Centre, Queen Mary University of London, who were integral to the development of the DNA methylation assays and processing of samples. Financial Support: This work was supported by Wellbeing of Women (D.J.C., grant number RTF318) and the Scottish Government’s Rural and Environment Science and Analytical Services Division including the Strategic Partnership for Animal Science Excellence (J.M.W., R.P.A., J.S.M. and C.L.A.).