C2H2 zinc finger genes of the Gli, Zic, KLF, SP, Wilms' tumour, Huckebein, Snail, Ovo, Spalt, Odd, Blimp-1, Fez and related gene families from Branchiostoma floridae

Single-cell Transcriptomic Studies Unveil Potential Nodes of the Notochord Gene Regulatory Network

Transcription factors (TFs) are DNA-binding proteins able to modulate the timing, location, and levels of gene expression by binding to regulatory DNA regions. Therefore, the repertoire of TFs present in the genome of a multicellular organism and the expression of variable constellations of TFs in different cellular cohorts determine the distinctive characteristics of developing tissues and organs. The information on tissue-specific assortments of TFs, their cross-regulatory interactions, and the genes/regulatory regions targeted by each TF is summarized in gene regulatory networks (GRNs), which provide genetic blueprints for the specification, development, and differentiation of multicellular structures. In this study, we review recent transcriptomic studies focused on the complement of TFs expressed in the notochord, a distinctive feature of all chordates. We analyzed notochord-specific datasets available from organisms representative of the three chordate subphyla, and highlighted lineage-specific variations in the suite of TFs expressed in their notochord. We framed the resulting findings within a provisional evolutionary scenario, which allows the formulation of hypotheses on the genetic/genomic changes that sculpted the structure and function of the notochord on an evolutionary scale.

Evolutionary origin of the chordate nervous system revealed by amphioxus developmental trajectories

The common ancestor of all vertebrates had a highly sophisticated nervous system, but questions remain about the evolution of vertebrate neural cell types. The amphioxus, a chordate that diverged before the origin of vertebrates, can inform vertebrate evolution. Here we develop and analyse a single-cell RNA-sequencing dataset from seven amphioxus embryo stages to understand chordate cell type evolution and to study vertebrate neural cell type origins. We identified many new amphioxus cell types, including homologues to the vertebrate hypothalamus and neurohypophysis, rooting the evolutionary origin of these structures. On the basis of ancestor-descendant reconstruction of cell trajectories of the amphioxus and other species, we inferred expression dynamics of transcription factor genes throughout embryogenesis and identified three ancient developmental routes forming chordate neurons. We characterized cell specification at the mechanistic level and generated mutant lines to examine the function of five key transcription factors involved in neural specification. Our results show three developmental origins for the vertebrate nervous system: an anterior FoxQ2-dependent mechanism that is deeply conserved in invertebrates, a less-conserved route leading to more posterior neurons in the vertebrate spinal cord and a mechanism for specifying neuromesoderm progenitors that is restricted to chordates. The evolution of neuromesoderm progenitors may have led to a dramatic shift in posterior neural and mesodermal cell fate decisions and the body elongation process in a stem chordate.

Amphioxus Gli knockout disrupts the development of left-right asymmetry but has limited impact on neural patterning

The Gli transcription factors are the primary mediators of Hedgehog (Hh) signaling. Vertebrate genomes contain multiple Gli paralogues with different functions downstream of Hh signal receipt, in part explaining the complexity of cellular responses to Hh that allow concentration-dependent target gene activation. Amphioxus is a chordate that split from the vertebrate lineage early in the evolution of chordates, before the genome duplications that occurred in early vertebrate evolution. It has a single Gli gene whose transcripts can be alternately spliced to yield two protein isoforms called GliS and GliL. We generated two knockout mutations in amphioxus Gli, one that affects the whole gene and a second that only affects GliL. Both knockouts showed major morphological and molecular defects in the development of left-right asymmetry, a phenotype that is similar but not identical to that previously found in Hh mutants. Hh signaling also patterns the amphioxus neural tube. Here, however, knockout of GliL showed no identifiable phenotype, while knockout of the full gene showed only small changes to the expression of one gene family, Olig. Other genes that were prominently affected by Hh knockout were not altered in expression in either knockout. Reasons for the differences between Hh and Gli knockouts in the pharynx and neural tube are discussed in the context of the likely different functions of amphioxus Gli isoforms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00195-w.

PpZAT5 suppresses the expression of a B-box gene PpBBX18 to inhibit anthocyanin biosynthesis in the fruit peel of red pear

