Coelimination and Survival in Gene Network Evolution: Dismantling the RA-Signaling in a Chordate

The ensemble of gene regulatory networks at mutation-selection balance

The evolution of diverse phenotypes both involves and is constrained by molecular interaction networks. When these networks influence patterns of expression, we refer to them as gene regulatory networks (GRNs). Here, we develop a model of GRN evolution analogous to work from quasi-species theory, which is itself essentially the mutation-selection balance model from classical population genetics extended to multiple loci. With this GRN model, we prove that-across a broad spectrum of selection pressures-the dynamics converge to a stationary distribution over GRNs. Next, we show from first principles how the frequency of GRNs at equilibrium is related to the topology of the genotype network, in particular, via a specific network centrality measure termed the eigenvector centrality. Finally, we determine the structural characteristics of GRNs that are favoured in response to a range of selective environments and mutational constraints. Our work connects GRN evolution to quasi-species theory-and thus to classical populations genetics-providing a mechanistic explanation for the observed distribution of GRNs evolving in response to various evolutionary forces, and shows how complex fitness landscapes can emerge from simple evolutionary rules.

Developmental biology of the larvacean Oikopleura dioica: Genome resources, functional screening, and imaging

The larvacean Oikopleura dioica is a cosmopolitan planktonic chordate and is closely related to vertebrates. It is characterized by a tadpole-shaped morphology with notochord flanked by muscle in the tail and brain on the dorsal side, a short life cycle of five days, a compact genome of approximately 56 Mb, a simple and transparent body with a small number of cells (~4000 in functional juveniles), invariant embryonic cell lineages, and fast development that ensures complete morphogenesis and organ formation 10 h after fertilization. With these features, this marine chordate is a promising and advantageous animal model in which genetic manipulation is feasible. In this review, we introduce relevant resources and modern techniques that have been developed: (1) Genome and transcriptomes. Oikopleura dioica has the smallest genome among non-parasitic metazoans. Its genome databases have been generated using three geographically distant O. dioica populations, and several intra-species sequence differences are becoming evident; (2) Functional genetic knockdown techniques. Comprehensive screening of genes is feasible using ovarian microinjection and double-strand DNA-induced gene knockdown; and (3) Live imaging of embryos and larvae. Application of these techniques has uncovered novel aspects of development, including meiotic cell arrest, left-right patterning, epidermal cell patterning, and mouth formation involving the connection of ectoderm and endoderm sheets. Oikopleura dioca has become very useful for developmental and evolutionary studies in chordates.

Origin and evolutionary landscape of Nr2f transcription factors across Metazoa

BACKGROUND

Nuclear Receptor Subfamily 2 Group F (Nr2f) orphan nuclear hormone transcription factors (TFs) are fundamental regulators of many developmental processes in invertebrates and vertebrates. Despite the importance of these TFs throughout metazoan development, previous work has not clearly outlined their evolutionary history.

RESULTS

We integrated molecular phylogeny with comparisons of intron/exon structure, domain architecture, and syntenic conservation to define critical evolutionary events that distinguish the Nr2f gene family in Metazoa. Our data indicate that a single ancestral eumetazoan Nr2f gene predated six main Bilateria subfamilies, which include single Nr2f homologs, here referred to as Nr2f1/2/5/6, that are present in invertebrate protostomes and deuterostomes, Nr2f1/2 homologs in agnathans, and Nr2f1, Nr2f2, Nr2f5, and Nr2f6 orthologs that are found in gnathostomes. Four cnidarian Nr2f1/2/5/6 and three agnathan Nr2f1/2 members are each due to independent expansions, while the vertebrate Nr2f1/Nr2f2 and Nr2f5/Nr2f6 members each form paralogous groups that arose from the established series of whole-genome duplications (WGDs). Nr2f6 members are the most divergent Nr2f subfamily in gnathostomes. Interestingly, in contrast to the other gnathostome Nr2f subfamilies, Nr2f5 has been independently lost in numerous vertebrate lineages. Furthermore, our analysis shows there are differential expansions and losses of Nr2f genes in teleosts following their additional rounds of WGDs.

CONCLUSION

Overall, our analysis of Nr2f gene evolution helps to reveal the origins and previously unrecognized relationships of this ancient TF family, which may allow for greater insights into the conservation of Nr2f functions that shape Metazoan body plans.

