Identification of a possible role of thymine DNA glycosylase (TDG) in epigenome maintenance

Thymine DNA glycosylase (TDG) was discovered as an enzyme capable of removing uracil (U) and thymine (T) from G/U and G/T mispairs, respectively. Owing to this ability, TDG was proposed to initiate restoration of C/G pairs at sites of cytosine or 5-methycytosine (5-meC) deamination. In addition to products of base deamination, the substrate spectrum of TDG covers a wide range of DNA base damages resulting from base oxidation and alkylation. TDG was also found to engage in physical and functional interactions with transcription factors, and more recent evidence supports additional interactions with the de novo DNA methyltransferases Dnmt3a and 3b in the context of gene transcription. Together with its biochemical properties, these observations suggest that TDG might be targeted to gene regulatory sequences as part of a macromolecular assembly to control their functional integrity. TDG may counteract the mutagenic effects of C and 5-meC deamination in CG-rich regions and/or be involved in the maintenance of CpG promoter methylation patterns. A tight regulation of CpG methylation at gene regulatory regions is critical for accurate gene expression, proper cellular differentiation and embryonic development. A somewhat surprising but in this context consistent finding was that, in contrast to other DNA glycosylases, TDG is essential for proper fetal development since a targeted knockout of the gene leads to embryonic lethality. To gain insights into the biological functions of TDG, we aimed to establish and apply biochemical fractionation procedures for high affinity purification and structural and functional characterization of TDG containing proteins complexes. The first part of the thesis was concerned with biochemical characterization of the protein interaction network of TDG in living mammalian cells. To this end, I applied different approaches allowing high affinity isolation of protein complexes from mammalian cells, such as the tandem affinity purification (TAP) method as well as immunoprecipitation of endogenous protein and of the TDGa isoform from TdgA overexpressing embryonic stem (ES) cells. These efforts, however, did not reveal any TDG interacting partners in subsequent mass spectrometry (MS) analyses. These results were surprising, as TDG was previously reported to interact with transcription factors and DNA methyltransferases. Remarkably, however, all previously identified protein interactors of TDG were discovered in screen with the respective partner proteins, and under conditions of simultaneous overexpression of both interacting proteins. The only proteins ever identified in screen with TDG were Sumo1 and Sumo3, which turned out to covalently modify the glycosylase. For this reason, we decided to pursue our search with classical cell fractionation experiments. We first did gel filtration experiments from total cell lysates and showed that TDG is indeed able to form distinct multiprotein complexes in undifferentiated mouse embryonic stem cells that may also contain the RNA helicase p68. Further subcellular fractionation experiments then revealed that TDG is present in all cell compartments, with a significant fraction of nuclear TDG being associated with chromatin, together with p68 and de novo DNA methyltransferases. Together with published findings, these results suggested that protein complexes containing TDG might act in a chromatin-associated context, at gene regulatory regions. The developmental phenotype of Tdg-/- knockout mice and the interactions of TDG with factors involved in developmental gene regulation (e.g. retinoic acid receptors RAR/RXR) implicate a function of TDG during early development and cell differentiation, at times governed by dynamic changes in gene expression, DNA methylation and histone modifications. Such changes have been studied using a well-established during in vitro differentiation of ES cells to lineage committed neuronal progenitors (NPs). We thus aimed to address the function of TDG as part of chromatin associated protein complexes during the process of retinoic acid induced differentiation of ES cells to NPs. In the second part of the thesis we made use of a this well-established in vitro differentiation system to examine the genome-wide localization of TDG to chromatin by TDG chromatin immunoprecipitation (ChIP) and to correlate TDG association to chromatin with gene expression and DNA methylation changes linked to cellular differentiation. TDG ChIP combined with high throughput sequencing showed that TDG associates with high preference to CpG islands in promoters of actively transcribed genes or genes poised for transcriptional activation. Such CpG rich sequences are normally unmethylated in mammalian genomes. Interestingly, we found TDG to localize to promoters of many genes controlling pluripotency (e.g. Oct4, Nanog) and developmental processes (e.g. Sfrp2, Tgfb2, Gata6), thus, supporting a function of TDG in cell differentiation and/or embryonic development. As different lines of circumstantial evidence have associated TDG with changes in CpG methylation following activation of hormone responsive gene promoters, we went on to further test genome-wide promoter methylation in Tdg+/- and Tdg-/- NPs making use of a combination of methylated DNA immunoprecipitation (MeDIP) and microarray technology. This showed that the loss of TDG does not affect global promoter DNA methylation. Nevertheless, there were a number of significant differences, suggesting that TDG might affect the CpG methylation pattern at some promoters. Also, owing to the limited resolution of the MeDIP method, however, we could not exclude an involvement of TDG in the control of DNA methylation of specific promoter CpGs. Additional bisulfite sequencing of promoters of TDG bound developmental genes (e.g. Sfrp2, Tgfb2) in NPs and differentiated mouse embryonic fibroblasts (MEFs) have indeed proved that loss of TDG affects local changes in DNA methylation at particular CpGs. Subsequent analysis of genome-wide gene expression profiles of ES cells and differentiated Tdg+/- and Tdg-/- NPs revealed that a limited number of genes (229) are differentially regulated in ES, whereas substantial differences in gene expression in were observed in NPs (1022 genes). This implicated a specific function of TDG in the regulation of cell differentiation triggered gene expression changes. Detailed analysis of the expression of the Pax6 gene, accurate regulation of which is essential for proper neuron development, showed that its promoter is bound by TDG and that its transcription is inappropriately regulated upon further differentiation of Tdg-/- NPs into the neuronal lineage. Whereas Tdg+/- NPs efficiently downregulated Pax6 (50x) and further differentiated into neuron-like cells, Tdg-/- NPs only partially downregulated Pax6 gene expression (6x) and underwent apoptosis at day 2 after plating in neuronal medium. This phenotype was complemented by expression of TDGa, clearly implicating TDG in the regulation of Pax6 expression during differentiation of ES cells to terminal neurons. We further observed misregulation of pluripotency genes (e.g. Oct4) regulated by TDG bound promoters during early differentiation of ES cells. In the absence of TDG, ES cells showed the tendency to enter spontaneous and/or RA induced differentiation, suggesting a role for TDG in the regulation of pluripotency. During RA induced differentiation we further observed the activation of the neuron specific gene Lrrtm2 exclusively in TDG proficient cells. In addition, ChIP experiments showed that transcription factors involved in the activation of the Lrrtm2 gene (e.g. COUP-TFI, RAR) are not recruited to the respective promoter in Tdg-/- cells, suggesting that TDG might act passively as a scaffold factor important for the recruitment of transcription factors to promoter regions. I set out to clarify the biological function of TDG by investigating its molecular interactions in mammalian cells. I found that TDG, as a DNA repair enzyme, associates tightly with chromatin, where it localizes with high preference to CpG island promoters of active genes and genes poised to be expressed. I also found that the loss of TDG causes misregulation of genes during cell differentiation and that this appears to be related to a function of TDG in establishing and/or maintaining CpG methylation pattern in gene regulatory sequences. These discoveries implicate a novel function of DNA repair, in the maintenance not only of the genome, but also the epigenome

