Solution structure of the K50 class homeodomain PITX2 bound to DNA and implications for mutations that cause Rieger syndrome

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.

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.

The interaction landscape between transcription factors and the nucleosome

Nucleosomes cover most of the genome and are thought to be displaced by transcription factors in regions that direct gene expression. However, the modes of interaction between transcription factors and nucleosomal DNA remain largely unknown. Here we systematically explore interactions between the nucleosome and 220 transcription factors representing diverse structural families. Consistent with earlier observations, we find that the majority of the studied transcription factors have less access to nucleosomal DNA than to free DNA. The motifs recovered from transcription factors bound to nucleosomal and free DNA are generally similar. However, steric hindrance and scaffolding by the nucleosome result in specific positioning and orientation of the motifs. Many transcription factors preferentially bind close to the end of nucleosomal DNA, or to periodic positions on the solvent-exposed side of the DNA. In addition, several transcription factors usually bind to nucleosomal DNA in a particular orientation. Some transcription factors specifically interact with DNA located at the dyad position at which only one DNA gyre is wound, whereas other transcription factors prefer sites spanning two DNA gyres and bind specifically to each of them. Our work reveals notable differences in the binding of transcription factors to free and nucleosomal DNA, and uncovers a diverse interaction landscape between transcription factors and the nucleosome.

Diversification of transcription factor-DNA interactions and the evolution of gene regulatory networks

Sequence-specific transcription factors (TFs) bind short DNA sequences in the genome to regulate the expression of target genes. In the last decade, numerous technical advances have enabled the determination of the DNA-binding specificities of many of these factors. Large-scale screens of many TFs enabled the creation of databases of TF DNA-binding specificities, typically represented as position weight matrices (PWMs). Although great progress has been made in determining and predicting binding specificities systematically, there are still many surprises to be found when studying a particular TF's interactions with DNA in detail. Paralogous TFs' binding specificities can differ in subtle ways, in a manner that is not immediately apparent from looking at their PWMs. These differences affect gene regulatory outputs and enable TFs to rewire transcriptional networks over evolutionary time. This review discusses recent observations made in the study of TF-DNA interactions that highlight the importance of continued in-depth analysis of TF-DNA interactions and their inherent complexity. This article is categorized under: Biological Mechanisms > Regulatory Biology.

A massively parallel reporter assay reveals context-dependent activity of homeodomain binding sites in vivo

Cone-rod homeobox (CRX) is a paired-like homeodomain transcription factor (TF) and a master regulator of photoreceptor development in vertebrates. The in vitro DNA binding preferences of CRX have been described in detail, but the degree to which in vitro binding affinity is correlated with in vivo enhancer activity is not known. In addition, paired-class homeodomain TFs can bind DNA cooperatively as both homodimers and heterodimers at inverted TAAT half-sites separated by 2 or 3 nucleotides. This dimeric configuration is thought to mediate target specificity, but whether monomeric and dimeric sites encode distinct levels of activity is not known. Here, we used a massively parallel reporter assay to determine how local sequence context shapes the regulatory activity of CRX binding sites in mouse photoreceptors. We assayed inactivating mutations in more than 1700 TF binding sites and found that dimeric CRX binding sites act as stronger enhancers than monomeric CRX binding sites. Furthermore, the activity of dimeric half-sites is cooperative, dependent on a strict 3-bp spacing, and tuned by the identity of the spacer nucleotides. Saturating single-nucleotide mutagenesis of 195 CRX binding sites showed that, on average, changes in TF binding site affinity are correlated with changes in regulatory activity, but this relationship is obscured when considering mutations across multiple cis-regulatory elements (CREs). Taken together, these results demonstrate that the activity of CRX binding sites is highly dependent on sequence context, providing insight into photoreceptor gene regulation and illustrating functional principles of homeodomain binding sites that may be conserved in other cell types.

Lysine Side-Chain Dynamics in the Binding Site of Homeodomain/DNA Complexes As Observed by NMR Relaxation Experiments and Molecular Dynamics Simulations

An important but poorly characterized contribution to the thermodynamics of protein-DNA interactions is the loss of entropy that occurs from restricting the conformational freedom of amino acid side chains. The effect of restricting the flexibility of several side chains at a protein-DNA interface may be comparable in many cases to the other factors that determine the binding thermodynamics and may, therefore, play a key role in dictating the binding affinity and/or specificity. Because the entropic contributions, including the presence and influence of side-chain dynamics, are especially difficult to estimate based on structural information, it is important to pursue experimental and theoretical studies that can provide direct information regarding these issues. We report on studies of a model system, the homeodomain/DNA complex, focusing on the Lys50 class of homeodomains where a key lysine residue in position 50 was shown previously to be critical for binding site specificity. NMR methodology was employed for determining the dynamics of lysine side-chain amino groups via 15N relaxation measurements in the Lys50-class homeodomains from the Drosophila protein Bicoid and the human protein Pitx2. In the case of Pitx2, complexes with both a consensus and a nonconsensus DNA binding site were examined. NMR-derived order parameters indicated moderate to substantial conformational freedom for the lysine NH3+ group in the complexes studied. To complement the experimental NMR measurements, molecular dynamics simulations were performed for the consensus complexes to gain further, detailed insights regarding the dynamics of the Lys50 side chain and other important residues in the protein-DNA interface.

Accurate prediction of functional, structural, and stability changes in PITX2 mutations using in silico bioinformatics algorithms

