Bispecific Forkhead Transcription Factor FoxN3 Recognizes Two Distinct Motifs with Different DNA Shapes

A role for pH dynamics regulating transcription factor DNA-binding selectivity

Intracellular pH (pHi) dynamics regulates diverse cell processes such as proliferation, dysplasia, and differentiation, often mediated by the protonation state of a functionally critical histidine residue in endogenous pH sensing proteins. How pHi dynamics can directly regulate gene expression or whether transcription factors can function as pH sensors has received limited attention. We tested the prediction that transcription factors with a histidine in their DNA-binding domain (DBD) that forms hydrogen bonds with nucleotides can have pH-regulated activity, which is relevant to more than 85 transcription factors in distinct families, including FOX, KLF, SOX, and MITF/Myc. Focusing on FOX family transcription factors, we use unbiased SELEX-seq to identify pH-dependent DNA-binding motif preferences and confirm pH-regulated binding affinities for FOXC2, FOXM1, and FOXN1 to a canonical FkhP DNA motif that are greater at pH 7.0 compared with pH 7.5 and for FOXN1 to a preferred FHL motif at higher pHi in cells. For FOXC2, we also find that greater activity for an FkhP motif at lower pH is dependent on a conserved histidine (His122) in the DBD. ChIP-seq and RNA-seq with FOXC2 also reveal pH-dependent differences in enriched promoter motifs. Our findings identify pH-regulated transcription factor-DNA binding selectivity with relevance to how pHi dynamics can regulate gene expression for myriad cell behaviours.

Clinical significance of FOXN3 expression in Indian breast cancer patients

Globally, breast cancer is the most common cancer to affect women. There are different molecular and pathological factors that are involved in the uncontrolled proliferation of breast cancer cells. FOXN3 that is member of Forkhead box family proteins is well recognized for having a crucial role in different biological processes and is reported to be dysregulated in various malignancies. The studies to evaluate the significance of the FOXN3 gene in progression of breast cancer are still under progress. We in the current study aim to examine the FOXN3 gene expression in Indian breast cancer patients and find its clinical relevance. FOXN3 expression analysis using RT-PCR, immunohistochemistry, and western blotting was performed in tumor and normal tissue collected from 142 sporadic breast cancer patients. To identify the genetic aberrations in FOXN3 gene automated DNA sequencing was done. FOXN3 expression study revealed the elevated expression of FOXN3 mRNA in 61.26% of the cases whereas FOXN3 protein was seen to be overexpressed in 59.15% cases. Further, it was found that the elevated expression of FOXN3 mRNA was significantly correlated with the post-menopausal (p = 0.003) status and positive lymph node status (p = 0.049) of the patients. The FOXN3 protein expression also exhibited the significant association with menopausal status (p = 0.008), lymph node status (p = 0.001) and clinical stage (p = 0.018) of the patients. However, we did not find any mutation in the DNA binding domain of the FOXN3 gene in the Indian breast cancer cases. Our findings indicates that overexpression of FOXN3 gene in Indian breast cancer cases can have a potential role in breast cancer progression especially in advanced clinical stages.

Mutations to transcription factor MAX allosterically increase DNA selectivity by altering folding and binding pathways

Understanding how proteins discriminate between preferred and non-preferred ligands ('selectivity') is essential for predicting biological function and a central goal of protein engineering efforts, yet the biophysical mechanisms underpinning selectivity remain poorly understood. Towards this end, we study how variants of the promiscuous transcription factor (TF) MAX (H. sapiens) alter DNA specificity and selectivity, yielding >1700 Kds and >500 rate constants in complex with multiple DNA sequences. Twenty-two of the 240 assayed MAX point mutations enhance selectivity, yet none of these mutations occur at residues that contact nucleotides in published structures. By applying thermodynamic and kinetic models to these results and previous observations for the highly similar yet far more selective TF Pho4 (S. cerevisiae), we find that these mutations enhance selectivity by altering partitioning between or affinity within conformations with different intrinsic selectivity, providing a mechanistic basis for allosteric modulation of ligand selectivity. These results highlight the importance of conformational heterogeneity in determining sequence selectivity and can guide future efforts to engineer selective proteins.

A role for pH dynamics regulating transcription factor DNA binding selectivity

