TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway.

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes.

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

TFAP2 paralogs facilitate chromatin access for MITF at pigmentation and cell proliferation genes.

In developing melanocytes and in melanoma cells, multiple paralogs of the Activating-enhancer-binding Protein 2 family of transcription factors (TFAP2) contribute to expression of genes encoding pigmentation regulators, but their interaction with Microphthalmia transcription factor (MITF), a master regulator of these cells, is unclear. Supporting the model that TFAP2 facilitates MITF's ability to activate expression of pigmentation genes, single-cell seq analysis of zebrafish embryos revealed that pigmentation genes are only expressed in the subset of mitfa-expressing cells that also express tfap2 paralogs. To test this model in SK-MEL-28 melanoma cells we deleted the two TFAP2 paralogs with highest expression, TFAP2A and TFAP2C, creating TFAP2 knockout (TFAP2-KO) cells. We then assessed gene expression, chromatin accessibility, binding of TFAP2A and of MITF, and the chromatin marks H3K27Ac and H3K27Me3 which are characteristic of active enhancers and silenced chromatin, respectively. Integrated analyses of these datasets indicate TFAP2 paralogs directly activate enhancers near genes enriched for roles in pigmentation and proliferation, and directly repress enhancers near genes enriched for roles in cell adhesion. Consistently, compared to WT cells, TFAP2-KO cells proliferate less and adhere to one another more. TFAP2 paralogs and MITF co-operatively activate a subset of enhancers, with the former necessary for MITF binding and chromatin accessibility. By contrast, TFAP2 paralogs and MITF do not appear to co-operatively inhibit enhancers. These studies reveal a mechanism by which TFAP2 profoundly influences the set of genes activated by MITF, and thereby the phenotype of pigment cells and melanoma cells.

Hypersensitivity of Zebrafish htr2b Mutant Embryos to Sertraline Indicates a Role for Serotonin Signaling in Cardiac Development.

ABSTRACT: Selective serotonin reuptake inhibitors (SSRIs) are antidepressants prescribed in 10% of pregnancies in the United States. Maternal use of SSRIs has been linked to an elevated rate of congenital heart defects, but the exact mechanism of pathogenesis is unknown. Previously, we have shown a decrease in cardiomyocyte proliferation, left ventricle size, and reduced cardiac expression of the serotonin receptor 5-HT 2B in offspring of mice exposed to the SSRI sertraline during pregnancy, relative to offspring of untreated mice. These results suggest that disruption of serotonin signaling leads to heart defects. Supporting this conclusion, we show here that zebrafish embryos exposed to sertraline develop with a smaller ventricle, reduced cardiomyocyte number, and lower cardiac expression of htr2b relative to untreated embryos. Moreover, zebrafish embryos homozygous for a nonsense mutation of htr2b ( htr2bsa16649 ) were sensitized to sertraline treatment relative to wild-type embryos. Specifically, the ventricle area was reduced in the homozygous htr2b mutants treated with sertraline compared with wild-type embryos treated with sertraline and homozygous htr2b mutants treated with vehicle control. Whereas long-term effects on left ventricle shortening fraction and stroke volume were observed by echocardiography in adult mice exposed to sertraline in utero, echocardiograms of adult zebrafish exposed to sertraline as embryos were normal. These results implicate the 5-HT 2B receptor functions in heart development and suggest zebrafish are a relevant animal model that can be used to investigate the connection between maternal SSRI use and elevated risk of congenital heart defects.

TFAP2 paralogs regulate melanocyte differentiation in parallel with MITF.

