Evolutionary conservation of otd/Otx2 transcription factor action: a genome-wide microarray analysis in Drosophila

O-GlcNAcylation regulates OTX2's proteostasis

O-GlcNAcylation is a key post-translational modification, playing a vital role in cell signaling during development, especially in the brain. In this study, we investigated the role of O-GlcNAcylation in regulating the homeobox protein OTX2, which contributes to various brain disorders, such as combined pituitary hormone deficiency, retinopathy, and medulloblastoma. Our research demonstrated that, under normal physiological conditions, the proteasome plays a pivotal role in breaking down endogenous OTX2. However, when the levels of OTX2 rise, it forms oligomers and/or aggregates that require macroautophagy for clearance. Intriguingly, we demonstrated that O-GlcNAcylation enhances the solubility of OTX2, thereby limiting the formation of these aggregates. Additionally, we unveiled an interaction between OTX2 and the chaperone protein CCT5 at the O-GlcNAc sites, suggesting a potential collaborative role in preventing OTX2 aggregation. Finally, our study demonstrated that while OTX2 physiologically promotes cell proliferation, an O-GlcNAc-depleted OTX2 is detrimental to cancer cells.

Cell-specific occupancy dynamics between the pioneer-like factor Opa/ZIC and Ocelliless/OTX regulate early head development in embryos

During development, embryonic patterning systems direct a set of initially uncommitted pluripotent cells to differentiate into a variety of cell types and tissues. A core network of transcription factors, such as Zelda/POU5F1, Odd-paired (Opa)/ZIC3 and Ocelliless (Oc)/OTX2, are conserved across animals. While Opa is essential for a second wave of zygotic activation after Zelda, it is unclear whether Opa drives head cell specification, in the Drosophila embryo. Our hypothesis is that Opa and Oc are interacting with distinct cis-regulatory regions for shaping cell fates in the embryonic head. Super-resolution microscopy and meta-analysis of single-cell RNAseq datasets show that opa’s and oc’s overlapping expression domains are dynamic in the head region, with both factors being simultaneously transcribed at the blastula stage. Additionally, analysis of single-embryo RNAseq data reveals a subgroup of Opa-bound genes to be Opa-independent in the cellularized embryo. Interrogation of these genes against Oc ChIPseq combined with in situ data, suggests that Opa is competing with Oc for the regulation of a subgroup of genes later in gastrulation. Specifically, we find that Oc binds to late, head-specific enhancers independently and activates them in a head-specific wave of zygotic transcription, suggesting distinct roles for Oc in the blastula and gastrula stages.

Grainyhead-like Transcription Factors in Craniofacial Development

Craniofacial development in vertebrates involves the coordinated growth, migration, and fusion of several facial prominences during embryogenesis, processes governed by strict genetic and molecular controls. A failure in any of the precise spatiotemporal sequences of events leading to prominence fusion often leads to anomalous facial, skull, and jaw formation-conditions termed craniofacial defects (CFDs). Affecting approximately 0.1% to 0.3% of live births, CFDs are a highly heterogeneous class of developmental anomalies, which are often underpinned by genetic mutations. Therefore, identifying novel disease-causing mutations in genes that regulate craniofacial development is a critical prerequisite to develop new preventive or therapeutic measures. The Grainyhead-like ( GRHL) transcription factors are one such gene family, performing evolutionarily conserved roles in craniofacial patterning. The antecedent member of this family, Drosophila grainyhead ( grh), is required for head skeleton development in fruit flies, loss or mutation of Grhl family members in mouse and zebrafish models leads to defects of both maxilla and mandible, and recently, mutations in human GRHL3 have been shown to cause or contribute to both syndromic (Van Der Woude syndrome) and nonsyndromic palatal clefts. In this review, we summarize the current knowledge regarding the craniofacial-specific function of the Grainyhead-like family in multiple model species, identify some of the major target genes regulated by the Grhl transcription factors in craniofacial patterning, and, by examining animal models, draw inferences as to how these data will inform the likely roles of GRHL factors in human CFDs comprising palatal clefting. By understanding the molecular networks regulated by Grhl2 and Grhl3 target genes in other systems, we can propose likely pathways that mediate the effects of these transcription factors in human palatogenesis.

Specific expression pattern of a novel Otx2 splicing variant during neural differentiation