Intracellular pH (pHi) dynamics regulates diverse cell processes such as proliferation, dysplasia, and differentiation, often mediated by the protonation state of a functionally critical histidine residue in endogenous pH sensing proteins. How pHi dynamics can directly regulate gene expression and whether transcription factors can function as pH sensors has received limited attention. We tested the prediction that transcription factors with a histidine in their DNA binding domain (DBD) that forms hydrogen bonds with nucleotides can have pH-regulated activity, which is relevant to more than 85 transcription factors in distinct families, including FOX, KLF, SOX and MITF/Myc. Focusing on FOX family transcription factors, we used unbiased SELEX-seq to identify pH-dependent DNA binding motif preferences, then confirm pH-regulated binding affinities for FOXC2, FOXM1, and FOXN1 to a canonical FkhP DNA motif that are 2.5 to 7.5 greater at pH 7.0 compared with pH 7.5. For FOXC2, we also find greater activity for an FkhP motif at lower pHi in cells and that pH-regulated binding and activity are dependent on a conserved histidine (His122) in the DBD. RNA-seq with FOXC2 also reveals pH-dependent differences in enriched promoter motifs. Our findings identify pH-regulated transcription factor-DNA binding selectivity with relevance to how pHi dynamics can regulate gene expression for myriad cell behaviours.

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Unraveling the role of the FHL family in cardiac diseases: Mechanisms, implications, and future directions

Heart diseases are a major cause of morbidity and mortality worldwide. Understanding the molecular mechanisms underlying these diseases is essential for the development of effective diagnostic and therapeutic strategies. The FHL family consists of five members: FHL1, FHL2, FHL3, FHL4, and FHL5/Act. These members exhibit different expression patterns in various tissues including the heart. FHL family proteins are implicated in cardiac remodeling, regulation of metabolic enzymes, and cardiac biomechanical stress perception. A large number of studies have explored the link between FHL family proteins and cardiac disease, skeletal muscle disease, and ovarian metabolism, but a comprehensive and in-depth understanding of the specific molecular mechanisms targeting FHL on cardiac disease is lacking. The aim of this review is to explore the structure and function of FHL family members, to comprehensively elucidate the mechanisms by which they regulate the heart, and to explore in depth the changes in FHL family members observed in different cardiac disorders, as well as the effects of mutations in FHL proteins on heart health.

Short tandem repeats bind transcription factors to tune eukaryotic gene expression

Short tandem repeats (STRs) are enriched in eukaryotic cis-regulatory elements and alter gene expression, yet how they regulate transcription remains unknown. We found that STRs modulate transcription factor (TF)-DNA affinities and apparent on-rates by about 70-fold by directly binding TF DNA-binding domains, with energetic impacts exceeding many consensus motif mutations. STRs maximize the number of weakly preferred microstates near target sites, thereby increasing TF density, with impacts well predicted by statistical mechanics. Confirming that STRs also affect TF binding in cells, neural networks trained only on in vivo occupancies predicted effects identical to those observed in vitro. Approximately 90% of TFs preferentially bound STRs that need not resemble known motifs, providing a cis-regulatory mechanism to target TFs to genomic sites.

Noncanonical binding of transcription factors: time to revisit specificity?

Transcription factors (TFs) are one of the most studied classes of DNA-binding proteins that have a direct functional impact on gene transcription and thus, on human physiology and disease. The mechanisms that TFs use for recognizing target DNA binding sites have been studied for nearly five decades, yet they remain poorly understood. It is classically assumed that a TF recognizes a specific sequence pattern, or motif, as its binding sites. However, recent studies are consistently finding examples of noncanonical binding, that is, TFs binding at sites that do not resemble their sequence motifs. Here we review the current literature on four major types of noncanonical TF binding, namely binding based on DNA shape readout, at Guanine-quadruplex structures, at repeat sequences, and bispecific binding. These examples point to a critical need for studies to unify our current observations, many of which are at odds with the "one TF, one motif" view, into a more comprehensive definition of the DNA-binding specificity of TFs.

DNA binding specificity of all four Saccharomyces cerevisiae forkhead transcription factors

