The Fox genes of Branchiostoma floridae

Transcriptional regulation analysis reveals the complexity of metamorphosis in the Pacific oyster (Crassostrea gigas)

Many marine invertebrate phyla are characterized by indirect development. These animals transit from planktonic larvae to benthic spats via settlement and metamorphosis, which contributes to their adaption to the marine environment. Studying the biological process of metamorphosis is, thus, key to understanding the origin and evolution of indirect development. Although numerous studies have been conducted on the relationship between metamorphosis and the marine environment, microorganisms, and neurohormones, little is known about gene regulation network (GRN) dynamics during metamorphosis. Metamorphosis-competent pediveligers of the Pacific oyster Crassostrea gigas were assayed in this study. By assaying gene expression patterns and open chromatin region changes of different samples of larvae and spats, the dynamics of molecular regulation during metamorphosis were examined. The results indicated significantly different gene regulation networks before, during and post-metamorphosis. Genes encoding membrane-integrated receptors and those related to the remodeling of the nervous system were upregulated before the initiation of metamorphosis. Massive biogenesis, e.g., of various enzymes and structural proteins, occurred during metamorphosis as inferred from the comprehensive upregulation of the protein synthesis system post epinephrine stimulation. Hierarchical downstream gene networks were then stimulated. Some transcription factors, including homeobox, basic helix-loop-helix and nuclear receptors, showed different temporal response patterns, suggesting a complex GRN during the transition stage. Nuclear receptors, as well as their retinoid X receptor partner, may participate in the GRN controlling oyster metamorphosis, indicating an ancient role of the nuclear receptor regulation system in animal metamorphosis.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00204-y.

The Fox Gene Repertoire in the Annelid Owenia fusiformis Reveals Multiple Expansions of the foxQ2 Class in Spiralia

Fox genes are a large and conserved family of transcription factors involved in many key biological processes, including embryogenesis and body patterning. Although the role of Fox genes has been studied in an array of model systems, comprehensive comparative studies in Spiralia-a large clade of invertebrate animals including molluscs and annelids-are scarce but much needed to better understand the evolutionary history of this gene family. Here, we reconstruct and functionally characterize the Fox gene complement in the annelid Owenia fusiformis, a slow evolving species and member of the sister group to all remaining annelids. The genome of O. fusiformis contains at least a single ortholog for 20 of the 22 Fox gene classes that are ancestral to Bilateria, including an ortholog of the recently discovered foxT class. Temporal and spatial expression dynamics reveal a conserved role of Fox genes in gut formation, mesoderm patterning, and apical organ and cilia formation in Annelida and Spiralia. Moreover, we uncover an ancestral expansion of foxQ2 genes in Spiralia, represented by 11 paralogs in O. fusiformis. Notably, although all foxQ2 copies have apical expression in O. fusiformis, they show variable spatial domains and staggered temporal activation, which suggest cooperation and sub-functionalization among foxQ2 genes for the development of apical fates in this annelid. Altogether, our study informs the evolution and developmental roles of Fox genes in Annelida and Spiralia generally, providing the basis to explore how regulatory changes in Fox gene expression might have contributed to developmental and morphological diversification in Spiralia.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

Phylogenetic analysis of forkhead transcription factors in the Panarthropoda

Fox genes encode transcription factors that contain a DNA binding domain, the forkhead domain, and are known from diverse animal species. The exact homology of the Fox genes of different species is debated and this makes inferences about the evolution of the Fox genes, and their duplications and losses difficult. We have performed phylogenetic analyses of the Fox gene complements of 32 panarthropod species. Our results confirm an ancestral complement of FoxA, FoxB, FoxC, FoxD, FoxF, FoxG, FoxJ1, FoxJ2/3, FoxK, FoxL1, FoxL2, FoxN1/4, FoxN2/3, FoxO, FoxP, and FoxQ2 in the Arthropoda, and additionally FoxH and FoxQ1 in the Panarthropoda (including tardigrades and onychophorans). We identify a novel Fox gene sub-family, that we designate as FoxT that includes two genes in Drosophila melanogaster, Circadianly Regulated Gene (Crg-1) and forkhead domain 3F (fd3F). In a very recent paper, the same new Fox gene sub-family was identified in insects (Lin et al. 2021). Our analysis confirms the presence of FoxT and shows that its members are present throughout Panarthropoda. We show that the hitherto unclassified gene CG32006 from the fly Drosophila melanogaster belongs to FoxJ1. We also detect gene losses: FoxE and FoxM were lost already in the panarthropod ancestor, whereas the loss of FoxH occurred in the arthropod ancestor. Finally, we find an ortholog of FoxQ1 in the bark scorpion Centruroides sculpturatus, confirmed not only by phylogenetic analysis, but also by forming an evolutionarily conserved gene cluster with FoxF, FoxC, and FoxL1. This suggests that FoxQ1 belongs to the ancestral Fox gene complement in panarthropods and also in chelicerates, but has been lost at the base of the mandibulate arthropods.

Amphioxus muscle transcriptomes reveal vertebrate-like myoblast fusion genes and a highly conserved role of insulin signalling in the metabolism of muscle

BACKGROUND

The formation and functioning of muscles are fundamental aspects of animal biology, and the evolution of 'muscle genes' is central to our understanding of this tissue. Feeding-fasting-refeeding experiments have been widely used to assess muscle cellular and metabolic responses to nutrition. Though these studies have focused on vertebrate models and only a few invertebrate systems, they have found similar processes are involved in muscle degradation and maintenance. Motivation for these studies stems from interest in diseases whose pathologies involve muscle atrophy, a symptom also triggered by fasting, as well as commercial interest in the muscle mass of animals kept for consumption. Experimentally modelling atrophy by manipulating nutritional state causes muscle mass to be depleted during starvation and replenished with refeeding so that the genetic mechanisms controlling muscle growth and degradation can be understood.

