Dissecting Gene Regulatory Networks Using Targeted Quantitative Proteomics

Gene regulatory networks control gene expression levels, and therefore play an essential role in mammalian development and function. Regulation of gene expression is the result of a complex interplay between DNA regulatory elements and their binding partners, known as transcription factors (TFs). Due to their vital role in development, intercellular signalling, cell cycle and disease development, elucidating the mechanisms by which TFs regulate gene expression is of crucial importance in the vast majority of biological processes. In particular, understanding how each TF contributes to the expression output of its respective target gene in space and time will help to elucidate how gene regulatory networks (GRNs) behave under different physiological or pathological conditions. Although extensive work has been accomplished in characterizing the key TFs involved in many biological processes, almost no quantitative information is currently available in the literature. To get a deep insight into the complex mechanisms underlying the regulation of gene expression, we need to acquire quantitative information, since TF abundance within the cell can be linked to their transcriptional capabilities. Such information would be of utmost importance to build accurate in silico quantitative DNA binding models that could predict and explain the particular properties of gene regulatory mechanisms. The quantification of TFs is a difficult task due their natural low abundance in cells, and their reliable detection is therefore very much dependent on the overall sensitivity of current technologies. In recent years, a new MS-based technology termed selected reaction monitoring (SRM) has gained popularity due to the targeted nature of its approach that allows the detection and quantification of proteins in complex samples with an exceptional sensitivity and specificity. I will show in this thesis, this approach is particularly well suited for targeting low abundant proteins such as TFs, which are otherwise difficult to identify with conventional shotgun proteomics experiments. Consequently, the main focus of my thesis research project entailed the development of an SRM-based platform aimed at quantifying TFs in absolute amounts based on in vitro protein expression during the terminal stage of adipogenesis, using the pre-adipocyte 3T3-L1 cell line. Interestingly, our initial efforts led to the creation of an atlas of TF-specific peptide data, which could be readily used for the design of quantitative assays. In the first phase, abundance measurements in terms of copies per cell were derived at precise differentiation time-points for two major adipogenic players, PPARγ and RXRα. In the second phase, we expanded the number of adipogenic TFs that can be monitored in one assay, allowing for the quantification of up to 10 TFs in one single, integrated SRM run. Such upscale increases the practical usefulness of the methodology while reducing the associated costs, and ultimately allows for non-negligible time-savings. The availability of absolute protein copy number data permitted us ultimately to examine the relationship between the number of genome-wide DNA binding events and TF molecules. We derived a quantitative DNA binding model that allowed the prediction of the number of PPARγ ChIP-seq binding events given its nuclear abundance, chromatin state, and DNA binding energetics. As such, we were able to explain the paradoxical observation of a significant increase in PPARγ binding sites despite a saturation in the number of PPARγ molecules. We thus demonstrate how TF abundance data can be modeled in conjunction with large-scale DNA occupancy and chromatin state data to further our understanding of gene regulatory mechanisms mediating cellular differentiation. We are now starting to build on our pioneering work to quantify in absolute terms key players of the entire, core adipogenic GRN, as such aiming to provide a quantitative explanation of the regulatory mechanisms at play during the terminal phase of adipocyte differentiation. Moreover, to increase the explicative power of our methodology and to alleviate the throughput limitation that comes with obtaining absolute protein measurements, we decided to perform copy-number estimates for a larger set of adipogenic TFs utilizing a modified version of our original approach. At the cost of a modest loss of accuracy, we are now aiming to develop a sensitive and robust methodology that will allow the quantification of entire GRNs at low cost and in a time-effective manner. This is consistent with the overall goal in life sciences or clinical research to improve our ability to accurately and reproducibly quantify entire pathways or biological networks to improve our systems understanding of biological processes.

The role of OCT4 and SOX2 in regulating chromatin accessibility and cell fate across the cell cycle in embryonic stem cells

Torgil Elias Friman

Precise spatiotemporal regulation of gene expression is essential for development and homeostasis of complex organisms. This is achieved in large part by sequence-specific transcription factors (TF) that bind to genomic regulatory elements to activate or repress transcription. DNA in eukaryotic nuclei is wrapped around histones into nucleosomes, which are inaccessible to most proteins. However, certain TFs can access their binding sites in the presence of histones and promote nucleosome eviction to increase chromatin accessibility, allowing for binding of other protein. These TFs, called pioneer factors, include OCT4 and SOX2, which are master regulators of embryonic stem (ES) cells. Dividing cells, including stem cells, must achieve proper gene regulation despite the progression of the cell cycle, which imposes several constraints. During S phase, the replication fork passes through the genome to make a copy of each chromosome, which leads to a loss of chromatin accessibility. Furthermore, the substantial chromatin compaction that occurs during mitosis to prepare for chromosome segregation is associated with eviction of most DNA-binding proteins and reduced accessibility at distal cis-regulatory elements including enhancers. How pioneer factors operate in the context of cell cycle progression is unclear. In this thesis, I will present work studying different features of OCT4 and SOX2, including how they interplay to regulate chromatin accessibility in ES cells, the impact of small endogenous fluctuations of their protein levels on cell fate and chromatin accessibility, and how their roles in regulating cell fate and chromatin accessibility is related to different stages of the cell cycle. Using live-cell imaging of fluorescently labeled proteins, we discovered that OCT4 and SOX2 were associated to mitotic chromatin. Using cell-cycle specific degradation strategies, we could show that the presence of SOX2 and OCT4 at the mitosis-G1 (M-G1) transition contributes to maintain the pluripotency of ES cells. We further showed that ES cells with high levels of OCT4 in G1, but not in S phase, are more prone to differentiate towards neuroectoderm and mesendoderm cell fates. Cells with high levels of OCT4 displayed increased chromatin accessibility at enhancers of differentiation genes, suggesting accessibility fluctuates with OCT4 levels and can prime cells for differentiation. To explore potential cell cycle-dependent regulation of chromatin accessibility, we degraded OCT4 at the M-G1 transition, which led to a reduction in accessibility at many OCT4-bound regions. When OCT4 returned later in the cell cycle, chromatin accessibility was recovered, although at many loci this recovery was incomplete, suggesting that OCT4 is required at the M-G1 transition to maintain these regions open. We then rapidly degraded OCT4 at other phases of the cell cycle. Surprisingly, loss of chromatin accessibility was highly similar at all cell cycle phases, showing that OCT4 is constantly required to maintain regulatory elements accessible. To explore the dynamic regulation of accessibility by OCT4, we performed time-course measurements of accessibility upon its degradation. This revealed that loss of accessibility was temporally correlated to the degradation of OCT4. These results suggest that chromatin accessibility is highly dynamic and continuously dependent on the presence of OCT4 across the cell cycle in ES cells.

