Analysis of homeodomain specificities allows the family-wide prediction of preferred recognition sites

LMX1B missense-perturbation of regulatory element footprints disrupts serotonergic forebrain axon arborization

Pathogenic coding mutations are prevalent in human neuronal transcription factors (TFs) but how they disrupt development is poorly understood. Lmx1b is a master transcriptional regulator of postmitotic Pet1 neurons that give rise to mature serotonin (5-HT) neurons; over two hundred pathogenic heterozygous mutations have been discovered in human LMX1B, yet their impact on brain development has not been investigated. Here, we developed mouse models with different LMX1B DNA-binding missense mutations. Missense heterozygosity broadly altered Pet1 neuron transcriptomes, but expression changes converged on axon and synapse genes. Missense heterozygosity effected highly specific deficits in the postnatal maturation of forebrain serotonin axon arbors, primarily in the hippocampus and motor cortex, which was associated with spatial memory defects. Digital genomic footprinting (DGF) revealed that missense heterozygosity caused complete loss of Lmx1b motif protection and chromatin accessibility at sites enriched for a distal active enhancer/active promoter histone signature and homeodomain binding motifs; at other bound Lmx1b motifs, varying levels of losses, gains, or no change in motif binding and accessibility were found. The spectrum of footprint changes was strongly associated with synapse and axon genes. Further, Lmx1b missense heterozygosity caused wide disruption of Lmx1b-dependent GRNs comprising diverse TFs expressed in Pet1 neurons. These findings reveal an unanticipated continuum of Lmx1b missense-forced perturbations on Pet1 neuron regulatory element TF binding and accessibility. Our work illustrates DGF's utility for gaining unique insight into how expressed TF missense mutations interfere with developing neuronal GRNs.

Predicting the DNA binding specificity of transcription factor mutants using family-level biophysically interpretable machine learning

Sequence-specific interactions of transcription factors (TFs) with genomic DNA underlie many cellular processes. High-throughput in vitro binding assays coupled with machine learning have made it possible to accurately define such molecular recognition in a biophysically interpretable way for hundreds of TFs across many structural families, providing new avenues for predicting how the sequence preference of a TF is impacted by disease-associated mutations in its DNA binding domain. We developed a method based on a reference-free tetrahedral representation of variation in base preference within a given structural family that can be used to accurately predict the effect of mutations in the protein sequence of the TF. Using the basic helix-loop-helix (bHLH) and homeodomain families as test cases, our results demonstrate the feasibility of accurately predicting the shifts (ΔΔΔG/RT) in binding free energy associated with TF mutants by leveraging high-quality DNA binding models for sets of homologous wild-type TFs.

Experimental approaches to investigate biophysical interactions between homeodomain transcription factors and DNA

Homeodomain transcription factors (TFs) bind to specific DNA sequences to regulate the expression of target genes. Structural work has provided insight into molecular identities and aided in unraveling structural features of these TFs. However, the detailed affinity and specificity by which these TFs bind to DNA sequences is still largely unknown. Qualitative methods, such as DNA footprinting, Electrophoretic Mobility Shift Assays (EMSAs), Systematic Evolution of Ligands by Exponential Enrichment (SELEX), Bacterial One Hybrid (B1H) systems, Surface Plasmon Resonance (SPR), and Protein Binding Microarrays (PBMs) have been widely used to investigate the biochemical characteristics of TF-DNA binding events. In addition to these qualitative methods, bioinformatic approaches have also assisted in TF binding site discovery. Here we discuss the advantages and limitations of these different approaches, as well as the benefits of utilizing more quantitative approaches, such as Mechanically Induced Trapping of Molecular Interactions (MITOMI), Microscale Thermophoresis (MST) and Isothermal Titration Calorimetry (ITC), in determining the biophysical basis of binding specificity of TF-DNA complexes and improving upon existing computational approaches aimed at affinity predictions.

Aberrant homeodomain-DNA cooperative dimerization underlies distinct developmental defects in two dominant CRX retinopathy models

Paired-class homeodomain (HD) transcription factors (TFs) play essential roles in vertebrate development, and their mutations are linked to human diseases. One unique feature of a paired-class HD is cooperative dimerization on specific palindrome DNA sequences. Yet, the functional significance of HD cooperative dimerization in animal development and its dysregulation in diseases remains elusive. Using the retinal TF cone-rod homeobox (CRX) as a model, we have studied how blindness-causing mutations in the paired HD, p.E80A and p.K88N, alter CRX's cooperative dimerization, leading to gene misexpression and photoreceptor developmental deficits in dominant manners. CRXE80A maintains binding at monomeric WT CRX motifs but is deficient in cooperative binding at dimeric motifs. CRXE80A's cooperativity defect impacts the exponential increase of photoreceptor gene expression in terminal differentiation and produces immature, nonfunctional photoreceptors in the Crx E80A retinas. CRXK88N is highly cooperative and localizes to ectopic genomic sites with strong enrichment of dimeric HD motifs. CRXK88N's altered biochemical properties disrupt CRX's ability to direct dynamic chromatin remodeling during development to activate photoreceptor differentiation programs and silence progenitor programs. Our study provides in vitro and in vivo molecular evidence that paired-class HD cooperative dimerization regulates neuronal development and that dysregulation of cooperative binding contributes to severe dominant blinding retinopathies.

Protein stability is determined by single-site bias rather than pairwise covariance

The biases revealed in protein sequence alignments have been shown to provide information related to protein structure, stability, and function. For example, sequence biases at individual positions can be used to design consensus proteins that are often more stable than naturally occurring counterparts. Likewise, correlations between pairs of residue can be used to predict protein structures. Recent work using Potts models show that together, single-site biases and pair correlations lead to improved predictions of protein fitness, activity, and stability. Here we use a Potts model to design groups of protein sequences with different amounts of single-site biases and pair correlations, and determine the thermodynamic stabilities of a representative set of sequences from each group. Surprisingly, sequences excluding pair correlations maximize stability, whereas sequences that maximize pair correlations are less stable, suggesting that pair correlations contribute to another aspect of protein fitness. Consistent with this interpretation, we find that for adenylate kinase, enzyme activity is greatly increased by maximizing pair correlations. The finding that elimination of covariant residue pairs increases protein stability suggests a route to enhance stability of designed proteins; indeed, this strategy produces hyperstable homeodomain and adenylate kinase proteins that retain significant activity.

Two coacting shadow enhancers regulate twin of eyeless expression during early Drosophila development