RESULTS

Using amphioxus, the earliest branching chordate lineage, we address the gap in previous work stemming from comparisons between distantly related vertebrate and invertebrate models. Our amphioxus feeding-fasting-refeeding muscle transcriptomes reveal a highly conserved myogenic program and that the pro-orthologues of many vertebrate myoblast fusion genes were present in the ancestral chordate, despite these invertebrate chordates having unfused mononucleate myocytes. We found that genes differentially expressed between fed and fasted amphioxus were orthologous to the genes that respond to nutritional state in vertebrates. This response is driven in a large part by the highly conserved IGF/Akt/FOXO pathway, where depleted nutrient levels result in activation of FOXO, a transcription factor with many autophagy-related gene targets.

CONCLUSION

Reconstruction of these gene networks and pathways in amphioxus muscle provides a key point of comparison between the distantly related groups assessed thus far, significantly refining the reconstruction of the ancestral state for chordate myoblast fusion genes and identifying the extensive role of duplicated genes in the IGF/Akt/FOXO pathway across animals. Our study elucidates the evolutionary trajectory of muscle genes as they relate to the increased complexity of vertebrate muscles and muscle development.

A Preliminary Single-Cell RNA-Seq Analysis of Embryonic Cells That Express Brachyury in the Amphioxus, Branchiostoma japonicum

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.

Analysis of Fox genes in Schmidtea mediterranea reveals new families and a conserved role of Smed-foxO in controlling cell death

The forkhead box (Fox) genes encode transcription factors that control several key aspects of development. Present in the ancestor of all eukaryotes, Fox genes underwent several duplications followed by loss and diversification events that gave rise to the current 25 families. However, few Fox members have been identified from the Lophotrochozoa clade, and specifically from planarians, which are a unique model for understanding development, due to the striking plasticity of the adult. The aim of this study was to identify and perform evolutionary and functional studies of the Fox genes of lophotrochozoan species and, specifically, of the planarian Schmidtea mediterranea. Generating a pipeline for identifying Forkhead domains and using phylogenetics allowed us the phylogenetic reconstruction of Fox genes. We corrected the annotation for misannotated genes and uncovered a new family, the QD, present in all metazoans. According to the new phylogeny, the 27 Fox genes found in Schmidtea mediterranea were classified into 12 families. In Platyhelminthes, family losses were accompanied by extensive gene diversification and the appearance of specific families, the A(P) and N(P). Among the newly identified planarian Fox genes, we found a single copy of foxO, which shows an evolutionary conserved role in controlling cell death.

The genetic factors of bilaterian evolution

The Cambrian explosion was a unique animal radiation ~540 million years ago that produced the full range of body plans across bilaterians. The genetic mechanisms underlying these events are unknown, leaving a fundamental question in evolutionary biology unanswered. Using large-scale comparative genomics and advanced orthology evaluation techniques, we identified 157 bilaterian-specific genes. They include the entire Nodal pathway, a key regulator of mesoderm development and left-right axis specification; components for nervous system development, including a suite of G-protein-coupled receptors that control physiology and behaviour, the Robo-Slit midline repulsion system, and the neurotrophin signalling system; a high number of zinc finger transcription factors; and novel factors that previously escaped attention. Contradicting the current view, our study reveals that genes with bilaterian origin are robustly associated with key features in extant bilaterians, suggesting a causal relationship.

Identification and expression profiles of Fox transcription factors in the Yesso scallop (Patinopecten yessoensis)

The forkhead box (Fox) gene family is a family of transcription factors that play important roles in a variety of biological processes in vertebrates, including early development and cell proliferation and differentiation. However, at present, studies on the mollusk Fox family are relatively lacking. In the present study, the Fox gene family of the Yesso scallop (Patinopecten yessoensis) was systematically identified. In addition, the expression profiles of the Fox gene family in early development and adult tissues were analyzed. The results showed that there were 26 Fox genes in P. yessoensis. Of the 26 genes, 24 belonged to 20 subfamilies. The Fox genes belonging to the I, Q1, R and S subfamilies were absent in P. yessoensis. The other 2 genes formed 2 independent clades with the Fox genes of other mollusks and protostomes. They might be new members of the Fox family and were named FoxY and FoxZ. P. yessoensis contained a FoxC-FoxL1 gene cluster similar in structure to that of Branchiostoma floridae, suggesting that the cluster might already exist in the ancestors of bilaterally symmetrical animals. The gene expression analysis of Fox showed that most of the genes were continuously expressed in multiple stages of early development, suggesting that Fox genes might be widely involved in the regulation of embryo and larval development of P. yessoensis. Nine Fox genes were specifically expressed in certain tissues, such as the nerve ganglia, foot, ovary, testis, and gills. For the 9 genes that were differentially expressed between the testis and ovary, their expression levels were analyzed during the 4 developmental stages of gonads. The results showed that FoxL2, FoxE and FoxY were highly expressed in the ovary during all developmental stages, while FoxZ was highly expressed in the testis during all developmental stages. The results suggested that these genes might play an important role in sex maintenance or gametogenesis. The present study could provide a reference for evolutionary and functional studies of the Fox family in metazoans.

cis-Regulatory analysis for later phase of anterior neuroectoderm-specific foxQ2 expression in sea urchin embryos

