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

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

Generation of human excitatory forebrain neurons by cooperative binding of proneural NGN2 and homeobox factor EMX1

Generation of defined neuronal subtypes from human pluripotent stem cells remains a challenge. The proneural factor NGN2 has been shown to overcome experimental variability observed by morphogen-guided differentiation and directly converts pluripotent stem cells into neurons, but their cellular heterogeneity has not been investigated yet. Here, we found that NGN2 reproducibly produces three different kinds of excitatory neurons characterized by partial coactivation of other neurotransmitter programs. We explored two principle approaches to achieve more precise specification: prepatterning the chromatin landscape that NGN2 is exposed to and combining NGN2 with region-specific transcription factors. Unexpectedly, the chromatin context of regionalized neural progenitors only mildly altered genomic NGN2 binding and its transcriptional response and did not affect neurotransmitter specification. In contrast, coexpression of region-specific homeobox factors such as EMX1 resulted in drastic redistribution of NGN2 including recruitment to homeobox targets and resulted in glutamatergic neurons with silenced nonglutamatergic programs. These results provide the molecular basis for a blueprint for improved strategies for generating a plethora of defined neuronal subpopulations from pluripotent stem cells for therapeutic or disease-modeling purposes.

Case report: A novel R246L mutation in the LMX1B homeodomain causes isolated nephropathy in a large Chinese family

BACKGROUND

Genetic factors contribute to chronic kidney disease (CKD) and end-stage renal disease (ESRD). Advances in genetic testing have enabled the identification of hereditary kidney diseases, including those caused by LMX1B mutations. LMX1B mutations can lead to nail-patella syndrome (NPS) or nail-patella-like renal disease (NPLRD) with only renal manifestations.

CASE PRESENTATION

The proband was a 13-year-old female who was diagnosed with nephrotic syndrome at the age of 6. Then she began intermittent hormone and drug therapy. When she was 13 years old, she was admitted to our hospital due to sudden chest tightness, which progressed to end-stage kidney disease (ESRD), requiring kidney replacement therapy. Whole-Exome Sequencing (WES) results suggest the presence of LMX1B gene mutation, c.737G > T, p.Arg246Leu. Tracing her family history, we found that her father, grandmother, uncle and 2 cousins all had hematuria, or proteinuria. In addition to the grandmother, a total of 9 members of the family performed WES. The members with kidney involved all carry the mutated gene. Healthy members did not have the mutated gene. It is characterized by co-segregation of genotype and phenotype. We followed the family for 9 year, the father developed ESRD at the age of 50 and started hemodialysis treatment. The rest patients had normal renal function. No extra-renal manifestations associated with NPS were found in any member of the family.

CONCLUSIONS

This study has successfully identified missense mutation, c.737G > T (p.Arg246Leu) in the homeodomain, which appears to be responsible for isolated nephropathy in the studied family. The arginine to leucine change at codon 246 likely disrupts the DNA-binding homeodomain of LMX1B. Previous research has documented 2 types of mutations at codon R246, namely R246Q and R246P, which are known to cause NPLRD. The newly discovered mutation, R246L, is likely to be another novel mutation associated with NPLRD, thus expanding the range of mutations at the crucial renal-critical codon 246 that contribute to the development of NPLRD. Furthermore, our findings suggest that any missense mutation occurring at the 246th amino acid position within the homeodomain of the LMX1B gene has the potential to lead to NPLRD.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Transcriptional Landscape of CUT-Class Homeobox Genes in Blastic Plasmacytoid Dendritic Cell Neoplasm

Homeobox genes encode developmental transcription factors regulating tissue-specific differentiation processes and drive cancerogenesis when deregulated. Dendritic cells (DCs) are myeloid immune cells occurring as two types, either conventional or plasmacytoid DCs. Recently, we showed that the expression of NKL-subclass homeobox gene VENTX is restricted to conventional DCs, regulating developmental genes. Here, we identified and investigated homeobox genes specifically expressed in plasmacytoid DCs (pDCs) and derived blastic plasmacytoid dendritic cell neoplasm (BPDCN). We analyzed gene expression data, performed RQ-PCR, protein analyses by Western blot and immuno-cytology, siRNA-mediated knockdown assays and subsequent RNA-sequencing and live-cell imaging. Screening of public gene expression data revealed restricted activity of the CUT-class homeobox gene CUX2 in pDCs. An extended analysis of this homeobox gene class in myelopoiesis showed that additional CUX2 activity was restricted to myeloid progenitors, while BPDCN patients aberrantly expressed ONECUT2, which remained silent in the complete myeloid compartment. ONECUT2 expressing BPDCN cell line CAL-1 served as a model to investigate its regulation and oncogenic activity. The ONECUT2 locus at 18q21 was duplicated and activated by IRF4, AUTS2 and TNF-signaling and repressed by BMP4-, TGFb- and IL13-signalling. Functional analyses of ONECUT2 revealed the inhibition of pDC differentiation and of CDKN1C and CASP1 expression, while SMAD3 and EPAS1 were activated. EPAS1 in turn enhanced survival under hypoxic conditions which thus may support dendritic tumor cells residing in hypoxic skin lesions. Collectively, we revealed physiological and aberrant activities of CUT-class homeobox genes in myelopoiesis including pDCs and in BPDCN, respectively. Our data may aid in the diagnosis of BPDCN patients and reveal novel therapeutic targets for this fatal malignancy.

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.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

Notch signaling and Bsh homeodomain activity are integrated to diversify Drosophila lamina neuron types

Notch signaling is an evolutionarily conserved pathway for specifying binary neuronal fates, yet how it specifies different fates in different contexts remains elusive. In our accompanying paper, using the Drosophila lamina neuron types (L1-L5) as a model, we show that the primary homeodomain transcription factor (HDTF) Bsh activates secondary HDTFs Ap (L4) and Pdm3 (L5) and specifies L4/L5 neuronal fates. Here we test the hypothesis that Notch signaling enables Bsh to differentially specify L4 and L5 fates. We show asymmetric Notch signaling between newborn L4 and L5 neurons, but they are not siblings; rather, Notch signaling in L4 is due to Delta expression in adjacent L1 neurons. While Notch signaling and Bsh expression are mutually independent, Notch is necessary and sufficient for Bsh to specify L4 fate over L5. The NotchON L4, compared to NotchOFF L5, has a distinct open chromatin landscape which allows Bsh to bind distinct genomic loci, leading to L4-specific identity gene transcription. We propose a novel model in which Notch signaling is integrated with the primary HDTF activity to diversify neuron types by directly or indirectly generating a distinct open chromatin landscape that constrains the pool of genes that a primary HDTF can activate.

Homeodomain proteins hierarchically specify neuronal diversity and synaptic connectivity