Cardiopharyngeal deconstruction and ancestral tunicate sessility

A central question in chordate evolution is the origin of sessility in adult ascidians, and whether the appendicularian complete free-living style represents a primitive or derived condition among tunicates¹. According to the 'a new heart for a new head' hypothesis, the evolution of the cardiopharyngeal gene regulatory network appears as a pivotal aspect to understand the evolution of the lifestyles of chordates^2-4. Here we show that appendicularians experienced massive ancestral losses of cardiopharyngeal genes and subfunctions, leading to the 'deconstruction' of two ancestral modules of the tunicate cardiopharyngeal gene regulatory network. In ascidians, these modules are related to early and late multipotency, which is involved in lineage cell-fate determination towards the first and second heart fields and siphon muscles. Our work shows that the deconstruction of the cardiopharyngeal gene regulatory network involved the regressive loss of the siphon muscle, supporting an evolutionary scenario in which ancestral tunicates had a sessile ascidian-like adult lifestyle. In agreement with this scenario, our findings also suggest that this deconstruction contributed to the acceleration of cardiogenesis and the redesign of the heart into an open-wide laminar structure in appendicularians as evolutionary adaptations during their transition to a complete pelagic free-living style upon the innovation of the food-filtering house⁵.

Nkx2-1 and FoxE regionalize glandular (mucus-producing) and thyroid-equivalent traits in the endostyle of the chordate Oikopleura dioica

The endostyle is a ventral pharyngeal organ used for internal filter feeding of basal chordates and is considered homologous to the follicular thyroid of vertebrates. It contains mucus-producing (glandular) and thyroid-equivalent regions organized along the dorsoventral (DV) axis. Although thyroid-related genes (Nkx2-1, FoxE, and thyroid peroxidase (TPO)) are known to be expressed in the endostyle, their roles in establishing regionalization within the organ have not been demonstrated. We report that Nkx2-1 and FoxE are essential for establishing DV axial identity in the endostyle of Oikopleura dioica. Genome and expression analyses showed von Willebrand factor-like (vWFL) and TPO/dual oxidase (Duox)/Nkx2-1/FoxE as orthologs of glandular and thyroid-related genes, respectively. Knockdown experiments showed that Nkx2-1 is necessary for the expression of glandular and thyroid-related genes, whereas FoxE is necessary only for thyroid-related genes. Moreover, Nkx2-1 expression is necessary for FoxE expression in larvae during organogenesis. The results demonstrate the essential roles of Nkx2-1 and FoxE in establishing regionalization in the endostyle, including (1) the Nkx2-1-dependent glandular region, and (2) the Nkx2-1/FoxE-dependent thyroid-equivalent region. DV axial regionalization may be responsible for organizing glandular and thyroid-equivalent traits of the pharynx along the DV axis.

Massive Gene Loss and Function Shuffling in Appendicularians Stretch the Boundaries of Chordate Wnt Family Evolution

Gene loss is a pervasive source of genetic variation that influences species evolvability, biodiversity and the innovation of evolutionary adaptations. To better understand the evolutionary patterns and impact of gene loss, here we investigate as a case study the evolution of the wingless (Wnt) family in the appendicularian tunicate Oikopleura dioica, an emergent EvoDevo model characterized by its proneness to lose genes among chordates. Genome survey and phylogenetic analyses reveal that only four of the thirteen Wnt subfamilies have survived in O. dioica—Wnt5, Wnt10, Wnt11, and Wnt16,—representing the minimal Wnt repertoire described in chordates. While the loss of Wnt4 and Wnt8 likely occurred in the last common ancestor of tunicates, representing therefore a synapomorphy of this subphylum, the rest of losses occurred during the evolution of appendicularians. This work provides the first complete Wnt developmental expression atlas in a tunicate and the first insights into the evolution of Wnt developmental functions in appendicularians. Our work highlights three main evolutionary patterns of gene loss: (1) conservation of ancestral Wnt expression domains not affected by gene losses; (2) function shuffling among Wnt paralogs accompanied by gene losses; and (3) extinction of Wnt expression in certain embryonic directly correlated with gene losses. Overall our work reveals that in contrast to “conservative” pattern of evolution of cephalochordates and vertebrates, O. dioica shows an even more radical “liberal” evolutionary pattern than that described ascidian tunicates, stretching the boundaries of the malleability of Wnt family evolution in chordates.