Mutations in the gene encoding transcription factor TFAP2A result in pigmentation anomalies in model organisms and premature hair graying in humans. However, the pleiotropic functions of TFAP2A and its redundantly-acting paralogs have made the precise contribution of TFAP2-type activity to melanocyte differentiation unclear. Defining this contribution may help to explain why TFAP2A expression is reduced in advanced-stage melanoma compared to benign nevi. To identify genes with TFAP2A-dependent expression in melanocytes, we profile zebrafish tissue and mouse melanocytes deficient in Tfap2a, and find that expression of a small subset of genes underlying pigmentation phenotypes is TFAP2A-dependent, including Dct, Mc1r, Mlph, and Pmel. We then conduct TFAP2A ChIP-seq in mouse and human melanocytes and find that a much larger subset of pigmentation genes is associated with active regulatory elements bound by TFAP2A. These elements are also frequently bound by MITF, which is considered the "master regulator" of melanocyte development. For example, the promoter of TRPM1 is bound by both TFAP2A and MITF, and we show that the activity of a minimal TRPM1 promoter is lost upon deletion of the TFAP2A binding sites. However, the expression of Trpm1 is not TFAP2A-dependent, implying that additional TFAP2 paralogs function redundantly to drive melanocyte differentiation, which is consistent with previous results from zebrafish. Paralogs Tfap2a and Tfap2b are both expressed in mouse melanocytes, and we show that mouse embryos with Wnt1-Cre-mediated deletion of Tfap2a and Tfap2b in the neural crest almost completely lack melanocytes but retain neural crest-derived sensory ganglia. These results suggest that TFAP2 paralogs, like MITF, are also necessary for induction of the melanocyte lineage. Finally, we observe a genetic interaction between tfap2a and mitfa in zebrafish, but find that artificially elevating expression of tfap2a does not increase levels of melanin in mitfa hypomorphic or loss-of-function mutants. Collectively, these results show that TFAP2 paralogs, operating alongside lineage-specific transcription factors such as MITF, directly regulate effectors of terminal differentiation in melanocytes. In addition, they suggest that TFAP2A activity, like MITF activity, has the potential to modulate the phenotype of melanoma cells.

Mutations in GDF11 and the extracellular antagonist, Follistatin, as a likely cause of Mendelian forms of orofacial clefting in humans.

Cleft lip with or without cleft palate (CL/P) is generally viewed as a complex trait with multiple genetic and environmental contributions. In 70% of cases, CL/P presents as an isolated feature and/or deemed nonsyndromic. In the remaining 30%, CL/P is associated with multisystem phenotypes or clinically recognizable syndromes, many with a monogenic basis. Here we report the identification, via exome sequencing, of likely pathogenic variants in two genes that encode interacting proteins previously only linked to orofacial clefting in mouse models. A variant in GDF11 (encoding growth differentiation factor 11), predicting a p.(Arg298Gln) substitution at the Furin protease cleavage site, was identified in one family that segregated with CL/P and both rib and vertebral hypersegmentation, mirroring that seen in Gdf11 knockout mice. In the second family in which CL/P was the only phenotype, a mutation in FST (encoding the GDF11 antagonist, Follistatin) was identified that is predicted to result in a p.(Cys56Tyr) substitution in the region that binds GDF11. Functional assays demonstrated a significant impact of the specific mutated amino acids on FST and GDF11 function and, together with embryonic expression data, provide strong evidence for the importance of GDF11 and Follistatin in the regulation of human orofacial development.

Irf6 directly regulates Klf17 in zebrafish periderm and Klf4 in murine oral epithelium, and dominant-negative KLF4 variants are present in patients with cleft lip and palate.

Non-syndromic (NS) cleft lip with or without cleft palate (CL/P) is a common disorder with a strong genetic underpinning. Genome-wide association studies have detected common variants associated with this disorder, but a large portion of the genetic risk for NSCL/P is conferred by unidentified rare sequence variants. Mutations in IRF6 (Interferon Regulatory Factor 6) and GRHL3 (Grainyhead-like 3) cause Van der Woude syndrome, which includes CL/P. Both genes encode members of a regulatory network governing periderm differentiation in model organisms. Here, we report that Krüppel-like factor 17 (Klf17), like Grhl3, acts downstream of Irf6 in this network in zebrafish periderm. Although Klf17 expression is absent from mammalian oral epithelium, a close homologue, Klf4, is expressed in this tissue and is required for the differentiation of epidermis. Chromosome configuration capture and reporter assays indicated that IRF6 directly regulates an oral-epithelium enhancer of KLF4. To test whether rare missense variants of KLF4 contribute risk for NSCL/P, we sequenced KLF4 in approximately 1000 NSCL/P cases and 300 controls. By one statistical test, missense variants of KLF4 as a group were enriched in cases versus controls. Moreover, two patient-derived KLF4 variants disrupted periderm differentiation upon forced expression in zebrafish embryos, suggesting that they have dominant-negative effect. These results indicate that rare NSCL/P risk variants can be found in members of the gene regulatory network governing periderm differentiation.

Identification of functional variants for cleft lip with or without cleft palate in or near PAX7, FGFR2, and NOG by targeted sequencing of GWAS loci.

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.

HOMER2, a stereociliary scaffolding protein, is essential for normal hearing in humans and mice.