How our brain generates diverse neuron types that assemble into precise neural circuits remains unclear. Using Drosophila lamina neuron types (L1-L5), we show that the primary homeodomain transcription factor (HDTF) brain-specific homeobox (Bsh) is initiated in progenitors and maintained in L4/L5 neurons to adulthood. Bsh activates secondary HDTFs Ap (L4) and Pdm3 (L5) and specifies L4/L5 neuronal fates while repressing the HDTF Zfh1 to prevent ectopic L1/L3 fates (control: L1-L5; Bsh-knockdown: L1-L3), thereby generating lamina neuronal diversity for normal visual sensitivity. Subsequently, in L4 neurons, Bsh and Ap function in a feed-forward loop to activate the synapse recognition molecule DIP-β, thereby bridging neuronal fate decision to synaptic connectivity. Expression of a Bsh:Dam, specifically in L4, reveals Bsh binding to the DIP-β locus and additional candidate L4 functional identity genes. We propose that HDTFs function hierarchically to coordinate neuronal molecular identity, circuit formation, and function. Hierarchical HDTFs may represent a conserved mechanism for linking neuronal diversity to circuit assembly and function.

A novel pathogenic variant c.262delA in PBX1 causing oligomeganephronia identified using whole-exome sequencing and a literature review

Oligomeganephronia (OMN) is a rare congenital renal hypoplasia reported more often in children than in adults. The diagnosis of OMN relies on renal biopsy and exhibits a significant reduction in the number of glomeruli and pronounced glomerular hypertrophy. Here, we report the case of an 8-year-old boy with recurrent proteinuria and abnormal external ears. A renal biopsy revealed large and rare glomeruli. The histological findings confirmed the diagnosis of OMN. Whole-exome sequencing of the patient revealed a new pathogenic variant in PBX1 (hg19, NM_002585, c.262delA, p.Thr88Glnfs*3). The PBX1 gene encodes a transcription factor whose pathogenic variants can result in congenital renal and urinary system anomalies, with or without hearing loss, abnormal ears, and developmental retardation (CAKUTED). This is the first report to detect PBX1 pathogenic variants in children with OMN, a novel phenotype of human PBX1 pathogenic variants. We performed functional prediction analyses of deletions in the corresponding structural domains. We summarized 27 cases of PBX1 single pathogenic variants reported between 2003 and 2023 in terms of truncating and missense pathogenic variants, which can deepen our understanding of the PBX1 structural domain and expand our knowledge of the PBX1 genotype and phenotype.

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.

Transcriptional control of embryonic and adult neural progenitor activity

Neural precursors generate neurons in the embryonic brain and in restricted niches of the adult brain in a process called neurogenesis. The precise control of cell proliferation and differentiation in time and space required for neurogenesis depends on sophisticated orchestration of gene transcription in neural precursor cells. Much progress has been made in understanding the transcriptional regulation of neurogenesis, which relies on dose- and context-dependent expression of specific transcription factors that regulate the maintenance and proliferation of neural progenitors, followed by their differentiation into lineage-specified cells. Here, we review some of the most widely studied neurogenic transcription factors in the embryonic cortex and neurogenic niches in the adult brain. We compare functions of these transcription factors in embryonic and adult neurogenesis, highlighting biochemical, developmental, and cell biological properties. Our goal is to present an overview of transcriptional regulation underlying neurogenesis in the developing cerebral cortex and in the adult brain.

Chromosomal Aberration t(14;17)(q32;q21) Simultaneously Activates HOXB5 and miR10a in Triple-Hit B-Cell Lymphoma

BCL2, BCL6 and MYC are major oncogenes in B-cell lymphoma. Their aberrant activation frequently occurs via chromosomal translocations which juxtapose light or heavy chain immunoglobulin (IG) genes to BCL2 and MYC or fuse diverse partner genes with BCL6. So-called double-hit lymphomas usually carry BCL2 and MYC rearrangements, while triple-hit lymphomas additionally bear BCL6-fusions. All these translocations are of diagnostic relevance and usually denote poor prognosis. Here, we genomically characterized classic follicular lymphoma (FL) cell line SC-1, thereby identifying t(14;18)(q32;q21) juxtaposing IGH and BCL2, t(8;14)(q24;q32) juxtaposing IGH and MYC, and t(3;3)(q25;q27) fusing MBNL1 to BCL6. In addition, we found that SC-1 carries a novel chromosomal rearrangement, t(14;17)(q32;q21), which, though present at establishment, has remained unreported until now. We further show that t(14;17)(q32;q21) juxtaposes IGH with the HOXB gene cluster at 17q21 and affect the oncogenic activation of both homeobox gene HOXB5 and neighboring micro-RNA gene miR10a. Moreover, we detected aberrant overexpression of HOXB5 in subsets of Burkitt lymphoma, FL, and multiple myeloma patients, confirming the clinical relevance of its deregulation. In SC-1, HOXB5 activation was additionally supported by co-expression of hematopoietic stem cell factor ZNF521, indicating an aberrant impact in cell differentiation. Functional investigations showed that HOXB5 represses the apoptotic driver BCL2L11 and promotes survival in the presence of etoposide, and that miR10a inhibits BCL6 and may thus play an oncogenic role in later stages of lymphomagenesis. Collectively, we characterize triple-hit B-cell line SC-1 and identify the aberrant expression of HOXB5 and miR10a, both novel oncogenes in B-cell lymphoma.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Rapid evolution of the embryonically-expressed homeobox gene LEUTX within primates

LEUTX is a homeodomain transcription factor expressed in the very early embryo with a function around embryonic genome activation. The LEUTX gene is found only in eutherian mammals, including humans, but unlike the majority of homeobox genes, the encoded amino acid sequence is very different between divergent mammalian species. However, whether dynamic evolution has also occurred between closely related mammalian species remains unclear. In this work, we perform a comparative genomics study of LEUTX within the primates, revealing dramatic evolutionary sequence change between closely related species. Positive selection has acted on sites in the LEUTX protein, including six sites within the homeodomain; this suggests that selection has driven changes in the set of downstream targets. Transfection into cell culture followed by transcriptomic analysis reveals small functional differences between human and marmoset LEUTX, suggesting rapid sequence evolution has fine-tuned the role of this homeodomain protein within the primates.

Unanticipated broad phylogeny of BEN DNA-binding domains revealed by structural homology searches

