A genetic model for interaction of the homeodomain recognition helix with DNA

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.

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.

Homeobox Gene Duplication and Divergence in Arachnids

Homeobox genes are key toolkit genes that regulate the development of metazoans and changes in their regulation and copy number have contributed to the evolution of phenotypic diversity. We recently identified a whole genome duplication (WGD) event that occurred in an ancestor of spiders and scorpions (Arachnopulmonata), and that many homeobox genes, including two Hox clusters, appear to have been retained in arachnopulmonates. To better understand the consequences of this ancient WGD and the evolution of arachnid homeobox genes, we have characterized and compared the homeobox repertoires in a range of arachnids. We found that many families and clusters of these genes are duplicated in all studied arachnopulmonates (Parasteatoda tepidariorum, Pholcus phalangioides, Centruroides sculpturatus, and Mesobuthus martensii) compared with nonarachnopulmonate arachnids (Phalangium opilio, Neobisium carcinoides, Hesperochernes sp., and Ixodes scapularis). To assess divergence in the roles of homeobox ohnologs, we analyzed the expression of P. tepidariorum homeobox genes during embryogenesis and found pervasive changes in the level and timing of their expression. Furthermore, we compared the spatial expression of a subset of P. tepidariorum ohnologs with their single copy orthologs in P. opilio embryos. We found evidence for likely subfunctionlization and neofunctionalization of these genes in the spider. Overall our results show a high level of retention of homeobox genes in spiders and scorpions post-WGD, which is likely to have made a major contribution to their developmental evolution and diversification through pervasive subfunctionlization and neofunctionalization, and paralleling the outcomes of WGD in vertebrates.

Brain-specific homeobox Bsx specifies identity of pineal gland between serially homologous photoreceptive organs in zebrafish

The pineal gland functioning as a photoreceptive organ in non-mammalian species is a serial homolog of the retina. Here we found that Brain-specific homeobox (Bsx) is a key regulator conferring individuality on the pineal gland between the two serially homologous photoreceptive organs in zebrafish. Bsx knock-down impaired the pineal development with reduced expression of exorh, the pineal-specific gene responsible for the photoreception, whereas it induced ectopic expression of rho, a retina-specific gene, in the pineal gland. Bsx remarkably transactivated the exorh promoter in combination with Otx5, but not with Crx, through its binding to distinct subtypes of PIRE, a DNA cis-element driving Crx/Otx-dependent pineal-specific gene expression. These results demonstrate that the identity of pineal photoreceptive neurons is determined by the combinatorial code of Bsx and Otx5, the former confers the pineal specificity at the tissue level and the latter determines the photoreceptor specificity at the cellular level.

Phenotypic Nonspecificity as the Result of Limited Specificity of Transcription Factor Function

Drosophila transcription factor (TF) function is phenotypically nonspecific. Phenotypic nonspecificity is defined as one phenotype being induced or rescued by multiple TFs. To explain this unexpected result, a hypothetical world of limited specificity is explored where all TFs have unique random distributions along the genome due to low information content of DNA sequence recognition and somewhat promiscuous cooperative interactions with other TFs. Transcription is an emergent property of these two conditions. From this model, explicit predictions are made. First, many more cases of TF nonspecificity are expected when examined. Second, the genetic analysis of regulatory sequences should uncover cis-element bypass and, third, genetic analysis of TF function should generally uncover differential pleiotropy. In addition, limited specificity provides evolutionary opportunity and explains the inefficiency of expression analysis in identifying genes required for biological processes.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

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.

The combination of sequence-specific and nonspecific DNA-binding modes of transcription factor SATB1

Transcription factor SATB1 (special AT-rich sequence binding protein 1) contains multiple DNA-binding domains (DBDs), i.e. two CUT-domain repeats (CUTr1 and CUTr2 from the N-terminus) and a homeodomain, and binds to the matrix attachment region (MAR) of DNA. Although CUTr1 and the homeodomain, but not CUTr2, are known to contribute to DNA binding, different research groups have not reached a consensus on which DBD is responsible for recognition of the target sequence in MAR, 5'-TAATA-3'. Here, we used isothermal titration calorimetry to demonstrate that CUTr1 has binding specificity to this motif, whereas the homeodomain shows affinity for a variety of DNAs without specificity. In line with nonspecific DNA-binding properties of the homeodomain, a mutation of the invariant Asn at position 51 of the homeodomain (typically in contact with the A base in a sequence-specific binding mode) did not affect the binding affinity significantly. The NMR analyses and computational modeling of the homeodomain, however, revealed the tertiary structure and DNA-binding mode that are typical of homeodomains capable of sequence-specific binding. We believe that the lack of highly conserved basic residues in the helix relevant to the base recognition loosens its fitting into the DNA groove and impairs the specific binding. The two DBDs, when fused in tandem, showed strong binding to DNA containing the 5'-TAATA-3' motif with an affinity constant >10(8) M(-1) and retained nonspecific binding activity. The combination of the sequence-specific and nonspecific DNA-binding modes of SATB1 should be advantageous in a search for target loci during transcriptional regulation.

Roles of microRNAs in cancers and development