Mutations in PITX2 have been implicated in several genetic disorders, particularly Axenfeld-Rieger syndrome. In order to determine the most reliable bioinformatics tools to assess the likely pathogenicity of PITX2 variants, the results of bioinformatics predictions were compared to the impact of variants on PITX2 structure and function. The MutPred, Provean, and PMUT bioinformatic tools were found to have the highest performance in predicting the pathogenicity effects of all 18 characterized missense variants in PITX2, all with sensitivity and specificity >93%. Applying these three programs to assess the likely pathogenicity of 13 previously uncharacterized PITX2 missense variants predicted 12/13 variants as deleterious, except A30V which was predicted as benign variant for all programs. Molecular modeling of the PITX2 homoedomain predicts that of the 31 known PITX2 variants, L54Q, F58L, V83F, V83L, W86C, W86S, and R91P alter PITX2's structure. In contrast, the remaining 24 variants are not predicted to change PITX2's structure. The results of molecular modeling, performed on all the PITX2 missense mutations located in the homeodomain, were compared with the findings of eight protein stability programs. CUPSAT was found to be the most reliable in predicting the effect of missense mutations on PITX2 stability. Our results showed that for PITX2, and likely other members of this homeodomain transcription factor family, MutPred, Provean, PMUT, molecular modeling, and CUPSAT can reliably be used to predict PITX2 missense variants pathogenicity.

Pitx2 in Embryonic and Adult Myogenesis

Skeletal muscle is a heterogeneous tissue that represents between 30 and 38% of the human body mass and has important functions in the organism, such as maintaining posture, locomotor impulse, or pulmonary ventilation. The genesis of skeletal muscle during embryonic development is a process controlled by an elaborate regulatory network combining the interplay of extrinsic and intrinsic regulatory mechanisms that transform myogenic precursor cells into functional muscle fibers through a finely tuned differentiation program. However, the capacity of generating muscle still remains once these fibers have matured. Adult myogenesis resembles many of the embryonic morphogenetic episodes and depends on the activation of satellite cells that have the potential to differentiate into new muscle fibers. Pitx2 is a member of the bicoid family of homeodomain transcription factors that play an important role in morphogenesis. In the last decade, Pitx2 has emerged as a key element involved in the fine-tuning mechanism that regulates skeletal-muscle development as well as the differentiation and cell fate of satellite cells in adult muscle. Here we present an integrative view of all aspects of embryonic and adult myogenesis in which Pitx2 is involved, from embryonic development to satellite-cell proliferation, fate specification, and differentiation. Those new Pitx2 functions on satellite-cell biology might open new perspectives to develop therapeutic strategies for muscular disorders.

Transcription factor family-specific DNA shape readout revealed by quantitative specificity models

Transcription factors (TFs) achieve DNA‐binding specificity through contacts with functional groups of bases (base readout) and readout of structural properties of the double helix (shape readout). Currently, it remains unclear whether DNA shape readout is utilized by only a few selected TF families, or whether this mechanism is used extensively by most TF families. We resequenced data from previously published HT‐SELEX experiments, the most extensive mammalian TF–DNA binding data available to date. Using these data, we demonstrated the contributions of DNA shape readout across diverse TF families and its importance in core motif‐flanking regions. Statistical machine‐learning models combined with feature‐selection techniques helped to reveal the nucleotide position‐dependent DNA shape readout in TF‐binding sites and the TF family‐specific position dependence. Based on these results, we proposed novel DNA shape logos to visualize the DNA shape preferences of TFs. Overall, this work suggests a way of obtaining mechanistic insights into TF–DNA binding without relying on experimentally solved all‐atom structures.

Regulatory logic driving stable levels of defective proventriculus expression during terminal photoreceptor specification in flies

How differential levels of gene expression are controlled in post-mitotic neurons is poorly understood. In the Drosophila retina, expression of the transcription factor Defective Proventriculus (Dve) at distinct cell type-specific levels is required for terminal differentiation of color- and motion-detecting photoreceptors. Here, we find that the activities of two cis-regulatory enhancers are coordinated to drive dve expression in the fly eye. Three transcription factors act on these enhancers to determine cell-type specificity. Negative autoregulation by Dve maintains expression from each enhancer at distinct homeostatic levels. One enhancer acts as an inducible backup ('dark' shadow enhancer) that is normally repressed but becomes active in the absence of the other enhancer. Thus, two enhancers integrate combinatorial transcription factor input, feedback and redundancy to generate cell type-specific levels of dve expression and stable photoreceptor fate. This regulatory logic may represent a general paradigm for how precise levels of gene expression are established and maintained in post-mitotic neurons.

Dynamic Local Polymorphisms in the Gbx1 Homeodomain Induced by DNA Binding

The Gastrulation Brain Homeobox 1 (Gbx1) gene encodes the Gbx1 homeodomain that targets TAATTA motifs in double-stranded DNA (dsDNA). Residues Glu17 and Arg52 in Gbx1 form a salt bridge, which is preserved in crystal structures and molecular dynamics simulations of homologous homeodomain-DNA complexes. In contrast, our nuclear magnetic resonance (NMR) studies show that DNA binding to Gbx1 induces dynamic local polymorphisms, which include breaking of the Glu17-Arg52 salt bridge. To study this interaction, we produced a variant with Glu17Arg and Arg52Glu mutations, which exhibited the same fold as the wild-type protein, but a 2-fold reduction in affinity for dsDNA. Analysis of the NMR structures of the Gbx1 homeodomain in the free form, the Gbx1[E17R,R52E] variant, and a Gbx1 homeodomain-DNA complex showed that stabilizing interactions of the Arg52 side chain with the DNA backbone are facilitated by transient breakage of the Glu17-Arg52 salt bridge in the DNA-bound Gbx1.

Mutations in zebrafish pitx2 model congenital malformations in Axenfeld-Rieger syndrome but do not disrupt left-right placement of visceral organs