Organization of protein sequences into domain families is a foundation for cataloging and investigating protein functions. However, long-standing strategies based on primary amino acid sequences are blind to the possibility that proteins with dissimilar sequences could have comparable tertiary structures. Building on our recent findings that in silico structural predictions of BEN family DNA-binding domains closely resemble their experimentally determined crystal structures, we exploited the AlphaFold2 database for comprehensive identification of BEN domains. Indeed, we identified numerous novel BEN domains, including members of new subfamilies. For example, while no BEN domain factors had previously been annotated in C. elegans, this species actually encodes multiple BEN proteins. These include key developmental timing genes of orphan domain status, sel-7 and lin-14, the latter being the central target of the founding miRNA lin-4. We also reveal that the domain of unknown function 4806 (DUF4806), which is widely distributed across metazoans, is structurally similar to BEN and comprises a new subtype. Surprisingly, we find that BEN domains resemble both metazoan and non-metazoan homeodomains in 3D conformation and preserve characteristic residues, indicating that despite their inability to be aligned by conventional methods, these DNA-binding modules are probably evolutionarily related. Finally, we broaden the application of structural homology searches by revealing novel human members of DUF3504, which exists on diverse proteins with presumed or known nuclear functions. Overall, our work strongly expands this recently identified family of transcription factors and illustrates the value of 3D structural predictions to annotate protein domains and interpret their functions.

The Role of IRX Homeobox Genes in Hematopoietic Progenitors and Leukemia

IRX genes are members of the TALE homeobox gene class and encode six related transcription factors (IRX1-IRX6) controlling development and cell differentiation of several tissues in humans. Classification of TALE homeobox gene expression patterns for the hematopoietic compartment, termed TALE-code, has revealed exclusive IRX1 activity in pro-B-cells and megakaryocyte erythroid progenitors (MEPs), highlighting its specific contribution to developmental processes at these early stages of hematopoietic lineage differentiation. Moreover, aberrant expression of IRX homeobox genes IRX1, IRX2, IRX3 and IRX5 has been detected in hematopoietic malignancies, including B-cell precursor acute lymphoblastic leukemia (BCP-ALL), T-cell ALL, and some subtypes of acute myeloid leukemia (AML). Expression analyses of patient samples and experimental studies using cell lines and mouse models have revealed oncogenic functions in cell differentiation arrest and upstream and downstream genes, thus, revealing normal and aberrant regulatory networks. These studies have shown how IRX genes play key roles in the development of both normal blood and immune cells, and hematopoietic malignancies. Understanding their biology serves to illuminate developmental gene regulation in the hematopoietic compartment, and may improve diagnostic classification of leukemias in the clinic and reveal new therapeutic targets and strategies.

Determination and Analysis of Hyper-Variable A Mating Types in Wild Strains of Lentinula edodes in Korea

The diversity of A mating type in wild strains of Lentinula edodes was extensively analyzed to characterize and utilize them for developing new cultivars. One hundred twenty-three A mating type alleles, including 67 newly discovered alleles, were identified from 106 wild strains collected for the past four decades in Korea. Based on previous studies and current findings, a total of 130 A mating type alleles have been found, 124 of which were discovered from wild strains, indicating the hyper-variability of A mating type alleles of L. edodes. About half of the A mating type alleles in wild strains were found in more than two strains, whereas the other half of the alleles were found in only one strain. About 90% of A mating type combinations in dikaryotic wild strains showed a single occurrence. Geographically, diverse A mating type alleles were intensively located in the central region of the Korean peninsula, whereas only allele A17 was observed throughout Korea. We also found the conservation of the TCCCAC motif in addition to the previously reported motifs, including ATTGT, ACAAT, and GCGGAG, in the intergenic regions of A mating loci. Sequence comparison among some alleles indicated that accumulated mutation and recombination would contribute to the diversification of A mating type alleles in L. edodes. Our data support the rapid evolution of A mating locus in L. edodes, and would help to understand the characteristics of A mating loci of wild strains in Korea and help to utilize them for developing new cultivars.

An example of parenchymal renal sparing in the context of complex malformations due to a novel mutation in the PBX1 gene

INTRODUCTION

PBX1 encodes the pre-B cell leukemia factor 1, a Three Amino acid Loop Extension (TALE) transcription factor crucial to regulate basic developmental processes. PBX1 loss-of-function variants have been initially described in association with renal malformations in both isolated and syndromic forms.

CASE REPORT

Herein, we report a male infant presenting multiple organ malformations (cleidosternal dysostosis, micrognathia, left lung hypoplasia, wide interatrial defect, pulmonary hypertension, total anomalous pulmonary venous return, intestinal malrotation) and carrying the heterozygous de novo c.868C > T (p.Arg290Trp) variant in PBX1. This novel variant affects the highly conserved homeodomain of the protein, leading to a non-conservative substitution and consequently altering its tridimensional structure and DNA-binding capacity.

CONCLUSION

So far, PBX1 has been reported in association with a broad spectrum of renal anomalies. However, given the role of this gene in many different developing processes, whole-exome sequencing can detect mutations in PBX1 even in patients with different phenotypes, not necessarily involving the renal primordium. This report presents a novel PBX1 variant with a predicted strong deleterious effect. The mutation leads to a non-conservative substitution in a very highly conserved domain of the protein, thus altering its tertiary structure and DNA-binding capacity.

The Hematopoietic TALE-Code Shows Normal Activity of IRX1 in Myeloid Progenitors and Reveals Ectopic Expression of IRX3 and IRX5 in Acute Myeloid Leukemia

Homeobox genes encode transcription factors that control basic developmental decisions. Knowledge of their hematopoietic activities casts light on normal and malignant immune cell development. Recently, we constructed the so-called lymphoid TALE-code that codifies expression patterns of all active TALE class homeobox genes in early hematopoiesis and lymphopoiesis. Here, we present the corresponding myeloid TALE-code to extend this gene signature, covering the entire hematopoietic system. The collective data showed expression patterns for eleven TALE homeobox genes and highlighted the exclusive expression of IRX1 in megakaryocyte-erythroid progenitors (MEPs), implicating this TALE class member in a specific myeloid differentiation process. Analysis of public profiling data from acute myeloid leukemia (AML) patients revealed aberrant activity of IRX1 in addition to IRX3 and IRX5, indicating an oncogenic role for these TALE homeobox genes when deregulated. Screening of RNA-seq data from 100 leukemia/lymphoma cell lines showed overexpression of IRX1, IRX3, and IRX5 in megakaryoblastic and myelomonocytic AML cell lines, chosen as suitable models for studying the regulation and function of these homeo-oncogenes. Genomic copy number analysis of IRX-positive cell lines demonstrated chromosomal amplification of the neighboring IRX3 and IRX5 genes at position 16q12 in MEGAL, underlying their overexpression in this cell line model. Comparative gene expression analysis of these cell lines revealed candidate upstream factors and target genes, namely the co-expression of GATA1 and GATA2 together with IRX1, and of BMP2 and HOXA10 with IRX3/IRX5. Subsequent knockdown and stimulation experiments in AML cell lines confirmed their activating impact in the corresponding IRX gene expression. Furthermore, we demonstrated that IRX1 activated KLF1 and TAL1, while IRX3 inhibited GATA1, GATA2, and FST. Accordingly, we propose that these regulatory relationships may represent major physiological and oncogenic activities of IRX factors in normal and malignant myeloid differentiation, respectively. Finally, the established myeloid TALE-code is a useful tool for evaluating TALE homeobox gene activities in AML.