MicroRNAs (miRNAs) are small noncoding RNA molecules thought to play an important role in regulating gene expression. Although knowledge of the biological functions of most miRNAs is still limited, miRNAs are thought to regulate the gene expression in various diseases and embryo development. In this chapter, the roles of miRNAs in human cancers are first examined from the viewpoint of up- and downregulation. Oncogenic miRNAs are involved in the overexpression/upregulation of cancers, whereas suppressive miRNAs are involved in the underexpression/downregulation of cancers. Statistical analysis of the positional nucleotide occurrence features of miRNAs revealed differences between the positional nucleotide occurrences of oncogenic and suppressive miRNAs. A miRNA gene-silencing score was then defined on the basis of the higher and lower levels of the statistical significances of positional nucleotides. Since the miRNA scores were closely related to miRNA frequencies, a method using the scores and nucleotide frequencies to distinguish whether a new miRNA is oncogenic or suppressive is proposed. This chapter also describes the roles of miRNAs in development. As miRNAs can act as cis-regulatory elements in the early embryonic development of Drosophila melanogaster, it is proposed that they mediate signal transduction between genes.

Study of base pair mutations in proline-rich homeodomain (PRH)-DNA complexes using molecular dynamics

Proline-rich homeodomain (PRH) is a regulatory protein controlling transcription and gene expression processes by binding to the specific sequence of DNA, especially to the sequence 5'-TAATNN-3'. The impact of base pair mutations on the binding between the PRH protein and DNA is investigated using molecular dynamics and free energy simulations to identify DNA sequences that form stable complexes with PRH. Three 20-ns molecular dynamics simulations (PRH-TAATTG, PRH-TAATTA and PRH-TAATGG complexes) in explicit solvent water were performed to investigate three complexes structurally. Structural analysis shows that the native TAATTG sequence forms a complex that is more stable than complexes with base pair mutations. It is also observed that upon mutation, the number and occupancy of the direct and water-mediated hydrogen bonds decrease. Free energy calculations performed with the thermodynamic integration method predict relative binding free energies of 0.64 and 2 kcal/mol for GC to AT and TA to GC mutations, respectively, suggesting that among the three DNA sequences, the PRH-TAATTG complex is more stable than the two mutated complexes. In addition, it is demonstrated that the stability of the PRH-TAATTA complex is greater than that of the PRH-TAATGG complex.

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.

Occludin is a direct target of thyroid transcription factor-1 (TTF-1/NKX2-1)

The thyroid transcription factor 1 gene (TTF-1 or NKX2-1) is essential to lung development; however, it is also a critical factor in lung cancer. TTF-1 is amplified in lung cancers, suggesting that it is a gain-of-function lung oncogene. Conversely, TTF-1 counters epithelial to mesenchymal transition in cell-based studies and inhibits progression of primary lung adenocarcinomas to metastases in an animal model of lung adenocarcinomas. The unifying theory regarding TTF-1 is that it exhibits both pro-oncogenic and anti-metastatic function depending on the cellular context. Occludin is the first discovered constituent of the epithelial tight junction; in recent years, a functional role of occludin as a tumor suppressor has begun to emerge. Here, we demonstrate that TTF-1 transactivated the expression of the epithelial tight junction molecules occludin (OCLN) and claudin-1 (CLDN1). We show that transcriptional activation occurred through a direct interaction of TTF-1 with the OCLN and CLDN1 promoters. Furthermore, in cells that lack TTF-1, exogenous TTF-1 expression dampened the inhibitory effect of TGF-β on occludin and claudin-1 content. Using cells derived from a genetically engineered mouse model of lung adenocarcinomas, we observed that silenced TTF-1 expression down-regulated occludin, which we supported with additional siRNA experiments. Finally, TTF-1 knockdown conferred human lung cancer cells resistance to anoikis, and expression of occludin restored cellular sensitivity to anoikis. Overexpression of occludin impeded migration and induced anoikis in lung carcinoma cells. Collectively, these data suggest that TTF-1 transcriptionally regulates occludin, which represents another avenue of TTF-1-mediated metastasis suppression.

Exploring the DNA-recognition potential of homeodomains

The recognition potential of most families of DNA-binding domains (DBDs) remains relatively unexplored. Homeodomains (HDs), like many other families of DBDs, display limited diversity in their preferred recognition sequences. To explore the recognition potential of HDs, we utilized a bacterial selection system to isolate HD variants, from a randomized library, that are compatible with each of the 64 possible 3' triplet sites (i.e., TAANNN). The majority of these selections yielded sets of HDs with overrepresented residues at specific recognition positions, implying the selection of specific binders. The DNA-binding specificity of 151 representative HD variants was subsequently characterized, identifying HDs that preferentially recognize 44 of these target sites. Many of these variants contain novel combinations of specificity determinants that are uncommon or absent in extant HDs. These novel determinants, when grafted into different HD backbones, produce a corresponding alteration in specificity. This information was used to create more explicit HD recognition models, which can inform the prediction of transcriptional regulatory networks for extant HDs or the engineering of HDs with novel DNA-recognition potential. The diversity of recovered HD recognition sequences raises important questions about the fitness barrier that restricts the evolution of alternate recognition modalities in natural systems.

HOXB4 can enhance the differentiation of embryonic stem cells by modulating the hematopoietic niche