BBX (B-box) proteins play a vital role in light-induced anthocyanin biosynthesis. PpBBX18 was an indispensable regulator for the induction of anthocyanin biosynthesis in the peel of red pear fruit (Pyrus pyrifolia Nakai.). However, the upstream regulation of BBX genes has not been well characterized. In this study, PpZAT5, a cysteine2/histidine2-type transcription factor, was discovered as the upstream negative regulator of PpBBX18. The results showed that PpZAT5 functions as a transcriptional repressor and directly binds to the CAAT motif of PpBBX18 and inhibits its expression. PpZAT5 expression was inhibited by light, which is converse to the expression pattern of anthocyanin-related structural genes. In addition, less anthocyanin accumulated in the PpZAT5-overexpressing pear calli than in the wild-type pear calli; on the contrary, more anthocyanin accumulated in PpZAT5-RNAi pear calli. Moreover, the crucial genes involved in light-induced anthocyanin biosynthesis were markedly down-regulated in the transcriptome of PpZAT5 overexpression pear calli compared to wild-type. In conclusion, our study indicates that PpBBX18 is negatively regulated by a C2H2-type transcriptional repressor, PpZAT5, which reduces anthocyanin content in pear. The present results demonstrate an upstream molecular mechanism of PpBBX18 and provide insights into light-induced anthocyanin biosynthesis.

Bioinformatics analysis of KLF2 as a potential prognostic factor in ccRCC and association with epithelial‑mesenchymal transition

Clear cell renal cell carcinoma (ccRCC) is a primary pathological subtype of RCC and has poor clinical outcome. Krüppel-like factors (KLFs), which are zinc-finger proteins, may be involved in ccRCC development and progression. KLFs belong to the zinc-finger family of DNA-binding transcription factors and regulate transcription of downstream target genes. KLFs are involved in cancer development. The present study aimed to investigate the role of KLFs in ccRCC prognosis. The Cancer Genome Atlas database and multifactorial analysis showed that KLFs were widely expressed in pan-cancers and KLF2 was an independent protective factor for ccRCC prognosis. Patients with low KLF2 expression had a low survival probability and expression of KLF2 was downregulated in patients with ccRCC with high pathological grade (II + III vs. I). In addition, western blot and reverse transcription-quantitative PCR revealed that KLF2 was expressed at low levels in ccRCC cell lines and overexpression of KLF2 inhibited cell migration. In addition, KLF2 expression was negatively correlated with methylation. KLF2 expression was elevated following treatment of ccRCC cells with DNA methyltransferase inhibitor. A prognostic risk index prediction model was constructed based on multiple Cox regression. The receiver operating characteristic curve was 0.780 (area under curve >0.5). Furthermore, Gene Ontology enrichment analysis showed that ‘cell adhesion’ and ‘junction’ were negatively correlated with KLF2 and that high-risk group exhibited significantly activated ‘epithelial-mesenchymal transition’. Western blot analysis showed that overexpression of KLF2 increased expression of E-cadherin, while decreasing levels of N-cadherin and vimentin. The present study highlighted the role of KLFs in ccRCC prognosis prediction and provides a research base for the search of validated prognostic biological markers for ccRCC.

Genome-Wide Identification and Expression Patterns of the C2H2-Zinc Finger Gene Family Related to Stress Responses and Catechins Accumulation in Camellia sinensis [L.] O. Kuntze

The C2H2-zinc finger protein (C2H2-ZFP) is essential for the regulation of plant development and widely responsive to diverse stresses including drought, cold and salt stress, further affecting the late flavonoid accumulation in higher plants. Tea is known as a popular beverage worldwide and its quality is greatly dependent on the physiological status and growing environment of the tea plant. To date, the understanding of C2H2-ZFP gene family in Camellia sinensis [L.] O. Kuntze is not yet available. In the present study, 134 CsC2H2-ZFP genes were identified and randomly distributed on 15 chromosomes. The CsC2H2-ZFP gene family was classified into four clades and gene structures and motif compositions of CsC2H2-ZFPs were similar within the same clade. Segmental duplication and negative selection were the main forces driving the expansion of the CsC2H2-ZFP gene family. Expression patterns suggested that CsC2H2-ZFPs were responsive to different stresses including drought, salt, cold and methyl jasmonate (MeJA) treatment. Specially, several C2H2-ZFPs showed a significant correlation with the catechins content and responded to the MeJA treatment, which might contribute to the tea quality and specialized astringent taste. This study will lay the foundations for further research of C2H2-type zinc finger proteins on the stress responses and quality-related metabolites accumulation in C. sinensis.

Krüppel-like factor/specificity protein evolution in the Spiralia and the implications for cephalopod visual system novelties

The cephalopod visual system is an exquisite example of convergence in biological complexity. However, we have little understanding of the genetic and molecular mechanisms underpinning its elaboration. The generation of new genetic material is considered a significant contributor to the evolution of biological novelty. We sought to understand if this mechanism may be contributing to cephalopod-specific visual system novelties. Specifically, we identified duplications in the Krüppel-like factor/specificity protein (KLF/SP) sub-family of C2H2 zinc-finger transcription factors in the squid Doryteuthis pealeii. We cloned and analysed gene expression of the KLF/SP family, including two paralogs of the DpSP6-9 gene. These duplicates showed overlapping expression domains but one paralog showed unique expression in the developing squid lens, suggesting a neofunctionalization of DpSP6-9a. To better understand this neofunctionalization, we performed a thorough phylogenetic analysis of SP6-9 orthologues in the Spiralia. We find multiple duplications and losses of the SP6-9 gene throughout spiralian lineages and at least one cephalopod-specific duplication. This work supports the hypothesis that gene duplication and neofunctionalization contribute to novel traits like the cephalopod image-forming eye and to the diversity found within Spiralia.