The Role of NKL Homeobox Genes in T-Cell Malignancies

Homeobox genes encode transcription factors controlling basic developmental processes. The homeodomain is encoded by the homeobox and mediates sequence-specific DNA binding and interaction with cofactors, thus operating as a basic regulatory platform. Similarities in their homeobox sequences serve to arrange these genes in classes and subclasses, including NKL homeobox genes. In accordance with their normal functions, deregulated homeobox genes contribute to carcinogenesis along with hematopoietic malignancies. We have recently described the physiological expression of eleven NKL homeobox genes in the course of hematopoiesis and termed this gene expression pattern NKL-code. Due to the developmental impact of NKL homeobox genes these data suggest a key role for their activity in the normal regulation of hematopoietic cell differentiation including T-cells. On the other hand, aberrant overexpression of NKL-code members or ectopical activation of non-code members has been frequently reported in lymphoid and myeloid leukemia/lymphoma, demonstrating their oncogenic impact in the hematopoietic compartment. Here, we provide an overview of the NKL-code in normal hematopoiesis and discuss the oncogenic role of deregulated NKL homeobox genes in T-cell malignancies.

Design and Characterization of a Cell-Penetrating Peptide Derived from the SOX2 Transcription Factor

SOX2 is an oncogenic transcription factor overexpressed in nearly half of the basal-like triple-negative breast cancers associated with very poor outcomes. Targeting and inhibiting SOX2 is clinically relevant as high SOX2 mRNA levels are positively correlated with decreased overall survival and progression-free survival in patients affected with breast cancer. Given its key role as a master regulator of cell proliferation, SOX2 represents an important scaffold for the engineering of dominant-negative synthetic DNA-binding domains (DBDs) that act by blocking or interfering with the oncogenic activity of the endogenous transcription factor in cancer cells. We have synthesized an interference peptide (iPep) encompassing a truncated 24 amino acid long C-terminus of SOX2 containing a potential SOX-specific nuclear localization sequence, and the determinants of the binding of SOX2 to the DNA and to its transcription factor binding partners. We found that the resulting peptide (SOX2-iPep) possessed intrinsic cell penetration and promising nuclear localization into breast cancer cells, and decreased cellular proliferation of SOX2 overexpressing cell lines. The novel SOX2-iPep was found to exhibit a random coil conformation predominantly in solution. Molecular dynamics simulations were used to characterize the interactions of both the SOX2 transcription factor and the SOX2-iPep with FGF4-enhancer DNA in the presence of the POU domain of the partner transcription factor OCT4. Predictions of the free energy of binding revealed that the iPep largely retained the binding affinity for DNA of parental SOX2. This work will enable the future engineering of novel dominant interference peptides to transport different therapeutic cargo molecules such as anti-cancer drugs into cells.

OTX2 Homeoprotein Functions in Adult Choroid Plexus

The choroid plexus is an important blood barrier that secretes cerebrospinal fluid, which essential for embryonic brain development and adult brain homeostasis. The OTX2 homeoprotein is a transcription factor that is critical for choroid plexus development and remains highly expressed in adult choroid plexus. Through RNA sequencing analyses of constitutive and conditional knockdown adult mouse models, we reveal putative functional roles for OTX2 in adult choroid plexus function, including cell signaling and adhesion, and show that OTX2 regulates the expression of factors that are secreted into the cerebrospinal fluid, notably transthyretin. We also show that Otx2 expression impacts choroid plexus immune and stress responses, and affects splicing, leading to changes in the mRNA isoforms of proteins that are implicated in the oxidative stress response and DNA repair. Through mass spectrometry analysis of OTX2 protein partners in the choroid plexus, and in known non-cell-autonomous target regions, such as the visual cortex and subventricular zone, we identify putative targets that are involved in cell adhesion, chromatin structure, and RNA processing. Thus, OTX2 retains important roles for regulating choroid plexus function and brain homeostasis throughout life.

Suboptimal Intermediates Underlie Evolution of the Bicoid Homeodomain

Changes in regulatory networks generate materials for evolution to create phenotypic diversity. For transcription networks, multiple studies have shown that alterations in binding sites of cis-regulatory elements correlate well with the gain or loss of specific features of the body plan. Less is known about alterations in the amino acid sequences of the transcription factors (TFs) that bind these elements. Here we study the evolution of Bicoid (Bcd), a homeodomain (HD) protein that is critical for anterior embryo patterning in Drosophila. The ancestor of Bcd (AncBcd) emerged after a duplication of a Zerknullt (Zen)-like ancestral protein (AncZB) in a suborder of flies. AncBcd diverged from AncZB, gaining novel transcriptional and translational activities. We focus on the evolution of the HD of AncBcd, which binds to DNA and RNA, and is comprised of four subdomains: an N-terminal arm (NT) and three helices; H1, H2, and Recognition Helix (RH). Using chimeras of subdomains and gene rescue assays in Drosophila, we show that robust patterning activity of the Bcd HD (high frequency rescue to adulthood) is achieved only when amino acid substitutions in three separate subdomains (NT, H1, and RH) are combined. Other combinations of subdomains also yield full rescue, but with lower penetrance, suggesting alternative suboptimal activities. Our results suggest a multistep pathway for the evolution of the Bcd HD that involved intermediate HD sequences with suboptimal activities, which constrained and enabled further evolutionary changes. They also demonstrate critical epistatic forces that contribute to the robust function of a DNA-binding domain.

Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms

Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.

The rax homeobox gene is mutated in the eyeless axolotl, Ambystoma mexicanum

BACKGROUND

Vertebrate eye formation requires coordinated inductive interactions between different embryonic tissue layers, first described in amphibians. A network of transcription factors and signaling molecules controls these steps, with mutations causing severe ocular, neuronal, and craniofacial defects. In eyeless mutant axolotls, eye morphogenesis arrests at the optic vesicle stage, before lens induction, and development of ventral forebrain structures is disrupted.

RESULTS

We identified a 5-bp deletion in the rax (retina and anterior neural fold homeobox) gene, which was tightly linked to the recessive eyeless (e) axolotl locus in an F2 cross. This frameshift mutation, in exon 2, truncates RAX protein within the homeodomain (P154fs35X). Quantitative RNA analysis shows that mutant and wild-type rax transcripts are equally abundant in E/e embryos. Translation appears to initiate from dual start codons, via leaky ribosome scanning, a conserved feature among gnathostome RAX proteins. Previous data show rax is expressed in the optic vesicle and diencephalon, deeply conserved among metazoans, and required for eye formation in other species.

CONCLUSION