Hematopoietic differentiation of embryonic stem cells (ESCs) in vitro has been used as a model to study early hematopoietic development, and it is well documented that hematopoietic differentiation can be enhanced by overexpression of HOXB4. HOXB4 is expressed in hematopoietic progenitor cells (HPCs) where it promotes self-renewal, but it is also expressed in the primitive streak of the gastrulating embryo. This led us to hypothesize that HOXB4 might modulate gene expression in prehematopoietic mesoderm and that this property might contribute to its prohematopoietic effect in differentiating ESCs. To test our hypothesis, we developed a conditionally activated HOXB4 expression system using the mutant estrogen receptor (ER(T2)) and showed that a pulse of HOXB4 prior to HPC emergence in differentiating ESCs led to an increase in hematopoietic differentiation. Expression profiling revealed an increase in the expression of genes associated with paraxial mesoderm that gives rise to the hematopoietic niche. Therefore, we considered that HOXB4 might modulate the formation of the hematopoietic niche as well as the production of hematopoietic cells per se. Cell mixing experiments supported this hypothesis demonstrating that HOXB4 activation can generate a paracrine as well as a cell autonomous effect on hematopoietic differentiation. We provide evidence to demonstrate that this activity is partly mediated by the secreted protein FRZB.

CRX ChIP-seq reveals the cis-regulatory architecture of mouse photoreceptors

Approximately 98% of mammalian DNA is noncoding, yet we understand relatively little about the function of this enigmatic portion of the genome. The cis-regulatory elements that control gene expression reside in noncoding regions and can be identified by mapping the binding sites of tissue-specific transcription factors. Cone-rod homeobox (CRX) is a key transcription factor in photoreceptor differentiation and survival, but its in vivo targets are largely unknown. Here, we used chromatin immunoprecipitation with massively parallel sequencing (ChIP-seq) on CRX to identify thousands of cis-regulatory regions around photoreceptor genes in adult mouse retina. CRX directly regulates downstream photoreceptor transcription factors and their target genes via a network of spatially distributed regulatory elements around each locus. CRX-bound regions act in a synergistic fashion to activate transcription and contain multiple CRX binding sites which interact in a spacing- and orientation-dependent manner to fine-tune transcript levels. CRX ChIP-seq was also performed on Nrl(-/-) retinas, which represent an enriched source of cone photoreceptors. Comparison with the wild-type ChIP-seq data set identified numerous rod- and cone-specific CRX-bound regions as well as many shared elements. Thus, CRX combinatorially orchestrates the transcriptional networks of both rods and cones by coordinating the expression of photoreceptor genes including most retinal disease genes. In addition, this study pinpoints thousands of noncoding regions of relevance to both Mendelian and complex retinal disease.

A cis-regulatory signature for chordate anterior neuroectodermal genes

One of the striking findings of comparative developmental genetics was that expression patterns of core transcription factors are extraordinarily conserved in bilaterians. However, it remains unclear whether cis-regulatory elements of their target genes also exhibit common signatures associated with conserved embryonic fields. To address this question, we focused on genes that are active in the anterior neuroectoderm and non-neural ectoderm of the ascidian Ciona intestinalis. Following the dissection of a prototypic anterior placodal enhancer, we searched all genomic conserved non-coding elements for duplicated motifs around genes showing anterior neuroectodermal expression. Strikingly, we identified an over-represented pentamer motif corresponding to the binding site of the homeodomain protein OTX, which plays a pivotal role in the anterior development of all bilaterian species. Using an in vivo reporter gene assay, we observed that 10 of 23 candidate cis-regulatory elements containing duplicated OTX motifs are active in the anterior neuroectoderm, thus showing that this cis-regulatory signature is predictive of neuroectodermal enhancers. These results show that a common cis-regulatory signature corresponding to K50-Paired homeodomain transcription factors is found in non-coding sequences flanking anterior neuroectodermal genes in chordate embryos. Thus, field-specific selector genes impose architectural constraints in the form of combinations of short tags on their target enhancers. This could account for the strong evolutionary conservation of the regulatory elements controlling field-specific selector genes responsible for body plan formation.

Regional identity in embryos is defined by a few specific transcription factors that activate a large number of target genes through binding to common tags in regulatory sequences. In chordates it is unclear if such tags can be identified in the cis-regulatory regions of regionally expressed genes. To address this question we focused on the anterior nervous system where Otx codes for a transcription factor that triggers expression of many other head-specific genes. We analyzed an element that is active in the region bordering the anterior nervous system in the marine invertebrate Ciona intestinalis. We found that the crucial binding sites have to be duplicated and close enough. One of the pairs is bound by OTX. We showed that anterior nervous system genes are often flanked by duplicated OTX binding sites. We confirmed by transgenic assays that about half of these genomic sequences are active and drive expression anteriorly. This study unravels a simple regulatory logic in the anterior enhancers. It indicates that although there are major changes in the organization of the binding sites at short evolutionary range, conserved expression patterns are partly generated by a duplicated organization of conserved binding sites for region-specific transcription factors.

cis-Regulatory sequences driving the expression of the Hbox12 homeobox-containing gene in the presumptive aboral ectoderm territory of the Paracentrotus lividus sea urchin embryo