Pitx2 is a conserved homeodomain transcription factor that has multiple functions during embryonic development. Mutations in human PITX2 cause autosomal dominant Axenfeld-Rieger syndrome (ARS), characterized by congenital eye and tooth malformations. Pitx2(-/-) knockout mouse models recapitulate aspects of ARS, but are embryonic lethal. To date, ARS treatments remain limited to managing individual symptoms due to an incomplete understanding of PITX2 function. In addition to regulating eye and tooth development, Pitx2 is a target of a conserved Nodal (TGFβ) signaling pathway that mediates left-right (LR) asymmetry of visceral organs. Based on its highly conserved asymmetric expression domain, the Nodal-Pitx2 axis has long been considered a common denominator of LR development in vertebrate embryos. However, functions of Pitx2 during asymmetric organ morphogenesis are not well understood. To gain new insight into Pitx2 function we used genome editing to create mutations in the zebrafish pitx2 gene. Mutations in the pitx2 homeodomain caused phenotypes reminiscent of ARS, including aberrant development of the cornea and anterior chamber of the eye and reduced or absent teeth. Intriguingly, LR asymmetric looping of the heart and gut was normal in pitx2 mutants. These results suggest conserved roles for Pitx2 in eye and tooth development and indicate Pitx2 is not required for asymmetric looping of zebrafish visceral organs. This work establishes zebrafish pitx2 mutants as a new animal model for investigating mechanisms underlying congenital malformations in ARS and high-throughput drug screening for ARS therapeutics. Additionally, pitx2 mutants present a unique opportunity to identify new genes involved in vertebrate LR patterning. We show Nodal signaling-independent of Pitx2-controls asymmetric expression of the fatty acid elongase elovl6 in zebrafish, pointing to a potential novel pathway during LR organogenesis.

Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

Structural effect of the L16Q, K50E, and R53P mutations on homeodomain of pituitary homeobox protein 2

The transcription factor pituitary homeobox protein 2 (PITX2) is involved in genetic control of development. Mutations in PITX2, most in the homeodomain, cause the autosomal-dominant disorder Rieger syndrome. The mutants L16Q, K50E and R53P destabilize the structure and disrupt DNA-binding activity. The biological functions of these mutants have been characterized but not the structural basis behind the loss of DNA-binding activity. We performed multiple molecular dynamics simulations at 37°C to investigate the structural and dynamic effects of the 3 PITX2 homeodomain mutants. Compared with the wild type (WT), the L16Q mutant induces a kink in the α3 helix, which is stabilized by the hydrogen bond of Q21-R59. The disruption in backbone hydrogen bonds of V47-N51 and W48-R52 leads to a kink formation in the α3 helix of K50E. The R53P mutant alters the relative orientation of helices, which is apparently stabilized by the formation of new hydrogen bonds of T38-Q11, T38-Q12, T38-R2, N39-R2, L40-Q1, L40-R2, and T41-Q4. The hydrophobic core residues F8, L13, L40 and V45 change their positions in all mutants to break the hydrophobic core. Thus, changes in helical orientations and hydrophobic core cause rearrangement of the DNA-binding surface and disrupt DNA-binding activity in the mutants. The structural and molecular dynamics properties of 3 PITX2 homeodomain mutants differ from those of the WT, especially in formation of a kink in the recognition helix, change in the packing of helices and disruption of the hydrophobic core. This structural basis for the loss of DNA-binding activity for these polymorphisms may help in understanding the effect of mutations on other homeodomains with other diseases.

Vsx2 controls eye organogenesis and retinal progenitor identity via homeodomain and non-homeodomain residues required for high affinity DNA binding

The homeodomain and adjacent CVC domain in the visual system homeobox (VSX) proteins are conserved from nematodes to humans. Humans with missense mutations in these regions of VSX2 have microphthalmia, suggesting both regions are critical for function. To assess this, we generated the corresponding mutations in mouse Vsx2. The homeodomain mutant protein lacked DNA binding activity and the knock-in mutant phenocopied the null mutant, ocular retardation J. The CVC mutant protein exhibited weakened DNA binding; and, although the corresponding knock-in allele was recessive, it unexpectedly caused the strongest phenotype, as indicated by severe microphthalmia and hyperpigmentation of the neural retina. This occurred through a cryptic transcriptional feedback loop involving the transcription factors Mitf and Otx1 and the Cdk inhibitor p27^Kip1. Our data suggest that the phenotypic severity of the CVC mutant depends on the weakened DNA binding activity elicited by the CVC mutation and a previously unknown protein interaction between Vsx2 and its regulatory target Mitf. Our data also suggest that an essential function of the CVC domain is to assist the homeodomain in high-affinity DNA binding, which is required for eye organogenesis and unhindered execution of the retinal progenitor program in mammals. Finally, the genetic and phenotypic behaviors of the CVC mutation suggest it has the characteristics of a recessive neomorph, a rare type of genetic allele.

Problems with the early development of the mammalian retina can cause congenital eye defects such as microphthalmia, in which the eye is dramatically smaller and functionally compromised. Severe microphthalmia is associated with mutations in the retinal-expressed visual system homeobox 2 (Vsx2) gene, but how Vsx2 controls retinal development, and ultimately eye formation, has remained unclear. We assessed the impact of two missense mutations, discovered in humans, on Vsx2 function and eye development in mice. One mutation altered a highly conserved residue of the homeodomain, and the other altered a highly conserved residue in the CVC domain, a region of unresolved function. Both mutations impacted the DNA binding properties of the protein, although to differing extents. Likewise, both mutations caused microphthalmia and disruptions in retinal development, also to differing extents and by distinct mechanisms. Our data suggest that Vsx2 acts as a gatekeeper of the retinal gene expression program by preventing the activation of interfering or competing gene expression programs. We propose that the evolutionary stable association between the VSX-class homeodomain and CVC domain set the stage for Vsx2 or its archetype to assume a gatekeeper function for retinal development and ultimately eye organogenesis.

Discovery, optimization and validation of an optimal DNA-binding sequence for the Six1 homeodomain transcription factor