The specification of anterior neuroectoderm is controlled by a highly conserved molecular mechanism in bilaterians. A forkhead family gene, foxQ2, is known to be one of the pivotal regulators, which is zygotically expressed in this region during embryogenesis of a broad range of bilaterians. However, what controls the expression of this essential factor has remained unclear to date. To reveal the regulatory mechanism of foxQ2, we performed cis-regulatory analysis of two foxQ2 genes, foxQ2a and foxQ2b, in a sea urchin Hemicentrotus pulcherrimus. In sea urchin embryos, foxQ2 is initially expressed in the entire animal hemisphere and subsequently shows narrower expression restricted to the anterior pole region. In this study, as a first step to understand the foxQ2 regulation, we focused on the later restricted expression and analyzed the upstream cis-regulatory sequences of foxQ2a and foxQ2b by using the constructs fused to short half-life green fluorescent protein. Based on deletion and mutation analyses of both foxQ2, we identified each of the five regulatory sequences, which were 4-9 bp long. Neither of the regulatory sequences contains any motifs for ectopic activation or spatial repression, suggesting that later mRNA localization is regulated in situ. We also suggest that the three amino acid loop extension-class homeobox gene Meis is involved in the maintenance of foxQ2b, the expression of which is dominant during embryogenesis.

Temporal and spatial expression of the Fox gene family in the Leech Helobdella austinensis

The Forkhead box (Fox) gene family is an evolutionarily ancient gene family named after the Drosophila melanogaster forkhead gene (fkh). Fox genes are highly conserved transcription factors critical for embryogenesis and carcinogenesis. In the current study, we report a whole-genome survey of Fox genes and their expression patterns in the leech Helobdella austienesis. Phylogenetic analysis suggests that some Fox genes of leeches are correlated with other Lophotrochozoa and vertebrate Fox genes. Here we have performed semiquantitative reverse transcription polymerase chain reaction and whole-mount in situ hybridization of Fox genes throughout the embryonic development of H. austinensis. We found that each one of the leech Fox genes (FoxA1, FoxA3, FoxC, FoxL2, FoxO1, and FoxO2) is expressed in a specific set of cells or tissue type. From Stages 9-11, Hau-FoxA1 was expressed in the foregut of the anterior region, and Hau-FoxL2 was expressed in mesodermal muscle fiber. Hau-FoxA3 was temporally expressed in the ventral neuroectoderm. At Stage 11, Hau-FoxC was expressed in the foregut. Hau-FoxO genes have a ubiquitous expression. Our results provide more insight on the evolutionary linkage and role of the Fox gene function in Bilateria.

Duplicated FoxA genes in the leech Helobdella: Insights into the evolution of direct development in clitellate annelids

BACKGROUND

As an adaptation to the land, the clitellate annelid had reorganized its embryogenesis to develop "directly" without the ancestral planktonic larval stage. To study the evolution of gut development in the directly developing clitellates, we characterized the expression pattern of the conserved gut gene, FoxA, in the embryonic development of the leech.

RESULTS

The leech has three FoxA paralogs. Hau-FoxA1 is first expressed in a subset of endoderm cells and then in the foregut and the midgut. Hau-FoxA2 is expressed in the stomodeum, which is secondarily derived from the anterior ectoderm in the clitellates rather than the tissue around the blastopore, the ancestral site of mouth formation in Phylum Annelida. Hau-FoxA3 is expressed during the morphogenesis of segmental ganglia from the ectodermal teloblast lineages, a clitellate-specific trait. Hau-FoxA1 and Hau-FoxA2 are also expressed during the morphogenesis of the leech-specific front sucker.

CONCLUSIONS

The expression patterns suggested that Hau-FoxA1 carries out most of the conserved function in the endoderm and gut development, while the other two duplicates appear to have evolved unique novel functions in the directly developing clitellate embryos. Therefore, neofunctionalization and co-option of FoxA might have made a significant contribution to the evolution of direct development in Clitellata. Developmental Dynamics 247:763-778, 2018. © 2018 Wiley Periodicals, Inc.

Unipotent progenitors contribute to the generation of sensory cell types in the nervous system of the cnidarian Nematostella vectensis

Nervous systems often consist of a large number of different types of neurons which are generated from neural stem and progenitor cells by a series of symmetric and asymmetric divisions. The origin and early evolution of these neural progenitor systems is not well understood. Here we use a cnidarian model organism, Nematostella vectensis, to gain insight into the generation of neural cell type diversity in a non-bilaterian animal. We identify NvFoxQ2d as a transcription factor that is expressed in a population of spatially restricted, proliferating ectodermal cells that are derived from NvSoxB(2)-expressing neural progenitor cells. Using a transgenic reporter line we show that the NvFoxQ2d cells undergo a terminal, symmetric division to generate a morphologically homogeneous population of putative sensory cells. The abundance of these cells, but not their proliferation status is affected by treatment with the γ-secretase inhibitor DAPT, suggesting regulation by Notch signalling. Our data suggest that intermediate progenitor cells and symmetric divisions contribute to the formation of the seemingly simple nervous system of a sea anemone.

A key role for foxQ2 in anterior head and central brain patterning in insects