Embryonic development is coordinated by networks of evolutionary conserved regulatory genes encoding transcription factors and components of cell signalling pathways. In the sea urchin embryo, a number of genes encoding transcription factors display territorial restricted expression. Among these, the zygotic Hbox12 homeobox gene is transiently transcribed in a limited number of cells of the animal-lateral half of the early Paracentrotus lividus embryo, whose descendants will constitute part of the ectoderm territory. To obtain insights on the regulation of Hbox12 expression, we have explored the cis-regulatory apparatus of the gene. In this paper, we show that the intergenic region of the tandem Hbox12 repeats drives GFP expression in the presumptive aboral ectoderm and that a 234 bp fragment, defined aboral ectoderm (AE) module, accounts for the restricted expression of the transgene. Within this module, a consensus sequence for a Sox factor and the binding of the Otx activator are both required for correct Hbox12 gene expression. Spatial restriction to the aboral ectoderm is achieved by a combination of different repressive sequence elements. Negative sequence elements necessary for repression in the endomesoderm map within the most upstream 60 bp region and nearby the Sox binding site. Strikingly, a Myb-like consensus is necessary for repression in the oral ectoderm, while down-regulation at the gastrula stage depends on a GA-rich region. These results suggest a role for Hbox12 in aboral ectoderm specification and represent our first attempt in the identification of the gene regulatory circuits involved in this process.

Functional analysis of human mutations in homeodomain transcription factor PITX3

BACKGROUND

The homeodomain-containing transcription factor PITX3 was shown to be essential for normal eye development in vertebrates. Human patients with point mutations in PITX3 demonstrate congenital cataracts along with anterior segment defects in some cases when one allele is affected and microphthalmia with brain malformations when both copies are mutated. The functional consequences of these human mutations remain unknown.

RESULTS

We studied the PITX3 mutant proteins S13N and G219fs to determine the type and severity of functional defects. Our results demonstrate alterations in DNA-binding profiles and/or transactivation activities and suggest a partial loss-of-function in both mutants with the G219fs form being more severely affected. No anomalies in cellular distribution and no dominant-negative effects were discovered for these mutants. Interestingly, the impairment of the G219fs activity varied between different ocular cell lines.

CONCLUSION

The G219fs mutation was found in multiple families affected with congenital cataracts along with anterior segment malformations in many members. Our data suggest that the presence/severity of anterior segment defects in families affected with G219fs may be determined by secondary factors that are expressed in the developing anterior segment structures and may modify the effect(s) of this mutation. The S13N mutant showed only minor alteration of transactivation ability and DNA binding pattern and may represent a rare polymorphism in the PITX3 gene. A possible contribution of this mutation to human disease needs to be further investigated.

SpGataE, a Strongylocentrotus purpuratus ortholog of mammalian Gata4/5/6: protein expression, interaction with putative target gene spec2a, and identification of friend of Gata factor SpFog1

In the sea urchin Strongylocentrotus purpuatus, SpGataE, an ortholog of the vertebrate zinc-finger transcription factors Gata4/5/6, occupies a key position in the gene regulatory network for endomesoderm specification. We have posited that in addition to regulating gene activity required for endomesoderm specification, SpGataE also represses the expression of the aboral ectoderm-specific spec2a gene in endomesoderm territories. Although the expression pattern of spgatae and its role in endomesoderm specification have been described in considerable detail, little is known about SpGataE protein accumulation and its interactions with target genes and coregulatory factors. Our purpose here was to gain further insight into the mechanisms by which SpGataE functions as a transcriptional regulator. To achieve this, we generated an anti-SpGataE antibody to determine the spatiotemporal expression pattern of SpGataE protein and establish whether it plays a role in repressing spec2a by binding to gata cis-regulatory elements within the endogenous spec2a enhancer. Because Gata proteins often associate with friend of Gata (Fog) coregulators, we identified an S. purpuratus fog ortholog, spfog1, and showed that SpGataE and SpFog1 physically interacted. Spfog1 transcripts were maximal by early blastula stage but continued thereafter to be expressed at low levels. Knockdown of spfog1 using antisense morpholino oligonucleotides did not produce notable effects on endomesoderm specification or spec2a enhancer activity, suggesting that SpGataE exerts these functions independently of SpFog1. In addition to providing new information on Gata and Fog proteins in sea urchins, the anti-SpGataE antibody developed here should be a useful reagent for future analysis of SpGataE function.

Conserved interactions of a compact highly active enhancer/promoter upstream of the rhodopsin kinase (GRK1) gene

Rhodopsin kinase (RK) is a conserved component of the light adaptation and recovery pathways shared among rod and cone photoreceptors of a variety of species. To gain insight into transcriptional mechanisms driving RK and potentially other genes of similar spatial profile, the components and the interactions of the highly compact enhancer/promoter region (E/P) upstream of the human RK gene were examined. Cross-species comparison outlined an active 49-bp widely shared E/P core as the major site of conservation in the entire 5' flanking sequence. The area consisted of a bicoid-type homeodomain recognition cassette and a unique T-rich module interacting with TATA-binding proteins. Homeodomain interactions involved primarily Crx and secondarily Otx2. Both strongly stimulated the E/P. In the absence of Crx, persistent E/P activity shifted from the outer retina to the inner to follow the Otx2 pattern. The spatial patterns were largely unaffected by the absence of rod transcription factors, Nrl and Nr2e3, and the RK transcriptional activity preceded the surge in rod-specific transcription. Conserved bicoid homeodomain factors thus appear to be the key factors governing localization of RK E/P activity in retina and photoreceptors.

In vitro analysis of DNA-protein interactions by proximity ligation

