spatzie: An R package for identifying significant transcription factor motif co-enrichment from enhancer-promoter interactions

Genomic interactions provide important context to our understanding of the state of the genome. One question is whether specific transcription factor interactions give rise to genome organization. We introduce spatzie, an R package and a website that implements statistical tests for significant transcription factor motif cooperativity between enhancer-promoter interactions. We conducted controlled experiments under realistic simulated data from ChIP-seq to confirm spatzie is capable of discovering co-enriched motif interactions even in noisy conditions. We then use spatzie to investigate cell type specific transcription factor cooperativity within recent human ChIA-PET enhancer-promoter interaction data. The method is available online at https://spatzie.mit.edu.

Landscape of allele-specific transcription factor binding in the human genome

Sequence variants in gene regulatory regions alter gene expression and contribute to phenotypes of individual cells and the whole organism, including disease susceptibility and progression. Single-nucleotide variants in enhancers or promoters may affect gene transcription by altering transcription factor binding sites. Differential transcription factor binding in heterozygous genomic loci provides a natural source of information on such regulatory variants. We present a novel approach to call the allele-specific transcription factor binding events at single-nucleotide variants in ChIP-Seq data, taking into account the joint contribution of aneuploidy and local copy number variation, that is estimated directly from variant calls. We have conducted a meta-analysis of more than 7 thousand ChIP-Seq experiments and assembled the database of allele-specific binding events listing more than half a million entries at nearly 270 thousand single-nucleotide polymorphisms for several hundred human transcription factors and cell types. These polymorphisms are enriched for associations with phenotypes of medical relevance and often overlap eQTLs, making candidates for causality by linking variants with molecular mechanisms. Specifically, there is a special class of switching sites, where different transcription factors preferably bind alternative alleles, thus revealing allele-specific rewiring of molecular circuitry.

The testis-specific transcription factor TCFL5 responds to A MYB to elaborate the male meiotic program in placental mammals

In male mice, the transcription factor (TF) A MYB initiates reprogramming of gene expression after spermatogonia enter meiosis. We report that A MYB activates Tcfl5, a testis-specific TF first produced in pachytene spermatocytes. Subsequently, A MYB and TCFL5 reciprocally reinforce their own transcription to establish an extensive circuit that regulates meiosis. TCFL5 promotes transcription of genes required for mRNA turnover, pachytene piRNA production, meiotic exit, and spermiogenesis. This transcriptional architecture is conserved in rhesus macaque, suggesting TCFL5 plays a central role in meiosis and spermiogenesis in placental mammals. Tcfl5em1/em1 mutants are sterile, and spermatogenesis arrests at the mid- or late-pachytene stage of meiosis.

DNA affinity purification sequencing and transcriptional profiling reveal new aspects of nitrogen regulation in a filamentous fungus

Sensing available nutrients and efficiently utilizing them is a challenge common to all organisms. The model filamentous fungus Neurospora crassa is capable of utilizing a variety of inorganic and organic nitrogen sources. Nitrogen utilization in N. crassa is regulated by a network of pathway-specific transcription factors that activate genes necessary to utilize specific nitrogen sources in combination with nitrogen catabolite repression regulatory proteins. We identified an uncharacterized pathway-specific transcription factor, amn-1, that is required for utilization of the nonpreferred nitrogen sources proline, branched-chain amino acids, and aromatic amino acids. AMN-1 also plays a role in regulating genes involved in responding to the simple sugar mannose, suggesting an integration of nitrogen and carbon metabolism. The utilization of nonpreferred nitrogen sources, which require metabolic processing before being used as a nitrogen source, is also regulated by the nitrogen catabolite regulator NIT-2. Using RNA sequencing combined with DNA affinity purification sequencing, we performed a survey of the role of NIT-2 and the pathway-specific transcription factors NIT-4 and AMN-1 in directly regulating genes involved in nitrogen utilization. Although previous studies suggested promoter binding by both a pathway-specific transcription factor and NIT-2 may be necessary for activation of nitrogen-responsive genes, our data show that pathway-specific transcription factors regulate genes involved in the catabolism of specific nitrogen sources, while NIT-2 regulates genes involved in utilization of all nonpreferred nitrogen sources, such as nitrogen transporters. Together, these transcription factors form a nutrient sensing network that allows N. crassa cells to regulate nitrogen utilization.