We cloned a new splicing variant of Otx2 gene, Otx2c. Otx2c lacks entire exon 4, most of the region encoding the homeodomain. More importantly, Otx2c harbors an early premature stop codon and bioinformatics analysis prefers it to be a non-protein coding RNA. In addition, this splicing variant is not simply a noise during mRNA processing, since it is mainly expressed in undifferentiated human embryonic stem cells but gradually decreased during differentiation. Therefore, we report here that a single pre-mRNA can generate both coding and non-coding RNAs through alternative splicing and this splicing activity is tightly regulated in different cell contexts.

A proline repeat polymorphism of the Frost gene of Drosophila melanogaster showing clinal variation but not associated with cold resistance

Genetic polymorphisms underlying adaptive shifts in thermal responses are poorly known even though studies are providing a detailed understanding of these responses at the cellular and physiological levels. The Frost gene of Drosophila melanogaster is a prime candidate for thermal adaptation; it is up-regulated under cold stress and knockdown of this gene influences cold resistance. Here we describe an amino-acid INDEL polymorphism in proline repeat number in the structural component of this gene. The two main repeats, accounting for more than 90% of alleles in eastern Australia, show a strong clinal pattern; the 6P allele was at a high frequency in tropical locations, and the 10P allele was common in temperate populations. However, the frequency of these alleles was not associated with three different assays of cold resistance. Adult transcription level of Frost was also unrelated to cold resistance as measured through post chill coma mobility. The functional significance of the proline repeat polymorphism therefore remains unclear despite its clinal pattern. The data also demonstrate the feasibility of using Roche/454 sequencing for establishing clinal patterns.

The role of eyg Pax gene in the development of the head vertex in Drosophila

The Drosophila head vertex is composed of three ocelli, stereotypic bristle patterns and characteristic cuticles. It is derived from the fusion of two eye-antenna discs. The head vertex primordium is located at the anterior-dorsal region of the eye disc. The orthodenticle (otd) homeobox gene is expressed in the primordium and is functionally required for its development and patterning. Here we show that the Pax gene eye gone (eyg) is expressed adjacent to the otd expression domain in the eye disc. otd is required and sufficient to repress eyg transcription, thereby preventing eyg from expressing in the head vertex primordium. In otd mutant, eyg expression is derepressed in the head vertex primordium and is a major negative effector to block head vertex development. Therefore, otd not only needs to induce downstream effector genes to execute the development and patterning of the head vertex development, but also needs to actively repress the negative regulator eyg. In addition, eyg is required for the development of the lateral bristles in the head vertex. So eyg plays both positive and negative roles in head vertex development.

Use with caution: developmental systems divergence and potential pitfalls of animal models

Transgenic animal models have played an important role in elucidating gene functions and the molecular basis development, physiology, behavior, and pathogenesis. Transgenic models have been so successful that they have become a standard tool in molecular genetics and biomedical studies and are being used to fulfill one of the main goals of the post-genomic era: to assign functions to each gene in the genome. However, the assumption that gene functions and genetic systems are conserved between models and humans is taken for granted, often in spite of evidence that gene functions and networks diverge during evolution. In this review, I discuss some mechanisms that generate functional divergence and highlight recent examples demonstrating that gene functions and regulatory networks diverge through time. These examples suggest that annotation of gene functions based solely on mutant phenotypes in animal models, as well as assumptions of conserved functions between species, can be wrong. Therefore, animal models of gene function and human disease may not provide appropriate information, particularly for rapidly evolving genes and systems.

Conserved chromosomal clustering of genes governed by chromatin regulators in Drosophila

BACKGROUND

The trithorax group (trxG) and Polycomb group (PcG) proteins are responsible for the maintenance of stable transcriptional patterns of many developmental regulators. They bind to specific regions of DNA and direct the post-translational modifications of histones, playing a role in the dynamics of chromatin structure.

RESULTS

We have performed genome-wide expression studies of trx and ash2 mutants in Drosophila melanogaster. Using computational analysis of our microarray data, we have identified 25 clusters of genes potentially regulated by TRX. Most of these clusters consist of genes that encode structural proteins involved in cuticle formation. This organization appears to be a distinctive feature of the regulatory networks of TRX and other chromatin regulators, since we have observed the same arrangement in clusters after experiments performed with ASH2, as well as in experiments performed by others with NURF, dMyc, and ASH1. We have also found many of these clusters to be significantly conserved in D. simulans, D. yakuba, D. pseudoobscura and partially in Anopheles gambiae.

CONCLUSION

The analysis of genes governed by chromatin regulators has led to the identification of clusters of functionally related genes conserved in other insect species, suggesting this chromosomal organization is biologically important. Moreover, our results indicate that TRX and other chromatin regulators may act globally on chromatin domains that contain transcriptionally co-regulated genes.