C2H2 Zinc Finger Proteins: Master Regulators of Abiotic Stress Responses in Plants

Abiotic stresses such as drought and salinity are major environmental factors that limit crop yields. Unraveling the molecular mechanisms underlying abiotic stress resistance is crucial for improving crop performance and increasing productivity under adverse environmental conditions. Zinc finger proteins, comprising one of the largest transcription factor families, are known for their finger-like structure and their ability to bind Zn2+. Zinc finger proteins are categorized into nine subfamilies based on their conserved Cys and His motifs, including the Cys2/His2-type (C2H2), C3H, C3HC4, C2HC5, C4HC3, C2HC, C4, C6, and C8 subfamilies. Over the past two decades, much progress has been made in understanding the roles of C2H2 zinc finger proteins in plant growth, development, and stress signal transduction. In this review, we focus on recent progress in elucidating the structures, functions, and classifications of plant C2H2 zinc finger proteins and their roles in abiotic stress responses.

Evolutionary History of GLIS Genes Illuminates Their Roles in Cell Reprograming and Ciliogenesis

The GLIS family transcription factors, GLIS1 and GLIS3, potentiate generation of induced pluripotent stem cells (iPSCs). In contrast, another GLIS family member, GLIS2, suppresses cell reprograming. To understand how these disparate roles arose, we examined evolutionary origins and genomic organization of GLIS genes. Comprehensive phylogenetic analysis shows that GLIS1 and GLIS3 originated during vertebrate whole genome duplication, whereas GLIS2 is a sister group to the GLIS1/3 and GLI families. This result is consistent with their opposing functions in cell reprograming. Glis1 evolved faster than Glis3, losing many protein-interacting motifs. This suggests that Glis1 acquired new functions under weakened evolutionary constraints. In fact, GLIS1 induces induced pluripotent stem cells more strongly. Transcriptomic data from various animal embryos demonstrate that glis1 is maternally expressed in some tetrapods, whereas vertebrate glis3 and invertebrate glis1/3 genes are rarely expressed in oocytes, suggesting that vertebrate (or tetrapod) Glis1 acquired a new expression domain and function as a maternal factor. Furthermore, comparative genomic analysis reveals that glis1/3 is part of a bilaterian-specific gene cluster, together with rfx3, ndc1, hspb11, and lrrc42. Because known functions of these genes are related to cilia formation and function, the last common ancestor of bilaterians may have acquired this cluster by shuffling gene order to establish more sophisticated epithelial tissues involving cilia. This evolutionary study highlights the significance of GLIS1/3 for cell reprograming, development, and diseases in ciliated organs such as lung, kidney, and pancreas.

Silencing of brain-expressed X-linked 2 (BEX2) promotes colorectal cancer metastasis through the Hedgehog signaling pathway

The incidence of colorectal cancer is increasing, and cancer metastasis is one of the major causes of poor outcomes. BEX2 has been reported to be involved in tumor development in several types of cancer, but its role in metastatic colorectal cancer remains largely undefined. Herein, we demonstrated that BEX2 knockout resulted in enhanced migratory and metastatic potential in colorectal cancer cells both in vitro and in vivo, and re-expression of BEX2 in knockout cells could reverse the enhanced migratory capacity. RNA-Seq results indicated that the hedgehog signaling pathway was activated after BEX2 knockout; moreover, the hedgehog signaling inhibitors, GANT61 and GDC-0449 could reverse the migratory enhancement of BEX2-/- colorectal cancer cells. We also demonstrated that the nuclear translocation of Zic2 after BEX2 silencing could activate the hedgehog signaling pathway, while Zic2 knockdown abrogated the migratory enhancement of BEX2-/- cells and inhibited the hedgehog signaling pathway. In summary, our findings suggest that BEX2 negatively modulates the hedgehog signaling pathway by retaining Zic2 in the cytoplasm in colorectal cancer cells, thereby inhibiting migration and metastasis of colorectal cancer cells.

Mutations in TFAP2B and previously unimplicated genes of the BMP, Wnt, and Hedgehog pathways in syndromic craniosynostosis