Quantifying the nucleotide preferences of DNA binding proteins is essential to understanding how transcription factors (TFs) interact with their targets in the genome. High-throughput in vitro binding assays have been used to identify the inherent DNA binding preferences of TFs in a controlled environment isolated from confounding factors such as genome accessibility, DNA methylation, and TF binding cooperativity. Unfortunately, many of the most common approaches for measuring binding preferences are not sensitive enough for the study of moderate-to-low affinity binding sites, and are unable to detect small-scale differences between closely related homologs. The Forkhead box (FOX) family of TFs is known to play a crucial role in regulating a variety of key processes from proliferation and development to tumor suppression and aging. By using the high-sequencing depth SELEX-seq approach to study all four FOX homologs in Saccharomyces cerevisiae, we have been able to precisely quantify the contribution and importance of nucleotide positions all along an extended binding site. Essential to this process was the alignment of our SELEX-seq reads to a set of candidate core sequences determined using a recently developed tool for the alignment of enriched k-mers and a newly developed approach for the reprioritization of candidate cores.

De novo distillation of thermodynamic affinity from deep learning regulatory sequence models of in vivo protein-DNA binding

Transcription factors (TF) are proteins that bind DNA in a sequence-specific manner to regulate gene transcription. Despite their unique intrinsic sequence preferences, in vivo genomic occupancy profiles of TFs differ across cellular contexts. Hence, deciphering the sequence determinants of TF binding, both intrinsic and context-specific, is essential to understand gene regulation and the impact of regulatory, non-coding genetic variation. Biophysical models trained on in vitro TF binding assays can estimate intrinsic affinity landscapes and predict occupancy based on TF concentration and affinity. However, these models cannot adequately explain context-specific, in vivo binding profiles. Conversely, deep learning models, trained on in vivo TF binding assays, effectively predict and explain genomic occupancy profiles as a function of complex regulatory sequence syntax, albeit without a clear biophysical interpretation. To reconcile these complementary models of in vitro and in vivo TF binding, we developed Affinity Distillation (AD), a method that extracts thermodynamic affinities de-novo from deep learning models of TF chromatin immunoprecipitation (ChIP) experiments by marginalizing away the influence of genomic sequence context. Applied to neural networks modeling diverse classes of yeast and mammalian TFs, AD predicts energetic impacts of sequence variation within and surrounding motifs on TF binding as measured by diverse in vitro assays with superior dynamic range and accuracy compared to motif-based methods. Furthermore, AD can accurately discern affinities of TF paralogs. Our results highlight thermodynamic affinity as a key determinant of in vivo binding, suggest that deep learning models of in vivo binding implicitly learn high-resolution affinity landscapes, and show that these affinities can be successfully distilled using AD. This new biophysical interpretation of deep learning models enables high-throughput in silico experiments to explore the influence of sequence context and variation on both intrinsic affinity and in vivo occupancy.

Genomic Analysis of Foxp3 Function in Regulatory T Cells

Regulatory T (Treg) cells are critical for tolerance to self-antigens and for preventing autoimmunity. Foxp3 has been identified as a Treg cell lineage-defining transcription factor controlling Treg cell differentiation and function. In this article, we review the current mechanistic and systemic understanding of Foxp3 function enabled by experimental and computational advances in high-throughput genomics.

Acetylation of Checkpoint suppressor 1 enhances its stability and promotes the progression of triple-negative breast cancer

Checkpoint suppressor 1 (CHES1), a transcriptional regulator, had been dysregulated in many types of malignancies including breast cancer, and its expression level is strongly associated with progression and prognosis of patients. However, the underlying regulatory mechanisms of CHES1 expression in the breast cancer and the effects of post-translational modifications (PTMs) on its functional performance remain to be fully investigated. Herein, we found that CHES1 had a high abundance in triple-negative breast cancer (TNBC) and its expression was tightly associated with malignant phenotype and poor outcomes of patients. Furthermore, we confirmed that CHES1 was an acetylated protein and its dynamic modification was mediated by p300 and HDAC1, and CHES1 acetylation enhanced its stability via decreasing its ubiquitination and degradation, which resulted in the high abundance of CHES1 in TNBC. RNA-seq and functional study revealed that CHES1 facilitated the activation of oncogenic genes and pathways leading to proliferation and metastasis of TNBC. Taken together, this research established a novel regulatory role of acetylation on the stability and activity of CHES1. The results demonstrate the significance of CHES1 acetylation and underlying mechanisms in the progression of TNBC, offering new potential candidate for molecular-targeted therapy in breast cancer.

Molecular basis for DNA recognition by the maternal pioneer transcription factor FoxH1