Balance between extracellular matrix production and macrophage survival by a Salmonella-specific SPI-2 encoded transcription factor

Cellulose is a major component of the Salmonella biofilm extracellular matrix and it is considered an antivirulence factor because it interferes with Salmonella survival inside macrophages and virulence in mice. Its synthesis is stimulated by CsgD, the master regulator of biofilm extracellular matrix formation in enterobacteria, which in turn is under the control of MlrA, a MerR-like transcription factor. In this work we identified a SPI-2 encoded Salmonella-specific transcription factor homolog to MlrA, MlrB, that represses transcription of its downstream gene, STM1389, also known as orf319, and of csgD inside host cells. MlrB is induced in laboratory media mimicking intracellular conditions and inside macrophages, and it is required for intramacrophage survival. An increased expression of csgD is observed in the absence of MlrB inside host cells. Interestingly, inactivation of the CsgD-controlled cellulose synthase coding-gene, bcsA, restored intramacrophage survival to rates comparable to wild type bacteria in the absence of MlrB. These data indicate that MlrB represses CsgD expression inside host cells and in consequence activation of the cellulose synthase. Our findings provide a novel link between biofilm formation and Salmonella virulence.

Reduced miR-184-3p expression occurring in Type 2 diabetic pancreatic islets protects β-cells from lipotoxic and proinflammatory apoptosis via a CRTC1-dependent mechanism

Loss of functional β-cell mass in Type 2 diabetes (T2D) involves molecular mechanisms including β-cell apoptosis, dysfunction, and/or dedifferentiation. MicroRNA miR-184-3p has been demonstrated to be involved in multiple β-cell functions including insulin secretion, proliferation and survival. However, downstream targets and upstream regulators of miR-184-3p have not yet been fully elucidated. Here, we showed that levels of miR-184-3p are reduced in human T2D pancreatic islets and that its reduction protected β-cells from lipotoxic- and inflammatory-induced apoptosis. Interestingly, CREB-Transcriptional Coactivator-1 (CRTC1) is a direct target of miR-184-3p and indeed its expression is upregulated in human T2D pancreatic islets. The downregulation of miR-184-3p in β-cells induced the upregulation of CRTC1 both at mRNA and protein level. Of note, miR-184-3p protection effect was dependent on CRTC1, since its silencing in human β-cells abrogates the protective mechanism exerted by miR-184-3p inhibition. Additionally, we found that the β-cell specific transcription factor NKX6.1, whose DNA binding sites were predicted to be present in human and mouse MIR184 gene promoter sequence, was reduced in T2D human pancreatic islets, in line with miR-184-3p downregulation, and was positively correlated with microRNA expression. Using chromatin immunoprecipitation analysis and mRNA silencing experiments, we demonstrated that NKX6.1 directly controls both human and murine miR-184 expression.In conclusion, we found that miR-184-3p expression is controlled by the β-cell specific transcription factor NKX6.1 and that miR-184-3p reduction protects β-cells from apoptosis through the upregulation of its target gene CRTC1.

Landscape of allele-specific transcription factor binding in the human genome

Sequence variants in gene regulatory regions alter gene expression and contribute to phenotypes of individual cells and the whole organism, including disease susceptibility and progression. Single-nucleotide variants in enhancers or promoters may affect gene transcription by altering transcription factor binding sites. Differential transcription factor binding in heterozygous genomic loci provides a natural source of information on such regulatory variants. We present a novel approach to call the allele-specific transcription factor binding events at single-nucleotide variants in ChIP-Seq data, taking into account the joint contribution of aneuploidy and local copy number variation, that is estimated directly from variant calls. We have conducted a meta-analysis of more than 7 thousand ChIP-Seq experiments and assembled the database of allele-specific binding events listing more than half a million entries at nearly 270 thousand single-nucleotide polymorphisms for several hundred human transcription factors and cell types. These polymorphisms are enriched for associations with phenotypes of medical relevance and often overlap eQTLs, making candidates for causality by linking variants with molecular mechanisms. Specifically, there is a special class of switching sites, where different transcription factors preferably bind alternative alleles, thus revealing allele-specific rewiring of molecular circuitry.