Anterior patterning of animals is based on a set of highly conserved transcription factors but the interactions within the protostome anterior gene regulatory network (aGRN) remain enigmatic. Here, we identify the red flour beetle Tribolium castaneum ortholog of foxQ2 (Tc-foxQ2) as a novel upstream component of the aGRN. It is required for the development of the labrum and higher order brain structures, namely the central complex and the mushroom bodies. We reveal Tc-foxQ2 interactions by RNAi and heat shock-mediated misexpression. Surprisingly, Tc-foxQ2 and Tc-six3 mutually activate each other, forming a novel regulatory module at the top of the aGRN. Comparisons of our results with those of sea urchins and cnidarians suggest that foxQ2 has acquired more upstream functions in the aGRN during protostome evolution. Our findings expand the knowledge on foxQ2 gene function to include essential roles in epidermal development and central brain patterning.

Divergence of ectodermal and mesodermal gene regulatory network linkages in early development of sea urchins

Developmental gene regulatory networks (GRNs) are assemblages of gene regulatory interactions that direct ontogeny of animal body plans. Studies of GRNs operating in the early development of euechinoid sea urchins have revealed that little appreciable change has occurred since their divergence ∼90 million years ago (mya). These observations suggest that strong conservation of GRN architecture was maintained in early development of the sea urchin lineage. Testing whether this holds for all sea urchins necessitates comparative analyses of echinoid taxa that diverged deeper in geological time. Recent studies highlighted extensive divergence of skeletogenic mesoderm specification in the sister clade of euechinoids, the cidaroids, suggesting that comparative analyses of cidaroid GRN architecture may confer a greater understanding of the evolutionary dynamics of developmental GRNs. Here I report spatiotemporal patterning of 55 regulatory genes and perturbation analyses of key regulatory genes involved in euechinoid oral-aboral patterning of nonskeletogenic mesodermal and ectodermal domains in early development of the cidaroid Eucidaris tribuloides These results indicate that developmental GRNs directing mesodermal and ectodermal specification have undergone marked alterations since the divergence of cidaroids and euechinoids. Notably, statistical and clustering analyses of echinoid temporal gene expression datasets indicate that regulation of mesodermal genes has diverged more markedly than regulation of ectodermal genes. Although research on indirect-developing euechinoid sea urchins suggests strong conservation of GRN circuitry during early embryogenesis, this study indicates that since the divergence of cidaroids and euechinoids, developmental GRNs have undergone significant, cell type-biased alterations.

Foxl2 and Its Relatives Are Evolutionary Conserved Players in Gonadal Sex Differentiation

Foxl2 is a member of the large family of Forkhead Box (Fox) domain transcription factors. It emerged during the last 15 years as a key player in ovarian differentiation and oogenesis in vertebrates and especially mammals. This review focuses on Foxl2 genes in light of recent findings on their evolution, expression, and implication in sex differentiation in animals in general. Homologs of Foxl2 and its paralog Foxl3 are found in all metazoans, but their gene evolution is complex, with multiple gains and losses following successive whole genome duplication events in vertebrates. This review aims to decipher the evolutionary forces that drove Foxl2/3 gene specialization through sub- and neo-functionalization during evolution. Expression data in metazoans suggests that Foxl2/3 progressively acquired a role in both somatic and germ cell gonad differentiation and that a certain degree of sub-functionalization occurred after its duplication in vertebrates. This generated a scenario where Foxl2 is predominantly expressed in ovarian somatic cells and Foxl3 in male germ cells. To support this hypothesis, we provide original results showing that in the pea aphid (insects) foxl2/3 is predominantly expressed in sexual females and showing that in bovine ovaries FOXL2 is specifically expressed in granulosa cells. Overall, current results suggest that Foxl2 and Foxl3 are evolutionarily conserved players involved in somatic and germinal differentiation of gonadal sex.

Molecular conservation of metazoan gut formation: evidence from expression of endomesoderm genes in Capitella teleta (Annelida)

BACKGROUND

Metazoan digestive systems develop from derivatives of ectoderm, endoderm and mesoderm, and vary in the relative contribution of each germ layer across taxa and between gut regions. In a small number of well-studied model systems, gene regulatory networks specify endoderm and mesoderm of the gut within a bipotential germ layer precursor, the endomesoderm. Few studies have examined expression of endomesoderm genes outside of those models, and thus, it is unknown whether molecular specification of gut formation is broadly conserved. In this study, we utilize a sequenced genome and comprehensive fate map to correlate the expression patterns of six transcription factors with embryonic germ layers and gut subregions during early development in Capitella teleta.

RESULTS

The genome of C. teleta contains the five core genes of the sea urchin endomesoderm specification network. Here, we extend a previous study and characterize expression patterns of three network orthologs and three additional genes by in situ hybridization during cleavage and gastrulation stages and during formation of distinct gut subregions. In cleavage stage embryos, Ct-otx, Ct-blimp1, Ct-bra and Ct-nkx2.1a are expressed in all four macromeres, the endoderm precursors. Ct-otx, Ct-blimp1, and Ct-nkx2.1a are also expressed in presumptive endoderm of gastrulae and later during midgut development. Additional gut-specific expression patterns include Ct-otx, Ct-bra, Ct-foxAB and Ct-gsc in oral ectoderm; Ct-otx, Ct-blimp1, Ct-bra and Ct-nkx2.1a in the foregut; and both Ct-bra and Ct-nkx2.1a in the hindgut.

CONCLUSIONS

Identification of core sea urchin endomesoderm genes in C. teleta indicates they are present in all three bilaterian superclades. Expression of Ct-otx, Ct-blimp1 and Ct-bra, combined with previously published Ct-foxA and Ct-gataB1 patterns, provide the most comprehensive comparison of these five orthologs from a single species within Spiralia. Each ortholog is likely involved in endoderm specification and midgut development, and several may be essential for establishment of the oral ectoderm, foregut and hindgut, including specification of ectodermal and mesodermal contributions. When the five core genes are compared across the Metazoa, their conserved expression patterns suggest that 'gut gene' networks evolved to specify distinct digestive system subregions, regardless of species-specific differences in gut architecture or germ layer contributions within each subregion.