Craniosynostosis (CS) is a frequent congenital malformation featuring premature fusion of cranial sutures; 15% of these children have syndromic disease, often due to rare mutations with large effect. While many genes causing Mendelian forms of syndromic CS have been identified, clinical sequencing often fails to identify a likely causative mutation. We performed whole-exome sequencing of 12 case-parent trios with previously negative genetic evaluations. The results identified likely pathogenic mutations in TFAP2B, KAT6A, GLI2, SOX11, CTNNA1, and GPC4 in these families, adding several loci to those known to cause syndromic CS. The findings have implications for determining risk of disease in subsequent offspring and demonstrate that unexplained syndromic CS cases are a particularly rich vein for discovery of CS loci.

Craniosynostosis (CS) is a frequent congenital anomaly featuring the premature fusion of 1 or more sutures of the cranial vault. Syndromic cases, featuring additional congenital anomalies, make up 15% of CS. While many genes underlying syndromic CS have been identified, the cause of many syndromic cases remains unknown. We performed exome sequencing of 12 syndromic CS cases and their parents, in whom previous genetic evaluations were unrevealing. Damaging de novo or transmitted loss of function (LOF) mutations were found in 8 genes that are highly intolerant to LOF mutation (P = 4.0 × 10⁻⁸); additionally, a rare damaging mutation in SOX11, which has a lower level of intolerance, was identified. Four probands had rare damaging mutations (2 de novo) in TFAP2B, a transcription factor that orchestrates neural crest cell migration and differentiation; this mutation burden is highly significant (P = 8.2 × 10⁻¹²). Three probands had rare damaging mutations in GLI2, SOX11, or GPC4, which function in the Hedgehog, BMP, and Wnt signaling pathways; other genes in these pathways have previously been implicated in syndromic CS. Similarly, damaging de novo mutations were identified in genes encoding the chromatin modifier KAT6A, and CTNNA1, encoding catenin α-1. These findings establish TFAP2B as a CS gene, have implications for assessing risk to subsequent children in these families, and provide evidence implicating other genes in syndromic CS. This high yield indicates the value of performing exome sequencing of syndromic CS patients when sequencing of known disease loci is unrevealing.

Molecular regionalization of the developing amphioxus neural tube challenges major partitions of the vertebrate brain

All vertebrate brains develop following a common Bauplan defined by anteroposterior (AP) and dorsoventral (DV) subdivisions, characterized by largely conserved differential expression of gene markers. However, it is still unclear how this Bauplan originated during evolution. We studied the relative expression of 48 genes with key roles in vertebrate neural patterning in a representative amphioxus embryonic stage. Unlike nonchordates, amphioxus develops its central nervous system (CNS) from a neural plate that is homologous to that of vertebrates, allowing direct topological comparisons. The resulting genoarchitectonic model revealed that the amphioxus incipient neural tube is unexpectedly complex, consisting of several AP and DV molecular partitions. Strikingly, comparison with vertebrates indicates that the vertebrate thalamus, pretectum, and midbrain domains jointly correspond to a single amphioxus region, which we termed Di-Mesencephalic primordium (DiMes). This suggests that these domains have a common developmental and evolutionary origin, as supported by functional experiments manipulating secondary organizers in zebrafish and mice.

Amphioxus Sp5 is a member of a conserved Specificity Protein complement and is modulated by Wnt/β-catenin signalling

A cluster of three Specificity Protein (Sp) genes (Sp1-4, Sp5 and Sp6-9) is thought to be ancestral in both chordates and the wider Eumetazoa. Sp5 and Sp6-9 gene groups are associated with embryonic growth zones, such as tailbuds, and are both Wnt/β-catenin signalling pathway members and targets. Currently, there are conflicting reports as to the number and identity of Sp genes in the cephalochordates, the sister group to the vertebrates and urochordates. We confirm the SP complement of Branchiostoma belcheri and Branchiostoma lanceolatum, as well as their genomic arrangement, protein domain structure and residue frequency. We assay Sp5 expression in B. lanceolatum embryos, and determine its response to pharmacologically increased β-catenin signalling. Branchiostoma possesses three Sp genes, located on the same genomic scaffold. Phylogenetic and domain structure analyses are consistent with their identification as SP1-4, SP5 and SP6-9, although SP1-4 contains a novel glutamine-rich N-terminal region. SP5 is expressed in axial mesoderm and neurectoderm, and marks the cerebral vesicle and presumptive pharynx. Early exposure to increased β-catenin caused ubiquitous SP5 expression in late gastrula, while later treatment at gastrula stages reduced SP5 expression in the posterior growth zone during axis elongation. Amphioxus possess a typical invertebrate eumetazoan SP complement, and SP5 expression in embryos is well conserved with vertebrate homologues. Its expression in the tailbud, a posterior growth zone, is consistent with expression seen in other bilaterians. Branchiostoma SP5 shows a dynamic response to Wnt/β-catenin signalling.

Bilaterian-like promoters in the highly compact Amphimedon queenslandica genome