Forkhead box H1 (FoxH1) is an essential maternal pioneer factor during embryonic development that binds to specific GG/GT-containing DNA target sequences. Here we have determined high-resolution structures of three FoxH1 proteins (from human, frog and fish species) and four DNAs to clarify the way in which FoxH1 binds to these sites. We found that the protein-DNA interactions extend to both the minor and major DNA grooves and are thus almost twice as extensive as those of other FOX family members. Moreover, we identified two specific amino acid changes in FoxH1 that allowed the recognition of GG/GT motifs. Consistent with the pioneer factor activity of FoxH1, we found that its affinity for nucleosomal DNA is even higher than for linear DNA fragments. The structures reported herein illustrate how FoxH1 binding to distinct DNA sites provides specificity and avoids cross-regulation by other FOX proteins that also operate during the maternal-zygotic transition and select canonical forkhead sites.

FOXL2 and FOXA1 cooperatively assemble on the TP53 promoter in alternative dimer configurations

Although both the p53 and forkhead box (FOX) family proteins are key transcription factors associated with cancer progression, their direct relationship is unknown. Here, we found that FOX family proteins bind to the non-canonical homotypic cluster of the p53 promoter region (TP53). Analysis of crystal structures of FOX proteins (FOXL2 and FOXA1) bound to the p53 homotypic cluster indicated that they interact with a 2:1 stoichiometry accommodated by FOX-induced DNA allostery. In particular, FOX proteins exhibited distinct dimerization patterns in recognition of the same p53-DNA; dimer formation of FOXA1 involved protein-protein interaction, but FOXL2 did not. Biochemical and biological functional analyses confirmed the cooperative binding of FOX proteins to the TP53 promoter for the transcriptional activation of TP53. In addition, up-regulation of TP53 was necessary for FOX proteins to exhibit anti-proliferative activity in cancer cells. These analyses reveal the presence of a discrete characteristic within FOX family proteins in which FOX proteins regulate the transcription activity of the p53 tumor suppressor via cooperative binding to the TP53 promoter in alternative dimer configurations.

The transcription factor FoxP3 can fold into two dimerization states with divergent implications for regulatory T cell function and immune homeostasis

FoxP3 is an essential transcription factor (TF) for immunologic homeostasis, but how it utilizes the common forkhead DNA-binding domain (DBD) to perform its unique function remains poorly understood. We here demonstrated that unlike other known forkhead TFs, FoxP3 formed a head-to-head dimer using a unique linker (Runx1-binding region [RBR]) preceding the forkhead domain. Head-to-head dimerization conferred distinct DNA-binding specificity and created a docking site for the cofactor Runx1. RBR was also important for proper folding of the forkhead domain, as truncation of RBR induced domain-swap dimerization of forkhead, which was previously considered the physiological form of FoxP3. Rather, swap-dimerization impaired FoxP3 function, as demonstrated with the disease-causing mutation R337Q, whereas a swap-suppressive mutation largely rescued R337Q-mediated functional impairment. Altogether, our findings suggest that FoxP3 can fold into two distinct dimerization states: head-to-head dimerization representing functional specialization of an ancient DBD and swap dimerization associated with impaired functions.

FOXN3 and GDNF Polymorphisms as Common Genetic Factors of Substance Use and Addictive Behaviors

Epidemiological and phenomenological studies suggest shared underpinnings between multiple addictive behaviors. The present genetic association study was conducted as part of the Psychological and Genetic Factors of Addictions study (n = 3003) and aimed to investigate genetic overlaps between different substance use, addictive, and other compulsive behaviors. Association analyses targeted 32 single-nucleotide polymorphisms, potentially addictive substances (alcohol, tobacco, cannabis, and other drugs), and potentially addictive or compulsive behaviors (internet use, gaming, social networking site use, gambling, exercise, hair-pulling, and eating). Analyses revealed 29 nominally significant associations, from which, nine survived an FDRbl correction. Four associations were observed between FOXN3 rs759364 and potentially addictive behaviors: rs759364 showed an association with the frequency of alcohol consumption and mean scores of scales assessing internet addiction, gaming disorder, and exercise addiction. Significant associations were found between GDNF rs1549250, rs2973033, CNR1 rs806380, DRD2/ANKK1 rs1800497 variants, and the “lifetime other drugs” variable. These suggested that genetic factors may contribute similarly to specific substance use and addictive behaviors. Specifically, FOXN3 rs759364 and GDNF rs1549250 and rs2973033 may constitute genetic risk factors for multiple addictive behaviors. Due to limitations (e.g., convenience sampling, lack of structured scales for substance use), further studies are needed. Functional correlates and mechanisms underlying these relationships should also be investigated.

