Genome control is operated by transcription factors (TFs) controlling their target genes by binding to promoters and enhancers. Conceptually, the interactions between TFs, their binding sites, and their functional targets are represented by gene regulatory networks (GRNs). Deciphering in vivo GRNs underlying organ development in an unbiased genome-wide setting involves identifying both functional TF-gene interactions and physical TF-DNA interactions. To reverse engineer the GRNs of eye development in Drosophila, we performed RNA-seq across 72 genetic perturbations and sorted cell types and inferred a coexpression network. Next, we derived direct TF-DNA interactions using computational motif inference, ultimately connecting 241 TFs to 5,632 direct target genes through 24,926 enhancers. Using this network, we found network motifs, cis-regulatory codes, and regulators of eye development. We validate the predicted target regions of Grainyhead by ChIP-seq and identify this factor as a general cofactor in the eye network, being bound to thousands of nucleosome-free regions.

The functionality of sensory neurons is defined by the expression of specific sensory receptor genes. During the development of the Drosophila larval eye, photoreceptor neurons (PRs) make a binary choice to express either the blue-sensitive Rhodopsin 5 (Rh5) or the green-sensitive Rhodopsin 6 (Rh6). Later during metamorphosis, ecdysone signaling induces a cell fate and sensory receptor switch: Rh5-PRs are re-programmed to express Rh6 and become the eyelet, a small group of extraretinal PRs involved in circadian entrainment. However, the genetic and molecular mechanisms of how the binary cell fate decisions are made and switched remain poorly understood. We show that interplay of two transcription factors Senseless (Sens) and Hazy control cell fate decisions, terminal differentiation of the larval eye and its transformation into eyelet. During initial differentiation, a pulse of Sens expression in primary precursors regulates their differentiation into Rh5-PRs and repression of an alternative Rh6-cell fate. Later, during the transformation of the larval eye into the adult eyelet, Sens serves as an anti-apoptotic factor in Rh5-PRs, which helps in promoting survival of Rh5-PRs during metamorphosis and is subsequently required for Rh6 expression. Comparably, during PR differentiation Hazy functions in initiation and maintenance of rhodopsin expression. Hazy represses Sens specifically in the Rh6-PRs, allowing them to die during metamorphosis. Our findings show that the same transcription factors regulate diverse aspects of larval and adult PR development at different stages and in a context-dependent manner.

In the Drosophila optic lobes, the medulla processes visual information coming from inner photoreceptors R7 and R8 and from lamina neurons. It contains approximately 40,000 neurons belonging to more than 70 different types. Here we describe how precise temporal patterning of neural progenitors generates these different neural types. Five transcription factors-Homothorax, Eyeless, Sloppy paired, Dichaete and Tailless-are sequentially expressed in a temporal cascade in each of the medulla neuroblasts as they age. Loss of Eyeless, Sloppy paired or Dichaete blocks further progression of the temporal sequence. We provide evidence that this temporal sequence in neuroblasts, together with Notch-dependent binary fate choice, controls the diversification of the neuronal progeny. Although a temporal sequence of transcription factors had been identified in Drosophila embryonic neuroblasts, our work illustrates the generality of this strategy, with different sequences of transcription factors being used in different contexts.

Exploring genes and pathways underlying intellectual disability (ID) provides insight into brain development and function, clarifying the complex puzzle of how cognition develops. As part of ongoing systematic studies to identify candidate ID genes, linkage analysis and next-generation sequencing revealed Zinc Finger and BTB Domain Containing 11 (ZBTB11) as a novel candidate ID gene. ZBTB11 encodes a little-studied transcription regulator, and the two identified missense variants in this study are predicted to disrupt canonical Zn2+-binding residues of its C2H2 zinc finger domain, leading to possible altered DNA binding. Using HEK293T cells transfected with wild-type and mutant GFP-ZBTB11 constructs, we found the ZBTB11 mutants being excluded from the nucleolus, where the wild-type recombinant protein is predominantly localized. Pathway analysis applied to ChIP-seq data deposited in the ENCODE database supports the localization of ZBTB11 in nucleoli, highlighting associated pathways such as ribosomal RNA synthesis, ribosomal assembly, RNA modification and stress sensing, and provides a direct link between subcellular ZBTB11 location and its function. Furthermore, given the report of prominent brain and spinal cord degeneration in a zebrafish Zbtb11 mutant, we investigated ZBTB11-ortholog knockdown in Drosophila melanogaster brain by targeting RNAi using the UAS/Gal4 system. The observed approximate reduction to a third of the mushroom body size-possibly through neuronal reduction or degeneration-may affect neuronal circuits in the brain that are required for adaptive behavior, specifying the role of this gene in the nervous system. In conclusion, we report two ID families segregating ZBTB11 biallelic mutations disrupting Zn2+-binding motifs and provide functional evidence linking ZBTB11 dysfunction to this phenotype.