The eyeless axolotl mutation is a null allele in the rax homeobox gene, with primary defects in neural ectoderm, including the retinal and hypothalamic primordia.

Nail-Patella syndrome with early onset end-stage renal disease in a child with a novel heterozygous missense mutation in the LMX1B homeodomain: A case report

Nail-Patella syndrome (NPS) is an inherited disease characterized by nail and skeletal anomalies, nephropathy and glaucoma. The diagnosis of NPS is based on clinical findings, including hypoplastic or absent patella, dystrophic nails, dysplasia of the elbows and iliac horns. However, the main determinant of NPS prognosis is nephropathy, which may range from asymptomatic proteinuria to end-stage renal disease. NPS is caused by heterozygous loss-of-function mutations in the LMX1B gene, which encodes the LIM homeodomain transcription factor LMX1B. LMX1B serves an essential role in the physiological development of dorsal-ventral limb structures, morphogenesis and function of podocytes, as well as in development of the anterior segments of the eyes, and in certain types of neurons. The present study aimed to identify the disease-causing mutation in a 2-year old girl with nephrotic syndrome that evolved rapidly to end-stage renal disease. The patient showed classical symptoms of NPS including dystrophic nails and an absence of the patellae. DNA sequence analysis identified a novel missense variant in exon 4 of LMX1B (c.709T>C, p.S237P); this substitution affected a conserved serine residue in the homeodomain of LMX1B and was predicted to be pathogenic. In silico modeling of the homeodomain revealed that the p.S237P mutation converted the A236-S237-F238 segment of α-helix 1 into a strand. It was hypothesized that this mutation affected binding of the transcription factor to its target DNA, thus abrogating transcription activation, which would explain the phenotype that manifested in the patient.

PBX1 promotes development of natural killer cells by binding directly to the Nfil3 promoter

The transcription factor nuclear factor interleukin-3-regulated protein (NFIL3, also called E4BP4) is crucial for commitment of natural killer (NK) cells from common lymphoid progenitors (CLPs). However, the identity of the factor that can regulate NFIL3 directly during the NK-cell development is not known. Here, we reveal that pre-B-cell leukemia transcription factor 1 (PBX1) can upregulate the NFIL3 expression directly. We used conditional knockout mice in which PBX1 in hematopoietic cells was specifically absent. The number of NK-committed progenitor pre-NKP cells and rNKP cells was reduced significantly in the absence of PBX1, which was consistent with NFIL3 deficiency. Also, the NFIL3 expression in NK cells was decreased if PBX1 was absent. We demonstrated that PBX1 was bound directly to the promoter of Nfil3 and facilitated transcription. Upon knockout of the binding site of PBX1 in the Nfil3 promoter, mice showed fewer NK-precursor cells and NK cells, just like that observed in Nfil3 knockout mice. Furthermore, asparagine N286 in the homeodomain of PBX1 controlled the binding of PBX1 to the Nfil3 promoter. Collectively, these findings demonstrate that the transcription factor PBX1 promotes the early development of NK cells by upregulating the Nfil3 expression directly.

Evolutionary Analysis of the Zinc Finger and Homeoboxes Family of Proteins Identifies Multiple Conserved Domains and a Common Early Chordate Ancestor

The Zinc Fingers and Homeoboxes (Zhx) proteins, Zhx1, Zhx2, and Zhx3, comprise a small family of proteins containing two amino-terminal C₂–H₂ zinc fingers and four or five carboxy-terminal homeodomains. These multiple homeodomains make Zhx proteins unusual because the majority of homeodomain-containing proteins contain a single homeodomain. Studies in cultured cells and mice suggest that Zhx proteins can function as positive or negative transcriptional regulators. Zhx2 regulates numerous hepatic genes, and all three Zhx proteins have been implicated in different cancers. Because Zhx proteins contain multiple predicted homeodomains, are associated with interesting physiological traits, and seem to be only present in the vertebrate lineage, we investigated the evolutionary history of this small family by comparing Zhx homologs from a wide range of chordates. This analysis indicates that the zinc finger motifs and homeodomains are highly similar among all Zhx proteins and also identifies additional Zhx-specific conserved regions, including a 13 amino acid amino-terminal motif that is nearly identical among all gnathostome Zhx proteins. We found single Zhx proteins in the sea lamprey (Petromyzon marinus) and in the nonvertebrate chordates sea squirt (Ciona intestinalis) and lancelet (Branchiostoma floridae); these Zhx proteins are most similar to gnathostome Zhx3. Based on our analyses, we propose that a duplication of the primordial Zhx gene gave rise to Zhx3 and the precursor to Zhx1 and Zhx2. A subsequent tandem duplication of this precursor generated Zhx1 and Zhx2 found in gnathostomes.

The Interaction of Natural Selection and GC Skew May Drive the Fast Evolution of a Sand Rat Homeobox Gene

Several processes can lead to strong GC skew in localized genomic regions. In most cases, GC skew should not affect conserved amino acids because natural selection will purge deleterious alleles. However, in the gerbil subfamily of rodents, several conserved genes have undergone radical alteration in association with strong GC skew. An extreme example concerns the highly conserved homeobox gene Pdx1, which is uniquely divergent and GC rich in the sand rat Psammomys obesus and close relatives. Here, we investigate the antagonistic interplay between very rare amino acid changes driven by GC skew and the force of natural selection. Using ectopic protein expression in cell culture, pulse-chase labeling, in vitro mutagenesis, and drug treatment, we compare properties of mouse and sand rat Pdx1 proteins. We find that amino acid change driven by GC skew resulted in altered protein stability, with a significantly longer protein half-life for sand rat Pdx1. Using a reversible inhibitor of the 26S proteasome, MG132, we find that sand rat and mouse Pdx1 are both degraded through the ubiquitin proteasome pathway. However, in vitro mutagenesis reveals this pathway operates through different amino acid residues. We propose that GC skew caused loss of a key ubiquitination site, conserved through vertebrate evolution, and that sand rat Pdx1 evolved or fixed a new ubiquitination site to compensate. Our results give molecular insight into the power of natural selection in the face of maladaptive changes driven by strong GC skew.

Intrinsic DNA Shape Accounts for Affinity Differences between Hox-Cofactor Binding Sites

Transcription factors bind to their binding sites over a wide range of affinities, yet how differences in affinity are encoded in DNA sequences is not well understood. Here, we report X-ray crystal structures of four heterodimers of the Hox protein AbdominalB bound with its cofactor Extradenticle to four target DNA molecules that differ in affinity by up to ∼20-fold. Remarkably, despite large differences in affinity, the overall structures are very similar in all four complexes. In contrast, the predicted shapes of the DNA binding sites (i.e., the intrinsic DNA shape) in the absence of bound protein are strikingly different from each other and correlate with affinity: binding sites that must change conformations upon protein binding have lower affinities than binding sites that have more optimal conformations prior to binding. Together, these observations suggest that intrinsic differences in DNA shape provide a robust mechanism for modulating affinity without affecting other protein-DNA interactions.