The Drosophila PAX6 homolog twin of eyeless (toy) sits at the pinnacle of the genetic pathway controlling eye development, the retinal determination network. Expression of toy in the embryo is first detectable at cellular blastoderm stage 5 in an anterior-dorsal band in the presumptive procephalic neuroectoderm, which gives rise to the primordia of the visual system and brain. Although several maternal and gap transcription factors that generate positional information in the embryo have been implicated in controlling toy, the regulation of toy expression in the early embryo is currently not well characterized. In this study, we adopt an integrated experimental approach utilizing bioinformatics, molecular genetic testing of putative enhancers in transgenic reporter gene assays and quantitative analysis of expression patterns in the early embryo, to identify 2 novel coacting enhancers at the toy gene. In addition, we apply mathematical modeling to dissect the regulatory landscape for toy. We demonstrate that relatively simple thermodynamic-based models, incorporating only 5 TF binding sites, can accurately predict gene expression from the 2 coacting enhancers and that the HUNCHBACK TF plays a critical regulatory role through a dual-modality function as an activator and repressor. Our analysis also reveals that the molecular architecture of the 2 enhancers is very different, indicating that the underlying regulatory logic they employ is distinct.

Ambivalent partnership of the Drosophila posterior class Hox protein Abdominal-B with Extradenticle and Homothorax

Hox proteins, a sub-group of the homeodomain (HD) transcription factor family, provide positional information for axial patterning in development and evolution. Hox protein functional specificity is reached, at least in part, through interactions with Pbc (Extradenticle (Exd) in Drosophila) and Meis/Prep (Homothorax (Hth) in Drosophila) proteins. Most of our current knowledge of Hox protein specificity stems from the study of anterior and central Hox proteins, identifying the molecular and structural bases for Hox/Pbc/Meis-Prep cooperative action. Posterior Hox class proteins, Abdominal-B (Abd-B) in Drosophila and Hox9-13 in vertebrates, have been comparatively less studied. They strongly diverge from anterior and central class Hox proteins, with a low degree of HD sequence conservation and the absence of a core canonical Pbc interaction motif. Here we explore how Abd-B function interface with that of Exd/Hth using several developmental contexts, studying mutual expression control, functional dependency and intrinsic protein requirements. Results identify cross-regulatory interactions setting relative expression and activity levels required for proper development. They also reveal organ-specific requirement and a binary functional interplay with Exd and Hth, either antagonistic, as previously reported, or synergistic. This highlights context specific use of Exd/Hth, and a similar context specific use of Abd-B intrinsic protein requirements.

LMX1B missense-perturbation of regulatory element footprints disrupts serotonergic forebrain axon arborization

Pathogenic coding mutations are prevalent in human neuronal transcription factors (TFs) but how they disrupt development is poorly understood. Lmx1b is a master transcriptional regulator of postmitotic Pet1 neurons that give rise to mature serotonin (5-HT) neurons; over two hundred pathogenic heterozygous mutations have been discovered in human LMX1B, yet their impact on brain development has not been investigated. Here, we developed mouse models with different LMX1B DNA-binding missense mutations. Missense heterozygosity broadly altered Pet1 neuron transcriptomes, but expression changes converged on axon and synapse genes. Missense heterozygosity effected highly specific deficits in the postnatal maturation of forebrain serotonin axon arbors, primarily in the hippocampus and motor cortex, which was associated with spatial memory defects. Digital genomic footprinting (DGF) revealed that missense heterozygosity caused complete loss of Lmx1b motif protection and chromatin accessibility at sites enriched for a distal active enhancer/active promoter histone signature and homeodomain binding motifs; at other bound Lmx1b motifs, varying levels of losses, gains or no change in motif binding and accessibility were found. The spectrum of footprint changes was strongly associated with synapse and axon genes. Further, Lmx1b missense heterozygosity caused wide disruption of Lmx1b-dependent GRNs comprising diverse TFs expressed in Pet1 neurons. These findings reveal an unanticipated continuum of Lmx1b missense-forced perturbations on Pet1 neuron regulatory element TF binding and accessibility. Our work illustrates the power of DGF for gaining unique insight into how TF missense mutations interfere with developing neuronal GRNs.

The functional evolution of collembolan Ubx on the regulation of abdominal appendage formation

Folsomia candida is a tiny soil-living arthropod belonging to the Collembola, which is an outgroup to Insecta. It resembles insects as having a pair of antennae and three pairs of thorax legs, while it also possesses three abdominal appendages: a ventral tube located in the first abdominal segment (A1), a retinaculum in A3, and a furca in A4. Collembolan Ubx and AbdA specify abdominal appendages, but they are unable to repress appendage marker gene Dll. The genetic basis of collembolan appendage formation and the mechanisms by which Ubx and AbdA regulate Dll transcription and appendage development remains unknown. In this study, we analysed the developmental transcriptomes of F. candida and identified candidate appendage formation genes, including Ubx (FcUbx). The expression data revealed the dominance of Dll over Ubx during the embryonic 3.5 and 4.5 days, suggesting that Ubx is deficient in suppressing Dll at early appendage formation stages. Furthermore, via electrophoretic mobility shift assays and dual luciferase assays, we found that the binding and repression capacity of FcUbx on Drosophila Dll resembles those of the longest isoform of Drosophila Ubx (DmUbx_Ib), while the regulatory mechanism of the C-terminus of FcUbx on Dll repression is similar to that of the crustacean Artemia franciscana Ubx (AfUbx), demonstrating that the function of collembolan Ubx is intermediate between that of Insecta and Crustacea. In summary, our study provides novel insights into collembolan appendage formation and sheds light on the functional evolution of Ubx. Additionally, we propose a model that collembolan Ubx regulates abdominal segments in a context-specific manner.

Mutational scanning of CRX classifies clinical variants and reveals biochemical properties of the transcriptional effector domain

The transcription factor (TF) cone-rod homeobox (CRX) is essential for the differentiation and maintenance of photoreceptor cell identity. Several human CRX variants cause degenerative retinopathies, but most are variants of uncertain significance. We performed a deep mutational scan (DMS) of nearly all possible single amino acid substitutions in CRX using a cell-based transcriptional reporter assay, curating a high-confidence list of nearly 2000 variants with altered transcriptional activity. In the structured homeodomain, activity scores closely aligned to a predicted structure and demonstrated position-specific constraints on amino acid substitution. In contrast, the intrinsically disordered transcriptional effector domain displayed a qualitatively different pattern of substitution effects, following compositional constraints without specific residue position requirements in the peptide chain. These compositional constraints were consistent with the acidic exposure model of transcriptional activation. We evaluated the performance of the DMS assay as a clinical variant classification tool using gold-standard classified human variants from ClinVar, identifying pathogenic variants with high specificity and moderate sensitivity. That this performance could be achieved using a synthetic reporter assay in a foreign cell type, even for a highly cell type-specific TF like CRX, suggests that this approach shows promise for DMS of other TFs that function in cell types that are not easily accessible. Together, the results of the CRX DMS identify molecular features of the CRX effector domain and demonstrate utility for integration into the clinical variant classification pipeline.