Evaluation of Analytical Methods to Study Aquifer Properties with Pumping Tests in Coastal Aquifers with Numerical Modelling (Motril-Salobreña Aquifer)

Two pumping tests were performed in the unconfined Motril-Salobreña detrital aquifer in a 250 m-deep well 300 m from the coastline containing both freshwater and saltwater. It is an artesian well as it is in the discharge zone of this coastal aquifer. The two observation wells where the drawdowns are measured record the influence of tidal fluctuations, and the well lithological columns reveal high vertical heterogeneity in the aquifer. The Theis and Cooper-Jacob approaches give average transmissivity (T) and storage coefficient (S) values of 1460 m2 /d and 0.027, respectively. Other analytical solutions, modified to be more accurate in the boundary conditions found in coastal aquifers, provide similar T values to those found with the Theis and Cooper-Jacob methods, but give very different S values or could not estimate them. Numerical modelling in a synthetic model was applied to analyse the sensitivity of the Theis and Cooper-Jacob approaches to the usual boundary conditions in coastal aquifers. The T and S values calculated from the numerical modelling drawdowns indicate that the regional flow, variable pumping flows, and tidal effect produce an error of under 10 % compared to results obtained with classic methods. Fluids of different density (freshwater and saltwater) cause an error of 20 % in estimating T and of over 100 % in calculating S. The factor most affecting T and S results in the pumping test interpretation is vertical heterogeneity in sediments, which can produce errors of over 100 % in both parameters.This research has been financed by Project CGL2012-32892 (Ministerio de Economía y Competitividad of Spain) and by the Research Group Sedimentary Geology and Groundwater (RNM-369) of the Junta de Andalucía