Resurrecting the role of transcription factor change in developmental evolution

A long-standing question in evolutionary and developmental biology concerns the relative contribution of cis-regulatory and protein changes to developmental evolution. Central to this argument is which mutations generate evolutionarily relevant phenotypic variation? A review of the growing body of evolutionary and developmental literature supports the notion that many developmentally relevant differences occur in the cis-regulatory regions of protein-coding genes, generally to the exclusion of changes in the protein-coding region of genes. However, accumulating experimental evidence demonstrates that many of the arguments against a role for proteins in the evolution of gene regulation, and the developmental evolution in general, are no longer supported and there is an increasing number of cases in which transcription factor protein changes have been demonstrated in evolution. Here, we review the evidence that cis-regulatory evolution is an important driver of phenotypic evolution and provide examples of protein-mediated developmental evolution. Finally, we present an argument that the evolution of proteins may play a more substantial, but thus far underestimated, role in developmental evolution.

Analysis of the Otd-dependent transcriptome supports the evolutionary conservation of CRX/OTX/OTD functions in flies and vertebrates

Homeobox transcription factors of the vertebrate CRX/OTX family play critical roles in photoreceptor neurons, the rostral brain and circadian processes. In mouse, the three related proteins, CRX, OTX1, and OTX2, fulfill these functions. In Drosophila, the single founding member of this gene family, called orthodenticle (otd), is required during embryonic brain and photoreceptor neuron development. We have used global gene expression analysis in late pupal heads to better characterize the post-embryonic functions of Otd in Drosophila. We have identified 61 genes that are differentially expressed between wild type and a viable eye-specific otd mutant allele. Among them, about one-third represent potentially direct targets of Otd based on their association with evolutionarily conserved Otd-binding sequences. The spectrum of biological functions associated with these gene targets establishes Otd as a critical regulator of photoreceptor morphology and phototransduction, as well as suggests its involvement in circadian processes. Together with the well-documented role of otd in embryonic patterning, this evidence shows that vertebrate and fly genes contribute to analogous biological processes, notwithstanding the significant divergence of the underlying genetic pathways. Our findings underscore the common evolutionary history of photoperception-based functions in vertebrates and invertebrates and support the view that a complex nervous system was already present in the last common ancestor of all bilateria.

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

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

De novo discovery of a tissue-specific gene regulatory module in a chordate

We engage the experimental and computational challenges of de novo regulatory module discovery in a complex and largely unstudied metazoan genome. Our analysis is based on the comprehensive characterization of regulatory elements of 20 muscle genes in the chordate, Ciona savignyi. Three independent types of data we generate contribute to the characterization of a muscle-specific regulatory module: (1) Positive elements (PEs), short sequences sufficient for strong muscle expression that are identified in a high-resolution in vivo analysis; (2) CisModules (CMs), candidate regulatory modules defined by clusters of overrepresented motifs predicted de novo; and (3) Conserved elements (CEs), short noncoding sequences of strong conservation between C. savignyi and C. intestinalis. We estimate the accuracy of the computational predictions by an analysis of the intersection of these data. As final biological validation of the discovered muscle regulatory module, we implement a novel algorithm to search the genome for instances of the module and identify seven novel enhancers.

Regulation of chemosensory and GABAergic motor neuron development by the C. elegans Aristaless/Arx homolog alr-1

Mutations in the highly conserved Aristaless-related homeodomain protein ARX have been shown to underlie multiple forms of X-linked mental retardation. Arx knockout mice exhibit thinner cerebral cortices because of decreased neural precursor proliferation, and also exhibit defects in the differentiation and migration of GABAergic interneurons. However, the role of ARX in the observed behavioral and developmental abnormalities is unclear. The regulatory functions of individual homeodomain proteins and the networks in which they act are frequently highly conserved across species, although these networks may be deployed in different developmental contexts. In Drosophila, aristaless mutants exhibit defects in the development of terminal appendages, and Aristaless has been shown to function with the LIM-homeodomain protein LIM1 to regulate leg development. Here, we describe the role of the Aristaless/Arx homolog alr-1 in C. elegans. We show that alr-1 acts in a pathway with the LIM1 ortholog lin-11 to regulate the development of a subset of chemosensory neurons. Moreover, we demonstrate that the differentiation of a GABAergic motoneuron subtype is affected in alr-1 mutants, suggesting parallels with ARX functions in vertebrates. Investigating ALR-1 functions in C. elegans may yield insights into the role of this important protein in neuronal development and the etiology of mental retardation.