DNA sequence and chromatin differentiate sequence-specific transcription factor binding in the human malaria parasite Plasmodium falciparum

Development of the malaria parasite, Plasmodium falciparum, is regulated by a limited number of sequence-specific transcription factors (TFs). However, the mechanisms by which these TFs recognize genome-wide binding sites is largely unknown. To address TF specificity, we investigated the binding of two TF subsets that either bind CACACA or GTGCAC DNA sequence motifs and further characterized two additional ApiAP2 TFs, PfAP2-G and PfAP2-EXP, which bind unique DNA motifs (GTAC and TGCATGCA). We also interrogated the impact of DNA sequence and chromatin context on P. falciparum TF binding by integrating high-throughput in vitro and in vivo binding assays, DNA shape predictions, epigenetic post-translational modifications, and chromatin accessibility. We found that DNA sequence context minimally impacts binding site selection for paralogous CACACA-binding TFs, while chromatin accessibility, epigenetic patterns, co-factor recruitment, and dimerization correlate with differential binding. In contrast, GTGCAC-binding TFs prefer different DNA sequence context in addition to chromatin dynamics. Finally, we determined that TFs that preferentially bind divergent DNA motifs may bind overlapping genomic regions due to low-affinity binding to other sequence motifs. Our results demonstrate that TF binding site selection relies on a combination of DNA sequence and chromatin features, thereby contributing to the complexity of P. falciparum gene regulatory mechanisms.

Multimodal Hox5 activity generates motor neuron diversity

Motor neurons (MNs) are the final output of circuits driving fundamental behaviors, such as respiration and locomotion. Hox proteins are essential in generating the MN diversity required for accomplishing these functions, but the transcriptional mechanisms that enable Hox paralogs to assign distinct MN subtype identities despite their promiscuous DNA binding motif are not well understood. Here we show that Hoxa5 modifies chromatin accessibility in all mouse spinal cervical MN subtypes and engages TALE co-factors to directly bind and regulate subtype-specific genes. We identify a paralog-specific interaction of Hoxa5 with the phrenic MN-specific transcription factor Scip and show that heterologous expression of Hoxa5 and Scip is sufficient to suppress limb-innervating MN identity. We also demonstrate that phrenic MN identity is stable after Hoxa5 downregulation and identify Klf proteins as potential regulators of phrenic MN maintenance. Our data identify multiple modes of Hoxa5 action that converge to induce and maintain MN identity.

Geometric deep learning of protein-DNA binding specificity

Predicting protein-DNA binding specificity is a challenging yet essential task for understanding gene regulation. Protein-DNA complexes usually exhibit binding to a selected DNA target site, whereas a protein binds, with varying degrees of binding specificity, to a wide range of DNA sequences. This information is not directly accessible in a single structure. Here, to access this information, we present Deep Predictor of Binding Specificity (DeepPBS), a geometric deep-learning model designed to predict binding specificity from protein-DNA structure. DeepPBS can be applied to experimental or predicted structures. Interpretable protein heavy atom importance scores for interface residues can be extracted. When aggregated at the protein residue level, these scores are validated through mutagenesis experiments. Applied to designed proteins targeting specific DNA sequences, DeepPBS was demonstrated to predict experimentally measured binding specificity. DeepPBS offers a foundation for machine-aided studies that advance our understanding of molecular interactions and guide experimental designs and synthetic biology.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Cooperative Gsx2-DNA binding requires DNA bending and a novel Gsx2 homeodomain interface

The conserved Gsx homeodomain (HD) transcription factors specify neural cell fates in animals from flies to mammals. Like many HD proteins, Gsx factors bind A/T-rich DNA sequences prompting the following question: How do HD factors that bind similar DNA sequences in vitro regulate specific target genes in vivo? Prior studies revealed that Gsx factors bind DNA both as a monomer on individual A/T-rich sites and as a cooperative homodimer to two sites spaced precisely 7 bp apart. However, the mechanistic basis for Gsx-DNA binding and cooperativity is poorly understood. Here, we used biochemical, biophysical, structural and modeling approaches to (i) show that Gsx factors are monomers in solution and require DNA for cooperative complex formation, (ii) define the affinity and thermodynamic binding parameters of Gsx2/DNA interactions, (iii) solve a high-resolution monomer/DNA structure that reveals that Gsx2 induces a 20° bend in DNA, (iv) identify a Gsx2 protein-protein interface required for cooperative DNA binding and (v) determine that flexible spacer DNA sequences enhance Gsx2 cooperativity on dimer sites. Altogether, our results provide a mechanistic basis for understanding the protein and DNA structural determinants that underlie cooperative DNA binding by Gsx factors.

Discerning the Role of DNA Sequence, Shape, and Flexibility in Recognition by Drosophila Transcription Factors

The precise spatial and temporal orchestration of gene expression is crucial for the ontogeny of an organism and is mainly governed by transcription factors (TFs). The mechanism of recognition of cognate sites amid millions of base pairs in the genome by TFs is still incompletely understood. In this study, we focus on DNA sequence composition, shape, and flexibility preferences of 28 quintessential TFs from Drosophila melanogaster that are critical to development and body patterning mechanisms. Our study finds that TFs exhibit distinct predilections for DNA shape, flexibility, and sequence compositions in the proximity of transcription factor binding sites (TFBSs). Notably, certain zinc finger proteins prefer GC-rich areas with less negative propeller twist, while homeodomains mainly seek AT-rich regions with a more negative propeller twist at their sites. Intriguingly, while numerous cofactors share similar binding site preferences and bind closer to each other in the genome, some cofactors that have different preferences bind farther apart. These findings shed light on TF DNA recognition and provide novel insights into possible cofactor binding and transcriptional regulation mechanisms.

The intrinsically disordered transactivation region of HOXA9 regulates its function by auto-inhibition of its DNA-binding activity

HOXA9 transcription factor is expressed in hematopoietic stem cells and is involved in the regulation of their differentiation and maturation to various blood cells. HOXA9 is linked to various leukemia and is a marker for poor prognosis of acute myeloid leukemia (AML). This protein has a conserved DNA-binding homeodomain and a transactivation domain. We show that this N-terminal transactivation domain is intrinsically disordered and inhibits DNA-binding by the homeodomain. Using NMR spectroscopy and molecular dynamics simulation, we show that the hexapeptide 197AANWLH202 in the disordered region transiently occludes the DNA-binding interface. The hexapeptide also forms a rigid segment, as determined by NMR dynamics, in an otherwise flexible disordered region. Interestingly, this hexapeptide is known to mediate the interaction of HOXA9 and its TALE partner proteins, such as PBX1, and help in cooperative DNA binding. Mutation of tryptophan to alanine in the hexapeptide abrogates the DNA-binding auto-inhibition. We propose that the disordered transactivation region plays a dual role in the regulation of HOXA9 function. In the absence of TALE partners, it inhibits DNA binding, and in the presence of TALE partners it interacts with the TALE protein and facilitates the cooperative DNA binding by the HOX-TALE complex.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

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.