Genome-wide identification, phylogeny, and gonadal expression of fox genes in Nile tilapia, Oreochromis niloticus

The fox genes play important roles in various biological processes, including sexual development. In the present study, we isolated 65 fox genes, belonging to 18 subfamilies named A-R, from Nile tilapia through genome-wide screening. Twenty-four of them have two or three (foxm1) copies. Furthermore, 16, 25, 68, and 45 fox members were isolated from nematodes, protochordates, teleosts, and tetrapods, respectively. Phylogenetic analyses indicated fox gene family had undergone three expansions parallel to the three rounds of genome duplication during evolution. We also analyzed the clustered fox genes and found that apparent linkage duplication existed in teleosts, which further supported fish-specific genome duplication hypothesis. In addition, species- and lineage-specific duplication is another reason for fox gene family expansion. Based on the four pairs of XX and XY gonadal transcriptome data from four critical developmental stages, we analyzed the expression profile of all fox genes and identified sexually dimorphic fox genes at each stage. All fox genes were detected in gonads, with 15 of them at the background expression level (total read per kb per million reads, RPKM < 10), 29 at moderate expression level (10 < total RPKM < 100), and 21 at high expression level (total RPKM > 100). There are 27, 24, 28, and 9 sexually dimorphic fox genes at 5, 30, 90, and 180 days after hatching (dah), respectively. foxq1a, foxf1, foxr1, and foxr1 were identified as the most differentially expressed genes at each stage. foxl2 was characterized as XX-dominant gene, while foxd5, foxi3, foxn3, foxj1a, foxj3b, and foxo6b were characterized as XY-dominant genes. qPCR and in situ hybridization of foxh1 and foxj1a were performed to confirm the expression profiles and to validate the transcriptome data. Our results suggest that fox genes might play important roles in sex determination and gonadal development in teleosts.

Genome-wide survey and expression analysis of the bHLH-PAS genes in the amphioxus Branchiostoma floridae reveal both conserved and diverged expression patterns between cephalochordates and vertebrates

Background

The bHLH-PAS transcription factors are found in both protostomes and deuterostomes. They are involved in many developmental and physiological processes, including regional differentiation of the central nervous system, tube-formation, hypoxia signaling, aromatic hydrocarbon sensing, and circadian rhythm regulation. To understand the evolution of these genes in chordates, we analyzed the bHLH-PAS genes of the basal chordate amphioxus (Branchiostoma floridae).

Results

From the amphioxus draft genome database, we identified ten bHLH-PAS genes, nine of which could be assigned to known orthologous families. The tenth bHLH-PAS gene could not be assigned confidently to any known bHLH family; however, phylogenetic analysis clustered this gene with arthropod Met family genes and two spiralian bHLH-PAS-containing sequences, suggesting that they may share the same ancestry. We examined temporal and spatial expression patterns of these bHLH-PAS genes in developing amphioxus embryos. We found that BfArnt, BfNcoa, BfSim, and BfHifα were expressed in the central nervous system in patterns similar to those of their vertebrate homologs, suggesting that their functions may be conserved. By contrast, the amphioxus BfAhr and BfNpas4 had expression patterns distinct from those in vertebrates. These results imply that there were changes in gene regulation after the divergence of cephalochordates and vertebrates.

Conclusions

We have identified ten bHLH-PAS genes from the amphioxus genome and determined the embryonic expression profiles for these genes. In addition to the nine currently recognized bHLH-PAS families, our survey suggests that the BfbHLHPAS-orphan gene along with arthropod Met genes and the newly identified spiralian bHLH-PAS-containing sequences represent an ancient group of genes that were lost in the vertebrate lineage. In a comparison with the expression patterns of the vertebrate bHLH-PAS paralogs, which are the result of whole-genome duplication, we found that although several members seem to retain conserved expression patterns during chordate evolution, many duplicated paralogs may have undergone subfunctionalization and neofunctionalization in the vertebrate lineage. In addition, our survey of amphioxus bHLH-PAS gene models from genome browser with experimentally verified cDNA sequences calls into question the accuracy of the current in silico gene annotation of the B. floridae genome.

The Fox/Forkhead transcription factor family of the hemichordate Saccoglossus kowalevskii

Background

The Fox gene family is a large family of transcription factors that arose early in organismal evolution dating back to at least the common ancestor of metazoans and fungi. They are key components of many gene regulatory networks essential for embryonic development. Although much is known about the role of Fox genes during vertebrate development, comprehensive comparative studies outside vertebrates are sparse. We have characterized the Fox transcription factor gene family from the genome of the enteropneust hemichordate Saccoglossus kowalevskii, including phylogenetic analysis, genomic organization, and expression analysis during early development. Hemichordates are a sister group to echinoderms, closely related to chordates and are a key group for tracing the evolution of gene regulatory mechanisms likely to have been important in the diversification of the deuterostome phyla.

Results

Of the 22 Fox gene families that were likely present in the last common ancestor of all deuterostomes, S. kowalevskii has a single ortholog of each group except FoxH, which we were unable to detect, and FoxQ2, which has three paralogs. A phylogenetic analysis of the FoxQ2 family identified an ancestral duplication in the FoxQ2 lineage at the base of the bilaterians. The expression analyses of all 23 Fox genes of S. kowalevskii provide insights into the evolution of components of the regulatory networks for the development of pharyngeal gill slits (foxC, foxL1, and foxI), mesoderm patterning (foxD, foxF, foxG), hindgut development (foxD, foxI), cilia formation (foxJ1), and patterning of the embryonic apical territory (foxQ2).