The regulatory systems underlying animal development must have evolved prior to the emergence of eumetazoans (cnidarians and bilaterians). Although representatives of earlier-branching animals - sponges ctenophores and placozoans - possess most of the developmental transcription factor families present in eumetazoans, the DNA regulatory elements that these transcription factors target remain uncharted. Here we characterise the core promoter sequences, U1 snRNP-binding sites (5' splice sites; 5'SSs) and polyadenylation sites (PASs) in the sponge Amphimedon queenslandica. Similar to unicellular opisthokonts, Amphimedon's genes are tightly packed in the genome and have small introns. In contrast, its genes possess metazoan-like core promoters populated with binding motifs previously deemed to be specific to vertebrates, including Nrf-1 and Krüppel-like elements. Also as in vertebrates, Amphimedon's PASs and 5'SSs are depleted downstream and upstream of transcription start sites, respectively, consistent with non-elongating transcripts being short-lived; PASs and 5'SSs are more evenly distributed in bidirectional promoters in Amphimedon. The presence of bilaterian-like regulatory DNAs in sponges is consistent with these being early and essential innovations of the metazoan gene regulatory repertoire.

The Role of WT1 in Embryonic Development and Normal Organ Homeostasis

The Wilms' tumor suppressor gene 1 (Wt1) is critically involved in a number of developmental processes in vertebrates, including cell differentiation, control of the epithelial/mesenchymal phenotype, proliferation, and apoptosis. Wt1 proteins act as transcriptional and post-transcriptional regulators, in mRNA splicing and in protein-protein interactions. Furthermore, Wt1 is involved in adult tissue homeostasis, kidney function, and cancer. For these reasons, Wt1 function has been extensively studied in a number of animal models to establish its spatiotemporal expression pattern and the developmental fate of the cells expressing this gene. In this chapter, we review the developmental anatomy of Wt1, collecting information about its dynamic expression in mesothelium, kidney, gonads, cardiovascular system, spleen, nervous system, lung, and liver. We also describe the adult expression of Wt1 in kidney podocytes, gonads, mesothelia, visceral adipose tissue, and a small fraction of bone marrow cells. We have reviewed the available animal models for Wt1-expressing cell lineage analysis, including direct Wt1 expression reporters and systems for permanent Wt1 lineage tracing, based on constitutive or inducible Cre recombinase expression under control of a Wt1 promoter. Finally we provide a number of laboratory protocols to be used with these animal models in order to assess reporter expression.

KLF/SP Transcription Factor Family Evolution: Expansion, Diversification, and Innovation in Eukaryotes

The Krüppel-like factor and specificity protein (KLF/SP) genes play key roles in critical biological processes including stem cell maintenance, cell proliferation, embryonic development, tissue differentiation, and metabolism and their dysregulation has been implicated in a number of human diseases and cancers. Although many KLF/SP genes have been characterized in a handful of bilaterian lineages, little is known about the KLF/SP gene family in nonbilaterians and virtually nothing is known outside the metazoans. Here, we analyze and discuss the origins and evolutionary history of the KLF/SP transcription factor family and associated transactivation/repression domains. We have identified and characterized the complete KLF/SP gene complement from the genomes of 48 species spanning the Eukarya. We have also examined the phylogenetic distribution of transactivation/repression domains associated with this gene family. We report that the origin of the KLF/SP gene family predates the divergence of the Metazoa. Furthermore, the expansion of the KLF/SP gene family is paralleled by diversification of transactivation domains via both acquisitions of pre-existing ancient domains as well as by the appearance of novel domains exclusive to this gene family and is strongly associated with the expansion of cell type complexity.

C2H2 zinc finger proteins of the SP/KLF, Wilms tumor, EGR, Huckebein, and Klumpfuss families in metazoans and beyond

Specificity proteins (SPs) and Krüppel-Like Factors (KLFs) are C2H2-type zinc finger transcription factors that play essential roles in differentiation, development, proliferation and cell death. SP/KLF proteins, similarly to Wilms tumor protein 1 (WT1), Early Growth Response (EGR), Huckebein, and Klumpfuss, prefer to bind GC-rich sequences such as GC-box and CACCC-box (GT-box). We searched various genomes and transcriptomes of metazoans and single-cell holozoans for members of these families. Seven groups of KLFs (KLFA-G) and three groups of SPs (SPA-C) were identified in the three lineages of Bilateria (Deuterostomia, Ecdysozoa, and Lophotrochozoa). The last ancestor of jawed vertebrates was inferred to have at least 18 KLFs (group A: KLF1/2/4/17, group B: KLF3/8/12; group C: KLF5/5l; group D: KLF6/7; group E: KLF9/13/16; group F: KLF10/KLF11; group G: KLF15/15l) and 10 SPs (group A: SP1/2/3/4; group B: SP5/5l; group C: SP6/7/8/9), since they were found in both cartilaginous and boned fishes. Placental mammals have added KLF14 (group E) and KLF18 (group A), and lost KLF5l (KLF5-like) and KLF15l (KLF15-like). Multiple KLF members were found in basal metazoans (Ctenophora, Porifera, Placozoa, and Cnidaria). Ctenophora has the least number of KLFs and no SPs, which could be attributed to its proposed sister group relationship to other metazoans or gene loss. While SP, EGR and Klumpfuss were only detected in metazoans, KLF, WT1, and Huckebein are present in nonmetazoan holozoans. Of the seven metazoan KLF groups, only KLFG, represented by KLF15 in human, was found in nonmetazoans. In addition, two nonmetazoan groups of KLFs are present in Choanoflagellatea and Filasterea. WT1 could be evolutionarily the earliest among these GC/GT-box-binding families due to its sole presence in Ichthyosporea.