Mutational scanning of CRX classifies clinical variants and reveals biochemical properties of the transcriptional effector domain

Cone-Rod Homeobox, encoded by CRX, is a transcription factor (TF) essential for the terminal differentiation and maintenance of mammalian photoreceptors. Structurally, CRX comprises an ordered DNA-binding homeodomain and an intrinsically disordered transcriptional effector domain. Although a handful of human variants in CRX have been shown to cause several different degenerative retinopathies with varying cone and rod predominance, as with most human disease genes the vast majority of observed CRX genetic variants are uncharacterized variants of uncertain significance (VUS). We performed a deep mutational scan (DMS) of nearly all possible single amino acid substitution variants in CRX, using an engineered cell-based transcriptional reporter assay. We measured the ability of each CRX missense variant to transactivate a synthetic fluorescent reporter construct in a pooled fluorescence-activated cell sorting assay and compared the activation strength of each variant to that of wild-type CRX to compute an activity score, identifying thousands of variants with altered transcriptional activity. We calculated a statistical confidence for each activity score derived from multiple independent measurements of each variant marked by unique sequence barcodes, curating a high-confidence list of nearly 2,000 variants with significantly altered transcriptional activity compared to wild-type CRX. We evaluated the performance of the DMS assay as a clinical variant classification tool using gold-standard classified human variants from ClinVar, and determined that activity scores could be used to identify pathogenic variants with high specificity. That this performance could be achieved using a synthetic reporter assay in a foreign cell type, even for a highly cell type-specific TF like CRX, suggests that this approach shows promise for DMS of other TFs that function in cell types that are not easily accessible. Per-position average activity scores closely aligned to a predicted structure of the ordered homeodomain and demonstrated position-specific residue requirements. The intrinsically disordered transcriptional effector domain, by contrast, displayed a qualitatively different pattern of substitution effects, following compositional constraints without specific residue position requirements in the peptide chain. The observed compositional constraints of the effector domain were consistent with the acidic exposure model of transcriptional activation. Together, the results of the CRX DMS identify molecular features of the CRX effector domain and demonstrate clinical utility for variant classification.

Aberrant homeodomain-DNA cooperative dimerization underlies distinct developmental defects in two dominant CRX retinopathy models

Paired-class homeodomain transcription factors (HD TFs) play essential roles in vertebrate development, and their mutations are linked to human diseases. One unique feature of paired-class HD is cooperative dimerization on specific palindrome DNA sequences. Yet, the functional significance of HD cooperative dimerization in animal development and its dysregulation in diseases remain elusive. Using the retinal TF Cone-rod Homeobox (CRX) as a model, we have studied how blindness-causing mutations in the paired HD, p.E80A and p.K88N, alter CRX's cooperative dimerization, lead to gene misexpression and photoreceptor developmental deficits in dominant manners. CRXE80A maintains binding at monomeric WT CRX motifs but is deficient in cooperative binding at dimeric motifs. CRXE80A's cooperativity defect impacts the exponential increase of photoreceptor gene expression in terminal differentiation and produces immature, non-functional photoreceptors in the CrxE80A retinas. CRXK88N is highly cooperative and localizes to ectopic genomic sites with strong enrichment of dimeric HD motifs. CRXK88N's altered biochemical properties disrupt CRX's ability to direct dynamic chromatin remodeling during development to activate photoreceptor differentiation programs and silence progenitor programs. Our study here provides in vitro and in vivo molecular evidence that paired-class HD cooperative dimerization regulates neuronal development and dysregulation of cooperative binding contributes to severe dominant blinding retinopathies.

CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development

Loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. Here, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. By engineering human induced pluripotent stem cells (hiPSCs) consisting of wild-type (WT), heterozygous (CDX2-Het), and homozygous null CDX2 (CDX2-KO) genotypes, differentiating these cells in a 2D gastruloid model, and subjecting these cells to single-nucleus RNA and ATAC sequencing, we identify several pathways that are dose-dependently regulated by CDX2 including VEGF and non-canonical WNT. snATAC-seq reveals that CDX2-Het cells retain a WT-like chromatin accessibility profile, suggesting accessibility alone is not sufficient to drive this variability in gene expression. Because the loss of CDX2 or TBXT phenocopy one another in vivo, we compared differentially expressed genes in our CDX2-KO to those from TBXT-KO hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and reveal pathways that may underlie the defects in vascular development and allantoic elongation seen in vivo.

Transcriptional precision in photoreceptor development and diseases - Lessons from 25 years of CRX research

The vertebrate retina is made up of six specialized neuronal cell types and one glia that are generated from a common retinal progenitor. The development of these distinct cell types is programmed by transcription factors that regulate the expression of specific genes essential for cell fate specification and differentiation. Because of the complex nature of transcriptional regulation, understanding transcription factor functions in development and disease is challenging. Research on the Cone-rod homeobox transcription factor CRX provides an excellent model to address these challenges. In this review, we reflect on 25 years of mammalian CRX research and discuss recent progress in elucidating the distinct pathogenic mechanisms of four CRX coding variant classes. We highlight how in vitro biochemical studies of CRX protein functions facilitate understanding CRX regulatory principles in animal models. We conclude with a brief discussion of the emerging systems biology approaches that could accelerate precision medicine for CRX-linked diseases and beyond.

Multimodal Hox5 activity generates motor neuron diversity

Motor neurons (MNs) are the final output of circuits driving fundamental behaviors, such as respiration and locomotion. Hox proteins are essential in generating the MN diversity required for accomplishing these functions, but the transcriptional mechanisms that enable Hox paralogs to assign distinct MN subtype identities despite their promiscuous DNA binding motif are not well understood. Here we show that Hoxa5 controls chromatin accessibility in all mouse spinal cervical MN subtypes and engages TALE co-factors to directly bind and regulate subtype-specific genes. We identify a paralog-specific interaction of Hoxa5 with the phrenic MN-specific transcription factor Scip and show that heterologous expression of Hoxa5 and Scip is sufficient to suppress limb-innervating MN identity. We also demonstrate that phrenic MN identity is stable after Hoxa5 downregulation and identify Klf proteins as potential regulators of phrenic MN maintenance. Our data identify multiple modes of Hoxa5 action that converge to induce and maintain MN identity.

CDX2 dose-dependently influences the gene regulatory network underlying human extraembryonic mesoderm development

