Asymmetric Conservation within Pairs of Co-Occurred Motifs Mediates Weak Direct Binding of Transcription Factors in ChIP-Seq Data.

(1) Background: Transcription factors (TFs) are main regulators of eukaryotic gene expression. The cooperative binding to genomic DNA of at least two TFs is the widespread mechanism of transcription regulation. Cooperating TFs can be revealed through the analysis of co-occurrence of their motifs. (2) Methods: We applied the motifs co-occurrence tool (MCOT) that predicted pairs of spaced or overlapped motifs (composite elements, CEs) for a single ChIP-seq dataset. We improved MCOT capability for the prediction of asymmetric CEs with one of the participating motifs possessing higher conservation than another does. (3) Results: Analysis of 119 ChIP-seq datasets for 45 human TFs revealed that almost for all families of TFs the co-occurrence with an overlap between motifs of target TFs and more conserved partner motifs was significantly higher than that for less conserved partner motifs. The asymmetry toward partner TFs was the most clear for partner motifs of TFs from the ETS (E26 Transformation Specific) family. (4) Conclusion: Co-occurrence with an overlap of less conserved motif of a target TF and more conserved motifs of partner TFs explained a substantial portion of ChIP-seq data lacking conserved motifs of target TFs. Among other TF families, conservative motifs of TFs from ETS family were the most prone to mediate interaction of target TFs with its weak motifs in ChIP-seq.

Exome sequencing identifies a de novo FOXA2 variant in a patient with syndromic diabetes.

OBJECTIVE: When diabetes is associated with congenital malformations, without autoimmune antibodies, a genetic cause is suspected. Here, we aimed to identify a defective gene that led to diabetes. RESEARCH DESIGN AND METHODS: We performed an exome analysis of an index case and his healthy parents. RESULTS: The child presented with childhood-onset diabetes, congenital hypopituitarism, cardiac malformation, and anal atresia. A DNA analysis revealed a heterozygous de novo pathogenic variant in the developmental transcription factor, forkhead box A2 (FOXA2). The mutation resided in the DNA-binding domain, which is highly conserved among species. Tridimensional molecular dynamics simulation modeling predicted an altered interaction between the mutated protein and DNA. CONCLUSIONS: A defect in the FOXA2 DNA-binding domain was associated with childhood-onset diabetes and multiple congenital anomalies, which reflected the pleiotropic nature of the gene. This report extends the recently described phenotype of neonatal hypoglycemia to later-onset diabetes. We suggest to include FOXA2 analysis for neonatal hypoglycemia and to implement a long-term follow-up, particularly for the risk of diabetes.

Structure of the Forkhead Domain of FOXA2 Bound to a Complete DNA Consensus Site.

FOXA2, a member of the forkhead family of transcription factors, plays essential roles in liver development and bile acid homeostasis. In this study, we report a 2.8 Å co-crystal structure of the FOXA2 DNA-binding domain (FOXA2-DBD) bound to a DNA duplex containing a forkhead consensus binding site (GTAAACA). The FOXA2-DBD adopts the canonical winged-helix fold, with helix H3 and wing 1 regions mainly mediating the DNA recognition. Although the wing 2 region was not defined in the structure, isothermal titration calorimetry assays suggested that this region was required for optimal DNA binding. Structure comparison with the FOXA3-DBD bound to DNA revealed more major groove contacts and fewer minor groove contacts in the FOXA2 structure than in the FOXA3 structure. Structure comparison with the FOXO1-DBD bound to DNA showed that different forkhead proteins could induce different DNA conformations upon binding to identical DNA sequences. Our findings provide the structural basis for FOXA2 protein binding to a consensus forkhead site and elucidate how members of the forkhead protein family bind different DNA sites.

Changes in Olfactory Receptor Expression Are Correlated With Odor Exposure During Early Development in the zebrafish (Danio rerio).

We have previously shown that exposure to phenyl ethyl alcohol (PEA) causes an increase in the expression of the transcription factor otx2 in the olfactory epithelium (OE) of juvenile zebrafish, and this change is correlated with the formation of an odor memory of PEA. Here, we show that the changes in otx2 expression are specific to ßPEA: exposure to αPEA did not affect otx2 expression. We identified 34 olfactory receptors (ORs) representing 16 families on 4 different chromosomes as candidates for direct regulation of OR expression via Otx2. Subsequent in silico analysis uncovered Hnf3b binding sites closely associated with Otx2 binding sites in the regions flanking the ORs. Analysis by quantitative polymerase chain reaction and RNA-seq of OR expression in developing zebrafish exposed to different isoforms of PEA showed that a subset of ORs containing both Otx2/Hnf3b binding sites were downregulated only in ßPEA-exposed juveniles and this change persisted through adult life. Localization of OR expression by in situ hybridization indicates the downregulation occurs at the level of RNA and not the number of cells expressing a given receptor. Finally, analysis of immediate early gene expression in the OE did not reveal changes in c-fos expression in response to either αPEA or ßPEA.