Conclusions

Comparisons of our results with data from echinoderms, chordates, and other bilaterians help to develop hypotheses about the developmental roles of Fox genes that likely characterized ancestral deuterostomes and bilaterians, and more recent clade-specific innovations.

Larval body patterning and apical organs are conserved in animal evolution

Background

Planktonic ciliated larvae are characteristic for the life cycle of marine invertebrates. Their most prominent feature is the apical organ harboring sensory cells and neurons of largely undetermined function. An elucidation of the relationships between various forms of primary larvae and apical organs is key to understanding the evolution of animal life cycles. These relationships have remained enigmatic due to the scarcity of comparative molecular data.

Results

To compare apical organs and larval body patterning, we have studied regionalization of the episphere, the upper hemisphere of the trochophore larva of the marine annelid Platynereis dumerilii. We examined the spatial distribution of transcription factors and of Wnt signaling components previously implicated in anterior neural development. Pharmacological activation of Wnt signaling with Gsk3β antagonists abolishes expression of apical markers, consistent with a repressive role of Wnt signaling in the specification of apical tissue. We refer to this Wnt-sensitive, six3- and foxq2-expressing part of the episphere as the ‘apical plate’. We also unraveled a molecular signature of the apical organ - devoid of six3 but expressing foxj, irx, nkx3 and hox - that is shared with other marine phyla including cnidarians. Finally, we characterized the cell types that form part of the apical organ by profiling by image registration, which allows parallel expression profiling of multiple cells. Besides the hox-expressing apical tuft cells, this revealed the presence of putative light- and mechanosensory as well as multiple peptidergic cell types that we compared to apical organ cell types of other animal phyla.

Conclusions

The similar formation of a six3+, foxq2+ apical plate, sensitive to Wnt activity and with an apical tuft in its six3-free center, is most parsimoniously explained by evolutionary conservation. We propose that a simple apical organ - comprising an apical tuft and a basal plexus innervated by sensory-neurosecretory apical plate cells - was present in the last common ancestors of cnidarians and bilaterians. One of its ancient functions would have been the control of metamorphosis. Various types of apical plate cells would then have subsequently been added to the apical organ in the divergent bilaterian lineages. Our findings support an ancient and common origin of primary ciliated larvae.

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.

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.

Evolution of new characters after whole genome duplications: insights from amphioxus

Additional copies of genes resulting from two whole genome duplications at the base of the vertebrates have been suggested as enabling the evolution of vertebrate-specific structures such as neural crest, a midbrain/hindbrain organizer and neurogenic placodes. These structures, however, did not evolve entirely de novo, but arose from tissues already present in an ancestral chordate. This review discusses the evolutionary history of co-option of old genes for new roles in vertebrate development as well as the relative contributions of changes in cis-regulation and in protein structure. Particular examples are the FoxD, FGF8/17/18 and Pax2/5/8 genes. Comparisons with invertebrate chordates (amphioxus and tunicates) paint a complex picture with co-option of genes into new structures occurring both after and before the whole genome duplications. In addition, while cis-regulatory changes are likely of primary importance in evolution of vertebrate-specific structures, changes in protein structure including alternative splicing are non-trivial.

A genome-wide view of transcription factor gene diversity in chordate evolution: less gene loss in amphioxus?

Previous studies of gene diversity in the homeobox superclass have shown that the Florida amphioxus Branchiostoma floridae has undergone remarkably little gene family loss. Here we use a combined BLAST and HMM search strategy to assess the family level diversity of four other transcription factor superclasses: the Paired/Pax genes, Tbx genes, Fox genes and Sox genes. We apply this across genomes from five chordate taxa, including B. floridae and Ciona intestinalis, plus two outgroup taxa. Our results show scattered gene family loss. However, as also found for homeobox genes, B. floridae has retained all ancient Pax, Tbx, Fox and Sox gene families that were present in the common ancestor of living chordates. We conclude that, at least in terms of transcription factor gene complexity, the genome of amphioxus has experienced remarkable stasis compared to the genomes of other chordates.

Sequencing and analysis of the Mediterranean amphioxus (Branchiostoma lanceolatum) transcriptome

Background

The basally divergent phylogenetic position of amphioxus (Cephalochordata), as well as its conserved morphology, development and genetics, make it the best proxy for the chordate ancestor. Particularly, studies using the amphioxus model help our understanding of vertebrate evolution and development. Thus, interest for the amphioxus model led to the characterization of both the transcriptome and complete genome sequence of the American species, Branchiostoma floridae. However, recent technical improvements allowing induction of spawning in the laboratory during the breeding season on a daily basis with the Mediterranean species Branchiostoma lanceolatum have encouraged European Evo-Devo researchers to adopt this species as a model even though no genomic or transcriptomic data have been available. To fill this need we used the pyrosequencing method to characterize the B. lanceolatum transcriptome and then compared our results with the published transcriptome of B. floridae.

Results

Starting with total RNA from nine different developmental stages of B. lanceolatum, a normalized cDNA library was constructed and sequenced on Roche GS FLX (Titanium mode). Around 1.4 million of reads were produced and assembled into 70,530 contigs (average length of 490 bp). Overall 37% of the assembled sequences were annotated by BlastX and their Gene Ontology terms were determined. These results were then compared to genomic and transcriptomic data of B. floridae to assess similarities and specificities of each species.