Proper regulation of gene dosage is critical for the development of the early embryo and the extraembryonic tissues that support it. Specifically, loss of Cdx2 in vivo leads to stunted development of the allantois, an extraembryonic mesoderm-derived structure critical for nutrient delivery and waste removal in the early embryo. In this study, we investigate how CDX2 dose-dependently influences the gene regulatory network underlying extraembryonic mesoderm development. We generate an allelic series for CDX2 in human induced pluripotent stem cells (hiPSCs) consisting of WT, heterozygous, and homozygous null CDX2 genotypes, differentiate these cells in a 2D gastruloid model, and subject these cells to multiomic single nucleus RNA and ATAC sequencing. We identify several genes that CDX2 dose-dependently regulate cytoskeletal integrity and adhesiveness in the extraembryonic mesoderm population, including regulators of the VEGF, canonical WNT, and non-canonical WNT signaling pathways. Despite these dose-dependent gene expression patterns, snATAC-seq reveals that heterozygous CDX2 expression is capable of inducing a WT-like chromatin accessibility profile, suggesting accessibility is not sufficient to drive gene expression when the CDX2 dosage is reduced. Finally, because the loss of CDX2 or TBXT phenocopy one another in vivo, we compare differentially expressed genes in our CDX2 knock-out model to those from TBXT knock-out hiPSCs differentiated in an analogous experiment. This comparison identifies several communally misregulated genes that are critical for cytoskeletal integrity and tissue permeability, including ANK3 and ANGPT1. Together, these results clarify how CDX2 dose-dependently regulates gene expression in the extraembryonic mesoderm and suggest these genes may underlie the defects in vascular development and allantoic elongation seen in the absence or reduction of CDX2 in vivo.

A Rapid and Simplified Method to Isolate Specific Regulators Based on Biotin-Avidin Binding Affinities in Crops

Transcription factors (TFs) are indispensable components of transcriptional regulatory pathways involved in crop growth and development. Herein, we developed a new method for the identification of upstream TFs specific to genes in crops based on the binding affinities of biotin and avidin. First, we constructed and verified the new biotin and avidin system (BAS) by a coprecipitation assay. Subsequently, the feasibility of DNA-based BAS (DBAS) was further proved by in vivo and in vitro assays. Furthermore, we cloned the promoter of rice OsNRT1.1B and the possible regulators were screened and identified. Additionally, partial candidates were validated by the electrophoresis mobility shift assay (EMSA), yeast one-hybrid, and luciferase activity assays. Remarkably, the results showed that the candidates PIP3 and PIP19 both responded to nitrate immediately and overexpression of PIP3 caused retard growth, which indicates that the candidates are functional and the new DBAS method is useful to isolate regulators in crops.

TALE transcription factors: Cofactors no more

Exd/PBX, Hth/MEIS and PREP proteins belong to the TALE (three-amino-acid loop extension) superclass of transcription factors (TFs) with an atypical homedomain (HD). Originally discovered as "cofactors" to HOX proteins, revisiting their traditional role in light of genome-wide experiments reveals a strong and reproducible pattern of HOX and TALE co-occupancy across diverse embryonic tissues. While confirming that TALE increases HOX specificity and selectivity in vivo, this wider outlook also reveals novel aspects of HOX:TALE collaboration, namely that HOX TFs generally require pre-bound TALE factors to access their functional binding sites in vivo. In contrast to the restricted expression domains of HOX TFs, TALE factors are largely ubiquitous, and PBX and PREP are expressed at the earliest developmental stages. PBX and MEIS control development of many organs and tissues and their dysregulation is associated with congenital disease and cancer. Accordingly, many instances of TALE cooperation with non HOX TFs have been documented in various systems. The model that emerges from these studies is that TALE TFs create a permissive chromatin platform that is selected by tissue-restricted TFs for binding. In turn, HOX and other tissue-restricted TFs selectively convert a ubiquitous pool of low affinity TALE binding events into high confidence, tissue-restricted binding events associated with transcriptional activation. As a result, TALE:TF complexes are associated with active chromatin and domain/lineage-specific gene activity. TALE ubiquitous expression and broad genomic occupancy, as well as the increasing examples of TALE tissue-specific partners, reveal a universal and obligatory role for TALE in the control of tissue and lineage-specific transcriptional programs, beyond their initial discovery as HOX co-factors.

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

DeepPBS: Geometric deep learning for interpretable prediction of protein-DNA binding specificity

Predicting specificity in protein-DNA interactions is a challenging yet essential task for understanding gene regulation. Here, we present Deep Predictor of Binding Specificity (DeepPBS), a geometric deep-learning model designed to predict binding specificity across protein families based on protein-DNA structures. The DeepPBS architecture allows investigation of different family-specific recognition patterns. DeepPBS can be applied to predicted structures, and can aid in the modeling of protein-DNA complexes. DeepPBS is interpretable and can be used to calculate protein heavy atom-level importance scores, demonstrated as a case-study on p53-DNA interface. When aggregated at the protein residue level, these scores conform well with alanine scanning mutagenesis experimental data. The inference time for DeepPBS is sufficiently fast for analyzing simulation trajectories, as demonstrated on a molecular-dynamics simulation of a Drosophila Hox-DNA tertiary complex with its cofactor. DeepPBS and its corresponding data resources offer a foundation for machine-aided protein-DNA interaction studies, guiding experimental choices and complex design, as well as advancing our understanding of molecular interactions.

Cooperative Gsx2-DNA Binding Requires DNA Bending and a Novel Gsx2 Homeodomain Interface

The conserved Gsx homeodomain (HD) transcription factors specify neural cell fates in animals from flies to mammals. Like many HD proteins, Gsx factors bind A/T-rich DNA sequences prompting the question - how do HD factors that bind similar DNA sequences in vitro regulate specific target genes in vivo? Prior studies revealed that Gsx factors bind DNA both as a monomer on individual A/T-rich sites and as a cooperative homodimer to two sites spaced precisely seven base pairs apart. However, the mechanistic basis for Gsx DNA binding and cooperativity are poorly understood. Here, we used biochemical, biophysical, structural, and modeling approaches to (1) show that Gsx factors are monomers in solution and require DNA for cooperative complex formation; (2) define the affinity and thermodynamic binding parameters of Gsx2/DNA interactions; (3) solve a high-resolution monomer/DNA structure that reveals Gsx2 induces a 20° bend in DNA; (4) identify a Gsx2 protein-protein interface required for cooperative DNA binding; and (5) determine that flexible spacer DNA sequences enhance Gsx2 cooperativity on dimer sites. Altogether, our results provide a mechanistic basis for understanding the protein and DNA structural determinants that underlie cooperative DNA binding by Gsx factors, thereby providing a deeper understanding of HD specificity.

Missense mutations in CRX homeodomain cause dominant retinopathies through two distinct mechanisms