Hereditary hearing loss is a clinically and genetically heterogeneous disorder. More than 80 genes have been implicated to date, and with the advent of targeted genomic enrichment and massively parallel sequencing (TGE+MPS) the rate of novel deafness-gene identification has accelerated. Here we report a family segregating post-lingual progressive autosomal dominant non-syndromic hearing loss (ADNSHL). After first excluding plausible variants in known deafness-causing genes using TGE+MPS, we completed whole exome sequencing in three hearing-impaired family members. Only a single variant, p.Arg185Pro in HOMER2, segregated with the hearing-loss phenotype in the extended family. This amino acid change alters a highly conserved residue in the coiled-coil domain of HOMER2 that is essential for protein multimerization and the HOMER2-CDC42 interaction. As a scaffolding protein, HOMER2 is involved in intracellular calcium homeostasis and cytoskeletal organization. Consistent with this function, we found robust expression in stereocilia of hair cells in the murine inner ear and observed that over-expression of mutant p.Pro185 HOMER2 mRNA causes anatomical changes of the inner ear and neuromasts in zebrafish embryos. Furthermore, mouse mutants homozygous for the targeted deletion of Homer2 present with early-onset rapidly progressive hearing loss. These data provide compelling evidence that HOMER2 is required for normal hearing and that its sequence alteration in humans leads to ADNSHL through a dominant-negative mode of action.

Zebrafish models of orofacial clefts.

Zebrafish is a model organism that affords experimental advantages toward investigating the normal function of genes associated with congenital birth defects. Here we summarize zebrafish studies of genes implicated in orofacial cleft (OFC). The most common use of zebrafish in this context has been to explore the normal function an OFC-associated gene product in craniofacial morphogenesis by inhibiting expression of its zebrafish ortholog. The most frequently deployed method has been to inject embryos with antisense morpholino oligonucleotides targeting the desired transcript. However, improvements in targeted mutagenesis strategies have led to widespread adoption of CRISPR/Cas9 technology. A second application of zebrafish has been for functional assays of gene variants found in OFC patients; such in vivo assays are valuable because the success of in silico methods for testing allele severity has been mixed. Finally, zebrafish have been used to test the tissue specificity of enhancers that harbor single nucleotide polymorphisms associated with risk for OFC. We review examples of each of these approaches in the context of genes that are implicated in syndromic and non-syndromic OFC. Developmental Dynamics 246:897-914, 2017. © 2017 Wiley Periodicals, Inc.

A single nucleotide polymorphism associated with isolated cleft lip and palate, thyroid cancer and hypothyroidism alters the activity of an oral epithelium and thyroid enhancer near FOXE1.

Three common diseases, isolated cleft lip and cleft palate (CLP), hypothyroidism and thyroid cancer all map to the FOXE1 locus, but causative variants have yet to be identified. In patients with CLP, the frequency of coding mutations in FOXE1 fails to account for the risk attributable to this locus, suggesting that the common risk alleles reside in nearby regulatory elements. Using a combination of zebrafish and mouse transgenesis, we screened 15 conserved non-coding sequences for enhancer activity, identifying three that regulate expression in a tissue specific pattern consistent with endogenous foxe1 expression. These three, located -82.4, -67.7 and +22.6 kb from the FOXE1 start codon, are all active in the oral epithelium or branchial arches. The -67.7 and +22.6 kb elements are also active in the developing heart, and the -67.7 kb element uniquely directs expression in the developing thyroid. Within the -67.7 kb element is the SNP rs7850258 that is associated with all three diseases. Quantitative reporter assays in oral epithelial and thyroid cell lines show that the rs7850258 allele (G) associated with CLP and hypothyroidism has significantly greater enhancer activity than the allele associated with thyroid cancer (A). Moreover, consistent with predicted transcription factor binding differences, the -67.7 kb element containing rs7850258 allele G is significantly more responsive to both MYC and ARNT than allele A. By demonstrating that this common non-coding variant alters FOXE1 expression, we have identified at least in part the functional basis for the genetic risk of these seemingly disparate disorders.

Dominant mutations in GRHL3 cause Van der Woude Syndrome and disrupt oral periderm development.

