Genome binding properties of Zic transcription factors underlie their changing functions during neuronal maturation

Background

The Zic family of transcription factors (TFs) promote both proliferation and maturation of cerebellar granule neurons (CGNs), raising the question of how a single, constitutively expressed TF family can support distinct developmental processes. Here we use an integrative experimental and bioinformatic approach to discover the regulatory relationship between Zic TF binding and changing programs of gene transcription during CGN differentiation.

Results

We first established a bioinformatic pipeline to integrate Zic ChIP-seq data from the developing mouse cerebellum with other genomic datasets from the same tissue. In newborn CGNs, Zic TF binding predominates at active enhancers that are co-bound by developmentally-regulated TFs including Atoh1, whereas in mature CGNs, Zic TF binding consolidates toward promoters where it co-localizes with activity-regulated TFs. We then performed CUT&RUN-seq in differentiating CGNs to define both the time course of developmental shifts in Zic TF binding and their relationship to gene expression. Mapping Zic TF binding sites to genes using chromatin looping, we identified the set of Zic target genes that have altered expression in RNA-seq from Zic1 or Zic2 knockdown CGNs.

Conclusion

Our data show that Zic TFs are required for both induction and repression of distinct, developmentally regulated target genes through a mechanism that is largely independent of changes in Zic TF binding. We suggest that the differential collaboration of Zic TFs with other TF families underlies the shift in their biological functions across CGN development.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

Multiomic Analysis of Neurons with Divergent Projection Patterns Identifies Novel Regulators of Axon Pathfinding

Axon pathfinding is a key step in neural circuits formation. However, the transcriptional mechanisms regulating its progression remain poorly understood. The binary decision of crossing or avoiding the midline taken by some neuronal axons during development represents a robust model to investigate the mechanisms that control the selection of axonal trajectories. Here, to identify novel regulators of axon guidance, this work compares the transcriptome and chromatin occupancy profiles of two neuronal subpopulations, ipsilateral (iRGC) and contralateral retinal ganglion cells (cRGC), with similar functions but divergent axon trajectories. These analyses retrieved a number of genes encoding for proteins not previously implicated in axon pathfinding. In vivo functional experiments confirm the implication of some of these candidates in axonal navigation. Among the candidate genes, γ-synuclein is identified as essential for inducing midline crossing. Footprint and luciferase assays demonstrate that this small-sized protein is regulated by the transcription factor (TF) Pou4f1 in cRGCs. It is also shown that Lhx2/9 are specifically expressed in iRGCs and control a program that partially overlaps with that regulated by Zic2, previously described as essential for iRGC specification. Overall, the analyses identify dozens of new molecules potentially involved in axon guidance and reveal the regulatory logic behind the selection of axonal trajectories.

Two distinct motifs for Zic-r.a drive specific gene expression in two cell lineages

Zic-r.a, a maternal transcription factor, specifies posterior fate in ascidian embryos. However, its direct target, Tbx6-r.b, does not contain typical Zic-r.a-binding sites in its regulatory region. Using an in vitro selection assay, we found that Zic-r.a binds to sites dissimilar to the canonical motif, by which it activates Tbx6-r.b in a sub-lineage of muscle cells. These sites with non-canonical motifs have weak affinity for Zic-r.a; therefore, it activates Tbx6-r.b only in cells expressing Zic-r.a abundantly. Meanwhile, we found that Zic-r.a expressed zygotically in late embryos activates neural genes through canonical sites. Because different zinc-finger domains of Zic-r.a are important for driving reporters with canonical and non-canonical sites, it is likely that the non-canonical motif is not a divergent version of the canonical motif. In other words, our data indicate that the non-canonical motif represents a motif distinct from the canonical motif. Thus, Zic-r.a recognizes two distinct motifs to activate two sets of genes at two timepoints in development. This article has an associated 'The people behind the papers' interview.

The regulatory pathway from genes directly activated by maternal factors to muscle structural genes in ascidian embryos

Striated muscle cells in the tail of ascidian tadpole larvae differentiate cell-autonomously. Although several key regulatory factors have been identified, the genetic regulatory pathway is not fully understood; comprehensive understanding of the regulatory pathway is essential for accurate modeling in order to deduce principles for gene regulatory network dynamics, and for comparative analysis on how ascidians have evolved the cell-autonomous gene regulatory mechanism. Here, we reveal regulatory interactions among three key regulatory factors, Zic-r.b, Tbx6-r.b and Mrf, and elucidate the mechanism by which these factors activate muscle structural genes. We reveal a cross-regulatory circuit among these regulatory factors, which maintains the expression of Tbx6-r.b and Mrf during gastrulation. Although these two factors combinatorially activate muscle structural genes in late-stage embryos, muscle structural genes are activated mainly by Tbx6-r.b before gastrulation. Time points when expression of muscle structural genes become first detectable are strongly correlated with the degree of Tbx6-r.b occupancy. Thus, the genetic pathway, starting with Tbx6-r.b and Zic-r.b, which are activated by maternal factors, and ending with expression of muscle structural genes, has been revealed.