Protein-binding DNA sequence elements encode a variety of regulated functions of genomes. Information about such elements is currently in a state of rapid growth, but improved methods are required to characterize the sequence specificity of DNA-binding proteins. We have established an in vitro method for specific and sensitive solution-phase analysis of interactions between proteins and nucleic acids in nuclear extracts, based on the proximity ligation assay. The reagent consumption is very low, and the excellent sensitivity of the assay enables analysis of as few as 1-10 cells. We show that our results are highly reproducible, quantitative, and in good agreement with both EMSA and predictions obtained by using a motif finding software. This assay can be a valuable tool to characterize in-depth the sequence specificity of DNA-binding proteins and to evaluate effects of polymorphisms in known transcription factor binding sites.

A 3' cis-regulatory region controls wingless expression in the Drosophila eye and leg primordia

The precise regulation of wingless (wg) expression in the Drosophila eye disc is key to control the anteroposterior and dorsoventral patterning of this disc. Here, we identify an eye disc-specific wg cis-regulatory element that functions as a regulatory rheostat. Pannier (Pnr), a transcription factor previously proposed to act as an upstream activator of wg, is sufficient to activate the eye disc enhancer but required for wg expression only in the peripodial epithelium of the disc. We propose that this regulation of wg by Pnr appeared associated to the development of the peripodial epithelium in higher dipterans and was added to an existing mechanism regulating the deployment of wingless in the dorsal region of the eye primordium. In addition, our analysis identifies a separate ventral disc enhancer that lies adjacent to the eye-specific one, and thus altogether, they define a 1-kb genomic region where disc-specific enhancers of the wg gene are located.

Nkx3.1 binds and negatively regulates the transcriptional activity of Sp-family members in prostate-derived cells

Nkx3.1 is a homeodomain-containing transcription factor that is expressed early in the development of the prostate gland and is believed to play an important role in the differentiation of prostatic epithelia. Loss of Nkx3.1 protein expression is often an early event in prostate tumorigenesis, and the abundance of Nkx3.1-negative epithelial cells increases with disease progression. In a number of systems, homeodomain proteins collaborate with zinc-finger-containing transcription factors to bind and regulate target genes. In the present paper, we report that Nkx3.1 collaborates with Sp-family members in the regulation of PSA (prostate-specific antigen) in prostate-derived cells. Nkx3.1 forms protein complexes with Sp proteins that are dependent on their respective DNA-binding domains and an N-terminal segment of Nkx3.1, and Nkx3.1 negatively regulates Sp-mediated transcription via Trichostatin A-sensitive and -insensitive mechanisms. A distal 1000 bp portion of the PSA promoter is required for transrepression by Nkx3.1, although Nkx3.1 DNA-binding activity is itself not required. We conclude that Nkx3.1 negatively regulates Sp-mediated transcription via the tethering of histone deacetylases and/or by inhibiting the association of Sp proteins with co-activators.

Studying statistical properties of regulatory DNA sequences, and their use in predicting regulatory regions in the eukaryotic genomes

There are no well-known properties in regulatory DNA analogous to those in coding sequences; their spatial location is not regular, the consensus regulatory elements are often degenerate and there are no understandable rules governing their evolution. This makes it difficult to recognize regulatory regions within genome. We review developments in the statistical characterization of regulatory regions and methods of their recognition in eukaryotic genomes.

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.

CHX10 targets a subset of photoreceptor genes

The homeobox gene CHX10 is required for retinal progenitor cell proliferation early in retinogenesis and subsequently for bipolar neuron differentiation. To clarify the molecular mechanisms employed by CHX10 we sought to identify its target genes. In a yeast one-hybrid assay Chx10 interacted with the Ret1 site of the photoreceptor-specific gene Rhodopsin. Gel shift assays using in vitro translated protein confirmed that CHX10 binds to Ret1, but not to the similar Rhodopsin sites Ret4 and BAT-1. Using retinal nuclear lysates, we observed interactions between Chx10 and additional photoreceptor-specific elements including the PCE-1 (Rod arrestin/S-antigen) and the Cone opsin locus control region (Red/green cone opsin). However, chromatin immunoprecipitation assays revealed that in vivo, Chx10 bound sites upstream of the Rod arrestin and Interphotoreceptor retinoid-binding protein genes but not Rhodopsin or Cone opsin. Thus, in a chromatin context, Chx10 associates with a specific subset of elements that it binds with comparable apparent affinity in vitro. Our data suggest that CHX10 may target these motifs to inhibit rod photoreceptor gene expression in bipolar cells.

Exo-Taq-Based Detection of DNA-Binding Protein for Homogeneous and Microarray Format

The study of DNA-protein interactions is of great importance to understand basic cellular processes such as transcription, replication and recombination. In this research, we developed a novel detection system for DNA-binding proteins (DBPs) involving the exonuclease (Exo) III and Taq DNA polymerase reactions. The system consists of three steps, as follows: the target DBP in the sample solution is incubated with probe DNA, and the probe is digested with Exo III and then extended with Taq using fluorescent dye-labeled dUTP as a substrate. The DBP protects the probe from digestion by Exo III. Therefore, only the DBP-bound probe allows the following extension. We examined this system using the lambda phage Cro repressor in a homogeneous format. The fluorescence image after gel electrophoresis showed a specific band. We also found that this system could be applied to the rapid and efficient detection of DBPs in stem and loop ds-DNA array formats. These results suggest that our method is useful as a new tool for analyzing DNA-protein interactions.