A genome-wide survey of genes encoding transcription factors in Japanese pearl oyster Pinctada fucata: II. Tbx, Fox, Ets, HMG, NFκB, bZIP, and C2H2 zinc fingers

To gain a better understanding of molluscan development and its relation to the evolution of their unique body plan, we performed a genomic survey of genes encoding transcription factors, such as Tbx, Fox, Ets, HMG, NFκB, bZIP, and C2H2 zinc finger proteins in the Japanese pearl oyster, Pinctada fucata. We annotated 133 transcription factor genes. Together with the orthologs of known deuterostome genes, we found several orphan genes in each class of transcription factor. Some possessed clear orthologs in other species of lophotrochozoans, while no counterpart genes were found in the deuterostomes or ecdysozoans. These observations suggest that a number of transcription factor genes are unique to lophotrochozoans, and thus additional research frontiers remain to be explored with regard to such transcription factors.

A new family of predicted Krüppel-like factor genes and pseudogenes in placental mammals

Krüppel-like factors (KLF) and specificity proteins (SP) constitute a family of zinc-finger-containing transcription factors that play important roles in a wide range of processes including differentiation and development of various tissues. The human genome possesses 17 KLF genes (KLF1-KLF17) and nine SP genes (SP1-SP9) with diverse functions. We used sequence similarity searches and gene synteny analysis to identify a new putative KLF gene/pseudogene named KLF18 that is present in most of the placental mammals with sequenced genomes. KLF18 is a chromosomal neighbor of the KLF17 gene and is likely a product of its duplication. Phylogenetic analyses revealed that mammalian predicted KLF18 proteins and KLF17 proteins experienced elevated rates of evolution and are grouped with KLF1/KLF2/KLF4 and non-mammalian KLF17. Predicted KLF18 proteins maintain conserved features in the zinc fingers of the SP/KLF family, while possessing repeats of a unique sequence motif in their N-terminal regions. No expression data have been reported for KLF18, suggesting that it either has highly restricted expression patterns and specialized functions, or could have become a pseudogene in extant placental mammals. Besides KLF18 genes/pseudogenes, we identified several KLF18-like genes such as Zfp352, Zfp352-like, and Zfp353 in the genomes of mouse and rat. These KLF18-like genes do not possess introns inside their coding regions, and gene expression data indicate that some of them may function in early embryonic development. They represent further expansions of KLF members in the murine lineage, most likely resulted from several events of retrotransposition and local gene duplication starting from an ancient spliced mRNA of KLF18.

GRG5/AES interacts with T-cell factor 4 (TCF4) and downregulates Wnt signaling in human cells and zebrafish embryos

Transcriptional control by TCF/LEF proteins is crucial in key developmental processes such as embryo polarity, tissue architecture and cell fate determination. TCFs associate with β-catenin to activate transcription in the presence of Wnt signaling, but in its absence act as repressors together with Groucho-family proteins (GRGs). TCF4 is critical in vertebrate intestinal epithelium, where TCF4-β-catenin complexes are necessary for the maintenance of a proliferative compartment, and their abnormal formation initiates tumorigenesis. However, the extent of TCF4-GRG complexes' roles in development and the mechanisms by which they repress transcription are not completely understood. Here we characterize the interaction between TCF4 and GRG5/AES, a Groucho family member whose functional relationship with TCFs has been controversial. We map the core GRG interaction region in TCF4 to a 111-amino acid fragment and show that, in contrast to other GRGs, GRG5/AES-binding specifically depends on a 4-amino acid motif (LVPQ) present only in TCF3 and some TCF4 isoforms. We further demonstrate that GRG5/AES represses Wnt-mediated transcription both in human cells and zebrafish embryos. Importantly, we provide the first evidence of an inherent repressive function of GRG5/AES in dorsal-ventral patterning during early zebrafish embryogenesis. These results improve our understanding of TCF-GRG interactions, have significant implications for models of transcriptional repression by TCF-GRG complexes, and lay the groundwork for in depth direct assessment of the potential role of Groucho-family proteins in both normal and abnormal development.