Mutations in interferon regulatory factor 6 (IRF6) account for ∼70% of cases of Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate. In 8 of 45 VWS-affected families lacking a mutation in IRF6, we found coding mutations in grainyhead-like 3 (GRHL3). According to a zebrafish-based assay, the disease-associated GRHL3 mutations abrogated periderm development and were consistent with a dominant-negative effect, in contrast to haploinsufficiency seen in most VWS cases caused by IRF6 mutations. In mouse, all embryos lacking Grhl3 exhibited abnormal oral periderm and 17% developed a cleft palate. Analysis of the oral phenotype of double heterozygote (Irf6(+/-);Grhl3(+/-)) murine embryos failed to detect epistasis between the two genes, suggesting that they function in separate but convergent pathways during palatogenesis. Taken together, our data demonstrated that mutations in two genes, IRF6 and GRHL3, can lead to nearly identical phenotypes of orofacial cleft. They supported the hypotheses that both genes are essential for the presence of a functional oral periderm and that failure of this process contributes to VWS.

Distinct requirements for wnt9a and irf6 in extension and integration mechanisms during zebrafish palate morphogenesis.

Development of the palate in vertebrates involves cranial neural crest migration, convergence of facial prominences and extension of the cartilaginous framework. Dysregulation of palatogenesis results in orofacial clefts, which represent the most common structural birth defects. Detailed analysis of zebrafish palatogenesis revealed distinct mechanisms of palatal morphogenesis: extension, proliferation and integration. We show that wnt9a is required for palatal extension, wherein the chondrocytes form a proliferative front, undergo morphological change and intercalate to form the ethmoid plate. Meanwhile, irf6 is required specifically for integration of facial prominences along a V-shaped seam. This work presents a mechanistic analysis of palate morphogenesis in a clinically relevant context.

Abnormal differentiation of dopaminergic neurons in zebrafish trpm7 mutant larvae impairs development of the motor pattern.

Transient receptor potential, melastatin-like 7 (Trpm7) is a combined ion channel and kinase implicated in the differentiation or function of many cell types. Early lethality in mice and frogs depleted of the corresponding gene impedes investigation of the functions of this protein particularly during later stages of development. By contrast, zebrafish trpm7 mutant larvae undergo early morphogenesis normally and thus do not have this limitation. The mutant larvae are characterized by multiple defects including melanocyte cell death, transient paralysis, and an ion imbalance that leads to the development of kidney stones. Here we report a requirement for Trpm7 in differentiation or function of dopaminergic neurons in vivo. First, trpm7 mutant larvae are hypomotile and fail to make a dopamine-dependent developmental transition in swim-bout length. Both of these deficits are partially rescued by the application of levodopa or dopamine. Second, histological analysis reveals that in trpm7 mutants a significant fraction of dopaminergic neurons lack expression of tyrosine hydroxylase, the rate-limiting enzyme in dopamine synthesis. Third, trpm7 mutants are unusually sensitive to the neurotoxin 1-methyl-4-phenylpyridinium, an oxidative stressor, and their motility is partially rescued by application of the iron chelator deferoxamine, an anti-oxidant. Finally, in SH-SY5Y cells, which model aspects of human dopaminergic neurons, forced expression of a channel-dead variant of TRPM7 causes cell death. In summary, a forward genetic screen in zebrafish has revealed that both melanocytes and dopaminergic neurons depend on the ion channel Trpm7. The mechanistic underpinning of this dependence requires further investigation.

Copy number variation analysis implicates the cell polarity gene glypican 5 as a human spina bifida candidate gene.

Neural tube defects (NTDs) are common birth defects of complex etiology. Family and population-based studies have confirmed a genetic component to NTDs. However, despite more than three decades of research, the genes involved in human NTDs remain largely unknown. We tested the hypothesis that rare copy number variants (CNVs), especially de novo germline CNVs, are a significant risk factor for NTDs. We used array-based comparative genomic hybridization (aCGH) to identify rare CNVs in 128 Caucasian and 61 Hispanic patients with non-syndromic lumbar-sacral myelomeningocele. We also performed aCGH analysis on the parents of affected individuals with rare CNVs where parental DNA was available (42 sets). Among the eight de novo CNVs that we identified, three generated copy number changes of entire genes. One large heterozygous deletion removed 27 genes, including PAX3, a known spina bifida-associated gene. A second CNV altered genes (PGPD8, ZC3H6) for which little is known regarding function or expression. A third heterozygous deletion removed GPC5 and part of GPC6, genes encoding glypicans. Glypicans are proteoglycans that modulate the activity of morphogens such as Sonic Hedgehog (SHH) and bone morphogenetic proteins (BMPs), both of which have been implicated in NTDs. Additionally, glypicans function in the planar cell polarity (PCP) pathway, and several PCP genes have been associated with NTDs. Here, we show that GPC5 orthologs are expressed in the neural tube, and that inhibiting their expression in frog and fish embryos results in NTDs. These results implicate GPC5 as a gene required for normal neural tube development.