Deregulated NKL Homeobox Genes in B-Cell Lymphoma

Recently, we have described physiological expression patterns of NKL homeobox genes in early hematopoiesis and in subsequent lymphopoiesis. We identified nine genes which constitute the so-called NKL-code. Aberrant overexpression of code-members or ectopically activated non-code NKL homeobox genes are described in T-cell leukemia and in T- and B-cell lymphoma, highlighting their oncogenic role in lymphoid malignancies. Here, we introduce the NKL-code in normal hematopoiesis and focus on deregulated NKL homeobox genes in B-cell lymphoma, including HLX, MSX1 and NKX2-2 in Hodgkin lymphoma; HLX, NKX2-1 and NKX6-3 in diffuse large B-cell lymphoma; and NKX2-3 in splenic marginal zone lymphoma. Thus, the roles of various members of the NKL homeobox gene subclass are considered in normal and pathological hematopoiesis in detail.

Development of a Molecular Marker Linked to the A4 Locus and the Structure of HD Genes in Pleurotus eryngii

Allelic differences in A and B mating-type loci are a prerequisite for the progression of mating in the genus Pleurotus eryngii; thus, the crossing is hampered by this biological barrier in inbreeding. Molecular markers linked to mating types of P. eryngii KNR2312 were investigated with randomly amplified polymorphic DNA to enhance crossing efficiency. An A4-linked sequence was identified and used to find the adjacent genomic region with the entire motif of the A locus from a contig sequenced by PacBio. The sequence-characterized amplified region marker 7-2₂₉₉ distinguished A4 mating-type monokaryons from KNR2312 and other strains. A BLAST search of flanked sequences revealed that the A4 locus had a general feature consisting of the putative HD1 and HD2 genes. Both putative HD transcription factors contain a homeodomain sequence and a nuclear localization sequence; however, valid dimerization motifs were found only in the HD1 protein. The ACAAT motif, which was reported to have relevance to sex determination, was found in the intergenic region. The SCAR marker could be applicable in the classification of mating types in the P. eryngii breeding program, and the A4 locus could be the basis for a multi-allele detection marker.

Drosophila Homeodomain-Interacting Protein Kinase (Hipk) Phosphorylates the Homeodomain Proteins Homeobrain, Empty Spiracles, and Muscle Segment Homeobox

The Drosophila homeodomain-interacting protein kinase (Hipk) is the fly representative of the well-conserved group of HIPKs in vertebrates. It was initially found through its characteristic interactions with homeodomain proteins. Hipk is involved in a variety of important developmental processes, such as the development of the eye or the nervous system. In the present study, we set Hipk and the Drosophila homeodomain proteins Homeobrain (Hbn), Empty spiracles (Ems), and Muscle segment homeobox (Msh) in an enzyme-substrate relationship. These homeoproteins are transcription factors that function during Drosophila neurogenesis and are, at least in part, conserved in vertebrates. We reveal a physical interaction between Hipk and the three homeodomain proteins in vivo using bimolecular fluorescence complementation (BiFC). In the course of in vitro phosphorylation analysis and subsequent mutational analysis we mapped several Hipk phosphorylation sites of Hbn, Ems, and Msh. The phosphorylation of Hbn, Ems, and Msh may provide further insight into the function of Hipk during development of the Drosophila nervous system.

The POU/Oct Transcription Factor Nubbin Controls the Balance of Intestinal Stem Cell Maintenance and Differentiation by Isoform-Specific Regulation

Drosophila POU/Oct transcription factors are required for many developmental processes, but their putative regulation of adult stem cell activity has not been investigated. Here, we show that Nubbin (Nub)/Pdm1, homologous to mammalian OCT1/POU2F1 and related to OCT4/POU5F1, is expressed in gut epithelium progenitor cells. We demonstrate that the nub-encoded protein isoforms, Nub-PB and Nub-PD, play opposite roles in the regulation of intestinal stem cell (ISC) maintenance and differentiation. Depletion of Nub-PB in progenitor cells increased ISC proliferation by derepression of escargot expression. Conversely, loss of Nub-PD reduced ISC proliferation, suggesting that this isoform is necessary for ISC maintenance, analogous to mammalian OCT4/POU5F1 functions. Furthermore, Nub-PB is required in enteroblasts to promote differentiation, and it acts as a tumor suppressor of Notch RNAi-driven hyperplasia. We suggest that a dynamic and well-tuned expression of Nub isoforms in progenitor cells is required for maintaining gut epithelium homeostasis.

•

Drosophila nubbin (nub) is expressed in adult midgut progenitor cells

•

The Nub-PB isoform drives differentiation and acts as a tumor suppressor

•

The Nub-PD isoform maintains intestinal stem cell proliferation

•

Nub-PD and Nub-PB regulate stem cell proliferation antagonistically

DIVARICATA AND RADIALIS INTERACTING FACTOR (DRIF) also interacts with WOX and KNOX proteins associated with wood formation in Populus trichocarpa

DIVARICATA AND RADIALIS INTERACTING FACTOR (DRIF) from snapdragon (Antirrhinum majus) is a MYB/SANT protein that interacts with related MYB/SANT proteins, RADIALIS and DIVARICATA, through its N-terminal MYB/SANT domain. In addition to the MYB/SANT domain, DRIF contains a C-terminal domain of unknown function (DUF3755). Here we describe novel protein-protein interactions involving a poplar (Populus trichocarpa) homolog of DRIF, PtrDRIF1. In addition to interacting with poplar homologs of RADIALIS (PtrRAD1) and DIVARICATA (PtrDIV4), PtrDRIF1 interacted with members of other families within the homeodomain-like superfamily, including PtrWOX13c, a WUSCHEL-RELATED HOMEOBOX protein, and PtrKNAT7, a KNOTTED1-LIKE HOMEOBOX protein. PtrRAD1 and PtrDIV4 interacted with the MYB/SANT-containing N-terminal portion of PtrDRIF1, whereas DUF3755 was both necessary and sufficient for interactions with PtrWOX13c and PtrKNAT7. Of the two MYB/SANT domains present in PtrDIV4, only the N-terminal MYB/SANT domain interacted with PtrDRIF1. GFP-PtrDRIF1 expressed alone or with PtrRAD1 localized to the cytoplasm, whereas co-expression of GFP-PtrDRIF1 with PtrDIV4, PtrWOX13c or PtrKNAT7 resulted in nuclear localization of GFP-PtrDRIF1. Modified yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) experiments using PtrDRIF1 as a bridge protein revealed that PtrDRIF1 simultaneously interacted with PtrRAD1 and PtrWOX13c, but could not form a heterotrimeric complex when PtrDIV4 was substituted for PtrRAD1. Moreover, a Y2H competition assay indicated that PtrKNAT7 inhibits the interaction between PtrRAD1 and PtrDRIF1. The discovery of an additional protein-protein interaction domain in DRIF proteins, DUF3755, and its ability to form heterodimers and heterotrimers involving MYB/SANT and wood-associated homeodomain proteins, implicates DRIF proteins as mediators of a broader array of processes than previously reported.