Bicoid cooperative DNA binding is critical for embryonic patterning in Drosophila

Cooperative interactions by DNA-binding proteins have been implicated in cell-fate decisions in a variety of organisms. To date, however, there are few examples in which the importance of such interactions has been explicitly tested in vivo. Here, we tested the importance of cooperative DNA binding by the Bicoid protein in establishing a pattern along the anterior-posterior axis of the early Drosophila embryo. We found that bicoid mutants specifically defective in cooperative DNA binding fail to direct proper development of the head and thorax, leading to embryonic lethality. The mutants did not faithfully stimulate transcription of downstream target genes such as hunchback (hb), giant, and Krüppel. Quantitative analysis of gene expression in vivo indicated that bcd cooperativity mutants were unable to accurately direct the extent to which hb is expressed along the anterior-posterior axis and displayed a reduced ability to generate sharp on/off transitions for hb gene expression. These failures in precise transcriptional control demonstrate the importance of cooperative DNA binding for embryonic patterning in vivo.

Fast and systematic genome-wide discovery of conserved regulatory elements using a non-alignment based approach

We describe a powerful new approach for discovering globally conserved regulatory elements between two genomes. The method is fast, simple and comprehensive, without requiring alignments. Its application to pairs of yeasts, worms, flies and mammals yields a large number of known and novel putative regulatory elements. Many of these are validated by independent biological observations, have spatial and/or orientation biases, are co-conserved with other elements and show surprising conservation across large phylogenetic distances.

Transcriptional activity of the paired-like homeodomain proteins CHX10 and VSX1

CHX10 and VSX1 are homeodomain (HD) proteins essential for normal retinal development. CHX10 is required first for retinal progenitor cell proliferation and later for bipolar cell differentiation, whereas VSX1 is important in the terminal differentiation of a subset of bipolar cells. Elucidating the transcriptional activity of CHX10 and VSX1 is required to understand how these factors control retinal development. We show that CHX10 and Vsx1 can function as transcriptional repressors. When tethered to a promoter by a heterologous LexA DNA-binding domain or its HD, CHX10 repressed multiple classes of activators in different immortalized cell lines. CHX10 blocked TATA-containing and TATA-less promoters, repressed at a distance, and inhibited a complex enhancer positioned upstream or downstream of the reporter gene, whereas retinoblastoma protein (RB) inhibited the downstream enhancer only. Interestingly, CHX10 mildly potentiated a subset of activators in chick neuronal cultures. Thus, CHX10 is both a versatile repressor and a context-specific weak activator. The CHX10 HD and CVC domains were sufficient for DNA binding and repression. VSX1 contains closely related homeo and CVC domains and, like CHX10, also repressed transcription. A VSX1 HD mutation, R166W, that impairs DNA binding and causes keratoconus in humans, hindered repressor function. Therefore, CHX10 and VSX1 may control retinal bipolar cell specification or differentiation by repressing genes required for the development of other cell types.

How to get ahead: the origin, evolution and function ofbicoid

In Drosophila, a Bcd protein gradient orchestrates patterning along the anteroposterior embryonic axis. However, studies of basal flies and other insects have revealed that bcd is a derived Hox3 gene found only in higher dipterans. To understand how bcd acquired its role in flies and how anteroposterior patterning mechanisms have evolved, I first review key features of bcd function in Drosophila: anterior localization and transcriptional and translation control of gene expression. I then discuss investigations of bcd in other higher dipterans that have provided insight into the evolution of regulatory interactions and the Bcd gradient. Finally, I review studies of Drosophila and other insects that address the evolution of bcd function and integration of bcd into ancestral regulatory mechanisms. I suggest further comparative studies may allow us to identify the intermediate steps in bcd evolution. This will make bcd a paradigm for the origin and evolution of genes and regulatory networks.

Structure, regulation, and function of micro1 in the sea urchin Hemicentrotus pulcherrimus

The animal-vegetal axis of sea urchin embryos is morphologically apparent at the 16-cell stage, when the mesomeres, macromeres, and micromeres align along it. At this stage, the micromere is the only autonomously specified blastomere that functions as a signaling center. We used a subtraction PCR survey to identify the homeobox gene micro1 as a micromere-specific gene. The micro1 gene is a representative of a novel family of paired-like class homeobox genes, along with PlHbox12 from Paracentrotus lividus and pmar1 from Strongylocentrotus purpuratus. In the present study, we showed that micro1 is a multicopy gene with six or more polymorphic loci, at least three of which are clustered in a 30-kb region of the genome. The micro1 gene is transiently expressed during early cleavage stages in the micromere. Recently, nuclear beta-catenin was shown to be essential for the specification of vegetal cell fates, including micromeres, and the temporal and spatial coincidence of micro1 expression with the nuclear entry of beta-catenin is highly suggestive. We demonstrated that micro1 is a direct target of beta-catenin. In addition, we showed that micro1 is necessary and sufficient for micromere specification. These observations on the structure, regulation, and function of micro1 lead to the conclusion that micro1 and pmar1 (and potentially PlHbox12) are orthologous.

A phage display selection of engrailed homeodomain mutants and the importance of residue Q50