Homeodomain transcription factors (HD TFs) are instrumental to vertebrate development. Mutations in HD TFs have been linked to human diseases, but their pathogenic mechanisms remain elusive. Here, we use Cone-Rod Homeobox (CRX) as a model to decipher the disease-causing mechanisms of two HD mutations, p.E80A and p.K88N, that produce severe dominant retinopathies. Through integrated analysis of molecular and functional evidence in vitro and in knock-in mouse models, we uncover two novel gain-of-function mechanisms: p.E80A increases CRX-mediated transactivation of canonical CRX target genes in developing photoreceptors; p.K88N alters CRX DNA-binding specificity resulting in binding at ectopic sites and severe perturbation of CRX target gene expression. Both mechanisms produce novel retinal morphological defects and hinder photoreceptor maturation distinct from loss-of-function models. This study reveals the distinct roles of E80 and K88 residues in CRX HD regulatory functions and emphasizes the importance of transcriptional precision in normal development.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance

PURPOSE

Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis.

METHODS

Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins.

RESULTS

Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis.

CONCLUSION

This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

A universal deep-learning model for zinc finger design enables transcription factor reprogramming

Cys2His2 zinc finger (ZF) domains engineered to bind specific target sequences in the genome provide an effective strategy for programmable regulation of gene expression, with many potential therapeutic applications. However, the structurally intricate engagement of ZF domains with DNA has made their design challenging. Here we describe the screening of 49 billion protein-DNA interactions and the development of a deep-learning model, ZFDesign, that solves ZF design for any genomic target. ZFDesign is a modern machine learning method that models global and target-specific differences induced by a range of library environments and specifically takes into account compatibility of neighboring fingers using a novel hierarchical transformer architecture. We demonstrate the versatility of designed ZFs as nucleases as well as activators and repressors by seamless reprogramming of human transcription factors. These factors could be used to upregulate an allele of haploinsufficiency, downregulate a gain-of-function mutation or test the consequence of regulation of a single gene as opposed to the many genes that a transcription factor would normally influence.

Prediction of cooperative homeodomain DNA binding sites from high-throughput-SELEX data

Homeodomain proteins constitute one of the largest families of metazoan transcription factors. Genetic studies have demonstrated that homeodomain proteins regulate many developmental processes. Yet, biochemical data reveal that most bind highly similar DNA sequences. Defining how homeodomain proteins achieve DNA binding specificity has therefore been a long-standing goal. Here, we developed a novel computational approach to predict cooperative dimeric binding of homeodomain proteins using High-Throughput (HT) SELEX data. Importantly, we found that 15 of 88 homeodomain factors form cooperative homodimer complexes on DNA sites with precise spacing requirements. Approximately one third of the paired-like homeodomain proteins cooperatively bind palindromic sequences spaced 3 bp apart, whereas other homeodomain proteins cooperatively bind sites with distinct orientation and spacing requirements. Combining structural models of a paired-like factor with our cooperativity predictions identified key amino acid differences that help differentiate between cooperative and non-cooperative factors. Finally, we confirmed predicted cooperative dimer sites in vivo using available genomic data for a subset of factors. These findings demonstrate how HT-SELEX data can be computationally mined to predict cooperativity. In addition, the binding site spacing requirements of select homeodomain proteins provide a mechanism by which seemingly similar AT-rich DNA sequences can preferentially recruit specific homeodomain factors.

Missense mutations in CRX homeodomain cause dominant retinopathies through two distinct mechanisms

Homeodomain transcription factors (HD TFs) are instrumental to vertebrate development. Mutations in HD TFs have been linked to human diseases, but their pathogenic mechanisms remain elusive. Here we use Cone-Rod Homeobox (CRX) as a model to decipher the disease-causing mechanisms of two HD mutations, p.E80A and p.K88N, that produce severe dominant retinopathies. Through integrated analysis of molecular and functional evidence in vitro and in knock-in mouse models, we uncover two novel gain-of-function mechanisms: p.E80A increases CRX-mediated transactivation of canonical CRX target genes in developing photoreceptors; p.K88N alters CRX DNA-binding specificity resulting in binding at ectopic sites and severe perturbation of CRX target gene expression. Both mechanisms produce novel retinal morphological defects and hinder photoreceptor maturation distinct from loss-of-function models. This study reveals the distinct roles of E80 and K88 residues in CRX HD regulatory functions and emphasizes the importance of transcriptional precision in normal development.

Intra-Domain Residue Coevolution in Transcription Factors Contributes to DNA Binding Specificity

Understanding the basis of the DNA-binding specificity of transcription factors (TFs) has been of long-standing interest. Despite extensive efforts to map millions of putative TF binding sequences, identifying the critical determinants for DNA binding specificity remains a major challenge. The coevolution of residues in proteins occurs due to a shared evolutionary history. However, it is unclear how coevolving residues in TFs contribute to DNA binding specificity. Here, we systematically collected publicly available data sets from multiple large-scale high-throughput TF-DNA interaction screening experiments for the major TF families with large numbers of TF members. These families included the Homeobox, HLH, bZIP_1, Ets, HMG_box, ZF-C4, and Zn_clus TFs. We detected TF subclass-determining sites (TSDSs) and showed that the TSDSs were more likely to coevolve with other TSDSs than with non-TSDSs, particularly for the Homeobox, HLH, Ets, bZIP_1, and HMG_box TF families. By in silico modeling, we showed that mutation of the highly coevolving residues could significantly reduce the stability of the TF-DNA complex. The distant residues from the DNA interface also contributed to TF-DNA binding activity. Overall, our study gave evidence that coevolved residues relate to transcriptional regulation and provided insights into the potential application of engineered DNA-binding domains and proteins. IMPORTANCE While unraveling DNA-binding specificity of TFs is the key to understanding the basis and molecular mechanism of gene expression regulation, identifying the critical determinants that contribute to DNA binding specificity remains a major challenge. In this study, we provided evidence showing that coevolving residues in TF domains contributed to DNA binding specificity. We demonstrated that the TSDSs were more likely to coevolve with other TSDSs than with non-TSDSs. Mutation of the coevolving residue pairs (CRPs) could significantly reduce the stability of THE TF-DNA complex, and even the distant residues from the DNA interface contribute to TF-DNA binding activity. Collectively, our study expands our knowledge of the interactions among coevolved residues in TFs, tertiary contacting, and functional importance in refined transcriptional regulation. Understanding the impact of coevolving residues in TFs will help understand the details of transcription of gene regulation and advance the application of engineered DNA-binding domains and protein.

The pioneering function of the hox transcription factors