Characterization of a fluvial aquifer at a range of depths and scales: the Triassic St Bees Sandstone Formation, Cumbria, UK

Fluvial sedimentary successions represent porous media that host groundwater and geothermal resources. Additionally, they overlie crystalline rocks hosting nuclear waste repositories in rift settings. The permeability characteristics of an arenaceous fluvial succession, the Triassic St Bees Sandstone Formation in England (UK), are described, from core-plug to well-test scale up to ~1 km depth. Within such lithified successions, dissolution associated with the circulation of meteoric water results in increased permeability (K~10−1–100 m/day) to depths of at least 150 m below ground level (BGL) in aquifer systems that are subject to rapid groundwater circulation. Thus, contaminant transport is likely to occur at relatively high rates. In a deeper investigation (> 150 m depth), where the aquifer has not been subjected to rapid groundwater circulation, well-test-scale hydraulic conductivity is lower, decreasing from K~10−2 m/day at 150–400 m BGL to 10−3 m/day down-dip at ~1 km BGL, where the pore fluid is hypersaline. Here, pore-scale permeability becomes progressively dominant with increasing lithostatic load. Notably, this work investigates a sandstone aquifer of fluvial origin at investigation depths consistent with highly enthalpy geothermal reservoirs (~0.7–1.1 km). At such depths, intergranular flow dominates in unfaulted areas with only minor contribution by bedding plane fractures. However, extensional faults represent preferential flow pathways, due to presence of high connective open fractures. Therefore, such faults may (1) drive nuclear waste contaminants towards the highly permeable shallow (< 150 m BGL) zone of the aquifer, and (2) influence fluid recovery in geothermal fields

Polaritonic Feshbach resonance

A Feshbach resonance occurs when the energy of two interacting free particles comes into resonance with a molecular bound state. When approaching this resonance, marked changes in the interaction strength between the particles can arise. Feshbach resonances provide a powerful tool for controlling the interactions in ultracold atomic gases, which can be switched from repulsive to attractive and have allowed a range of many-body quantum physics effects to be explored. Here we demonstrate a Feshbach resonance based on the polariton spinor interactions in a semiconductor microcavity. By tuning the energy of two polaritons with anti-parallel spins across the biexciton bound state energy, we show an enhancement of attractive interactions and a prompt change to repulsive interactions. A mean-field two-channel model quantitatively reproduces the experimental results. This observation paves the way for a new tool for tuning polariton interactions and to move forward into quantum correlated polariton physics

Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow.

The hundreds of cichlid fish species in Lake Malawi constitute the most extensive recent vertebrate adaptive radiation. Here we characterize its genomic diversity by sequencing 134 individuals covering 73 species across all major lineages. The average sequence divergence between species pairs is only 0.1-0.25%. These divergence values overlap diversity within species, with 82% of heterozygosity shared between species. Phylogenetic analyses suggest that diversification initially proceeded by serial branching from a generalist Astatotilapia-like ancestor. However, no single species tree adequately represents all species relationships, with evidence for substantial gene flow at multiple times. Common signatures of selection on visual and oxygen transport genes shared by distantly related deep-water species point to both adaptive introgression and independent selection. These findings enhance our understanding of genomic processes underlying rapid species diversification, and provide a platform for future genetic analysis of the Malawi radiation

An integrative approach for a network based meta-analysis of viral RNAi screens.