Novel insights from a multiomics dissection of the hayflick limit

The process wherein dividing cells exhaust proliferative capacity and enter into replicative senescence has become a prominent model for cellular aging in vitro. Despite decades of study, this cellular state is not fully understood in culture and even much less so during aging. Here, we revisit Leonard Hayflick’s original observation of replicative senescence in WI-38 human lung fibroblasts equipped with a battery of modern techniques including RNA-seq, single-cell RNA-seq, proteomics, metabolomics, and ATAC-seq. We find evidence that the transition to a senescent state manifests early, increases gradually, and corresponds to a concomitant global increase in DNA accessibility in nucleolar and lamin associated domains. Furthermore, we demonstrate that senescent WI-38 cells acquire a striking resemblance to myofibroblasts in a process similar to the epithelial to mesenchymal transition (EMT) that is regulated by t YAP1/TEAD1 and TGF-β2. Lastly, we show that verteporfin inhibition of YAP1/TEAD1 activity in aged WI-38 cells robustly attenuates this gene expression program.

Regulation of retinal amacrine cell generation by miR-216b and Foxn3

The mammalian retina contains a complex mixture of different types of neurons. We find that microRNA miR-216b is preferentially expressed in postmitotic retinal amacrine cells in the mouse retina, and expression of miR-216a/b and miR-217 in retina depend in part on Ptf1a, a transcription factor required for amacrine cell differentiation. Surprisingly, ectopic expression of miR-216b directed the formation of additional amacrine cells and reduced bipolar neurons in the developing retina. We identify the Foxn3 mRNA as a retinal target of miR-216b by Argonaute PAR-CLIP and reporter analysis. Inhibition of Foxn3, a transcription factor, in the postnatal developing retina by RNAi increased the formation of amacrine cells and reduced bipolar cell formation. Foxn3 disruption by CRISPR in embryonic retinal explants also increased amacrine cell formation, whereas Foxn3 overexpression inhibited amacrine cell formation prior to Ptf1a expression. Co-expression of Foxn3 partially reversed the effects of ectopic miR-216b on retinal cell formation. Our results identify Foxn3 as a novel regulator of interneuron formation in the developing retina and suggest that miR-216b likely regulates Foxn3 and other genes in amacrine cells.

Motif models proposing independent and interdependent impacts of nucleotides are related to high and low affinity transcription factor binding sites in Arabidopsis

Position weight matrix (PWM) is the traditional motif model representing the transcription factor (TF) binding sites. It proposes that the positions contribute independently to TFs binding affinity, although this hypothesis does not fit the data perfectly. This explains why PWM hits are missing in a substantial fraction of ChIP-seq peaks. To study various modes of the direct binding of plant TFs, we compiled the benchmark collection of 111 ChIP-seq datasets for Arabidopsis thaliana, and applied the traditional PWM, and two alternative motif models BaMM and SiteGA, proposing the dependencies of the positions. The variation in the stringency of the recognition thresholds for the models proposed that the hits of PWM, BaMM, and SiteGA models are associated with the sites of high/medium, any, and low affinity, respectively. At the medium recognition threshold, about 60% of ChIP-seq peaks contain PWM hits consisting of conserved core consensuses, while BaMM and SiteGA provide hits for an additional 15% of peaks in which a weaker core consensus is compensated through intra-motif dependencies. The presence/absence of these dependencies in the motifs of alternative/traditional models was confirmed by the dependency logo DepLogo visualizing the position-wise partitioning of the alignments of predicted sites. We exemplify the detailed analysis of ChIP-seq profiles for plant TFs CCA1, MYC2, and SEP3. Gene ontology (GO) enrichment analysis revealed that among the three motif models, the SiteGA had the highest portions of genes with the significantly enriched GO terms among all predicted genes. We showed that both alternative motif models provide for traditional PWM greater extensions in predicted sites for TFs MYC2/SEP3 with condition/tissue specific functions, compared to those for TF CCA1 with housekeeping functions. Overall, the combined application of standard and alternative motif models is beneficial to detect various modes of the direct TF-DNA interactions in the maximal portion of ChIP-seq loci.

Mechanistic Insights into the Preference for Tandem Binding Sites in DNA Recognition by FOXM1