Intragenic spatial patterns of codon usage bias in prokaryotic and eukaryotic genomes

To study the roles of translational accuracy, translational efficiency, and the Hill-Robertson effect in codon usage bias, we studied the intragenic spatial distribution of synonymous codon usage bias in four prokaryotic (Escherichia coli, Bacillus subtilis, Sulfolobus tokodaii, and Thermotoga maritima) and two eukaryotic (Saccharomyces cerevisiae and Drosophila melanogaster) genomes. We generated supersequences at each codon position across genes in a genome and computed the overall bias at each codon position. By quantitatively evaluating the trend of spatial patterns using isotonic regression, we show that in yeast and prokaryotic genomes, codon usage bias increases along translational direction, which is consistent with purifying selection against nonsense errors. Fruit fly genes show a nearly symmetric M-shaped spatial pattern of codon usage bias, with less bias in the middle and both ends. The low codon usage bias in the middle region is best explained by interference (the Hill-Robertson effect) between selections at different codon positions. In both yeast and fruit fly, spatial patterns of codon usage bias are characteristically different from patterns of GC-content variations. Effect of expression level on the strength of codon usage bias is more conspicuous than its effect on the shape of the spatial distribution.

Transcriptional signature of an adult brain tumor in Drosophila

Background

Mutations and gene expression alterations in brain tumors have been extensively investigated, however the causes of brain tumorigenesis are largely unknown. Animal models are necessary to correlate altered transcriptional activity and tumor phenotype and to better understand how these alterations cause malignant growth. In order to gain insights into the in vivo transcriptional activity associated with a brain tumor, we carried out genome-wide microarray expression analyses of an adult brain tumor in Drosophila caused by homozygous mutation in the tumor suppressor gene brain tumor (brat).

Results

Two independent genome-wide gene expression studies using two different oligonucleotide microarray platforms were used to compare the transcriptome of adult wildtype flies with mutants displaying the adult brat^k06028 mutant brain tumor. Cross-validation and stringent statistical criteria identified a core transcriptional signature of brat^k06028 neoplastic tissue. We find significant expression level changes for 321 annotated genes associated with the adult neoplastic brat^k06028 tissue indicating elevated and aberrant metabolic and cell cycle activity, upregulation of the basal transcriptional machinery, as well as elevated and aberrant activity of ribosome synthesis and translation control. One fifth of these genes show homology to known mammalian genes involved in cancer formation.

Conclusion

Our results identify for the first time the genome-wide transcriptional alterations associated with an adult brain tumor in Drosophila and reveal insights into the possible mechanisms of tumor formation caused by homozygous mutation of the translational repressor brat.

Inter-species differences of co-expression of neighboring genes in eukaryotic genomes

BACKGROUND

There is increasing evidence that gene order within the eukaryotic genome is not random. In yeast and worm, adjacent or neighboring genes tend to be co-expressed. Clustering of co-expressed genes has been found in humans, worm and fruit flies. However, in mice and rats, an effect of chromosomal distance (CD) on co-expression has not been investigated yet. Also, no cross-species comparison has been made so far. We analyzed the effect of CD as well as normalized distance (ND) using expression data in six eukaryotic species: yeast, fruit fly, worm, rat, mouse and human.

RESULTS

We analyzed 24 sets of expression data from the six species. Highly co-expressed pairs were sorted into bins of equal sized intervals of CD, and a co-expression rate (CoER) in each bin was calculated. In all datasets, a higher CoER was obtained in a short CD range than a long distance range. These results show that across all studied species, there was a consistent effect of CD on co-expression. However, the results using the ND show more diversity. Intra- and inter-species comparisons of CoER reveal that there are significant differences in the co-expression rates of neighboring genes among the species. A pair-wise BLAST analysis finds 8-30 % of the highly co-expressed pairs are duplicated genes.

CONCLUSION

We confirmed that in the six eukaryotic species, there was a consistent tendency that neighboring genes are likely to be co-expressed. Results of pair-wised BLAST indicate a significant effect of non-duplicated pairs on co-expression. A comparison of CD and ND suggests the dominant effect of CD.

The Drosophila melanogaster genome