Species-specific strategies underlying conserved functions of metabolic transcription factors.

The winged helix protein FOXA2 and the nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ) are highly conserved, regionally expressed transcription factors (TFs) that regulate networks of genes controlling complex metabolic functions. Cistrome analysis for Foxa2 in mouse liver and PPARγ in mouse adipocytes has previously produced consensus-binding sites that are nearly identical to those used by the corresponding TFs in human cells. We report here that, despite the conservation of the canonical binding motif, the great majority of binding regions for FOXA2 in human liver and for PPARγ in human adipocytes are not in the orthologous locations corresponding to the mouse genome, and vice versa. Of note, TF binding can be absent in one species despite sequence conservation, including motifs that do support binding in the other species, demonstrating a major limitation of in silico binding site prediction. Whereas only approximately 10% of binding sites are conserved, gene-centric analysis reveals that about 50% of genes with nearby TF occupancy are shared across species for both hepatic FOXA2 and adipocyte PPARγ. Remarkably, for both TFs, many of the shared genes function in tissue-specific metabolic pathways, whereas species-unique genes fail to show enrichment for these pathways. Nonetheless, the species-unique genes, like the shared genes, showed the expected transcriptional regulation by the TFs in loss-of-function experiments. Thus, species-specific strategies underlie the biological functions of metabolic TFs that are highly conserved across mammalian species. Analysis of factor binding in multiple species may be necessary to distinguish apparent species-unique noise and reveal functionally relevant information.

Regulated Fox-2 isoform expression mediates protein 4.1R splicing during erythroid differentiation.

A regulated splicing event in protein 4.1R pre-mRNA-the inclusion of exon 16-encoding peptides for spectrin-actin binding-occurs in late erythroid differentiation. We defined the functional significance of an intronic splicing enhancer, UGCAUG, and its cognate splicing factor, mFox2A, on exon 16 splicing during differentiation. UGCAUG displays cell-type-specific splicing regulation in a test neutral reporter and has a dose-dependent enhancing effect. Erythroid cells express 2 UGCAUG-binding mFox-2 isoforms, an erythroid differentiation-inducible mFox-2A and a commonly expressed mFox-2F. When overexpressed, both enhanced internal exon splicing in an UGCAUG-dependent manner, with mFox-2A exerting a much stronger effect than mFox-2F. A significant reciprocal increase in mFox-2A and decrease in mFox-2F occurred during erythroid differentiation and correlated with exon 16 inclusion. Furthermore, isoform-specific expression reduction reversed mFox-2A-enhancing activity, but not that of mFox-2F on exon 16 inclusion. Our results suggest that an erythroid differentiation-inducible mFox-2A isoform is a critical regulator of the differentiation-specific exon 16 splicing switch, and that its up-regulation in late erythroid differentiation is vital for exon 16 splicing.

The Foxa family of transcription factors in development and metabolism.

The Foxa subfamily of winged helix/forkhead box (Fox) transcription factors has been the subject of genetic and biochemical study for over 15 years. During this time its three members, Foxa1, Foxa2 and Foxa3, have been found to play important roles in multiple stages of mammalian life, beginning with early development, continuing during organogenesis, and finally in metabolism and homeostasis in the adult. Foxa2 is required for the formation of the node and notochord, and in its absence severe defects in gastrulation, neural tube patterning, and gut morphogenesis result in embryonic lethality. Foxa1 and Foxa2 cooperate to establish competence in foregut endoderm and are required for normal development of endoderm-derived organs such as the liver, pancreas, lungs, and prostate. In post-natal life, members of the Foxa family control glucose metabolism through the regulation of multiple target genes in the liver, pancreas, and adipose tissue. Insight into the unique molecular basis of Foxa function has been obtained from recent genetic and genomic data, which identify the Foxa proteins as 'pioneer factors' whose binding to promoters and enhancers enable chromatin access for other tissue-specific transcription factors.