Mutants of engrailed homeodomain (HD) that retain DNA-binding activity were isolated using a phage display selection. This selection was used to enrich for active DNA-binding clones from a complex library consisting of over a billion members. A more focused library of mutant homeodomains consisting of all possible amino acid combinations at two DNA-contacting residues (I47 and Q50) was constructed and screened for members capable of binding tightly and specifically to the engrailed consensus sequence, TAATTA. The isolated mutants largely recapitulated the distribution of amino acids found at these positions in natural homeodomains thus validating the in vitro selection conditions. In particular, the unequivocal advantage enjoyed by glutamine at residue 50 is surprising in light of reports that minimize the importance of this residue. Here, the subtle contributions of residue Q50 are demonstrated to play a functionally important role in specific recognition of DNA. These results highlight the complex subtlety of protein-DNA interactions, underscoring the value of the first reported in vitro selection of a homeodomain.

Otx2 regulates the extent, identity and fate of neuronal progenitor domains in the ventral midbrain

The specification of distinct neuronal cell-types is controlled by inducing signals whose interpretation in distinct areas along the central nervous system provides neuronal progenitors with a precise and typical expression code of transcription factors. To gain insights into this process, we investigated the role of Otx2 in the specification of identity and fate of neuronal progenitors in the ventral midbrain. To achieve this, Otx2 was inactivated by Cre recombinase under the transcriptional control of En1. Lack of Otx2 in the ventrolateral and posterior midbrain results in a dorsal expansion of Shh expression and in a dorsal and anterior rotation of the midbrain-hindbrain boundary and Fgf8 expression. Indeed, in this mutant correct positioning of the ventral site of midbrain-hindbrain boundary and Fgf8 expression are efficiently controlled by Otx1 function, thus allowing the study of the identity and fate of neuronal progenitors of the ventral midbrain in the absence of Otx2. Our results suggest that Otx2 acts in two ways: by repressing Nkx2.2 in the ventral midbrain and maintaining the Nkx6.1-expressing domain through dorsal antagonism on Shh. Failure of this control affects the identity code and fate of midbrain progenitors, which exhibit features in common with neuronal precursors of the rostral hindbrain even though the midbrain retains its regional identity and these neuronal precursors are rostral to Fgf8 expression. Dopaminergic neurons are greatly reduced in number, red nucleus precursors disappear from the ventral midbrain where a relevant number of serotonergic neurons are generated. These results indicate that Otx2 is an essential regulator of the identity, extent and fate of neuronal progenitor domains in the ventral midbrain and provide novel insights into the mechanisms by which neuronal diversity is generated in the central nervous system.

A HOX gene mutation in a family with isolated congenital vertical talus and Charcot-Marie-Tooth disease

Congenital vertical talus (CVT), also known as "rocker-bottom foot" deformity, is a dislocation of the talonavicular joint, with rigid dorsal dislocation of the navicular over the neck of the talus. This condition is usually associated with multiple other congenital deformities and only rarely is an isolated deformity. The reported familial cases are consistent with an autosomal dominant mode of inheritance with incomplete penetrance. In contrast, Charcot-Marie-Tooth disease (CMT) is thought to be a completely distinct heterogeneous group of disorders, with foot abnormalities that typically develop a high-arched "claw foot" appearance later in life. In the present study, DNA was isolated from 36 members of a single upstate (northern) New York white family of Italian descent in which both CVT and CMT were segregating. Whole-genome linkage analysis with Affymetrix GeneChip Mapping 10K Array defined a 7-Mb critical region on chromosome 2q31, which led to candidate-gene sequencing of six HOX genes and detection of a single missense mutation, M319K (956T-->A), in the HOXD10 gene. In the study family, this mutation was fully penetrant and exhibited significant evidence of linkage (LOD 6.33; theta =0), and it very likely accounts for both CVT and CMT in heterozygotes.

Dissecting the Engrailed homeodomain-DNA interaction by phage-displayed shotgun scanning

Phage-displayed alanine shotgun scanning was used to dissect contributions by engrailed homedomain (En-HD) residues 17 through 46, which indirectly influence recognition of DNA. The relative contributions of such indirect contacts, quantified by shotgun scanning, highlight previously unexplored En-HD residues. Two motifs dominate En-HD function in this region. First, two surface-exposed aromatic residues (F20 and Y25) bracket the hydrophobic core. Second, two sets of turn-forming residues are highlighted, including carboxamide-requiring residues E22/N23 and a leucine/isoleucine splint. The En-HD hydrophobic core exhibits a surprising degree of malleability, as demonstrated by homolog shotgun scanning. Most selectants from in vitro shotgun scanning mirror the consensus human homeodomain sequence. Thus, natural evolution and in vitro selection use similar selection criteria: affinity, specificity, and stability. However, homolog shotgun scanning identified mutations capable of improving the affinity and specificity of En-HD.

From sequence to structure and back again: approaches for predicting protein-DNA binding

Gene regulation in higher organisms is achieved by a complex network of transcription factors (TFs). Modulating gene expression and exploring gene function are major aims in molecular biology. Furthermore, the identification of putative target genes for a certain TF serve as powerful tools for specific targeting of rational drugs.

Detecting the short and variable transcription factor binding sites (TFBSs) in genomic DNA is an intriguing challenge for computational and structural biologists. Fast and reliable computational methods for predicting TFBSs on a whole-genome scale offer several advantages compared to the current experimental methods that are rather laborious and slow. Two main approaches are being explored, advanced sequence-based algorithms and structure-based methods.