Evolutionary aspects and enzymology of metazoan carotenoid cleavage oxygenases

The carotenoids are terpenoid fat-soluble pigments produced by plants, algae, and several bacteria and fungi. They are ubiquitous components of animal diets. Carotenoid cleavage oxygenase (CCO) superfamily members are involved in carotenoid metabolism and are present in all kingdoms of life. Throughout the animal kingdom, carotenoid oxygenases are widely distributed and they are completely absent only in two unicellular organisms, Monosiga and Leishmania. Mammals have three paralogs 15,15'-β-carotene oxygenase (BCO1), 9',10'-β-carotene oxygenase (BCO2) and RPE65. The first two enzymes are classical carotenoid oxygenases: they cleave carbon‑carbon double bonds and incorporate two atoms of oxygen in the substrate at the site of cleavage. The third, RPE65, is an unusual family member, it is the retinoid isomerohydrolase in the visual cycle that converts all-trans-retinyl ester into 11-cis-retinol. Here we discuss evolutionary aspects of the carotenoid cleavage oxygenase superfamily and their enzymology to deduce what insight we can obtain from their evolutionary conservation.

The Cis-Regulatory Code for Kelch-like 21/30 Specific Expression in Ciona robusta Sensory Organs

The tunicate Ciona robusta is an emerging model system to study the evolution of the nervous system. Due to their small embryos and compact genomes, tunicates, like Ciona robusta, have great potential to comprehend genetic circuitry underlying cell specific gene repertoire, among different neuronal cells. Their simple larvae possess a sensory vesicle comprising two pigmented sensory organs, the ocellus and the otolith. We focused here on Klhl21/30, a gene belonging to Kelch family, that, in Ciona robusta, starts to be expressed in pigmented cell precursors, becoming specifically maintained in the otolith precursor during embryogenesis. Evolutionary analyses demonstrated the conservation of Klhl21/30 in all the chordates. Cis-regulatory analyses and CRISPR/Cas9 mutagenesis of potential upstream factors, revealed that Klhl21/30 expression is controlled by the combined action of three transcription factors, Mitf, Dmrt, and Msx, which are downstream of FGF signaling. The central role of Mitf is consistent with its function as a fundamental regulator of vertebrate pigment cell development. Moreover, our results unraveled a new function for Dmrt and Msx as transcriptional co-activators in the context of the Ciona otolith.

Rimbp, a New Marker for the Nervous System of the Tunicate Ciona robusta

Establishment of presynaptic mechanisms by proteins that regulate neurotransmitter release in the presynaptic active zone is considered a fundamental step in animal evolution. Rab3 interacting molecule-binding proteins (Rimbps) are crucial components of the presynaptic active zone and key players in calcium homeostasis. Although Rimbp involvement in these dynamics has been described in distantly related models such as fly and human, the role of this family in most invertebrates remains obscure. To fill this gap, we defined the evolutionary history of Rimbp family in animals, from sponges to mammals. We report, for the first time, the expression of the two isoforms of the unique Rimbp family member in Ciona robusta in distinct domains of the larval nervous system. We identify intronic enhancers that are able to drive expression in different nervous system territories partially corresponding to Rimbp endogenous expression. The analysis of gene expression patterns and the identification of regulatory elements of Rimbp will positively impact our understanding of this family of genes in the context of Ciona embryogenesis.

A genome database for a Japanese population of the larvacean Oikopleura dioica