Conclusion

We obtained a high-quality amphioxus (B. lanceolatum) reference transcriptome using a high throughput sequencing approach. We found that 83% of the predicted genes in the B. floridae complete genome sequence are also found in the B. lanceolatum transcriptome, while only 41% were found in the B. floridae transcriptome obtained with traditional Sanger based sequencing. Therefore, given the high degree of sequence conservation between different amphioxus species, this set of ESTs may now be used as the reference transcriptome for the Branchiostoma genus.

EST and transcriptome analysis of cephalochordate amphioxus--past, present and future

The cephalochordates, commonly known as amphioxus or lancelets, are now considered the most basal chordate group, and the studies of these organisms therefore offer important insights into various levels of evolutionary biology. In the past two decades, the investigation of amphioxus developmental biology has provided key knowledge for understanding the basic patterning mechanisms of chordates. Comparative genome studies of vertebrates and amphioxus have uncovered clear evidence supporting the hypothesis of two-round whole-genome duplication thought to have occurred early in vertebrate evolution and have shed light on the evolution of morphological novelties in the complex vertebrate body plan. Complementary to the amphioxus genome-sequencing project, a large collection of expressed sequence tags (ESTs) has been generated for amphioxus in recent years; this valuable collection represents a rich resource for gene discovery, expression profiling and molecular developmental studies in the amphioxus model. Here, we review previous EST analyses and available cDNA resources in amphioxus and discuss their value for use in evolutionary and developmental studies. We also discuss the potential advantages of applying high-throughput, next-generation sequencing (NGS) technologies to the field of amphioxus research.

Nodal-dependent mesendoderm specification requires the combinatorial activities of FoxH1 and Eomesodermin

Vertebrate mesendoderm specification requires the Nodal signaling pathway and its transcriptional effector FoxH1. However, loss of FoxH1 in several species does not reliably cause the full range of loss-of-Nodal phenotypes, indicating that Nodal signals through additional transcription factors during early development. We investigated the FoxH1-dependent and -independent roles of Nodal signaling during mesendoderm patterning using a novel recessive zebrafish FoxH1 mutation called midway, which produces a C-terminally truncated FoxH1 protein lacking the Smad-interaction domain but retaining DNA–binding capability. Using a combination of gel shift assays, Nodal overexpression experiments, and genetic epistasis analyses, we demonstrate that midway more accurately represents a complete loss of FoxH1-dependent Nodal signaling than the existing zebrafish FoxH1 mutant schmalspur. Maternal-zygotic midway mutants lack notochords, in agreement with FoxH1 loss in other organisms, but retain near wild-type expression of markers of endoderm and various nonaxial mesoderm fates, including paraxial and intermediate mesoderm and blood precursors. We found that the activity of the T-box transcription factor Eomesodermin accounts for specification of these tissues in midway embryos. Inhibition of Eomesodermin in midway mutants severely reduces the specification of these tissues and effectively phenocopies the defects seen upon complete loss of Nodal signaling. Our results indicate that the specific combinations of transcription factors available for signal transduction play critical and separable roles in determining Nodal pathway output during mesendoderm patterning. Our findings also offer novel insights into the co-evolution of the Nodal signaling pathway, the notochord specification program, and the chordate branch of the deuterostome family of animals.

Multiple signaling pathways function combinatorially to form and pattern the primary tissue layers of almost all organisms, by interacting with each other and by utilizing different pathway components to perform specific roles. Here we investigated the combinatorial aspects of the Nodal signaling pathway, which is essential for proper induction of mesoderm and endoderm in vertebrates. We identified a new mutation in the zebrafish FoxH1 gene, which encodes a Nodal pathway transcription factor, a protein that responds to Nodal signals to carry out the pathway's cellular functions by regulating target gene expression. Using this mutation, we determined that FoxH1 acts in a combinatorial fashion with two other transcription factors, called Mixer and Eomesodermin, to carry out all roles of the Nodal pathway during early development. Through genetic manipulation, we were able to identify the discrete functions regulated by different combinations of these three transcription factors. Our results indicate that the availability of specific Nodal-responsive transcription factors dictates the functions of the Nodal pathway in specific areas of the developing embryo. Our work also provides evidence that the FoxH1 family of transcription factors evolved concomitantly, and perhaps causally, with the chordate branch of animals, to which all vertebrates including humans belong.

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.

Alternative splicing and gene duplication in the evolution of the FoxP gene subfamily

The FoxP gene subfamily of transcription factors is defined by its characteristic 110 amino acid long DNA-binding forkhead domain and plays essential roles in vertebrate biology. Its four members, FoxP1–P4, have been extensively characterized functionally. FoxP1, FoxP2, and FoxP4 are involved in lung, heart, gut, and central nervous system (CNS) development. FoxP3 is necessary and sufficient for the specification of regulatory T cells (Tregs) of the adaptive immune system.