The Six1 transcription factor is a homeodomain protein involved in controlling gene expression during embryonic development. Six1 establishes gene expression profiles that enable skeletal myogenesis and nephrogenesis, among others. While several homeodomain factors have been extensively characterized with regards to their DNA-binding properties, relatively little is known of the properties of Six1. We have used the genomic binding profile of Six1 during the myogenic differentiation of myoblasts to obtain a better understanding of its preferences for recognizing certain DNA sequences. DNA sequence analyses on our genomic binding dataset, combined with biochemical characterization using binding assays, reveal that Six1 has a much broader DNA-binding sequence spectrum than had been previously determined. Moreover, using a position weight matrix optimization algorithm, we generated a highly sensitive and specific matrix that can be used to predict novel Six1-binding sites with highest accuracy. Furthermore, our results support the idea of a mode of DNA recognition by this factor where Six1 itself is sufficient for sequence discrimination, and where Six1 domains outside of its homeodomain contribute to binding site selection. Together, our results provide new light on the properties of this important transcription factor, and will enable more accurate modeling of Six1 function in bioinformatic studies.

¹H, ¹³C and ¹⁵N chemical shift assignments for the human Pitx2 homeodomain in complex with a 22-base hairpin DNA

As one of the most important eukaryotic DNA-binding motifs, the homeodomain has been identified in over one thousand proteins. Transcription factors containing homeodomains play critical roles in diverse biological processes and are often implicated in human genetic disorders. The human Pitx2 protein is a central downstream transcription factor in both the TGF-β and Wnt/β-catenin signaling pathways, and involved in the regulation of cell differentiation, cell pattern formation and tissue maintenance. To better understand how Pitx2 regulates its target genes, we are investigating the mechanisms that confer DNA-binding specificity on its homeodomain. Here, we report complete (1)H, (13)C and (15)N protein chemical shift assignments of the human Pitx2 homeodomain in complex with a non-consensus DNA sequence.

Structural and biophysical insights into the ligand-free Pitx2 homeodomain and a ring dermoid of the cornea inducing homeodomain mutant

The homeodomain-containing transcription factor Pitx2 (pituitary homeobox protein 2) is present in many developing embryonic tissues, including the heart. Its homeodomain is responsible for the recognition and binding to target DNA sequences and thus constitutes a major functional unit in the Pitx2 protein. Nuclear magnetic resonance techniques were employed to determine the solution structure of the native Pitx2 homeodomain and a R24H mutant that causes autosomal dominantly inherited ring dermoid of the cornea syndrome. The structures reveal that both isoforms possess the canonical homeodomain fold. However, the R24H mutation results in a 2-fold increase in DNA binding affinity and a 5 °C decrease in thermal stability, while changing the dynamic environment of the homeodomain only locally. When introduced into full-length Pitx2c, the mutation results in an only 25% loss of transactivation activity. Our data correlate well with clinical observations suggesting a milder deficiency for the R24H mutation compared to those of other Pitx2 homeodomain mutations.

A zebrafish model of axenfeld-rieger syndrome reveals that pitx2 regulation by retinoic acid is essential for ocular and craniofacial development

PURPOSE

The homeobox transcription factor PITX2 is a known regulator of mammalian ocular development, and human PITX2 mutations are associated with Axenfeld-Rieger syndrome (ARS). However, the treatment of patients with ARS remains mostly supportive and palliative.

METHODS

The authors used molecular genetic, pharmacologic, and embryologic techniques to study the biology of ARS in a zebrafish model that uses transgenes to mark neural crest and muscle cells in the head.

RESULTS

The authors demonstrated in vivo that pitx2 is a key downstream target of retinoic acid (RA) in craniofacial development, and this pathway is required for coordinating neural crest, mesoderm, and ocular development. pitx2a knockdown using morpholino oligonucleotides disrupts jaw and pharyngeal arch formation and recapitulates ocular characteristics of ARS, including corneal and iris stroma maldevelopment. These phenotypes could be rescued with human PITX2A mRNA, demonstrating the specificity of the knockdown and evolutionary conservation of pitx2a function. Expression of the ARS dominant negative human PITX2A K50E allele also caused ARS-like phenotypes. Similarly, inhibition of RA synthesis in the developing eye (genetic or pharmacologic) disrupted craniofacial and ocular development, and human PITX2A mRNA partially rescued these defects.

CONCLUSIONS

RA regulation of pitx2 is essential for coordinating interactions among neural crest, mesoderm, and developing eye. The marked evolutionary conservation of Pitx2 function in eye and craniofacial development makes zebrafish a potentially powerful model of ARS, amenable to in vivo experimentation and development of potential therapies.

Transcription factor heterogeneity and epiblast pluripotency

Stem cells are defined by the simultaneous possession of the seemingly incongruent properties of self-renewal and multi-lineage differentiation potential. To maintain a stem cell population, these opposing forces must be balanced. Transcription factors that function to direct pluripotent cell identity are not all equally distributed throughout the pluripotent cell population. While Oct4 levels are relatively homogeneous, other transcription factors, such as Nanog, are more heterogeneously expressed. Moreover, Oct4 positive cells fluctuate between states of high Nanog expression associated with a high probability of self-renewal and low Nanog expression associated with an increased propensity to differentiate. As embryonic stem (ES) cells transit to the more developmentally advanced epiblast stem cell (EpiSC) state, the levels of pluripotency transcription factors are modulated. Such modulations are blunted in cells that overexpress Nanog and this may underlie the resistance of Nanog-overexpressing cells to transit to an EpiSC state. Interestingly, increasing the levels of Nanog in EpiSC can facilitate reversion to the ES cell state. Together these observations suggest that Nanog lies close to the top of the hierarchy of pluripotent transcription factor regulation.

Novel NOBOX loss-of-function mutations account for 6.2% of cases in a large primary ovarian insufficiency cohort