Sonic hedgehog signaling in the developing CNS where it has been and where it is going

Sonic Hedgehog (Shh) is one of three mammalian orthologs of the Hedgehog (Hh) family of secreted proteins first identified for their role in patterning the Drosophila embryo. In this review, we will highlight some of the outstanding questions regarding how Shh signaling controls embryonic development. We will mainly consider its role in the developing mammalian central nervous system (CNS) where the pathway plays a critical role in orchestrating the specification of distinct cell fates within ventral regions, a process of exquisite complexity that is necessary for the proper wiring and hence function of the mature system. Embryonic development is a process that plays out in both the spatial and the temporal dimensions, and it is becoming increasingly clear that our understanding of Shh signaling in the CNS is grounded in an appreciation for the dynamic nature of this process. In addition, any consideration of Hh signaling must by necessity include a consideration of data from many different model organisms and systems. In many cases, the extent to which insights gained from these studies are applicable to the CNS remains to be determined, yet they provide a strong framework in which to explore its role in CNS development. We will also discuss how Shh controls cell fate diversification through the regulation of patterned target gene expression in the spinal cord, a region where our understanding of the morphogenetic action of graded Shh signaling is perhaps the furthest advanced.

Ancient homeobox gene loss and the evolution of chordate brain and pharynx development: deductions from amphioxus gene expression

Homeobox genes encode a large superclass of transcription factors with widespread roles in animal development. Within chordates there are over 100 homeobox genes in the invertebrate cephalochordate amphioxus and over 200 in humans. Set against this general trend of increasing gene number in vertebrate evolution, some ancient homeobox genes that were present in the last common ancestor of chordates have been lost from vertebrates. Here, we describe the embryonic expression of four amphioxus descendants of these genes--AmphiNedxa, AmphiNedxb, AmphiMsxlx and AmphiNKx7. All four genes are expressed with a striking asymmetry about the left-right axis in the pharyngeal region of neurula embryos, mirroring the pronounced asymmetry of amphioxus embryogenesis. AmphiMsxlx and AmphiNKx7 are also transiently expressed in an anterior neural tube region destined to become the cerebral vesicle. These findings suggest significant rewiring of developmental gene regulatory networks occurred during chordate evolution, coincident with homeobox gene loss. We propose that loss of otherwise widely conserved genes is possible when these genes function in a confined role in development that is subsequently lost or significantly modified during evolution. In the case of these homeobox genes, we propose that this has occurred in relation to the evolution of the chordate pharynx and brain.

Expression and phylogenetic analysis of the zic gene family in the evolution and development of metazoans

Background

zic genes are members of the gli/glis/nkl/zic super-family of C2H2 zinc finger (ZF) transcription factors. Homologs of the zic family have been implicated in patterning neural and mesodermal tissues in bilaterians. Prior to this study, the origin of the metazoan zic gene family was unknown and expression of zic gene homologs during the development of early branching metazoans had not been investigated.

Results

Phylogenetic analyses of novel zic candidate genes identified a definitive zic homolog in the placozoan Trichoplax adhaerens, two gli/glis/nkl-like genes in the ctenophore Mnemiopsis leidyi, confirmed the presence of three gli/glis/nkl-like genes in Porifera, and confirmed the five previously identified zic genes in the cnidarian Nematostella vectensis. In the cnidarian N. vectensis, zic homologs are expressed in ectoderm and the gastrodermis (a bifunctional endomesoderm), in presumptive and developing tentacles, and in oral and sensory apical tuft ectoderm. The Capitella teleta zic homolog (Ct-zic) is detectable in a subset of the developing nervous system, the foregut, and the mesoderm associated with the segmentally repeated chaetae. Lastly, expression of gli and glis homologs in Mnemiopsis. leidyi is detected exclusively in neural cells in floor of the apical organ.

Conclusions

Based on our analyses, we propose that the zic gene family arose in the common ancestor of the Placozoa, Cnidaria and Bilateria from a gli/glis/nkl-like gene and that both ZOC and ZF-NC domains evolved prior to cnidarian-bilaterian divergence. We also conclude that zic expression in neural ectoderm and developing neurons is pervasive throughout the Metazoa and likely evolved from neural expression of an ancestral gli/glis/nkl/zic gene. zic expression in bilaterian mesoderm may be related to the expression in the gastrodermis of a cnidarian-bilaterian common ancestor.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background

The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians.

Results

In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern.

Conclusions

The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

Origin and genetic evolution of the vertebrate skeleton