In Drosophila melanogaster, in silico predictions identify one unique FoxP subfamily gene member (CG16899) with no described function. We characterized this gene and established that it generates by alternative splicing two isoforms that differ in the forkhead DNA-binding domain. In D. melanogaster, both isoforms are expressed in the embryonic CNS, but in hemocytes, only isoform A is expressed, hinting to a putative modulation through alternative splicing of FoxP1 function in immunity and/or other hemocyte-dependent processes. Furthermore, we show that in vertebrates, this novel alternative splicing pattern is conserved for FoxP1. In mice, this new FoxP1 isoform is expressed in brain, liver, heart, testes, thymus, and macrophages (equivalent in function to hemocytes). This alternative splicing pattern has arisen at the base of the Bilateria, probably through exon tandem duplication. Moreover, our phylogenetic analysis suggests that in vertebrates, FoxP1 is more related to the FoxP gene ancestral form and the other three paralogues, originated through serial duplications, which only retained one of the alternative exons. Also, the newly described isoform differs from the other in amino acids critical for DNA-binding specificity. The integrity of its fold is maintained, but the molecule has lost the direct hydrogen bonding to DNA bases leading to a putatively lower specificity and possibly affinity toward DNA.

With the present comparative study, through the integration of experimental and in silico studies of the FoxP gene subfamily across the animal kingdom, we establish a new model for the FoxP gene in invertebrates and for the vertebrate FoxP1 paralogue. Furthermore, we present a scenario for the structural evolution of this gene class and reveal new previously unsuspected levels of regulation for FoxP1 in the vertebrate system.

Clustered Fox genes in lophotrochozoans and the evolution of the bilaterian Fox gene cluster

FoxC, FoxF, FoxL1 and FoxQ1 genes have been shown to be clustered in some animal genomes, with mesendodermal expression hypothesised as a selective force maintaining cluster integrity. Hypotheses are, however, constrained by a lack of data from the Lophotrochozoa. Here we characterise members of the FoxC, FoxF, FoxL1 and FoxQ1 families from the annelid Capitella teleta and the molluscs Lottia gigantea and Patella vulgata. We cloned FoxC, FoxF, FoxL1 and FoxQ1 genes from C. teleta, and FoxC, FoxF and FoxL1 genes from P. vulgata, and established their expression during development. We also examined their genomic organisation in C. teleta and L. gigantea, and investigated local syntenic relationships. Our results show mesodermal and anterior gut expression is a common feature of these genes in lophotrochozoans. In L. gigantea FoxC, FoxF and FoxL1 are closely linked, while in C. teleta Ct-foxC and Ct-foxL1 are closely linked, with Ct-foxF and Ct-foxQ1 on different scaffolds. Adjacent to these genes there is limited evidence of local synteny. This demonstrates conservation of genomic organisation and expression of these genes can be traced in all three bilaterian Superphyla. These data are evaluated against competing theories for the long-term maintenance of gene clusters.

Evolutionary genomics of the Fox genes: origin of gene families and the ancestry of gene clusters

Over the past decade genomic approaches have begun to revolutionise the study of animal diversity. In particular, genome sequencing programmes have spread beyond the traditional model species to encompass an increasing diversity of animals from many different phyla, as well as unicellular eukaryotes that are closely related to the animals. Whole genome sequences allow researchers to establish, with reasonable confidence, the full complement of any particular family of genes in a genome. Comparison of gene complements from appropriate genomes can reveal the evolutionary history of gene families, indicating when both gene diversification and gene loss have occurred. More than that, however, assembled genomes allow the genomic environment in which individual genes are found to be analysed and compared between species. This can reveal how gene diversification occurred. Here, we focus on the Fox genes, drawing from multiple animal genomes to develop an evolutionary framework explaining the timing and mechanism of origin of the diversity of animal Fox genes. Ancient linkages between genes are a prominent feature of the Fox genes, depicting a history of gene clusters, some of which may be relevant to understanding Fox gene function.

It's a long way from amphioxus: descendants of the earliest chordate

The origin of chordates and the consequent genesis of vertebrates were major events in natural history. The amphioxus (lancelet) is now recognised as the closest extant relative to the stem chordate and is the only living invertebrate that retains a vertebrate-like development and body plan through its lifespan, despite more than 500 million years of independent evolution from the stem vertebrate. The inspiring data coming from its recently sequenced genome confirms that amphioxus has a prototypical chordate genome with respect to gene content and structure, and even chromosomal organisation. Pushed by joint efforts of amphioxus researchers, amphioxus is now entering a new era, namely its maturation as a laboratory model, through the availability of a large amount of molecular data and the advent of experimental manipulation of the embryo. These two facts may well serve to illuminate the hidden secrets of the genetic changes that generated, among other vertebrates, ourselves.

A cDNA resource for the cephalochordate amphioxus Branchiostoma floridae

Cephalochordates are the basal invertebrate chordates within the phylum Chordata. They are widely used as a model system for research in evolutionary developmental biology (EvoDevo) to understand the basic patterning mechanisms for the chordate body plan and the origin of vertebrates. Recently, the genome of the cephalochordate Branchiostoma floridae was sequenced, which further brings this organism to the front for comparative genomic studies. In this paper, we report the generation of large-scale 5'- and 3'-expressed sequence tags (ESTs) from B. floridae and the complementary deoxyribonucleic acid (cDNA) resource for this species. Both 5'- and 3'-ESTs were sequenced for approximately 140,000 cDNA clones derived from five developmental stages, and the cDNA clones were subsequently grouped into independent clusters using 3'-EST sequences. We identified 21,229 cDNA clusters, and each corresponds to a unique transcript species from B. floridae. We then chose 24,020 cDNA clones representing all of these 21,229 clusters to generate the "Branchiostoma floridae Gene Collection Release 1." We also constructed a database with a searchable interface for this EST dataset and the related information on "Branchiostoma floridae Gene Collection Release 1." This set of cDNA clones along with our cDNA database will serve as an important resource for future research in this basal chordate. This Gene Collection and the original 140,000 individual cDNA clones are available to the research community upon request.