FOXM1 is an essential proliferation-associated transcription factor that controls the activation of a number of cell cycle regulatory genes. Unlike other forkhead box (FOX) transcription factors, FOXM1 has been shown to prefer binding tandem regulatory DNA sites. However, the underlying reason for such preference is not clear. Here, we showed that the tandem DNA motif, named DIV2, is widely distributed in the promoter region of FOXM1 target genes. The binding of FOXM1 on the DIV2 site differs dramatically from other sites, which is in a highly cooperative fashion, with a much enhanced thermal stability and can be clearly detected by EMSA. The crystal structure of FOXM1 in complex with the DIV2 DNA reveals that the cooperative binding is likely to be driven by intermolecular protein-protein interactions (PPIs). Further half-site spacer insertion assays showed that FOXM1 can bind another site, DIV0, in a similar manner to the DIV2 site. Given the high occurrence of the DIV2 and DIV0 sites in FOXM1 target genes, our results suggest that FOXM1 prefers tandem DNA sites to enable cooperative DNA recognition, and such binding characteristics may further confer its specificity during transcriptional regulation.

Toward a mechanistic understanding of DNA binding by forkhead transcription factors and its perturbation by pathogenic mutations

Forkhead box (FOX) proteins are an evolutionarily conserved family of transcription factors that play numerous regulatory roles in eukaryotes during developmental and adult life. Dysfunction of FOX proteins has been implicated in a variety of human diseases, including cancer, neurodevelopment disorders and genetic diseases. The FOX family members share a highly conserved DNA-binding domain (DBD), which is essential for DNA recognition, binding and function. Since the first FOX structure was resolved in 1993, >30 FOX structures have been reported to date. It is clear now that the structure and DNA recognition mechanisms vary among FOX members; however, a systematic review on this aspect is lacking. In this manuscript, we present an overview of the mechanisms by which FOX transcription factors bind DNA, including protein structures, DNA binding properties and disease-causing mutations. This review should enable a better understanding of FOX family transcription factors for basic researchers and clinicians.

Application of alternative de novo motif recognition models for analysis of structural heterogeneity of transcription factor binding sites: a case study of FOXA2 binding sites

The most popular model for the search of ChIP-seq data for transcription factor binding sites (TFBS) is the positional weight matrix (PWM). However, this model does not take into account dependencies between nucleotide occurrences in different site positions. Currently, two recently proposed models, BaMM and InMoDe, can do as much. However, application of these models was usually limited only to comparing their recognition accuracies with that of PWMs, while none of the analyses of the co-prediction and relative positioning of hits of different models in peaks has yet been performed. To close this gap, we propose the pipeline called MultiDeNA. This pipeline includes stages of model training, assessing their recognition accuracy, scanning ChIP-seq peaks and their classification based on scan results. We applied our pipeline to 22 ChIP-seq datasets of TF FOXA2 and considered PWM, dinucleotide PWM (diPWM), BaMM and InMoDe models. The combination of these four models allowed a significant increase in the fraction of recognized peaks compared to that for the sole PWM model: the increase was 26.3 %. The BaMM model provided the main contribution to the recognition of sites. Although the major fraction of predicted peaks contained TFBS of different models with coincided positions, the medians of the fraction of peaks containing the predictions of sole models were 1.08, 0.49, 4.15 and 1.73 % for PWM, diPWM, BaMM and InMoDe, respectively. Thus, FOXA2 BSs were not fully described by only a sole model, which indicates theirs heterogeneity. We assume that the BaMM model is the most successful in describing the structure of the FOXA2 BS in ChIP-seq datasets under study.

Systematic identification of non-canonical transcription factor motifs

Sequence-specific transcription factors (TFs) recognize motifs of related nucleotide sequences at their DNA binding sites. Upon binding at these sites, TFs regulate critical molecular processes such as gene expression. It is widely assumed that a TF recognizes a single “canonical” motif, although recent studies have identified additional “non-canonical” motifs for some TFs. A comprehensive approach to identify non-canonical DNA binding motifs and the functional importance of those motifs’ matches in the human genome is necessary for fully understanding the mechanisms of TF-regulated molecular processes in human cells. To address this need, we developed a statistical pipeline for in vitro HT-SELEX data that identifies and characterizes the distributions of non-canonical TF motifs in a stringent manner. Analyzing ~170 human TFs’ HT-SELEX data, we found non-canonical motifs for 19 TFs (11%). These non-canonical motifs occur independently of the TFs’ canonical motifs. Non-canonical motif occurrences in the human genome show similar evolutionary conservation to canonical motif occurrences, explain TF binding in locations without canonical motifs, and occur within gene promoters and epigenetically marked regulatory sequences in human cell lines and tissues. Our approach and collection of non-canonical motifs expand current understanding of functionally relevant DNA binding sites for human TFs.