Ascidian Zic Genes

Ascidians are tunicates, which constitute the sister group of vertebrates. The ascidian genome contains two Zic genes, called Zic-r.a (also called Macho-1) and Zic-r.b (ZicL). The latter is a multi-copy gene, and the precise copy number has not yet been determined. Zic-r.a is maternally expressed, and soon after fertilization Zic-r.a mRNA is localized in the posterior pole of the zygote. Zic-r.a protein is translated there and is involved in specification of posterior fate; in particular it is important for specification of muscle fate. Zic-r.a is also expressed zygotically in neural cells of the tailbud stage. On the other hand, Zic-r.b is first expressed in marginal cells of the vegetal hemisphere of 32-cell embryos and then in neural cells that contribute to the central nervous system during gastrulation. Zic-r.b is required first for specification of mesodermal tissues and then for specification of the central nervous system. Their upstream and downstream genetic pathways have been studied extensively by functional assays, which include gene knockdown and chromatin immunoprecipitation assays. Thus, ascidian Zic genes play central roles in specification of mesodermal and neural fates.

Lophotrochozoan Zic Genes

Lophotrochozoa is a sister taxon of Ecdysozoa in the Protostomia that includes mollusks, annelids, brachiopods, and platyhelminths. Recent studies have clarified the structure, expression, and roles of lophotrochozoan Zic family genes. Zic genes in oligochaete annelid Tubifex tubifex (freshwater sludge worm) and polychaete annelid Capitella teleta (bristle worm) are commonly expressed in a subset of developing brain and mesoderm derivatives. The latter includes the naïve mesoderm and the associated chaetal sacs in each body segment, although the segmentation processes differ between the two species. Furthermore, in brachiopod Terebratalia transversa (lamp shell), Zic is expressed in the anterior ectodermal domains and mesodermal derivatives, including those associated with the chaetal sacs. This result suggests the common involvement of Zic genes in the development of chaetae, a lophotrochozoan novelty acquired in the course of evolution. In addition, the highly simplified lophotrochozoan Dicyema acuticephalum (dicyemid mesozoan, a cephalopod endoparasite), which lost its gut, nervous system, and muscles during evolution, expresses its Zic genes in hermaphroditic gonads, highlighting the role of Zic genes in germ cell development. The role of Zic in head regeneration was revealed in studies on platyhelminth Schmidtea mediterranea (freshwater planarian). Planarian Zic expression was induced in a subpopulation of neoblasts that includes adult pluripotent stem cells. It is needed for head regeneration and production of an anterior signaling center. Suppression of Wnt-β-catenin signaling underlies Zic-mediated head regeneration, reminiscent of Wnt-β-catenin suppression by vertebrate Zic genes. Taken together, studies on the lophotrochozoan Zic genes are essential to understanding not only the roles of these genes in body plan evolution but also the molecular mechanism underlying adult stem cell regulation.

Cnidarian Zic Genes