The aim of this review is to outline the computational and experimental methods currently being applied in the field of protein-DNA interactions. With a focus on the former, the current state of the art in modeling these interactions is discussed. Surveying sequence and structure-based methods for predicting TFBSs, we conclude that in order to achieve a sound and specific method applicable on genomic sequences it is desirable and important to bring these two approaches together.

The human tyrosine hydroxylase gene promoter

13.329 kilobases of the single copy human tyrosine hydroxylase (hTH) gene were isolated from a genomic library. The 5' flanking 11 kilobases fused to the reporter green fluorescent protein (GFP) drove high level expression in TH+ cells of the substantia nigra of embryonic and adult transgenic mice as determined by double label fluorescence microscopy. To provide a basis for future analysis of polymorphisms and structure-function studies, the previously unreported distal 10.5 kilobases of the hTH promoter were sequenced with an average coverage of 20-fold, the remainder with 4-fold coverage. Sequence features identified included four perfect matches to the bicoid binding element (BBE, consensus: BBTAATCYV) all of which exhibited specific binding by electrophoretic mobility shift assay (EMSA). Comparison to published sequences of mouse and rat TH promoters revealed five areas of exceptional homology shared by these species in the upstream TH promoter region -2 kb to -9 kb relative to the transcription start site. Within these conserved regions (CRs I-V), potential recognition sites for NR4A2 (Nurr1), HNF-3beta, HOXA4, and HOXA5 were shared across human, mouse, and rat TH promoters.

DNA microarrays with unimolecular hairpin double-stranded DNA probes: fabrication and exploration of sequence-specific DNA/protein interactions

We have fabricated double-stranded DNA (dsDNA) microarrays containing unimolecular hairpin dsDNA probes immobilized on glass slides. The unimolecular hairpin dsDNA microarrays were manufactured by four steps: Firstly, synthesizing single-stranded DNA (ssDNA) oligonucleotides with two reverse-complementary sequences at 3' hydroxyl end and an overhang sequence at 5' amino end. Secondly, microspotting ssDNA on glutaraldehyde-derived glass slide to form ssDNA microarrays. Thirdly, annealing two reverse-complementary sequences to form hairpin primer at 3' end of immobilized ssDNA and thus to create partial-dsDNA microarray. Fourthly, enzymatically extending hairpin primer to convert partial-dsDNA microarrays into complete-dsDNA microarray. The excellent efficiency and high accuracy of the enzymatic synthesis were demonstrated by incorporation of fluorescently labeled dUTPs in Klenow extension and digestion of dsDNA microarrays with restriction endonuclease. The accessibility and specificity of the DNA-binding proteins binding to dsDNA microarrays were verified by binding Cy3-labeled NF-kappaB to dsDNA microarrays. The dsDNA microarrays have great potential to provide a high-throughput platform for investigation of sequence-specific DNA/protein interactions involved in gene expression regulation, restriction and so on.

Stepwise induced fit in the pico- to nanosecond time scale governs the complexation of the even-skipped transcriptional repressor homeodomain to DNA

Induced fit effects in the complex of a DNA decamer with two even-skipped transcriptional repressor homeodomain molecules were investigated by means of molecular dynamics simulations. Dynamics of these effects are found to be in the time scale from pico- to nanoseconds. First steps are made by the fast-moving DNA backbone phosphates, which upon binding change their B(I)/B(II) substate distribution. Further rearrangements in the DNA double helix induced upon complexation, like bending of the helix axis, changes of the minor groove width, and of different helical parameters, are slower and occur within a few nanoseconds. The flexibility of the DNA, especially of its backbone, seems thereby to play an important role for specific DNA ligand recognition.

Fabrication of Unimolecular Double-stranded DNA Microarrays on Solid Surfaces for Probing DNA-Protein/Drug Interactions

We present a novel method for fabricating unimolecular double-stranded DNA microarrays on solid surfaces, which were used to probe sequence-specific DNA/protein interactions. For manufacturing the unimolecular double-stranded DNA microarrays, two kinds of special single-stranded oligonucleotides, constant oligonucleotide and target oligonucleotide, were chemically synthesized. The constant oligonucleotides with internal aminated dT were used to capture and immobilize the target oligonucleotides onto the solid surface, and also to provide a primer for later enzymatic extension reactions, while target oligonucleotides took the role of harbouring DNA-binding sites of DNA-binding proteins. The variant target oligonucleotides were annealed and ligated with the constant oligonucleotides to form the new unimolecular oligonucleotides for microspotting. The prepared unimolecular oligonucleotides were microspotted on aldehyde-derivatized glass slides to make partial-dsDNA microarrays. Finally, the partial-dsDNA microarrays were converted into a unimolecular complete-dsDNA microarray by a DNA polymerase extension reaction. The efficiency and accuracy of the polymerase synthesis were demonstrated by the fluorescent-labeled dUTP incorporation in the enzymatic extension reaction and the restriction endonuclease digestion of the fabricated unimolecular complete-dsDNA microarray. The accessibility and specificity of the sequence-specific DNA-binding proteins binding to the immobilized unimolecular dsDNA probes were demonstrated by the binding of Cy3 labeled NF-κB (p50·p50) to the unimolecular dsDNA microarray. This unimolecular dsDNA microarray provides a general technique for high-throughput DNA-protein or DNA-drugs interactions.