The larvacean Oikopleura dioica is a planktonic chordate and is a tunicate that belongs to the closest relatives to vertebrates. Its simple and transparent body, invariant embryonic cell lineages, and short life cycle of 5 days make it a promising model organism for the study of developmental biology. The genome browser OikoBase was established in 2013 using Norwegian O. dioica. However, genome information for other populations is not available, even though many researchers have studied local populations. In the present study, we sequenced using Illumina and PacBio RSII technologies the genome of O. dioica from a southwestern Japanese population that was cultured in our laboratory for 3 years. The genome of Japanese O. dioica was assembled into 576 scaffold sequences with a total length and N50 length of 56.6 and 1.5 Mb, respectively. A total of 18,743 gene models (transcript models) were predicted in the genome assembly, named OSKA2016. In addition, 19,277 non-redundant transcripts were assembled using RNA-seq data. The OSKA2016 has global sequence similarity of only 86.5% when compared with the OikoBase, highlighting the sequence difference between the two far distant O. dioica populations on the globe. The genome assembly, transcript assembly, and transcript models were incorporated into ANISEED (https://www.aniseed.cnrs.fr/) for genome browsing and BLAST searches. Mapping of reads obtained from male- or female-specific genome libraries yielded male-specific scaffolds in the OSKA2016 and revealed that over 2.6 Mb of sequence were included in the male-specific Y-region. The genome and transcriptome resources from two distinct populations will be useful datasets for developmental biology, evolutionary biology, and molecular ecology using this model organism.

A chordate species lacking Nodal utilizes calcium oscillation and Bmp for left-right patterning

Larvaceans are chordates with a tadpole-like morphology. In contrast to most chordates of which early embryonic morphology is bilaterally symmetric and the left-right (L-R) axis is specified by the Nodal pathway later on, invariant L-R asymmetry emerges in four-cell embryos of larvaceans. The asymmetric cell arrangements exist through development of the tailbud. The tail thus twists 90° in a counterclockwise direction relative to the trunk, and the tail nerve cord localizes on the left side. Here, we demonstrate that larvacean embryos have nonconventional L-R asymmetries: 1) L- and R-cells of the two-cell embryo had remarkably asymmetric cell fates; 2) Ca²⁺ oscillation occurred through embryogenesis; 3) Nodal, an evolutionarily conserved left-determining gene, was absent in the genome; and 4) bone morphogenetic protein gene (Bmp) homolog Bmp.a showed right-sided expression in the tailbud and larvae. We also showed that Ca²⁺ oscillation is required for Bmp.a expression, and that BMP signaling suppresses ectopic expression of neural genes. These results indicate that there is a chordate species lacking Nodal that utilizes Ca²⁺ oscillation and Bmp.a for embryonic L-R patterning. The right-side Bmp.a expression may have arisen via cooption of conventional BMP signaling in order to restrict neural gene expression on the left side.

The evolutionary landscape of the Rab family in chordates

Intracellular traffic amongst organelles represents a key feature for eukaryotes and is orchestrated principally by members of Rab family, the largest within Ras superfamily. Given that variations in Rab repertoire have been fundamental in animal diversification, we provided the most exhaustive survey regarding the Rab toolkit of chordates. Our findings reveal the existence of 42 metazoan conserved subfamilies exhibiting a univocal intron/exon structure preserved from cnidarians to vertebrates. Since the current view does not capture the Rab complexity, we propose a new Rab family classification in three distinct monophyletic clades. The Rab complement of chordates shows a dramatic diversification due to genome duplications and independent gene duplications and losses with sharp differences amongst cephalochordates, tunicates and gnathostome vertebrates. Strikingly, the analysis of the domain architecture of this family highlighted the existence of chimeric calcium-binding Rabs, which are animal novelties characterized by a complex evolutionary history in gnathostomes and whose role in cellular metabolism is obscure. This work provides novel insights in the knowledge of Rab family: our hypothesis is that chordates represent a hotspot of Rab variability, with many events of gene gains and losses impacting intracellular traffic capabilities. Our results help to elucidate the role of Rab members in the transport amongst endomembranes and shed light on intracellular traffic routes in vertebrates. Then, since the predominant role of Rabs in the molecular communication between different cellular districts, this study paves to way to comprehend inherited or acquired human disorders provoked by dysfunctions in Rab genes.

Genes lost during the transition from land to water in cetaceans highlight genomic changes associated with aquatic adaptations

The transition from land to water in whales and dolphins (cetaceans) was accompanied by remarkable adaptations. To reveal genomic changes that occurred during this transition, we screened for protein-coding genes that were inactivated in the ancestral cetacean lineage. We found 85 gene losses. Some of these were likely beneficial for cetaceans, for example, by reducing the risk of thrombus formation during diving (F12 and KLKB1), erroneous DNA damage repair (POLM), and oxidative stress-induced lung inflammation (MAP3K19). Additional gene losses may reflect other diving-related adaptations, such as enhanced vasoconstriction during the diving response (mediated by SLC6A18) and altered pulmonary surfactant composition (SEC14L3), while loss of SLC4A9 relates to a reduced need for saliva. Last, loss of melatonin synthesis and receptor genes (AANAT, ASMT, and MTNR1A/B) may have been a precondition for adopting unihemispheric sleep. Our findings suggest that some genes lost in ancestral cetaceans were likely involved in adapting to a fully aquatic lifestyle.

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of “gene vestiges” in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