To understand the ancestral and evolved roles of zic homologs, it is important to reconstruct the putative roles of ancient zic homologs in the animal phylogeny. Most studies of zic genes have been conducted in model systems that are members of the bilaterian phylum. However, two additional phyla have zic homologs encoded in their genomes. The three animal phyla that contain zic homologs all share a common ancestor and collectively are termed the parahoxozoans (cnidarians (corals, sea anemones, and jellyfish), placozoans (Trichoplax adhaerens), and bilaterians (chordates, insects, nematodes, annelids, echinoderms, etc.). In this chapter we briefly discuss our understanding of zic genes in the parahoxozoans with a particular focus on how expression of cnidarian zic homologs in the medusozoan Hydra vulgaris and the anthozoan Nematostella vectensis informs our understanding of the putative ancestral roles zic homologs played in the cnidarian-bilaterian common ancestor.

Zic Genes in Nematodes: A Role in Nervous System Development and Wnt Signaling

Transcription factors of the Zic family play important roles during animal development, and their misregulation has been implicated in several human diseases. Zic proteins are present in nematodes, and their function has been mostly studied in the model organism C. elegans. C. elegans possesses only one Zic family member, REF-2. Functional studies have shown that this factor plays a key role during the development of the nervous system, epidermis, and excretory system. In addition, they have revealed that the C. elegans Zic protein acts as an atypical mediator of the Wnt/β-catenin pathway. In other animals including vertebrates, Zic factors are also regulators of nervous system development and modulators of Wnt signaling, suggesting that these are evolutionary ancient functions of Zic proteins.

Zebrafish Rfx4 controls dorsal and ventral midline formation in the neural tube

BACKGROUND

Rfx winged-helix transcription factors, best known as key regulators of core ciliogenesis, also play ciliogenesis-independent roles during neural development. Mammalian Rfx4 controls neural tube morphogenesis via both mechanisms.

RESULTS

We set out to identify conserved aspects of rfx4 gene function during vertebrate development and to establish a new genetic model in which to analyze these mechanisms further. To this end, we have generated frame-shift alleles in the zebrafish rfx4 locus using CRISPR/Cas9 mutagenesis. Using RNAseq-based transcriptome analysis, in situ hybridization and immunostaining we identified a requirement for zebrafish rfx4 in the forming midlines of the caudal neural tube. These functions are mediated, least in part, through transcriptional regulation of several zic genes in the dorsal hindbrain and of foxa2 in the ventral hindbrain and spinal cord (floor plate).

CONCLUSIONS

The midline patterning functions of rfx4 are conserved, because rfx4 regulates transcription of foxa2 and zic2 in zebrafish and in mouse. In contrast, zebrafish rfx4 function is dispensable for forebrain morphogenesis, while mouse rfx4 is required for normal formation of forebrain ventricles in a ciliogenesis-dependent manner. Collectively, this report identifies conserved aspects of rfx4 function and establishes a robust new genetic model for in-depth dissection of these mechanisms. Developmental Dynamics 247:650-659, 2018. © 2017 Wiley Periodicals, Inc.

Odd-Paired: The Drosophila Zic Gene

Zinc finger in the cerebellum (Zic) proteins are a family of transcription factors with multiple roles during development, particularly in neural tissues. The founding member of the Zic family is the Drosophila odd-paired (opa) gene. The Opa protein has a DNA binding domain containing five Cys2His2-type zinc fingers and has been shown to act as a sequence-specific DNA binding protein. Opa has significant homology to mammalian Zic1, Zic2, and Zic3 within the zinc finger domain and in two other conserved regions outside that domain. opa was initially identified as a pair-rule gene, part of the hierarchy of genes that establish the segmental body plan of the early Drosophila embryo. However, its wide expression pattern during embryogenesis indicates it plays additional roles. Embryos deficient in opa die before hatching with aberrant segmentation but also with defects in larval midgut formation. Post-embryonically, opa plays important roles in adult head development and circadian rhythm. Based on extensive neural expression, opa is predicted to be involved in many aspects of neural development and behavior, like other proteins of the Zic family. Consensus DNA binding sites have been identified for Opa and have been shown to activate transcription in vivo. However, there is evidence Opa may serve as a transcriptional regulator in the absence of direct DNA binding, as has been seen for other Zic proteins.

Antagonism between β-catenin and Gata.a sequentially segregates the germ layers of ascidian embryos

Many animal embryos use nuclear β-catenin (nβ-catenin) during the segregation of endomesoderm (or endoderm) from ectoderm. This mechanism is thus likely to be evolutionarily ancient. In the ascidian embryo, nβ-catenin reiteratively drives binary fate decisions between ectoderm and endomesoderm at the 16-cell stage, and then between endoderm and margin (mesoderm and caudal neural) at the 32-cell stage. At the 16-cell stage, nβ-catenin activates endomesoderm genes in the vegetal hemisphere. At the same time, nβ-catenin suppresses the DNA-binding activity of a maternal transcription factor, Gata.a, through a physical interaction, and Gata.a thereby activates its target genes only in the ectodermal lineage. In the present study, we found that this antagonism between nβ-catenin and Gata.a also operates during the binary fate switch at the 32-cell stage. Namely, in marginal cells where nβ-catenin is absent, Gata.a directly activates its target, Zic-r.b (ZicL), to specify the marginal cell lineages. Thus, the antagonistic action between nβ-catenin and Gata.a is involved in two consecutive stages of germ layer segregation in ascidian embryos.

Odd-paired controls frequency doubling in Drosophila segmentation by altering the pair-rule gene regulatory network

The Drosophila embryo transiently exhibits a double-segment periodicity, defined by the expression of seven 'pair-rule' genes, each in a pattern of seven stripes. At gastrulation, interactions between the pair-rule genes lead to frequency doubling and the patterning of 14 parasegment boundaries. In contrast to earlier stages of Drosophila anteroposterior patterning, this transition is not well understood. By carefully analysing the spatiotemporal dynamics of pair-rule gene expression, we demonstrate that frequency-doubling is precipitated by multiple coordinated changes to the network of regulatory interactions between the pair-rule genes. We identify the broadly expressed but temporally patterned transcription factor, Odd-paired (Opa/Zic), as the cause of these changes, and show that the patterning of the even-numbered parasegment boundaries relies on Opa-dependent regulatory interactions. Our findings indicate that the pair-rule gene regulatory network has a temporally modulated topology, permitting the pair-rule genes to play stage-specific patterning roles.

How targets select activation or repression in response to Wnt

In metazoans, the Wnt signaling pathway plays a key role in the regulation of binary decisions during development. During this process different sets of target genes are activated in cells where the Wnt pathway is active (classic target genes) versus cells where the pathway is inactive (opposite target genes). While the mechanism of transcriptional activation is well understood for classic target genes, how opposite target genes are activated in the absence of Wnt remains poorly characterized. Here we discuss how the key transcriptional mediator of the Wnt pathway, the TCF family member POP-1, regulates opposite target genes during C. elegans development. We examine recent findings suggesting that the direction of the transcriptional output (activation or repression) can be determined by the way TCF is recruited and physically interacts with its target gene.

Induction and dorsal restriction of Paired-box 3 (Pax3) gene expression in the caudal neuroectoderm is mediated by integration of multiple pathways on a short neural crest enhancer

Pax3 encodes a paired-box transcription factor with key roles in neural crest and neural tube ontogenesis. Robust control of Pax3 neural expression is ensured by two redundant sets of cis-regulatory modules (CRMs) that integrate anterior-posterior (such as Wnt-βCatenin signaling) as well as dorsal-ventral (such as Shh-Gli signaling) instructive cues. In previous work, we sought to characterize the Wnt-mediated regulation of Pax3 expression and identified the Cdx transcription factors (Cdx1/2/4) as critical intermediates in this process. We identified the neural crest enhancer-2 (NCE2) from the 5'-flanking region of Pax3 as a Cdx-dependent CRM that recapitulates the restricted expression of Pax3 in the mouse caudal neuroectoderm. While this is consistent with a key role in relaying the inductive signal from posteriorizing Wnt ligands, the broad expression of Cdx proteins in the tailbud region is not consistent with the restricted activity of NCE2. This implies that other positive and/or negative inputs are required and, here, we report a novel role for the transcription factor Zic2 in this regulation. Our data strongly suggests that Zic2 is involved in the induction (as a direct Pax3NCE2 activator and Cdx neural cofactor) as well as the maintenance of Pax3 dorsal restriction (as a target of the ventral Shh repressive input). We also provide evidence that the inductive Cdx-Zic2 interaction is integrated on NCE2 with a positive input from the neural-specific transcription factor Sox2. Altogether, our data provide important mechanistic insights into the coordinated integration of different signaling pathways on a short Pax3 CRM.

Stem cell-dependent formation of a functional anterior regeneration pole in planarians requires Zic and Forkhead transcription factors

Planarians can regenerate their head within days. This process depends on the direction of adult stem cells to wound sites and the orchestration of their progenitors to commit to appropriate lineages and to arrange into patterned tissues. We identified a zinc finger transcription factor, Smed-ZicA, as a downstream target of Smed-FoxD, a Forkhead transcription factor required for head regeneration. Smed-zicA and Smed-FoxD are co-expressed with the Wnt inhibitor notum and the Activin inhibitor follistatin in a cluster of cells at the anterior-most tip of the regenerating head - the anterior regeneration pole - and in surrounding stem cell progeny. Depletion of Smed-zicA and Smed-FoxD by RNAi abolishes notum and follistatin expression at the pole and inhibits head formation downstream of initial polarity decisions. We suggest a model in which ZicA and FoxD transcription factors synergize to control the formation of Notum- and Follistatin-producing anterior pole cells. Pole formation might constitute an early step in regeneration, resulting in a signaling center that orchestrates cellular events in the growing tissue.

Pax3 synergizes with Gli2 and Zic1 in transactivating the Myf5 epaxial somite enhancer

Both Glis, the downstream effectors of hedgehog signaling, and Zic transcription factors are required for Myf5 expression in the epaxial somite. Here we demonstrate a novel synergistic interaction between members of both families and Pax3, a paired-domain transcription factor that is essential for both myogenesis and neural crest development. We show that Pax3 synergizes with both Gli2 and Zic1 in transactivating the Myf5 epaxial somite (ES) enhancer in concert with the Myf5 promoter. This synergy is dependent on conserved functional domains of the proteins, as well as on a novel homeodomain motif in the Myf5 promoter and the essential Gli motif in the ES enhancer. Importantly, overexpression of Zic1 and Pax3 in the 10T1/2 mesodermal cell model results in enrichment of these factors at the endogenous Myf5 locus and induction of Myf5 expression. In our previous work, we showed that by enhancing nuclear translocation of Gli factors, Zics provide spatiotemporal patterning for Gli family members in the epaxial induction of Myf5 expression. Our current study indicates a complementary mechanism in which association with DNA-bound Pax3 strengthens the ability of both Zic1 and Gli2 to transactivate Myf5 in the epaxial somite.