Mechanism of forkhead transcription factors binding to a novel palindromic DNA site

Forkhead transcription factors bind a canonical consensus DNA motif, RYAAAYA (R = A/G, Y = C/T), as a monomer. However, the molecular mechanisms by which forkhead transcription factors bind DNA as a dimer are not well understood. In this study, we show that FOXO1 recognizes a palindromic DNA element DIV2, and mediates transcriptional regulation. The crystal structure of FOXO1/DIV2 reveals that the FOXO1 DNA binding domain (DBD) binds the DIV2 site as a homodimer. The wing1 region of FOXO1 mediates the dimerization, which enhances FOXO1 DNA binding affinity and complex stability. Further biochemical assays show that FOXO3, FOXM1 and FOXI1 also bind the DIV2 site as homodimer, while FOXC2 can only bind this site as a monomer. Our structural, biochemical and bioinformatics analyses not only provide a novel mechanism by which FOXO1 binds DNA as a homodimer, but also shed light on the target selection of forkhead transcription factors.

Structural Basis for DNA Recognition by FOXG1 and the Characterization of Disease-causing FOXG1 Mutations

Forkhead box G1 (FOXG1) is a transcription factor mainly expressed in the brain that plays a critical role in the development and regionalization of the forebrain. Aberrant expression of FOXG1 has implications in FOXG1 syndrome, a serious neurodevelopmental disorder. Here, we report the crystal structure of the FOXG1 DNA-binding domain (DBD) in complex with the forkhead consensus DNA site DBE2 at the resolution of 1.6 Å. FOXG1-DBD adopts a typical winged helix fold. Compared to those of other FOX-DBD/DBE2 structures, the N terminus, H3 helix and wing2 region of FOXG1-DBD exhibit differences in DNA recognition. The FOXG1-DBD wing2 region adopts a unique architecture composed of two β-strands that differs from all other known FOX-DBD wing2 folds. Mutation assays revealed that the disease-causing mutations within the FOXG1-DBD affect DNA binding, protein thermal stability, or both. Our report provides initial insight into how FOXG1 binds DNA and sheds light on how disease-causing mutations in FOXG1-DBD affect its DNA-binding ability.

FOXN3 suppresses the growth and invasion of papillary thyroid cancer through the inactivation of Wnt/β-catenin pathway

Forkhead box N3 (FOXN3) is a subtype of FOX family that has been demonstrated to be implicated in several cancers. However, the role of FOXN3 in papillary thyroid carcinoma (PTC) and its mechanisms have not yet been investigated. Our results showed that FOXN3 was markedly down regulated in PTC tissues and cell lines. Overexpression of FOXN3 suppressed the proliferation, colony formation, migration, and invasion in PTC cells. Overexpression of FOXN3 also prevented EMT process in PTC cells, as shown by the increased E-cadherin expression level and decreased expression levels of N-cadherin and vimentin. In addition, overexpression of FOXN3 inhibited tumor growth of PTC in vivo. Furthermore, overexpression of FOXN3 caused significant decreases in expression levels of β-catenin, c-Myc, and cyclin D1. Additionally, activation of Wnt/β-catenin pathway reversed the effects of FOXN3 on PTC cells. In conclusion, these findings indicated that FOXN3 exerted a tumor suppressive activity in PTC, which was mediated by Wnt/β-catenin pathway.

The crystal structure of human forkhead box N1 in complex with DNA reveals the structural basis for forkhead box family specificity