Ever since the discovery that the Hox family of transcription factors establish morphological diversity in the developing embryo, major efforts have been directed towards understanding Hox-dependent patterning. This has led to important discoveries, notably on the mechanisms underlying the collinear expression of Hox genes and Hox binding specificity. More recently, several studies have provided evidence that Hox factors have the capacity to bind their targets in an inaccessible chromatin context and trigger the switch to an accessible, transcriptional permissive, chromatin state. In this review, we provide an overview of the evidences supporting that Hox factors behave as pioneer factors and discuss the potential mechanisms implicated in Hox pioneer activity as well as the significance of this functional property in Hox-dependent patterning.

Hox function and specificity – A tissue centric view

Since their discovery, the Hox genes, with their incredible power to reprogram the identity of complete body regions, a phenomenon called homeosis, have captured the fascination of many biologists. Recent research has provided new insights into the function of Hox proteins in different germ layers and the mechanisms they employ to control tissue morphogenesis. We focus in this review on the ectoderm and mesoderm to highlight new findings and discuss them with regards to established concepts of Hox target gene regulation. Furthermore, we highlight the molecular mechanisms involved the transcriptional repression of specific groups of Hox target genes, and summarize the role of Hox mediated gene silencing in tissue development. Finally, we reflect on recent findings identifying a large number of tissue-specific Hox interactor partners, which open up new avenues and directions towards a better understanding of Hox function and specificity in different tissues.

Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation

Sequence-specific recognition of transcription factor (TF) binding motifs in the target site of DNA over the vast amount of non-target DNA is of primary importance for the transcriptional regulation of gene expression by the TFs. Binding of TFs to the target site of DNA relies not only on the direct contact formation but also on the structural and conformational features of DNA. Recognition of DNA structural features or shape readout by proteins is an important factor in the context of TF-DNA interaction. Based on the atomistic molecular simulation, here we report the sequence-dependent unique structural features, solvation, and ion-binding properties of biologically relevant AT- and GC-rich human TF binding motifs in DNA. Counterion and water distribution around the motif is found to be sensitive to the motif sequence, which is accompanied with the DNA shape features. The motif sequence affects the electrostatic potential along the grooves, and cytosine methylation alters the DNA shape features. Characteristic solvation properties of TF binding motif DNA fragments infer that an ionic environment and hydration influences are essential to describe TF-DNA interactions.

Singular value decomposition of protein sequences as a method to visualize sequence and residue space

Singular value decomposition (SVD) of multiple sequence alignments (MSAs) is an important and rigorous method to identify subgroups of sequences within the MSA, and to extract consensus and covariance sequence features that define the alignment and distinguish the subgroups. This information can be correlated to structure, function, stability, and taxonomy. However, the mathematics of SVD is unfamiliar to many in the field of protein science. Here, we attempt to present an intuitive yet comprehensive description of SVD analysis of MSAs. We begin by describing the underlying mathematics of SVD in a way that is both rigorous and accessible. Next, we use SVD to analyze sequences generated with a simplified model in which the extent of sequence conservation and covariance between different positions is controlled, to show how conservation and covariance produce features in the decomposed coordinate system. We then use SVD to analyze alignments of two protein families, the homeodomain and the Ras superfamilies. Both families show clear evidence of sequence clustering when projected into singular value space. We use k-means clustering to group MSA sequences into specific clusters, show how the residues that distinguish these clusters can be identified, and show how these clusters can be related to taxonomy and function. We end by providing a description a set of Python scripts that can be used for SVD analysis of MSAs, displaying results, and identifying and analyzing sequence clusters. These scripts are freely available on GitHub.

Learning probabilistic protein-DNA recognition codes from DNA-binding specificities using structural mappings

Knowledge of how proteins interact with DNA is essential for understanding gene regulation. Although DNA-binding specificities for thousands of transcription factors (TFs) have been determined, the specific amino acid-base interactions comprising their structural interfaces are largely unknown. This lack of resolution hampers attempts to leverage these data in order to predict specificities for uncharacterized TFs or TFs mutated in disease. Here we introduce recognition code learning via automated mapping of protein-DNA structural interfaces (rCLAMPS), a probabilistic approach that uses DNA-binding specificities for TFs from the same structural family to simultaneously infer both which nucleotide positions are contacted by particular amino acids within the TF as well as a recognition code that relates each base-contacting amino acid to nucleotide preferences at the DNA positions it contacts. We apply rCLAMPS to homeodomains, the second largest family of TFs in metazoans and show that it learns a highly effective recognition code that can predict de novo DNA-binding specificities for TFs. Furthermore, we show that the inferred amino acid-nucleotide contacts reveal whether and how nucleotide preferences at individual binding site positions are altered by mutations within TFs. Our approach is an important step toward automatically uncovering the determinants of protein-DNA specificity from large compendia of DNA-binding specificities and inferring the altered functionalities of TFs mutated in disease.

HOX genes in stem cells: Maintaining cellular identity and regulation of differentiation

Stem cells display a unique cell type within the body that has the capacity to self-renew and differentiate into specialized cell types. Compared to pluripotent stem cells, adult stem cells (ASC) such as mesenchymal stem cells (MSCs) and hematopoietic stem cells (HSCs) exhibit restricted differentiation capabilities that are limited to cell types typically found in the tissue of origin, which implicates that there must be a certain code or priming determined by the tissue of origin. HOX genes, a subset of homeobox genes encoding transcription factors that are generally repressed in undifferentiated pluripotent stem cells, emerged here as master regulators of cell identity and cell fate during embryogenesis, and in maintaining this positional identity throughout life as well as specifying various regional properties of respective tissues. Concurrently, intricate molecular circuits regulated by diverse stem cell-typical signaling pathways, balance stem cell maintenance, proliferation and differentiation. However, it still needs to be unraveled how stem cell-related signaling pathways establish and regulate ASC-specific HOX expression pattern with different temporal-spatial topography, known as the HOX code. This comprehensive review therefore summarizes the current knowledge of specific ASC-related HOX expression patterns and how these were integrated into stem cell-related signaling pathways. Understanding the mechanism of HOX gene regulation in stem cells may provide new ways to manipulate stem cell fate and function leading to improved and new approaches in the field of regenerative medicine.

Specificity of the Hox member Deformed is determined by transcription factor levels and binding site affinities

Hox proteins have similar binding specificities in vitro, yet they control different morphologies in vivo. This paradox has been partially solved with the identification of Hox low-affinity binding sites. However, anterior Hox proteins are more promiscuous than posterior Hox proteins, raising the question how anterior Hox proteins achieve specificity. We use the AP2x enhancer, which is activated in the maxillary head segment by the Hox TF Deformed (Dfd). This enhancer lacks canonical Dfd-Exd sites but contains several predicted low-affinity sites. Unexpectedly, these sites are strongly bound by Dfd-Exd complexes and their conversion into optimal Dfd-Exd sites results only in a modest increase in binding strength. These small variations in affinity change the sensitivity of the enhancer to different Dfd levels, resulting in perturbed AP-2 expression and maxillary morphogenesis. Thus, Hox-regulated morphogenesis seems to result from the co-evolution of Hox binding affinity and Hox dosage for precise target gene regulation.

