Inferring biological networks from genome-wide transcriptional and fitness data

In the last 15 years, the increased use of high throughput biology techniques such as genome-wide gene expression profiling, fitness profiling and protein interactomics has led to the generation of an extraordinary amount of data. The abundance of such diverse data has proven to be an essential foundation for understanding the complexities of molecular mechanisms and underlying pathways within a biological system. This thesis demonstrates the capabilities and applications of using biological networks to extrapolate biological information from the wealth of data available in the yeast species Saccharomyces cerevisiae and Schizosaccharomyces pombe. This study marks the first time a mutual information based network inference approach has been applied to a set of specific genome-wide expression and fitness compendia. In particular, this work has generated hypotheses in S. pombe that have led to a deeper understanding of the relationship between ribosomal proteins and energy metabolism, a recently discovered pathway termed riboneogenesis. Experimental validation of this hypothesis has led to new theories on the role of energy metabolism enzymes in controlling ribosome biogenesis in S. pombe, including the novel finding that fructose-1, 6-bisphosphatase (FBP1) may have roles in both gluconeogenesis and riboneogenesis. This thesis also demonstrates how the use of multi-level data allows for comprehensive insight into nuclear functions of the S. pombe nonsense-mediated mRNA decay protein, UPF1. This study provides substantial evidence demonstrating the role of UPF1 in DNA replication. The applicability of fitness data in identifying targets of metal and metalloid toxicity in S. cerevisiae has also been investigated

The N-terminus of Stag1 is required to repress the 2C program by maintaining rRNA expression and nucleolar integrity

Our understanding of how STAG proteins contribute to cell identity and disease have largely been studied from the perspective of chromosome topology and protein-coding gene expression. Here, we show that STAG1 is the dominant paralog in mouse embryonic stem cells (mESCs) and is required for pluripotency. mESCs express a wide diversity of naturally occurring Stag1 isoforms, resulting in complex regulation of both the levels of STAG paralogs and the proportion of their unique terminal ends. Skewing the balance of these isoforms impacts cell identity. We define a novel role for STAG1, in particular its N-terminus, in regulating repeat expression, nucleolar integrity, and repression of the two-cell (2C) state to maintain mESC identity. Our results move beyond protein-coding gene regulation via chromatin loops to new roles for STAG1 in nucleolar structure and function, and offer fresh perspectives on how STAG proteins, known to be cancer targets, contribute to cell identity and disease

Enhancer accessibility and CTCF occupancy underlie asymmetric TAD architecture and cell type specific genome topology.

Cohesin and CTCF are master regulators of genome topology. How these ubiquitous proteins contribute to cell-type specific genome structure is poorly understood. Here, we explore quantitative aspects of topologically associated domains (TAD) between pluripotent embryonic stem cells (ESC) and lineage-committed cells. ESCs exhibit permissive topological configurations which manifest themselves as increased inter- TAD interactions, weaker intra-TAD interactions, and a unique intra-TAD connectivity whereby one border makes pervasive interactions throughout the domain. Such 'stripe' domains are associated with both poised and active chromatin landscapes and transcription is not a key determinant of their structure. By tracking the developmental dynamics of stripe domains, we show that stripe formation is linked to the functional state of the cell through cohesin loading at lineage-specific enhancers and developmental control of CTCF binding site occupancy. We propose that the unique topological configuration of stripe domains represents a permissive landscape facilitating both productive and opportunistic gene regulation and is important for cellular identity

Enhancer accessibility and CTCF occupancy underlie asymmetric TAD architecture and cell type specific genome topology.

Cohesin and CTCF are master regulators of genome topology. How these ubiquitous proteins contribute to cell-type specific genome structure is poorly understood. Here, we explore quantitative aspects of topologically associated domains (TAD) between pluripotent embryonic stem cells (ESC) and lineage-committed cells. ESCs exhibit permissive topological configurations which manifest themselves as increased inter- TAD interactions, weaker intra-TAD interactions, and a unique intra-TAD connectivity whereby one border makes pervasive interactions throughout the domain. Such 'stripe' domains are associated with both poised and active chromatin landscapes and transcription is not a key determinant of their structure. By tracking the developmental dynamics of stripe domains, we show that stripe formation is linked to the functional state of the cell through cohesin loading at lineage-specific enhancers and developmental control of CTCF binding site occupancy. We propose that the unique topological configuration of stripe domains represents a permissive landscape facilitating both productive and opportunistic gene regulation and is important for cellular identity

Cohesin-independent STAG proteins interact with RNA and localise to R-loops to promote complex loading