Drosophila's importance as a model organism made it an obvious choice to be among the first genomes sequenced, and the Release 1 sequence of the euchromatic portion of the genome was published in March 2000. This accomplishment demonstrated that a whole genome shotgun (WGS) strategy could produce a reliable metazoan genome sequence. Despite the attention to sequencing methods, the nucleotide sequence is just the starting point for genome-wide analyses; at a minimum, the genome sequence must be interpreted using expressed sequence tag (EST) and complementary DNA (cDNA) evidence and computational tools to identify genes and predict the structures of their RNA and protein products. The functions of these products and the manner in which their expression and activities are controlled must then be assessed-a much more challenging task with no clear endpoint that requires a wide variety of experimental and computational methods. We first review the current state of the Drosophila melanogaster genome sequence and its structural annotation and then briefly summarize some promising approaches that are being taken to achieve an initial functional annotation.

Otx/otd Homeobox Genes Specify Distinct Sensory Neuron Identities in C. elegans

The mechanisms by which the diverse functional identities of neurons are generated are poorly understood. C. elegans responds to thermal and chemical stimuli using 12 types of sensory neurons. The Otx/otd homolog ttx-1 specifies the identities of the AFD thermosensory neurons. We show here that ceh-36 and ceh-37, the remaining two Otx-like genes in the C. elegans genome, specify the identities of AWC, ASE, and AWB chemosensory neurons, defining a role for this gene family in sensory neuron specification. All C. elegans Otx genes and rat Otx1 can substitute for ceh-37 and ceh-36, but only ceh-37 functionally substitutes for ttx-1. Functional substitution in the AWB neurons is mediated by activation of the same downstream target lim-4 by different Otx genes. Misexpression experiments indicate that although the specific identity adopted upon expression of an Otx gene may be constrained by the cellular context, individual Otx genes preferentially promote distinct neuronal identities.

Analysis of the eye developmental pathway in Drosophila using DNA microarrays

Pax-6 genes encode evolutionarily conserved transcription factors capable of activating the gene-expression program required to build an eye. When ectopically expressed in Drosophila imaginal discs, Pax-6 genes induce the eye formation on the corresponding appendages of the adult fly. We used two different Drosophila full-genome DNA microarrays to compare gene expression in wild-type leg discs versus leg discs where eyeless, one of the two Drosophila Pax-6 genes, was ectopically expressed. We validated these data by analyzing the endogenous expression of selected genes in eye discs and identified 371 genes that are expressed in the eye imaginal discs and up-regulated when an eye morphogenetic field is ectopically induced in the leg discs. These genes mainly encode transcription factors involved in photoreceptor specification, signal transducers, cell adhesion molecules, and proteins involved in cell division. As expected, genes already known to act downstream of eyeless during eye development were identified, together with a group of genes that were not yet associated with eye formation.

Gliogenesis inDrosophila: genome-wide analysis of downstream genes ofglial cells missingin the embryonic nervous system

In Drosophila, the glial cells missing (gcm) gene encodes a transcription factor that controls the determination of glial versus neuronal fate. In gcm mutants, presumptive glial cells are transformed into neurons and, conversely, when gcm is ectopically misexpressed, presumptive neurons become glia. Although gcm is thought to initiate glial cell development through its action on downstream genes that execute the glial differentiation program, little is known about the identity of these genes. To identify gcm downstream genes in a comprehensive manner, we used genome-wide oligonucleotide arrays to analyze differential gene expression in wild-type embryos versus embryos in which gcm is misexpressed throughout the neuroectoderm. Transcripts were analyzed at two defined temporal windows during embryogenesis. During the first period of initial gcm action on determination of glial cell precursors, over 400 genes were differentially regulated. Among these are numerous genes that encode other transcription factors, which underscores the master regulatory role of gcm in gliogenesis. During a second later period, when glial cells had already differentiated, over 1200 genes were differentially regulated. Most of these genes, including many genes for chromatin remodeling factors and cell cycle regulators, were not differentially expressed at the early stage, indicating that the genetic control of glial fate determination is largely different from that involved in maintenance of differentiated cells. At both stages, glial-specific genes were upregulated and neuron-specific genes were downregulated, supporting a model whereby gcm promotes glial development by activating glial genes, while simultaneously repressing neuronal genes. In addition, at both stages, numerous genes that were not previously known to be involved in glial development were differentially regulated and, thus, identified as potential new downstream targets of gcm. For a subset of the differentially regulated genes, tissue-specific in vivo expression data were obtained that confirmed the transcript profiling results. This first genome-wide analysis of gene expression events downstream of a key developmental transcription factor presents a novel level of insight into the repertoire of genes that initiate and maintain cell fate choices in CNS development.