Despite the central role of Hox genes in controlling morphogenesis, the DNA binding of different Hox members is relatively similar. Here they show that specificity of Hox member Dfd relies on a precise balance of transcription factors and binding site affinities.

Hox genes collaborate with helix-loop-helix factor Grainyhead to promote neuroblast apoptosis along the anterior-posterior axis of the Drosophila larval central nervous system

Hox genes code for a family of a homeodomain (HD) containing transcription factors that use TALE-HD containing factors Pbx/Exd and Meis/Hth to specify the development of the anterior-posterior (AP) axis of an organism. However, the absence of TALE-HD containing factors from specific tissues emphasizes the need to identify and validate new Hox cofactors. In Drosophila central nervous system (CNS), Hox execute segment-specific apoptosis of neural stem cells (neuroblasts-NBs) and neurons. In abdominal segments of larval CNS, Hox gene Abdominal-A (AbdA) mediates NB apoptosis with the help of Exd and bHLH factor Grainyhead (Grh) using a 717 bp apoptotic enhancer. In this study, we show that this enhancer is critical for abdominal NB apoptosis and relies on two separable set of DNA binding motifs responsible for its initiation and maintenance. Our results also show that AbdA and Grh interact through their highly conserved DNA binding domains, and the DNA binding specificity of AbdA-HD is important for it to interact with Grh and essential for it to execute NB apoptosis in CNS. We also establish that Grh is required for Hox-dependent NB apoptosis in Labial and Sex Combs Reduced (Scr) expressing regions of the CNS, and it can physically interact with all the Hox proteins in vitro. Our biochemical and functional data collectively support the idea that Grh can function as a Hox cofactor and help them carry out their in vivo roles during development.

Transcription factor paralogs orchestrate alternative gene regulatory networks by context-dependent cooperation with multiple cofactors

In eukaryotes, members of transcription factor families often exhibit similar DNA binding properties in vitro, yet orchestrate paralog-specific gene regulatory networks in vivo. The serially homologous first (T1) and third (T3) thoracic legs of Drosophila, which are specified by the Hox proteins Scr and Ubx, respectively, offer a unique opportunity to address this paradox in vivo. Genome-wide analyses using epitope-tagged alleles of both Hox loci in the T1 and T3 leg imaginal discs, the precursors to the adult legs and ventral body regions, show that ~8% of Hox binding is paralog-specific. Binding specificity is mediated by interactions with distinct cofactors in different domains: the Hox cofactor Exd acts in the proximal domain and is necessary for Scr to bind many of its paralog-specific targets, while in the distal leg domain, the homeodomain protein Distal-less (Dll) enhances Scr binding to a different subset of loci. These findings reveal how Hox paralogs, and perhaps paralogs of other transcription factor families, orchestrate alternative downstream gene regulatory networks with the help of multiple, context-specific cofactors.

Diversification and Functional Evolution of HOX Proteins

Gene duplication and divergence is a major contributor to the generation of morphological diversity and the emergence of novel features in vertebrates during evolution. The availability of sequenced genomes has facilitated our understanding of the evolution of genes and regulatory elements. However, progress in understanding conservation and divergence in the function of proteins has been slow and mainly assessed by comparing protein sequences in combination with in vitro analyses. These approaches help to classify proteins into different families and sub-families, such as distinct types of transcription factors, but how protein function varies within a gene family is less well understood. Some studies have explored the functional evolution of closely related proteins and important insights have begun to emerge. In this review, we will provide a general overview of gene duplication and functional divergence and then focus on the functional evolution of HOX proteins to illustrate evolutionary changes underlying diversification and their role in animal evolution.

Genetic Regulation of Vertebrate Forebrain Development by Homeobox Genes

Forebrain development in vertebrates is regulated by transcription factors encoded by homeobox, bHLH and forkhead gene families throughout the progressive and overlapping stages of neural induction and patterning, regional specification and generation of neurons and glia from central nervous system (CNS) progenitor cells. Moreover, cell fate decisions, differentiation and migration of these committed CNS progenitors are controlled by the gene regulatory networks that are regulated by various homeodomain-containing transcription factors, including but not limited to those of the Pax (paired), Nkx, Otx (orthodenticle), Gsx/Gsh (genetic screened), and Dlx (distal-less) homeobox gene families. This comprehensive review outlines the integral role of key homeobox transcription factors and their target genes on forebrain development, focused primarily on the telencephalon. Furthermore, links of these transcription factors to human diseases, such as neurodevelopmental disorders and brain tumors are provided.

Single‐cell RNA sequencing of motoneurons identifies regulators of synaptic wiring in Drosophila embryos

The correct wiring of neuronal circuits is one of the most complex processes in development, since axons form highly specific connections out of a vast number of possibilities. Circuit structure is genetically determined in vertebrates and invertebrates, but the mechanisms guiding each axon to precisely innervate a unique pre‐specified target cell are poorly understood. We investigated Drosophila embryonic motoneurons using single‐cell genomics, imaging, and genetics. We show that a cell‐specific combination of homeodomain transcription factors and downstream immunoglobulin domain proteins is expressed in individual cells and plays an important role in determining cell‐specific connections between differentiated motoneurons and target muscles. We provide genetic evidence for a functional role of five homeodomain transcription factors and four immunoglobulins in the neuromuscular wiring. Knockdown and ectopic expression of these homeodomain transcription factors induces cell‐specific synaptic wiring defects that are partly phenocopied by genetic modulations of their immunoglobulin targets. Taken together, our data suggest that homeodomain transcription factor and immunoglobulin molecule expression could be directly linked and function as a crucial determinant of neuronal circuit structure.

The power of the (imperfect) palindrome: Sequence-specific roles of palindromic motifs in gene regulation

In human languages, a palindrome reads the same forward as backward (e.g., 'madam'). In regulatory DNA, a palindrome is an inverted sequence repeat that allows a transcription factor to bind as a homodimer or as a heterodimer with another type of transcription factor. Regulatory palindromes are typically imperfect, that is, the repeated sequences differ in at least one base pair, but the functional significance of this asymmetry remains poorly understood. Here, we review the use of imperfect palindromes in Drosophila photoreceptor differentiation and mammalian steroid receptor signaling. Moreover, we discuss mechanistic explanations for the predominance of imperfect palindromes over perfect palindromes in these two gene regulatory contexts. Lastly, we propose to elucidate whether specific imperfectly palindromic variants have specific regulatory functions in steroid receptor signaling and whether such variants can help predict transcriptional outcomes as well as the response of individual patients to drug treatments.