Most studies of cohesin function consider the Stromalin Antigen (STAG/SA) proteins as core complex members given their ubiquitous interaction with the cohesin ring. Here, we provide functional data to support the notion that the SA subunit is not a mere passenger in this structure, but instead plays a key role in the localization of cohesin to diverse biological processes and promotes loading of the complex at these sites. We show that in cells acutely depleted for RAD21, SA proteins remain bound to chromatin, cluster in 3D and interact with CTCF, as well as with a wide range of RNA binding proteins involved in multiple RNA processing mechanisms. Accordingly, SA proteins interact with RNA and are localised to R-loops where they contribute to R-loop regulation. Our results place SA1 within R-loop domains upstream of the cohesin complex and reveal a role for SA1 in cohesin loading which is independent of NIPBL, the canonical cohesin loader. We propose that SA1 takes advantage of structural R-loop platforms to link cohesin loading and chromatin structure with diverse functions. Since SA proteins are pan-cancer targets, and R-loops play an increasingly prevalent role in cancer biology, our results have important implications for the mechanistic understanding of SA proteins in cancer and disease

Cohesin-independent STAG proteins interact with RNA and R-loops and promote complex loading

Most studies of cohesin function consider the Stromalin Antigen (STAG/SA) proteins as core complex members given their ubiquitous interaction with the cohesin ring. Here, we provide functional data to support the notion that the SA subunit is not a mere passenger in this structure, but instead plays a key role in the localization of cohesin to diverse biological processes and promotes loading of the complex at these sites. We show that in cells acutely depleted for RAD21, SA proteins remain bound to chromatin, cluster in 3D and interact with CTCF, as well as with a wide range of RNA binding proteins involved in multiple RNA processing mechanisms. Accordingly, SA proteins interact with RNA, RNA binding proteins and R-loops, even in the absence of cohesin. Our results place SA1 on chromatin upstream of the cohesin ring and reveal a role for SA1 in cohesin loading which is independent of NIPBL, the canonical cohesin loader. We propose that SA1 takes advantage of structural R-loop platforms to link cohesin loading and chromatin structure with diverse functions. Since SA proteins are pan-cancer targets, and R-loops play an increasingly prevalent role in cancer biology, our results have important implications for the mechanistic understanding of SA proteins in cancer and disease

Genome-wide chromosomal association of Upf1 is linked to Pol II transcription in Schizosaccharomyces pombe

Although the RNA helicase Upf1 has hitherto been examined mostly in relation to its cytoplasmic role in nonsense mediated mRNA decay (NMD), here we report high-throughput ChIP data indicating genome-wide association of Upf1 with active genes in Schizosaccharomyces pombe. This association is RNase sensitive, correlates with Pol II transcription and mRNA expression levels. Changes in Pol II occupancy were detected in a Upf1 deficient (upf1∆) strain, prevalently at genes showing a high Upf1 relative to Pol II association in wild-type. Additionally, an increased Ser2 Pol II signal was detected at all highly transcribed genes examined by ChIP-qPCR. Furthermore, upf1cells are hypersensitive to the transcription elongation inhibitor 6-azauracil. A significant proportion of the genes associated with Upf1 in wild-type conditions are also mis-regulated in upf1. These data envisage that by operating on the nascent transcript, Upf1 might influence Pol II phosphorylation and transcription

Ribosomal proteins' association with transcription sites peaks at tRNA genes in Schizosaccharomyces pombe

Ribosomal proteins (RPs) are essential components of ribosomes, but several RPs are also present at transcription sites of eukaryotic chromosomes. Here, we report a genome-wide ChIP-on-chip analysis of the association of three representative 60S RPs with sites in the Schizosaccharomyces pombe chromosomes. All three proteins tend to bind at the same subset of coding and noncoding loci. The data demonstrate selective RNA-dependent interactions between RPs and many transcription sites and suggest that the RPs bind as components of a preassembled multiprotein complex, perhaps 60S or pre-60S subunits. These findings further indicate that the presence of RPs complexes at transcription sites might be a general feature of eukaryotic cells and functionally important. Unexpectedly, the RPs' chromosomal association is highest at centromeres and tRNA genes—the RPs were found at 167 of the 171 tRNA genes assayed. These findings raise the intriguing possibility that RP complexes are involved in tRNA biogenesis and possibly centromere functions

A novel intergenic enhancer that regulates Bdnf expression in developing cortical neurons

Summary: Brain-derived neurotrophic factor (BDNF) promotes neuronal differentiation and survival and is implicated in the pathogenesis of many neurological disorders. Here, we identified a novel intergenic enhancer located 170 kb from the Bdnf gene, which promotes the expression of Bdnf transcript variants during mouse neuronal differentiation and activity. Following Bdnf activation, enhancer-promoter contacts increase, and the region moves away from the repressive nuclear periphery. Bdnf enhancer activity is necessary for neuronal clustering and dendritogenesis in vitro, and for cortical development in vivo. Our findings provide the first evidence of a regulatory mechanism whereby the activation of a distal enhancer promotes Bdnf expression during brain development