Uncovering the interaction between genomes and the environment is a principal challenge of modern genomics and preventive medicine. While theoretical models are well defined, little is known of the G × E interactions in humans. We used an integrative approach to comprehensively assess the interactions between 1.6 million data points, encompassing a range of environmental exposures, health, and gene expression levels, coupled with whole-genome genetic variation. From ∼1000 individuals of a founder population in Quebec, we reveal a substantial impact of the environment on the transcriptome and clinical endophenotypes, overpowering that of genetic ancestry. Air pollution impacts gene expression and pathways affecting cardio-metabolic and respiratory traits, when controlling for genetic ancestry. Finally, we capture four expression quantitative trait loci that interact with the environment (air pollution). Our findings demonstrate how the local environment directly affects disease risk phenotypes and that genetic variation, including less common variants, can modulate individual’s response to environmental challenges.

Controlling the transcriptional program is essential to maintain the identity and the biological functions of a cell. The Mediator and Cohesin complexes have been established as central cofactors controlling the transcriptional program in normal cells. However, the distribution, recruitment and importance of these complexes in cancer cells have not been fully investigated. Here we show that FOXA and master transcription factors are part of the core transcriptional regulatory circuitry of cancer cells and are essential to recruit M ediator and Cohesin. Indeed, Mediator and Cohesin occupied the enhancer and promoter regions of actively transcribed genes and maintained the proliferation and colony forming potential. Through integration of publically available ChIP-Seq datasets, we predicted the core transcriptional regulatory circuitry of each cancer cell. Unexpectedly, for all cells investigated, the pioneer transcription factors FOXA1 and/or FOXA2 were identified in addition to cell-specific master transcription factors. Loss of both types of transcription factors phenocopied the loss of Mediator and Cohesin. Lastly, the master and pioneer transcription factors were essential to recruit Mediator and Cohesin to regulatory regions of actively transcribed genes. Our study proposes that maintenance of the cancer cell state is dependent on recruitment of Mediator and Cohesin through FOXA and master transcription factors.

In addition to its role in sister chromatid cohesion, genome stability and integrity, the cohesin complex is involved in gene transcription. Mutations in core cohesin subunits SMC1A, SMC3 and RAD21, or their regulators NIPBL and HDAC8, cause Cornelia de Lange syndrome (CdLS). Recent evidence reveals that gene expression dysregulation could be the underlying mechanism for CdLS. These findings raise intriguing questions regarding the potential role of cohesin-mediated transcriptional control and pathogenesis. Here, we identified numerous dysregulated genes occupied by cohesin by combining the transcriptome of CdLS cell lines carrying mutations in SMC1A gene and ChIP-Seq data. Genome-wide analyses show that genes changing in expression are enriched for cohesin-binding. In addition, our results indicate that mutant cohesin impairs both RNA polymerase II (Pol II) transcription initiation at promoters and elongation in the gene body. These findings highlight the pivotal role of cohesin in transcriptional regulation and provide an explanation for the typical gene dysregulation observed in CdLS patients.

The binding of PRDM9 to chromatin is a key step in the induction of DNA double-strand breaks associated with meiotic recombination hotspots; it is normally expressed solely in germ cells. We interrogated 1879 cancer samples in 39 different cancer types and found that PRDM9 is unexpectedly expressed in 20% of these tumors even after stringent gene homology correction. The expression levels of PRDM9 in tumors are significantly higher than those found in healthy neighboring tissues and in healthy nongerm tissue databases. Recurrently mutated regions located within 5 Mb of the PRDM9 loci, as well as differentially expressed genes in meiotic pathways, correlate with PRDM9 expression. In samples with aberrant PRDM9 expression, structural variant breakpoints frequently neighbor the DNA motif recognized by PRDM9, and there is an enrichment of structural variants at sites of known meiotic PRDM9 activity. This study is the first to provide evidence of an association between aberrant expression of the meiosis-specific gene PRDM9 with genomic instability in cancer.