Primary ovarian insufficiency (POI) is a disorder associated with female infertility, which affects approximately 1% of women under 40 years of age. A genetic component has been suggested as one possible cause of the majority of cases of nonsyndromic forms. Newborn Ovary Homeobox (NOBOX) is an ovary-specific gene, playing a critical role in ovary in mice, as its absence leads to sterility mimicking a POI. In this study, we sequenced NOBOX in a cohort of 178 women with idiopathic POI. Among 19 identified variations, we described one nonsense (c.907C>T/p.R303X) and four missense (c.271G>T/p.G91W, c.349C>T/p.R117W, c.1025G>C/p.S342T, and c.1048G>T/p.V350L) NOBOX heterozygous mutations in 12 patients. We reproduced each of the five mutations and tested their effects on the signaling activity in transfected cells. We demonstrated that these mutations compromised the ability of the proteins to bind to and transactivate the well-known growth differentiation factor 9 (GDF9) promoter. The pattern of our findings suggests that the genetic mechanism in humans responsible for POI in women involves haploinsufficiency rather than dominant negative gene action. The identification, characterization, and the very high 6.2% prevalence of these new mutations in POI patients suggest considering NOBOX as the first autosomal candidate gene involved in this syndrome.

Effect of CpG methylation on DNA binding protein: molecular dynamics simulations of the homeodomain PITX2 bound to the methylated DNA

A large number of studies have argued that aberrant CpG methylation is associated with some human cancers. One possible mechanism of the cancer caused by CpG methylation is the gene repression, which is a binding-inhibition of the sequence-specific transcription factors bound to specific DNA-binding sites. Exploring the effects of CpG methylation on the structure and the thermodynamic property of DNA-binding transcription factors will help to gain an insight into how CpG methylation affects the repression of gene transcription in cancer. We have performed molecular dynamics (MD) simulations and free energy calculations of the protein bound to the native or the methylated DNA, in which the solution structure of the K50-class homeodomain PITX2 bound to DNA was used as a template. The simulation results indicate that the methylated CpG located at the DNA major groove can enhance the protein-DNA interactions, and the residue side-chains near the methylated CpG pair appear to have an unusually high affinity with DNA. The structural analysis and calculated energy values demonstrate that the binding-induced structural changes were further encouraged as the CpG methylation upon the complexation. Moreover, the CpG methylation may reduce the unfavorable effect of the conformational entropy and increase the electrostatic contribution to the binding free energy of DNA-PITX2. The changes in specific binding sites and the excessive affinity between DNA and protein caused by the CpG methylation could affect the gene transcriptional activity.

The Ets factor Etv1 interacts with Tpit protein for pituitary pro-opiomelanocortin (POMC) gene transcription

Pro-opiomelanocortin (POMC) is expressed in two lineages of the pituitary, the anterior lobe corticotrophs and the intermediate lobe melanotrophs. POMC expression in these two lineages is highly dependent on the cell-restricted transcription factor Tpit. As Tpit intervenes relatively late in differentiation of those lineages, we have been searching for other transcription factors that may participate in their gene expression program. On the basis of similarity with the Tpit expression profile, we identified Ets variant gene 1 (Etv1/Er81) as a putative POMC transcription factor. Using Etv1-lacZ knockin mice, we describe preferential Etv1 expression in pituitary POMC cells and also in posterior lobe pituicytes. We further show that Etv1 enhances POMC transcription on its own and in synergy with Tpit. The Ets-binding site located within the Tpit/Pitx regulatory element is necessary for Etv1 activity in POMC-expressing AtT-20 cells but dispensable for synergy with Tpit. Etv1 and Tpit interact together in coimmunoprecipitation experiments. Furthermore, Etv1 is present at the POMC promoter, and siRNA-mediated knockdown of Etv1 in AtT-20 cells produces a significant decrease in POMC expression. Etv1 knockout pituitaries show normal POMC cell distribution and normal POMC mRNA abundance, suggesting compensation by other factors. The coordinate expression of Etv1 with POMC cell differentiation and its interaction with the highly cell-restricted Tpit factor indicate that Etv1 participates in a combinatorial code for pituitary cell-specific gene expression.

1H, 13C and 15N chemical shift assignments for the human Pitx2 homeodomain and a R24H homeodomain mutant

The homeodomain is one of the most important eukaryotic DNA-binding motifs and has been identified in over one thousand proteins. Homeodomain proteins play critical roles in diverse biological processes, including cell differentiation and cell pattern formation. The human Pitx2 homeodomain binds several different DNA sequences and is a pivotal component of both the TGF-β and Wnt/β-catenin signaling pathways. As the recognition of specific DNA sequences represents an essential biochemical function of all DNA-binding proteins, we have chosen the Pitx2 homeodomain model to investigate the mechanisms that convey biological specificity in these protein-DNA interactions. Here, we report complete chemical shift assignments of the human Pitx2 homeodomain and the R24H mutation that induces ring dermoid of the cornea syndrome.

A novel loss-of-function mutation in OTX2 in a patient with anophthalmia and isolated growth hormone deficiency

Heterozygous mutations of the gene encoding transcription factor OTX2 were recently shown to be responsible for ocular as well as pituitary abnormalities. Here, we describe a patient with unilateral anophthalmia and short stature. Endocrine evaluation of the hypothalamic-pituitary axis revealed isolated growth hormone deficiency (IGHD) with small anterior pituitary gland, invisible stalk, ectopic posterior lobe, and right anophthalmia on brain magnetic resonance imaging. DNA was analyzed for mutations in the HESX1, SOX2, and OTX2 genes. Molecular analysis yielded a novel heterozygous OTX2 mutation (c.270A>T, p.R90S) within the homeodomain. Functional analysis revealed that the mutation inhibited both the DNA binding and transactivation activities of the protein. This novel loss-of-function mutation is associated with anophthalmia and IGHD in a patient of Sephardic Jewish descent. We recommend that patients with GH deficiency and ocular malformation in whom genetic analysis for classic transcription factor genes (PROP1, POU1F1, HESX1, and LHX4) failed to identify alterations should be checked for the presence of mutations in the OTX2 gene.