Novel Tfap2-mediated control of soxE expression facilitated the evolutionary emergence of the neural crest.

Gene duplication has been proposed to drive the evolution of novel morphologies. After gene duplication, it is unclear whether changes in the resulting paralogs' coding-regions, or in their cis-regulatory elements, contribute most significantly to the assembly of novel gene regulatory networks. The Transcription Factor Activator Protein 2 (Tfap2) was duplicated in the chordate lineage and is essential for development of the neural crest, a tissue that emerged with vertebrates. Using a tfap2-depleted zebrafish background, we test the ability of available gnathostome, agnathan, cephalochordate and insect tfap2 paralogs to drive neural crest development. With the exception of tfap2d (lamprey and zebrafish), all are able to do so. Together with expression analyses, these results indicate that sub-functionalization has occurred among Tfap2 paralogs, but that neo-functionalization of the Tfap2 protein did not drive the emergence of the neural crest. We investigate whether acquisition of novel target genes for Tfap2 might have done so. We show that in neural crest cells Tfap2 directly activates expression of sox10, which encodes a transcription factor essential for neural crest development. The appearance of this regulatory interaction is likely to have coincided with that of the neural crest, because AP2 and SoxE are not co-expressed in amphioxus, and because neural crest enhancers are not detected proximal to amphioxus soxE. We find that sox10 has limited ability to restore the neural crest in Tfap2-deficient embryos. Together, these results show that mutations resulting in novel Tfap2-mediated regulation of sox10 and other targets contributed to the evolution of the neural crest.

A mutation in the Srrm4 gene causes alternative splicing defects and deafness in the Bronx waltzer mouse.

Sensory hair cells are essential for hearing and balance. Their development from epithelial precursors has been extensively characterized with respect to transcriptional regulation, but not in terms of posttranscriptional influences. Here we report on the identification and functional characterization of an alternative-splicing regulator whose inactivation is responsible for defective hair-cell development, deafness, and impaired balance in the spontaneous mutant Bronx waltzer (bv) mouse. We used positional cloning and transgenic rescue to locate the bv mutation to the splicing factor-encoding gene Ser/Arg repetitive matrix 4 (Srrm4). Transcriptome-wide analysis of pre-mRNA splicing in the sensory patches of embryonic inner ears revealed that specific alternative exons were skipped at abnormally high rates in the bv mice. Minigene experiments in a heterologous expression system confirmed that these skipped exons require Srrm4 for inclusion into the mature mRNA. Sequence analysis and mutagenesis experiments showed that the affected transcripts share a novel motif that is necessary for the Srrm4-dependent alternative splicing. Functional annotations and protein-protein interaction data indicated that the encoded proteins cluster in the secretion and neurotransmission pathways. In addition, the splicing of a few transcriptional regulators was found to be Srrm4 dependent, and several of the genes known to be targeted by these regulators were expressed at reduced levels in the bv mice. Although Srrm4 expression was detected in neural tissues as well as hair cells, analyses of the bv mouse cerebellum and neocortex failed to detect splicing defects. Our data suggest that Srrm4 function is critical in the hearing and balance organs, but not in all neural tissues. Srrm4 is the first alternative-splicing regulator to be associated with hearing, and the analysis of bv mice provides exon-level insights into hair-cell development.

Gene regulatory evolution and the origin of macroevolutionary novelties: insights from the neural crest.

The appearance of novel anatomic structures during evolution is driven by changes to the networks of transcription factors, signaling pathways, and downstream effector genes controlling development. The nature of the changes to these developmental gene regulatory networks (GRNs) is poorly understood. A striking test case is the evolution of the GRN controlling development of the neural crest (NC). NC cells emerge from the neural plate border (NPB) and contribute to multiple adult structures. While all chordates have a NPB, only in vertebrates do NPB cells express all the genes constituting the neural crest GRN (NC-GRN). Interestingly, invertebrate chordates express orthologs of NC-GRN components in other tissues, revealing that during vertebrate evolution new regulatory connections emerged between transcription factors primitively expressed in the NPB and genes primitively expressed in other tissues. Such interactions could have evolved by two mechanisms. First, transcription factors primitively expressed in the NPB may have evolved new DNA and/or cofactor binding properties (protein neofunctionalization). Alternately, cis-regulatory elements driving NPB expression may have evolved near genes primitively expressed in other tissues (cis-regulatory neofunctionalization). Here we discuss how gene duplication can, in principle, promote either form of neofunctionalization. We review recent published examples of interspecies gene-swap, or regulatory-element-swap, experiments that test both models. Such experiments have yielded little evidence to support the importance of protein neofunctionalization in the emergence of the NC-GRN, but do support the importance of novel cis-regulatory elements in this process. The NC-GRN is an excellent model for the study of gene regulatory and macroevolutionary innovation.