ORTHOSCOPE: An Automatic Web Tool for Phylogenetically Inferring Bilaterian Orthogroups with User-Selected Taxa

Identification of orthologous or paralogous relationships of coding genes is fundamental to all aspects of comparative genomics. For accurate identification of orthologs among deeply diversified bilaterian lineages, precise estimation of gene trees is indispensable, given the complicated histories of genes over millions of years. By estimating gene trees, orthologs can be identified as members of an orthogroup, a set of genes descended from a single gene in the last common ancestor of all the species being considered. In addition to comparisons with a given species tree, purposeful taxonomic sampling increases the accuracy of gene tree estimation and orthogroup identification. Although some major phylogenetic relationships of bilaterians are gradually being unraveled, the scattering of published genomic data among separate web databases is becoming a significant hindrance to identification of orthogroups with appropriate taxonomic sampling. By integrating more than 250 metazoan gene models predicted in genome projects, we developed a web tool called ORTHOSCOPE to identify orthogroups of specific protein-coding genes within major bilaterian lineages. ORTHOSCOPE allows users to employ several sequences of a specific molecule and broadly accepted nodes included in a user-specified species tree as queries and to evaluate the reliability of estimated orthogroups based on topologies and node support values of estimated gene trees. A test analysis using data from 36 bilaterians was accomplished within 140 s. ORTHOSCOPE results can be used to evaluate orthologs identified by other stand-alone programs using genome-scale data. ORTHOSCOPE is freely available at https://www.orthoscope.jp or https://github.com/jun-inoue/orthoscope (last accessed December 28, 2018).

The Singularity of Cetacea Behavior Parallels the Complete Inactivation of Melatonin Gene Modules

Melatonin, the hormone of darkness, is a peculiar molecule found in most living organisms. Emerging as a potent broad-spectrum antioxidant, melatonin was repurposed into extra roles such as the modulation of circadian and seasonal rhythmicity, affecting numerous aspects of physiology and behaviour, including sleep entrainment and locomotor activity. Interestingly, the pineal gland—the melatonin synthesising organ in vertebrates—was suggested to be absent or rudimentary in some mammalian lineages, including Cetacea. In Cetacea, pineal regression is paralleled by their unique bio-rhythmicity, as illustrated by the unihemispheric sleeping behaviour and long-term vigilance. Here, we examined the genes responsible for melatonin synthesis (Aanat and Asmt) and signalling (Mtnr1a and Mtnr1b) in 12 toothed and baleen whale genomes. Based on an ample genomic comparison, we deduce that melatonin-related gene modules are eroded in Cetacea.

Oikopleura dioica: An Emergent Chordate Model to Study the Impact of Gene Loss on the Evolution of the Mechanisms of Development

The urochordate Oikopleura dioica is emerging as a nonclassical animal model in the field of evolutionary developmental biology (a.k.a. evo-devo) especially attractive for investigating the impact of gene loss on the evolution of mechanisms of development. This is because this organism fulfills the requirements of an animal model (i.e., has a simple and accessible morphology, a short generation time and life span, and affordable culture in the laboratory and amenable experimental manipulation), but also because O. dioica occupies a key phylogenetic position to understand the diversification and origin of our own phylum, the chordates. During its evolution, O. dioica genome has suffered a drastic process of compaction, becoming the smallest known chordate genome, a process that has been accompanied by exacerbating amount of gene losses. Interestingly, however, despite the extensive gene losses, including entire regulatory pathways essential for the embryonic development of other chordates, O. dioica retains the typical chordate body plan. This unexpected situation led to the formulation of the so-called inverse paradox of evo-devo, that is, when a genetic diversity is able to maintain a phenotypic unity. This chapter reviews the biological features of O. dioica as a model animal, along with the current data on the evolution of its genes and genome. We pay special attention to the numerous examples of gene losses that have taken place during the evolution of this unique animal model, which is helping us to understand to which the limits of evo-devo can be pushed off.