Microphthalmia in Texel sheep is associated with a missense mutation in the paired-like homeodomain 3 (PITX3) gene

Microphthalmia in sheep is an autosomal recessive inherited congenital anomaly found within the Texel breed. It is characterized by extremely small or absent eyes and affected lambs are absolutely blind. For the first time, we use a genome-wide ovine SNP array for positional cloning of a Mendelian trait in sheep. Genotyping 23 cases and 23 controls using Illumina's OvineSNP50 BeadChip allowed us to localize the causative mutation for microphthalmia to a 2.4 Mb interval on sheep chromosome 22 by association and homozygosity mapping. The PITX3 gene is located within this interval and encodes a homeodomain-containing transcription factor involved in vertebrate lens formation. An abnormal development of the lens vesicle was shown to be the primary event in ovine microphthalmia. Therefore, we considered PITX3 a positional and functional candidate gene. An ovine BAC clone was sequenced, and after full-length cDNA cloning the PITX3 gene was annotated. Here we show that the ovine microphthalmia phenotype is perfectly associated with a missense mutation (c.338G>C, p.R113P) in the evolutionary conserved homeodomain of PITX3. Selection against this candidate causative mutation can now be used to eliminate microphthalmia from Texel sheep in production systems. Furthermore, the identification of a naturally occurring PITX3 mutation offers the opportunity to use the Texel as a genetically characterized large animal model for human microphthalmia.

Structural basis for DNA recognition by the human PAX3 homeodomain

The transcription regulatory protein PAX3 binds to cognate DNA sequences through two DNA-binding domains, a paired domain and a homeodomain, and has important functions during neurogenesis and myogenesis. In humans, mutations in the PAX3 gene cause Waardenburg syndrome, whereas a chromosomal translocation that generates a PAX3-FOXO1 fusion gene is associated with the development of alveolar rhabdomyosarcoma. We have determined the crystal structure of the human PAX3 homeodomain in complex with a palindromic DNA containing two inverted TAATC sequences at 1.95 A resolution. Two homeodomains bind to DNA as a symmetric dimer, inducing a 3 degrees bend in the DNA helix. The N-terminal arm of the homeodomain inserts into the minor groove and makes direct and water-mediated interactions with bases and the sugar-phosphate backbone. The recognition helix fits directly into the major groove, and an elaborate network of structurally conserved water molecules mediates the majority of protein-DNA interactions. The structure elucidates the role of serine 50 in selection of the CG sequence immediately 3' of the TAAT motif by PAX class homeodomains and provides insights into the molecular mechanisms by which certain Waardenburg syndrome-associated missense mutations could destabilize the fold of the PAX3 homeodomain whereas others could affect its interaction with DNA.

Transcriptional regulation and stabilization of left-right neuronal identity in C. elegans

At discrete points in development, transient signals are transformed into long-lasting cell fates. For example, the asymmetric identities of two Caenorhabditis elegans olfactory neurons called AWC(ON) and AWC(OFF) are specified by an embryonic signaling pathway, but maintained throughout the life of an animal. Here we show that the DNA-binding protein NSY-7 acts to convert a transient, partially differentiated state into a stable AWC(ON) identity. Expression of an AWC(ON) marker is initiated in nsy-7 loss-of-function mutants, but subsequently lost, so that most adult animals have two AWC(OFF) neurons and no AWC(ON) neurons. nsy-7 encodes a protein with distant similarity to a homeodomain. It is expressed in AWC(ON), and is an early transcriptional target of the embryonic signaling pathway that specifies AWC(ON) and AWC(OFF); its expression anticipates future AWC asymmetry. The NSY-7 protein binds a specific optimal DNA sequence that was identified through a complete biochemical survey of 8-mer DNA sequences. This sequence is present in the promoter of an AWC(OFF) marker and essential for its asymmetric expression. An 11-base-pair (bp) sequence required for AWC(OFF) expression has two activities: One region activates expression in both AWCs, and the overlapping NSY-7-binding site inhibits expression in AWC(ON). Our results suggest that NSY-7 responds to transient embryonic signaling by repressing AWC(OFF) genes in AWC(ON), thus acting as a transcriptional selector for a randomly specified neuronal identity.

Conservation of mechanisms mediating gonadotrophin-releasing hormone 1 stimulation of human luteinizing hormone beta subunit transcription

Gonadotrophin-releasing hormone (GNRH1) regulates pituitary luteinizing hormone (LH). Previous studies have delineated a mechanism for GNRH1-induced LHbeta subunit gene (Lhb) transcription, the rate-limiting step in LH production. GNRH1 induces expression of early growth response 1 (EGR1), which interacts with steroidogenic factor 1 (SF1) and paired-like homeodomain transcription factor 1 (PITX1) to regulate Lhb promoter activity. Though the cis-elements for these factors are conserved across species, regulation of human LHB transcription has not been thoroughly investigated. We therefore characterized LHB transcriptional regulation by GNRH1 using promoter-reporter analyses in LbetaT2 cells. GNRH1 stimulated LHB transcription via an extracellular signal-regulated kinase 1/2 pathway. EGR1 bound to two binding sites on the LHB promoter and this binding was increased by GNRH1. Mutation of either site or knockdown of endogenous EGR1 decreased basal and/or GNRH1-regulated promoter activity. The human LHB promoter also contains low and high affinity SF1 binding sites. Mutation of these elements or depletion of endogenous SF1 impaired basal and ligand-induced transcription. Knockdown of PITX1 or PITX2 isoforms impaired GNRH1 induction, and endogenous PITX1 bound to the candidate PITX binding site on the LHB promoter. Thus, the mechanism described for GNRH1 regulation of Lhb in other species is largely conserved for human LHB. We also uncover a previously unappreciated role for PITX2 isoforms in this process.