Changing expressions: a hypothesis for the origin of the vascular plant life cycle

Plant life cycles underwent fundamental changes during the initial colonization of the land in the Early Palaeozoic, shaping the direction of evolution. Fossils reveal unanticipated diversity, including new variants of meiotic cell division and leafless gametophytes with mycorrhizal-like symbioses, rhizoids, vascular tissues and stomata. Exceptional fossils from the 407-Ma Rhynie chert (Scotland) play a key role in unlocking this diversity. These fossils are reviewed against progress in our understanding of the plant tree of life and recent advances in developmental genetics. Combining data from different sources sheds light on a switch in life cycle that gave rise to the vascular plants. One crucial step was the establishment of a free-living sporophyte from one that was an obligate matrotroph borne on the gametophyte. It is proposed that this difficult evolutionary transition was achieved through expansion of gene expression primarily from the gametophyte to the sporophyte, establishing a now extinct life cycle variant that was more isomorphic than heteromorphic. These changes also linked for the first time in one developmental system rhizoids, vascular tissues and stomata, putting in place the critical components that regulate transpiration and forming a physiological platform of primary importance to the diversification of vascular plants.This article is part of a discussion meeting issue 'The Rhynie cherts: our earliest terrestrial ecosystem revisited'.

A TALE of shrimps: Genome-wide survey of homeobox genes in 120 species from diverse crustacean taxa

The homeodomain-containing proteins are an important group of transcription factors found in most eukaryotes including animals, plants and fungi. Homeobox genes are responsible for a wide range of critical developmental and physiological processes, ranging from embryonic development, innate immune homeostasis to whole-body regeneration. With continued fascination on this key class of proteins by developmental and evolutionary biologists, multiple efforts have thus far focused on the identification and characterization of homeobox orthologs from key model organisms in attempts to infer their evolutionary origin and how this underpins the evolution of complex body plans. Despite their importance, the genetic complement of homeobox genes has yet been described in one of the most valuable groups of animals representing economically important food crops. With crustacean aquaculture being a growing industry worldwide, it is clear that systematic and cross-species identification of crustacean homeobox orthologs is necessary in order to harness this genetic circuitry for the improvement of aquaculture sustainability. Using publicly available transcriptome data sets, we identified a total of 4183 putative homeobox genes from 120 crustacean species that include food crop species, such as lobsters, shrimps, crayfish and crabs. Additionally, we identified 717 homeobox orthologs from 6 other non-crustacean arthropods, which include the scorpion, deer tick, mosquitoes and centipede. This high confidence set of homeobox genes will now serve as a key resource to the broader community for future functional and comparative genomics studies.

Islet Coordinately Regulates Motor Axon Guidance and Dendrite Targeting through the Frazzled/DCC Receptor

Motor neuron axon targeting in the periphery is correlated with the positions of motor neuron inputs in the CNS, but how these processes are coordinated to form a myotopic map remains poorly understood. We show that the LIM homeodomain factor Islet (Isl) controls targeting of both axons and dendrites in Drosophila motor neurons through regulation of the Frazzled (Fra)/DCC receptor. Isl is required for fra expression in ventrally projecting motor neurons, and isl and fra mutants have similar axon guidance defects. Single-cell labeling indicates that isl and fra are also required for dendrite targeting in a subset of motor neurons. Finally, overexpression of Fra rescues axon and dendrite targeting defects in isl mutants. These results indicate that Fra acts downstream of Isl in both the periphery and the CNS, demonstrating how a single regulatory relationship is used in multiple cellular compartments to coordinate neural circuit wiring.

Hox functional diversity: Novel insights from flexible motif folding and plastic protein interaction

How the formidable diversity of forms emerges from developmental and evolutionary processes is one of the most fascinating questions in biology. The homeodomain-containing Hox proteins were recognized early on as major actors in diversifying animal body plans. The molecular mechanisms underlying how this transcription factor family controls a large array of context- and cell-specific biological functions is, however, still poorly understood. Clues to functional diversity have emerged from studies exploring how Hox protein activity is controlled through interactions with PBC class proteins, also evolutionary conserved HD-containing proteins. Recent structural data and molecular dynamic simulations add further mechanistic insights into Hox protein mode of action, suggesting that flexible folding of protein motifs allows for plastic protein interaction. As we discuss in this review, these findings define a novel type of Hox-PBC interaction, weak and dynamic instead of strong and static, hence providing novel clues to understanding Hox transcriptional specificity and diversity.

Advances in Understanding Mating Type Gene Organization in the Mushroom-Forming Fungus Flammulina velutipes

The initiation of sexual development in the important edible and medicinal mushroom Flammulina velutipes is controlled by special genes at two different, independent, mating type (MAT) loci: HD and PR. We expanded our understanding of the F. velutipes mating type system by analyzing the MAT loci from a series of strains. The HD locus of F. velutipes houses homeodomain genes (Hd genes) on two separated locations: sublocus HD-a and HD-b. The HD-b subloci contained strain-specific Hd1/Hd2 gene pairs, and crosses between strains with different HD-b subloci indicated a role in mating. The function of the HD-a sublocus remained undecided. Many, but not all strains contained the same conserved Hd2 gene at the HD-a sublocus. The HD locus usually segregated as a whole, though we did detect one new HD locus with a HD-a sublocus from one parental strain, and a HD-b sublocus from the other. The PR locus of F. velutipes contained pheromone receptor (STE3) and pheromone precursor (Pp) genes at two locations, sublocus PR-a and PR-b. PR-a and PR-b both contained sets of strain-specific STE3 and Pp genes, indicating a role in mating. PR-a and PR-b cosegregated in our experiments. However, the identification of additional strains with identical PR-a, yet different PR-b subloci, demonstrated that PR subloci can recombine within the PR locus. In conclusion, at least three of the four MAT subloci seem to participate in mating, and new HD and PR loci can be generated through intralocus recombination in F. velutipes.

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.

Regulation of Silk Genes by Hox and Homeodomain Proteins in the Terminal Differentiated Silk Gland of the Silkworm Bombyx mori