Forkhead box N1 (FOXN1) is a member of the forkhead box family of transcription factors and plays an important role in thymic epithelial cell differentiation and development. FOXN1 mutations in humans and mice give rise to the "nude" phenotype, which is marked by athymia. FOXN1 belongs to a subset of the FOX family that recognizes an alternative forkhead-like (FHL) consensus sequence (GACGC) that is different from the more widely recognized forkhead (FKH) sequence RYAAAYA (where R is purine, and Y is pyrimidine). Here, we present the FOXN1 structure in complex with DNA containing an FHL motif at 1.6 Å resolution, in which the DNA sequence is recognized by a mixture of direct and water-mediated contacts provided by residues in an α-helix inserted in the DNA major groove (the recognition helix). Comparisons with the structure of other FOX family members revealed that the FKH and FHL DNA sequences are bound in two distinct modes, with partially different registers for the protein DNA contacts. We identified a single alternative rotamer within the recognition helix itself as an important determinant of DNA specificity and found protein sequence features in the recognition helix that could be used to predict the specificity of other FOX family members. Finally, we demonstrate that the C-terminal region of FOXN1 is required for high-affinity DNA binding and that FOXN1 has a significantly reduced affinity for DNA that contains 5'-methylcytosine, which may have implications for the role of FOXN1 in thymic involution.

Mechanistic insights into transcription factor cooperativity and its impact on protein-phenotype interactions

Recent high-throughput transcription factor (TF) binding assays revealed that TF cooperativity is a widespread phenomenon. However, a global mechanistic and functional understanding of TF cooperativity is still lacking. To address this, here we introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data. We identify DNA shape as driver for cooperativity, with a particularly strong effect for Forkhead-Ets pairs. Follow-up experiments reveal a local shape preference at the Ets-DNA-Forkhead interface and decreased cooperativity upon loss of the interaction. Additionally, we discover many functional associations for cooperatively bound TFs. Examination of the link between FOXO1:ETV6 and lymphomas reveals that their joint expression levels improve patient clinical outcome stratification. Altogether, our results demonstrate that inter-family cooperative TF binding is driven by position-specific DNA readout mechanisms, which provides an additional regulatory layer for downstream biological functions.

Although transcription factor (TF) cooperativity is widespread, a global mechanistic understanding of the role of TF cooperativity is still lacking. Here the authors introduce a statistical learning framework that provides structural insight into TF cooperativity and its functional consequences based on next generation sequencing data and provide mechanistic insights into TF cooperativity and its impact on protein-phenotype interactions.

FOXN3 controls liver glucose metabolism by regulating gluconeogenic substrate selection

The FOXN3 gene locus is associated with fasting blood glucose levels in non-diabetic human population genetic studies. The blood glucose-modifying variation within this gene regulates the abundance of both FOXN3 protein and transcript in primary human hepatocytes, with the hyperglycemia risk allele causing increases in both FOXN3 protein and transcript. Using transgenic and knock-out zebrafish models, we showed previously that FOXN3 is a transcriptional repressor that regulates fasting blood glucose by altering liver gene expression of MYC, a  master transcriptional regulator of glucose utilization, and by modulating pancreatic α cell mass and function through an unknown mechanism. Since homozygous Foxn3 null mice die perinatally, and heterozygous carries of the null allele are smaller than wild-type siblings, we examine the metabolic effects of decreasing mouse liver Foxn3 expression in adult life, performing dynamic endocrine tests not feasible in adult zebrafish. Fasting glucose, glucagon, and insulin; and dynamic responses to glucose, insulin, pyruvate, glutamine, and glucagon were measured. Gluconeogenic and amino acid catabolic gene expression was examined in livers, as well. Knocking down liver Foxn3 expression via transduction with adeno-associated virus serotype 8 particles encoding a short hairpin RNA targeting Fonx3 decreases fasting glucose and increases Myc expression, without altering fasting glucagon or fasting insulin. Liver Foxn3 knock-down confers increases glucose tolerance, has no effect on insulin tolerance or response to glucagon challenge, blunts pyruvate and glutamine tolerance, and modulates expression of amino acid transporters and catabolic enzymes. We conclude that liver Foxn3 regulates substrate selection for gluconeogenesis.

Recent Advances in Understanding FOXN3 in Breast Cancer, and Other Malignancies

FOXN3 (forkhead box N3; CHES1: check point suppressor 1) belongs to the forkhead box (FOX) protein family. FOXN3 displays transcriptional inhibitory activity, and is involved in cell cycle regulation and tumorigenesis. FOXN3 is a tumor suppresser and alterations in FOXN3 are found in of a variety of cancers including melanoma, osteosarcoma, and hepatocellular carcinoma. While the roles of FOXN3 role in some cancers have been explored, its role in breast cancer remains unclear. Here we describe current state of knowledge of FOXN3 functions, and focus on its roles (known and potential) in breast cancer.