The six family of homeobox genes in development and cancer

The homeobox gene superfamily encodes transcription factors that act as master regulators of development through their ability to activate or repress a diverse range of downstream target genes. Numerous families exist within the homeobox gene superfamily, and are classified on the basis of conservation of their homeodomains as well as additional motifs that contribute to DNA binding and to interactions with other proteins. Members of one such family, the Six family, form a transcriptional complex with Eya and Dach proteins, and together these proteins make up part of the retinal determination network first identified in Drosophila. This network is highly conserved in both invertebrate and vertebrate species, where it influences the development of numerous organs in addition to the eye, primarily through regulation of cell proliferation, survival, migration, and invasion. Mutations in Six, Eya, and Dach genes have been identified in a variety of human genetic disorders, demonstrating their critical role in human development. In addition, aberrant expression of Six, Eya, and Dach occurs in numerous human tumors, and Six1, in particular, plays a causal role both in tumor initiation and in metastasis. Emerging evidence for the importance of Six family members and their cofactors in numerous human tumors suggests that targeting of this complex may be a novel and powerful means to inhibit both tumor growth and progression.

Asymmetric lower-limb malformations in individuals with homeobox PITX1 gene mutation

Clubfoot is one of the most common severe musculoskeletal birth defects, with a worldwide incidence of 1 in 1000 live births. In the present study, we describe a five-generation family with asymmetric right-sided predominant idiopathic clubfoot segregating as an autosomal-dominant condition with incomplete penetrance. Other lower-limb malformations, including patellar hypoplasia, oblique talus, tibial hemimelia, developmental hip dysplasia, and preaxial polydactyly, were also present in some family members. Genome-wide linkage analysis with Affymetrix GeneChip Mapping 10K mapping data from 13 members of this family revealed a multipoint LOD(max) of 3.31 on chromosome 5q31. A single missense mutation (c.388G-->A) was identified in PITX1, a bicoid-related homeodomain transcription factor critical for hindlimb development, and segregated with disease in this family. The PITX1 E130K mutation is located in the highly conserved homeodomain and reduces the ability of PITX1 to transactivate a luciferase reporter. The PITX1 E130K mutation also suppresses wild-type PITX1 activity in a dose-dependent manner, suggesting dominant-negative effects on transcription. The propensity for right-sided involvement in tibial hemimelia and clubfoot suggests that PITX1, or pathways involving PITX1, may be involved in their etiology. Implication of a gene involved in early limb development in clubfoot pathogenesis also suggests additional pathways for future investigations of idiopathic clubfoot etiology in humans.

FIZ1 is part of the regulatory protein complex on active photoreceptor-specific gene promoters in vivo

Background

FIZ1 (Flt-3 Interacting Zinc-finger) is a broadly expressed protein of unknown function. We reported previously that in the mammalian retina, FIZ1 interacts with NRL (Neural-Retina Leucine-zipper), an essential transcriptional activator of rod photoreceptor-specific genes. The concentration of FIZ1 in the retina increases during photoreceptor terminal maturation, when two key transcription factors NRL and CRX (Cone-Rod Homeobox) become detectable on the promoters of photoreceptor-specific genes (i.e. Rhodopsin, Pde6b). To determine if FIZ1 is involved in regulating CRX-mediated transcriptional activation, we examined FIZ1 subcellular location in mouse neural retina, its ability to interact with CRX, and its association with CRX/NRL target genes.

Results

FIZ1 is present in the nucleus of adult photoreceptors as well as other retinal neurons as shown by transmission electron microscopy with nano-gold labeling. FIZ1 and CRX were co-precipitated from retinal nuclear extracts with antibodies to either protein. Chromatin immunoprecipitation (ChIP) assays revealed that FIZ1 is part of the protein complex on several rod and cone gene promoters, within photoreceptor cells of the mouse retina. FIZ1 complexes with CRX or NRL on known NRL- and CRX-responsive elements, as shown by electrophoretic mobility shift assays with FIZ1 antibody. FIZ1 can directly bind to CRX, as demonstrated using yeast two-hybrid and GST pull-down assays. Co-transfection assays demonstrated that FIZ1 increases CRX-mediated activation of Opsin test promoters. Quantitative ChIP analysis revealed an increased association of FIZ1 with the Rhodopsin promoter in adult (P-25) neural retina versus immature (P-3) neural retina. The quantity of transcriptionally active RNA Polymerase-II within the Rhodopsin gene (Rho) was significantly increased in the adult neural retina, compared to the immature retina.

Conclusion

FIZ1 directly interacts with CRX to enhance CRX's transactivation activity for target genes. Developmentally, in neural retina tissue, the increased association of FIZ1 with CRX target genes corresponds to an increased association of transcriptionally active Pol-II within the Rho gene. Together with previous findings, our results suggest that FIZ1 may act as a transcriptional co-regulator of photoreceptor-specific genes, recruited by at least two photoreceptor-specific transcription factors, CRX and NRL. Further studies are underway to elucidate the exact role of FIZ1 in photoreceptor gene expression, development and maintenance.

Paired-like homeodomain transcription factors 1 and 2 regulate follicle-stimulating hormone beta-subunit transcription through a conserved cis-element