EPFL2019

Absolute quantification of transcription factors during cellular differentiation using multiplexed targeted proteomics

Bart Deplancke, Paola Alejandra Gilardoni, Irina Krier, Marc Moniatte, Felix Naef, Sunil Kumar Raghav, Adrian Schmid, Jovan Simicevic, Benjamin Zoller

The cellular abundance of transcription factors (TFs) is an important determinant of their regulatory activities. Deriving TF copy numbers is therefore crucial to understanding how these proteins control gene expression. We describe a sensitive selected reaction monitoring-based mass spectrometry assay that allowed us to determine the copy numbers of up to ten proteins simultaneously. We applied this approach to profile the absolute levels of key TFs, including PPARγ and RXRα, during terminal differentiation of mouse 3T3-L1 pre-adipocytes. Our analyses revealed that individual TF abundance differs dramatically (from ∼250 to >300,000 copies per nucleus) and that their dynamic range during differentiation can vary up to fivefold. We also formulated a DNA binding model for PPARγ based on TF copy number, binding energetics and local chromatin state. This model explains the increase in PPARγ binding sites during the final differentiation stage that occurs despite a concurrent saturation in PPARγ copy number.

Nature Publishing Group2013

Development and application of experimental tools to elucidate the gene regulatory network underlying adipogenesis

Carine Delattre-Gubelmann

Excess fat mass in obese individuals, which is characterized by an increase of white adipocyte cell size and number, significantly raises the risk of developing metabolic syndrome symptoms as well as cancer. A detailed understanding of the transcriptional mechanisms mediating white adipocyte differentiation (i.e., adipogenesis) is required to counteract the growing obesity epidemics. In this thesis, we have implemented two combinatory approaches to detect new protein-DNA interactions in the mouse, one of which targets the adipogenic gene regulatory network (GRN) that underlies the white adipocyte differentiation process. As a first approach to shed light on GRNs, we have developed and validated a powerful approach that combines a high-throughput yeast one-hybrid-based (Y1H) system to detect transcription factors (TFs)-DNA element interactions, together with a novel microfluidics-based technique, enabling the identification and characterization of individual binding sites for detected TFs within the respective regulatory sequences at unprecedented throughput and resolution. Using well-described regulatory elements as well as an uncharacterized enhancer which can play also a role in adipogenesis, we show that our approach in a first phase uncovers many known as well as novel TF-DNA interactions, of which a large proportion were independently validated. In addition, we further demonstrate how our approach can be used in a second phase to characterize detected interactions by enabling the generation of a relative DNA occupancy landscape over the length of the respective DNA element. Furthermore, we provide evidence that this system enables the detection of novel interactions that may have in vivo relevance. Given the strict requirements for direct and well-characterized protein-DNA interaction data to generate and model gene regulatory networks, the presented two-tiered approach should be of great value to enhance our understanding of the regulatory principles underlying mouse gene expression. To identify transcription factors involved in white adipocyte differentiation, we overexpressed TFs in preadipocytes and quantified their effect on differentiation. To this end, we have created a mouse open-reading frame resource consisting of more than 750 fully sequence-verified TF clones. We transferred these TF clones to a Tet-On lentiviral vector and used them to infect white preadipocytes. Based on microscopy-assisted quantification of lipids after induction of differentiation, we identified putatively novel TF candidates that strongly enhanced differentiation based on the increased lipid accumulation, indicating a potential role in adipogenesis. Gene expression and ChIP-seq experiments have been performed for the three top candidates (ZEB1, ZFP30 and ZFP277) to further study their function within the adipogenic regulatory network. To assess whether the effect of these TFs on adipogenesis could be reproduced in vivo, murine adipose stromal-vascular fraction containing adipocyte precursor cells were implanted into mice fed a high-fat diet, after each of the three TF candidates were knockdown or overexpressed in these cells. This experiment and the RNA-seq analyses are currently being processed while this thesis is being written. Furthermore, due to the absence of quantitative real-time PCR (qPCR) primer designing tools that would allow for gene-specific expression profiling in a high-throughput fashion, we design a user-friendly platform: GETPrime. This platform uniquely combines and automates several features critical to optimal qPCR primer design for all annotated Homo sapiens, Mus musculus, Caenorhabditis elegans, Drosophila melanogaster and Danio rerio genes in the Ensembl database. GETPrime primers have been extensively validated experimentally, demonstrating their high quality as well as high transcript specificity in complex samples.