Differentiation of zebrafish melanophores depends on transcription factors AP2 alpha and AP2 epsilon.

A model of the gene-regulatory-network (GRN), governing growth, survival, and differentiation of melanocytes, has emerged from studies of mouse coat color mutants and melanoma cell lines. In this model, Transcription Factor Activator Protein 2 alpha (TFAP2A) contributes to melanocyte development by activating expression of the gene encoding the receptor tyrosine kinase Kit. Next, ligand-bound Kit stimulates a pathway activating transcription factor Microphthalmia (Mitf), which promotes differentiation and survival of melanocytes by activating expression of Tyrosinase family members, Bcl2, and other genes. The model predicts that in both Tfap2a and Kit null mutants there will be a phenotype of reduced melanocytes and that, because Tfap2a acts upstream of Kit, this phenotype will be more severe, or at least as severe as, in Tfap2a null mutants in comparison to Kit null mutants. Unexpectedly, this is not the case in zebrafish or mouse. Because many Tfap2 family members have identical DNA-binding specificity, we reasoned that another Tfap2 family member may work redundantly with Tfap2a in promoting Kit expression. We report that tfap2e is expressed in melanoblasts and melanophores in zebrafish embryos and that its orthologue, TFAP2E, is expressed in human melanocytes. We provide evidence that Tfap2e functions redundantly with Tfap2a to maintain kita expression in zebrafish embryonic melanophores. Further, we show that, in contrast to in kita mutants where embryonic melanophores appear to differentiate normally, in tfap2a/e doubly-deficient embryonic melanophores are small and under-melanized, although they retain expression of mitfa. Interestingly, forcing expression of mitfa in tfap2a/e doubly-deficient embryos partially restores melanophore differentiation. These findings reveal that Tfap2 activity, mediated redundantly by Tfap2a and Tfap2e, promotes melanophore differentiation in parallel with Mitf by an effector other than Kit. This work illustrates how analysis of single-gene mutants may fail to identify steps in a GRN that are affected by the redundant activity of related proteins.

Identification of early requirements for preplacodal ectoderm and sensory organ development.

Preplacodal ectoderm arises near the end of gastrulation as a narrow band of cells surrounding the anterior neural plate. This domain later resolves into discrete cranial placodes that, together with neural crest, produce paired sensory structures of the head. Unlike the better-characterized neural crest, little is known about early regulation of preplacodal development. Classical models of ectodermal patterning posit that preplacodal identity is specified by readout of a discrete level of Bmp signaling along a DV gradient. More recent studies indicate that Bmp-antagonists are critical for promoting preplacodal development. However, it is unclear whether Bmp-antagonists establish the proper level of Bmp signaling within a morphogen gradient or, alternatively, block Bmp altogether. To begin addressing these issues, we treated zebrafish embryos with a pharmacological inhibitor of Bmp, sometimes combined with heat shock-induction of Chordin and dominant-negative Bmp receptor, to fully block Bmp signaling at various developmental stages. We find that preplacodal development occurs in two phases with opposing Bmp requirements. Initially, Bmp is required before gastrulation to co-induce four transcription factors, Tfap2a, Tfap2c, Foxi1, and Gata3, which establish preplacodal competence throughout the nonneural ectoderm. Subsequently, Bmp must be fully blocked in late gastrulation by dorsally expressed Bmp-antagonists, together with dorsally expressed Fgf and Pdgf, to specify preplacodal identity within competent cells abutting the neural plate. Localized ventral misexpression of Fgf8 and Chordin can activate ectopic preplacodal development anywhere within the zone of competence, whereas dorsal misexpression of one or more competence factors can activate ectopic preplacodal development in the neural plate. Conversely, morpholino-knockdown of competence factors specifically ablates preplacodal development. Our work supports a relatively simple two-step model that traces regulation of preplacodal development to late blastula stage, resolves two distinct phases of Bmp dependence, and identifies the main factors required for preplacodal competence and specification.