The current understanding of the origin and evolution of the genetic cassette for the vertebrate skeletal system is reviewed. Molecular phylogenetic analyses of fibrillar collagen genes, which encode the main component of both cartilage and mineralized bone, suggest that genome duplications in vertebrate ancestors were essential for producing distinct collagen fibers for cartilage and mineralized bone. Several data Indicate co-expression of the ancestral copy of fibrillar collagen with the SoxE and Runx transcription factors. Therefore, the genetic cassette may have already existed in protochordate ancestors, and may operate in the development of the pharyngeal gill skeleton. Accompanied by genome duplications in vertebrate ancestors, this genetic cassette may have also been duplicated and co-opted for cartilage and bone. Subsequently, the genetic cassette for cartilage recruited novel genetic material via domain shuffling. Aggrecan, acquired by means of domain shuffling, performs an essential role in cartilage as a shock absorber. In contrast, the cassette for bone recruited new genetic material produced by tandem duplication of the SPARC/osteonectin genes. Some of the duplicated copies of SPARC/osteonectin became secretory Cabinding phosphoproteins (SCPPs) performing a central role in mineralization by regulating the calcium phosphate concentration. Comparative genome analysis revealed similar molecular evolutionary histories for the genetic cassettes for cartilage and bone, namely duplication of the ancestral genetic cassette and recruitment of novel genetic material.

A clustered set of three Sp-family genes is ancestral in the Metazoa: evidence from sequence analysis, protein domain structure, developmental expression patterns and chromosomal location

Background

The Sp-family of transcription factors are evolutionarily conserved zinc finger proteins present in many animal species. The orthology of the Sp genes in different animals is unclear and their evolutionary history is therefore controversially discussed. This is especially the case for the Sp gene buttonhead (btd) which plays a key role in head development in Drosophila melanogaster, and has been proposed to have originated by a recent gene duplication. The purpose of the presented study was to trace orthologs of btd in other insects and reconstruct the evolutionary history of the Sp genes within the metazoa.

Results

We isolated Sp genes from representatives of a holometabolous insect (Tribolium castaneum), a hemimetabolous insect (Oncopeltus fasciatus), primitively wingless hexapods (Folsomia candida and Thermobia domestica), and an amphipod crustacean (Parhyale hawaienis). We supplemented this data set with data from fully sequenced animal genomes. We performed phylogenetic sequence analysis with the result that all Sp factors fall into three monophyletic clades. These clades are also supported by protein domain structure, gene expression, and chromosomal location. We show that clear orthologs of the D. melanogaster btd gene are present even in the basal insects, and that the Sp5-related genes in the genome sequence of several deuterostomes and the basal metazoans Trichoplax adhaerens and Nematostella vectensis are also orthologs of btd.

Conclusions

All available data provide strong evidence for an ancestral cluster of three Sp-family genes as well as synteny of this Sp cluster and the Hox cluster. The ancestral Sp gene cluster already contained a Sp5/btd ortholog, which strongly suggests that btd is not the result of a recent gene duplication, but directly traces back to an ancestral gene already present in the metazoan ancestor.

Insights into the evolution of the snail superfamily from metazoan wide molecular phylogenies and expression data in annelids

BACKGROUND

An important issue concerning the evolution of duplicated genes is to understand why paralogous genes are retained in a genome even though the most likely fate for a redundant duplicated gene is nonfunctionalization and thereby its elimination. Here we study a complex superfamily generated by gene duplications, the snail related genes that play key roles during animal development. We investigate the evolutionary history of these genes by genomic, phylogenetic, and expression data studies.

RESULTS

We systematically retrieved the full complement of snail related genes in several sequenced genomes. Through phylogenetic analysis, we found that the snail superfamily is composed of three ancestral families, snail, scratchA and scratchB. Analyses of the organization of the encoded proteins point out specific molecular signatures, indicative of functional specificities for Snail, ScratchA and ScratchB proteins. We also report the presence of two snail genes in the annelid Platynereis dumerilii, which have distinct expression patterns in the developing mesoderm, nervous system, and foregut. The combined expression of these two genes is identical to that of two independently duplicated snail genes in another annelid, Capitella spI, but different aspects of the expression patterns are differentially shared among paralogs of Platynereis and Capitella.

CONCLUSION

Our study indicates that the snail and scratchB families have expanded through multiple independent gene duplications in the different bilaterian lineages, and highlights potential functional diversifications of Snail and ScratchB proteins following duplications, as, in several instances, paralogous proteins in a given species show different domain organizations. Comparisons of the expression pattern domains of the two Platynereis and Capitella snail paralogs provide evidence for independent subfunctionalization events which have occurred in these two species. We propose that the snail related genes may be especially prone to subfunctionalization, and this would explain why the snail superfamily underwent so many independent duplications leading to maintenance of functional paralogs.