Diatom bloom-derived biotoxins cause aberrant development and gene expression in the appendicularian chordate Oikopleura dioica

Investigating environmental hazards than could affect appendicularians is of prime ecological interest because they are among the most abundant components of the mesozooplankton. This work shows that embryo development of the appendicularian Oikopleura dioica is compromised by diatom bloom-derived biotoxins, even at concentrations in the same range as those measured after blooms. Developmental gene expression analysis of biotoxin-treated embryos uncovers an aberrant golf ball-like phenotype affecting morphogenesis, midline convergence, and tail elongation. Biotoxins induce a rapid upregulation of defensome genes, and considerable delay and silencing of zygotic transcription of developmental genes. Upon a possible future intensification of blooms associated with ocean warming and acidification, our work puts an alert on the potential impact that an increase of biotoxins may have on marine food webs, and points to defensome genes as molecular biosensors that marine ecologists could use to monitor the genetic stress of natural populations exposed to microalgal blooms.

De novo draft assembly of the Botrylloides leachii genome provides further insight into tunicate evolution

Tunicates are marine invertebrates that compose the closest phylogenetic group to the vertebrates. These chordates present a particularly diverse range of regenerative abilities and life-history strategies. Consequently, tunicates provide an extraordinary perspective into the emergence and diversity of these traits. Here we describe the genome sequencing, annotation and analysis of the Stolidobranchian Botrylloides leachii. We have produced a high-quality 159 Mb assembly, 82% of the predicted 194 Mb genome. Analysing genome size, gene number, repetitive elements, orthologs clustering and gene ontology terms show that B. leachii has a genomic architecture similar to that of most solitary tunicates, while other recently sequenced colonial ascidians have undergone genome expansion. In addition, ortholog clustering has identified groups of candidate genes for the study of colonialism and whole-body regeneration. By analysing the structure and composition of conserved gene linkages, we observed examples of cluster breaks and gene dispersions, suggesting that several lineage-specific genome rearrangements occurred during tunicate evolution. We also found lineage-specific gene gain and loss within conserved cell-signalling pathways. Such examples of genetic changes within conserved cell-signalling pathways commonly associated with regeneration and development that may underlie some of the diverse regenerative abilities observed in tunicates. Overall, these results provide a novel resource for the study of tunicates and of colonial ascidians.

Reporter Analyses Reveal Redundant Enhancers that Confer Robustness on Cis-Regulatory Mechanisms

Reporter analyses of Hox1 and Brachyury (Bra) genes have revealed examples of redundant enhancers that provide regulatory robustness. Retinoic acid (RA) activates through an RA-response element the transcription of Hox1 in the nerve cord of the ascidian Ciona intestinalis. We also found a weak RA-independent neural enhancer within the second intron of Hox1. The Hox1 gene in the larvacean Oikopleura dioica is also expressed in the nerve cord. The O. dioica genome, however, does not contain the RA receptor-encoding gene, and the expression of Hox1 has become independent of RA. We have found that the upstream sequence of the O. dioica Hox1 was able to activate reporter gene expression in the nerve cord of the C. intestinalis embryo, suggesting that an RA-independent regulatory system in the nerve cord might be common in larvaceans and ascidians. This RA-independent redundant regulatory system may have facilitated the Oikopleura ancestor losing RA signaling without an apparent impact on Hox1 expression domains. On the other hand, vertebrate Bra is expressed in the ventral mesoderm and notochord, whereas its ascidian ortholog is exclusively expressed in the notochord. Fibroblast growth factor (FGF) induces Bra in the ventral mesoderm in vertebrates, whereas it induces Bra in the notochord in ascidians. Disruption of the FGF signal does not completely silence Bra expression in ascidians, suggesting that FGF-dependent and independent enhancers might comprise a redundant regulatory system in ascidians. The existence of redundant enhancers, therefore, provides regulatory robustness that may facilitate the acquisition of new expression domains.

Metallothioneins of the urochordate Oikopleura dioica have Cys-rich tandem repeats, large size and cadmium-binding preference

Oikopleura dioica has the longest metallothionein described so far, made of repeats generated by a modular and step-wise evolution.