The silk gland of the silkworm Bombyx mori is a long tubular organ that is divided into several subparts along its anteroposterior (AP) axis. As a trait of terminal differentiation of the silk gland, several silk protein genes are expressed with unique regional specificities. Most of the Hox and some of the homeobox genes are also expressed in the differentiated silk gland with regional specificities. The expression patterns of Hox genes in the silk gland roughly correspond to those in embryogenesis showing "colinearity". The central Hox class protein Antennapedia (Antp) directly regulates the expression of several middle silk gland-specific silk genes, whereas the Lin-1/Isl-1/Mec3 (LIM)-homeodomain transcriptional factor Arrowhead (Awh) regulates the expression of posterior silk gland-specific genes for silk fiber proteins. We summarize our results and discuss the usefulness of the silk gland of Bombyx mori for analyzing the function of Hox genes. Further analyses of the regulatory mechanisms underlying the region-specific expression of silk genes will provide novel insights into the molecular bases for target-gene selection and regulation by Hox and homeodomain proteins.

Postmitotic diversification of olfactory neuron types is mediated by differential activities of the HMG-box transcription factor SOX-2

Diversification of neuron classes is essential for functions of the olfactory system, but the underlying mechanisms that generate individual olfactory neuron types are only beginning to be understood. Here we describe a role of the highly conserved HMG-box transcription factor SOX-2 in postmitotic specification and alternative differentiation of the Caenorhabditis elegans AWC and AWB olfactory neurons. We show that SOX-2 partners with different transcription factors to diversify postmitotic olfactory cell types. SOX-2 functions cooperatively with the OTX/OTD transcription factor CEH-36 to specify an AWC "ground state," and functions with the LIM homeodomain factor LIM-4 to suppress this ground state and drive an AWB identity instead. Our findings provide novel insights into combinatorial codes that drive terminal differentiation programs in the nervous system and reveal a biological function of the deeply conserved Sox2 protein that goes beyond its well-known role in stem cell biology.

Loading of PAX3 to Mitotic Chromosomes Is Mediated by Arginine Methylation and Associated with Waardenburg Syndrome

PAX3 is a transcription factor critical to gene regulation in mammalian development. Mutations in PAX3 are associated with Waardenburg syndrome (WS), but the mechanism of how mutant PAX3 proteins cause WS remains unclear. Here, we found that PAX3 loads on mitotic chromosomes using its homeodomain. PAX3 WS mutants with mutations in homeodomain lose the ability to bind mitotic chromosomes. Moreover, loading of PAX3 on mitotic chromosomes requires arginine methylation, which is regulated by methyltransferase PRMT5 and demethylase JMJD6. Mutant PAX3 proteins that lose mitotic chromosome localization block cell proliferation and normal development of zebrafish. These results reveal the molecular mechanism of PAX3s loading on mitotic chromosomes and the importance of this localization pattern in normal development. Our findings suggest that PAX3 WS mutants interfere with the normal functions of PAX3 in a dominant negative manner, which is important to the understanding of the pathogenesis of Waardenburg syndrome.

Controlled expression of Drosophila homeobox loci using the Hostile takeover system

BACKGROUND

Hostile takeover (Hto) is a Drosophila protein trapping system that allows the investigator to both induce a gene and tag its product. The Hto transposon carries a GAL4-regulated promoter expressing an exon encoding a FLAG-mCherry tag. Upon expression, the Hto exon can splice to a downstream genomic exon, generating a fusion transcript and tagged protein.

RESULTS

Using rough-eye phenotypic screens, Hto inserts were recovered at eight homeobox or Pax loci: cut, Drgx/CG34340, Pox neuro, araucan, shaven/D-Pax2, Zn finger homeodomain 2, Sex combs reduced (Scr), and the abdominal-A region. The collection yields diverse misexpression phenotypes. Ectopic Drgx was found to alter the cytoskeleton and cell adhesion in ovary follicle cells. Hto expression of cut, araucan, or shaven gives phenotypes similar to those of the corresponding UAS-cDNA constructs. The cut and Pox neuro phenotypes are suppressed by the corresponding RNAi constructs. The Scr and abdominal-A inserts do not make fusion proteins, but may act by chromatin- or RNA-based mechanisms.

CONCLUSIONS

Hto can effectively express tagged homeodomain proteins from their endogenous loci; the Minos vector allows inserts to be obtained even in transposon cold-spots. Hto screens may recover homeobox genes at high rates because they are particularly sensitive to misexpression.

Reinforcing the egg-timer: recruitment of novel lophotrochozoa homeobox genes to early and late development in the pacific oyster

The metazoan superclade Lophotrochozoa includes mollusks, annelids, and several other animal phyla. It is reasonable to assume that this organismal diversity may be traced, in part, to changes in developmentally important genes, such as the homeobox genes. Although most comparative studies have focussed on ancient homeobox gene families conserved across bilaterians, there are also "novel" homeobox genes that have arisen more recently in evolution, presumably by duplication followed by radical divergence and functional change. We classify 136 homeobox genes in the genome sequence of the Pacific oyster, Crassostrea gigas. The genome shows an unusually low degree of homeobox gene clustering, with disruption of the NK, Hox, and ParaHox gene clusters. Among the oyster genes, 31 do not fall into ancient metazoan or bilaterian homeobox gene families; we deduce that they originated in the lophotrochozoan clade. We compared eight lophotrochozoan genomes to trace the pattern of homeobox gene evolution across this clade, allowing us to define 19 new lophotrochozoan-specific clades within the ANTP, PRD, TALE, ZF, SIX, and CUT classes. Using transcriptome data, we compared temporal expression of each homeobox gene in oyster development, and discovered that the lophotrochozoan-specific homeobox genes have peak expression either in early development (egg to gastrula) or in late development (after the trochophore larval stage), but rarely in between. This finding is consistent with the egg-timer, hourglass or phylotypic stage model of developmental evolution, in which there is a conserved central phase of development, but more evolutionarily labile early and late phases.

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain-DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1-DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein-DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein-DNA complexes. The order of stability of protein-DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein-DNA complexes. Among specific protein-DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.

Meis: New friends of Pax

The generation of neuronal diversity in the mammalian brain is a multistep process, beginning with the regional patterning of neural stem- and progenitor cell domains, the commitment of these cells toward a general neuronal fate, followed by the selection of a particular neuronal subtype and the differentiation of postmitotic neurons. Each of these steps as well as the transitions between them require precisely controlled changes in transcriptional programs. Although a large number of transcription factors are known to regulate neurogenesis in the embryonic and adult central nervous system, the sheer number of neuronal cell types in the brain and the complexity of the cellular processes that accompany their production suggest that transcription factors act cooperatively to control individual steps in neurogenesis. In fact, combinatorial regulation by sets of transcription factors has emerged as a versatile mode to control cell fate specification. Here, I discuss our recent finding that members of the MEIS-subfamily of TALE-transcription factors, originally identified as HOX cofactors in non-neural tissues, function in concert with PAX-proteins in the regulation of cell fate specification and neuronal differentiation in the embryonic and adult brain.