Paired-like homeodomain transcription factors (PITX) regulate the activity of pituitary hormone-encoding genes. Here, we examined mechanisms through which the family of PITX proteins control murine FSH beta-subunit (Fshb) transcription. We observed that endogenous PITX1 and PITX2 isoforms from murine LbetaT2 gonadotrope cells could bind a highly conserved proximal cis-element. Transfection of PITX1 or PITX2C in heterologous cells stimulated both murine and human Fshb/FSHB promoter-reporter activities, and in both cases, mutation of the critical cis-element abrogated these effects. In homologous LbetaT2 cells, the same mutation decreased basal reporter activity and greatly reduced activin A-stimulated transcription from murine and human promoter-reporters. Transfecting dominant-negative forms of PITX1 or PITX2C or knocking down PITX1 or -2 expression by RNA interference in LbetaT2 cells inhibited murine Fshb transcription, confirming roles for endogenous PITX proteins. Both PITX1 and PITX2C interacted with Smad3 (an effector of the activin signaling cascade in these cells) in coprecipitation experiments, and the PITX binding site mutation greatly inhibited Smad2/3/4-stimulated Fshb transcription. In summary, both PITX1 and PITX2C regulate murine and human Fshb/FSHB transcription through a conserved cis-element in the proximal promoter. Furthermore, the data indicate both common and distinct mechanisms of PITX1 and PITX2C action.

Crystal structure and DNA binding of the homeodomain of the stem cell transcription factor Nanog

The transcription factor Nanog is an upstream regulator in early mammalian development and a key determinant of pluripotency in embryonic stem cells. Nanog binds to promoter elements of hundreds of target genes and regulates their expression by an as yet unknown mechanism. Here, we report the crystal structure of the murine Nanog homeodomain (HD) and analysis of its interaction with a DNA element derived from the Tcf3 promoter. Two Nanog amino acid pairs, unique among HD sequences, appear to affect the mechanism of nonspecific DNA recognition as well as maintain the integrity of the structural scaffold. To assess selective DNA recognition by Nanog, we performed electrophoretic mobility shift assays using a panel of modified DNA binding sites and found that Nanog HD preferentially binds the TAAT(G/T)(G/T) motif. A series of rational mutagenesis experiments probing the role of six variant residues of Nanog on its DNA binding function establish their role in affecting binding affinity but not binding specificity. Together, the structural and functional evidence establish Nanog as a distant member of a Q50-type HD despite having considerable variation at the sequence level.

Molecular and developmental mechanisms of anterior segment dysgenesis

Anterior segment dysgenesis (ASD) is a failure of the normal development of the tissues of the anterior segment of the eye. It leads to anomalies in the structure of the mature anterior segment, associated with an increased risk of glaucoma and corneal opacity. Several different gene mutations have been identified underlying these anomalies with the majority of ASD genes encoding transcriptional regulators. In this review, the role of the ASD genes, PITX2 and FOXC1, is considered in relation to the embryology of the anterior segment, the biochemical function of these proteins, and their role in development and disease aetiology. The emerging view is that these genes act in concert to specify a population of mesenchymal progenitor cells, mainly of neural crest origin, as they migrate anteriorly around the embryonic optic cup. These same genes then regulate mesenchymal cell differentiation to give rise to distinct anterior segment tissues. Development appears critically sensitive to gene dosage, and variation in the normal level of transcription factor activity causes a range of anterior segment anomalies. Interplay between PITX2 and FOXC1 in the development of different anterior segment tissues may partly explain the phenotypic variability and the genetic heterogeneity characteristic of ASD.

The structures that underlie normal reproductive function

The mechanisms and physiology of reproductive function have fascinated scientists throughout time. Recent cellular and molecular level structural studies have provided unprecedented insights into reproductive systems and signaling networks. This 'cutting edge' editorial provides a recent example in each of these areas, namely, the anatomical integrity of the follicle, the molecular structure of activin with its binding partners and the molecular regulation of inhibin. These three examples of structure informing function help explain reproductive health and may provide solutions to reproductive disease.

Molecular dissection reveals decreased activity and not dominant negative effect in human OTX2 mutants

The paired-type homeodomain transcription factor Otx2 is essential for forebrain and eye development. Severe ocular malformations in humans have recently been associated with heterozygous OTX2 mutations. To document the molecular defects in human mutants, Otx2 structural characterization was carried out. A collection of deletion and point mutants was created to perform transactivation, DNA binding, and subcellular localization analyses. Transactivation was ascribed to both N- and C-termini of the protein, and DNA binding to the minimal homeodomain, where critical amino acid residues were identified. Acute nuclear localization appeared controlled by a nuclear localization sequence located within the homeodomain which acts in conjunction with a novel nuclear retention domain that we unraveled located in the central part of the protein. This region, which is poorly conserved among Otx proteins, was also endowed with dominant negative activity suggesting that it might confer unique properties to Otx2. Molecular diagnostic of human mutant OTX2 proteins discriminates hypomorphic and loss of function mutations from other mutations that may not be relevant to ocular pathology.

The solution structure of the native K50 Bicoid homeodomain bound to the consensus TAATCC DNA-binding site

The solution structure of the homeodomain of the Drosophila morphogenic protein Bicoid (Bcd) complexed with a TAATCC DNA site is described. Bicoid is the only known protein that uses a homeodomain to regulate translation, as well as transcription, by binding to both RNA and DNA during early Drosophila development; in addition, the Bcd homeodomain can recognize an array of different DNA sites. The dual functionality and broad recognition capabilities signify that the Bcd homeodomain may possess unique structural/dynamic properties. Bicoid is the founding member of the K50 class of homeodomain proteins, containing a lysine residue at the critical 50th position (K50) of the homeodomain sequence, a residue required for DNA and RNA recognition; Bcd also has an arginine residue at the 54th position (R54), which is essential for RNA recognition. Bcd is the only known homeodomain with the K50/R54 combination of residues. The Bcd structure indicates that this homeodomain conforms to the conserved topology of the homeodomain motif, but exhibits a significant variation from other homeodomain structures at the end of helix 1. A key result is the observation that the side-chains of the DNA-contacting residues K50, N51 and R54 all show strong signs of flexibility in the protein-DNA interface. This finding is supportive of the adaptive-recognition theory of protein-DNA interactions.