BACKGROUND: Big data is becoming ubiquitous in biology, and poses significant challenges in data analysis and interpretation. RNAi screening has become a workhorse of functional genomics, and has been applied, for example, to identify host factors involved in infection for a panel of different viruses. However, the analysis of data resulting from such screens is difficult, with often low overlap between hit lists, even when comparing screens targeting the same virus. This makes it a major challenge to select interesting candidates for further detailed, mechanistic experimental characterization. RESULTS: To address this problem we propose an integrative bioinformatics pipeline that allows for a network based meta-analysis of viral high-throughput RNAi screens. Initially, we collate a human protein interaction network from various public repositories, which is then subjected to unsupervised clustering to determine functional modules. Modules that are significantly enriched with host dependency factors (HDFs) and/or host restriction factors (HRFs) are then filtered based on network topology and semantic similarity measures. Modules passing all these criteria are finally interpreted for their biological significance using enrichment analysis, and interesting candidate genes can be selected from the modules. CONCLUSIONS: We apply our approach to seven screens targeting three different viruses, and compare results with other published meta-analyses of viral RNAi screens. We recover key hit genes, and identify additional candidates from the screens. While we demonstrate the application of the approach using viral RNAi data, the method is generally applicable to identify underlying mechanisms from hit lists derived from high-throughput experimental data, and to select a small number of most promising genes for further mechanistic studies. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13015-015-0035-7) contains supplementary material, which is available to authorized users

Monitoring the Nutrient Composition of Food Prepared Out-of-Home in the United Kingdom: Database Development and Case Study.

BACKGROUND: Hand transcribing nutrient composition data from websites requires extensive human resources and is prone to error. As a result, there are limited nutrient composition data on food prepared out of the home in the United Kingdom. Such data are crucial for understanding and monitoring the out-of-home food environment, which aids policy making. Automated data collection from publicly available sources offers a potential low-resource solution to address this gap. OBJECTIVE: In this paper, we describe the first UK longitudinal nutritional database of food prepared out of the home, MenuTracker. As large chains will be required to display calorie information on their UK menus from April 2022, we also aimed to identify which chains reported their nutritional information online in November 2021. In a case study to demonstrate the utility of MenuTracker, we estimated the proportions of menu items exceeding recommended energy and nutrient intake (eg, &gt;600 kcal per meal). METHODS: We have collated nutrient composition data of menu items sold by large chain restaurants quarterly since March 2021. Large chains were defined as those with 250 employees or more (those covered by the new calorie labeling policy) or belonging to the top 100 restaurants based on sales volume. We developed scripts in Python to automate the data collection process from business websites. Various techniques were used to harvest web data and extract data from nutritional tables in PDF format. RESULTS: Automated Python programs reduced approximately 85% of manual work, totaling 500 hours saved for each wave of data collection. As of January 2022, MenuTracker has 76,405 records from 88 large out-of-home food chains at 4 different time points (ie, March, June, September, and December) in 2021. In constructing the database, we found that one-quarter (24.5%, 256/1043) of large chains, which are likely to be subject to the United Kingdom's calorie menu labeling regulations, provided their nutritional information online in November 2021. Across these chains, 24.7% (16,391/66,295) of menu items exceeded the UK government's recommendation of a maximum of 600 kcal for a single meal. Comparable figures were 46.4% (29,411/63,416) for saturated fat, 34.7% (21,964/63,388) for total fat, 17.6% (11,260/64,051) for carbohydrates, 17.8% (11,434/64,059) for sugar, and 35.2% (22,588/64,086) for salt. Furthermore, 0.7% to 7.1% of the menu items exceeded the maximum daily recommended intake for these nutrients. CONCLUSIONS: MenuTracker is a valuable resource that harnesses the power of data science techniques to use publicly available data online. Researchers, policy makers, and consumers can use MenuTracker to understand and assess foods available from out-of-home food outlets. The methods used in development are available online and can be used to establish similar databases elsewhere.This work was funded by UKRI grant MC_UU_00006/7. For the purpose of open access, the author has applied a Creative Commons Attribution (CC BY) licence to any Author Accepted Manuscript version arising. YH is supported through a Gates Cambridge Scholarship. DT is supported by a PhD studentship awarded by the National Institute for Health Research (NIHR), School for Public Health Research (Grant No. PD‐SPH‐2015‐10025). No funders had any role in the study design; collection, analysis and interpretation of data; the writing of the manuscript; or the decision to submit the manuscript for publication