Laboratory study of Fritillaria lifecycle reveals key morphogenetic events leading to genus-specific anatomy.

A fascinating variety of adult body plans can be found in the Tunicates, the closest existing relatives of vertebrates. A distinctive feature of the larvacean class of pelagic tunicates is the presence of a highly specialized surface epithelium that produces a cellulose test, the "larvacean house". While substantial differences exist between the anatomy of larvacean families, most of the ontogeny is derived from the observations of a single genus, Oikopleura. We present the first study of Fritillaria development based on the observation of individuals reproduced in the laboratory. Like the other small epipelagic species Oikopleura dioica, the larvae of Fritillaria borealis grow rapidly in the laboratory, and they acquire the adult form within a day. We could show that major morphological differences exhibited by Fritillaria and Oikopleura adults originate from a key developmental stage during larval organogenesis. Here, the surface epithelium progressively retracts from the posterior digestive organs of Fritillaria larvae, and it establishes house-producing territories around the pharynx. Our results show that the divergence between larvacean genera was associated with a profound rearrangement of the mechanisms controlling the differentiation of the larval ectoderm.

Conservative route to genome compaction in a miniature annelid.

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

Evolution of the U2 Spliceosome for Processing Numerous and Highly Diverse Non-canonical Introns in the Chordate Fritillaria borealis.

An overwhelming majority of eukaryotic introns have GT/AG ends, whose identities play a critical role for their recognition and removal by the U2 spliceosome, a well-conserved complex of protein and RNAs. Introns with other splice sites exist at very low frequencies in various genomes, and some of them are processed by the U12 spliceosome. Here, we show that, in the chordate Fritillaria borealis, the majority of old introns have been lost and replaced by introns with highly divergent splice sites. The new introns of F. borealis are exceptionally diverse, though more frequently AG/AC or AG/AT, and features of thousands of them support an origin from transposons. They cannot be processed in human cells, but their splicing is rescued by mutating terminal dinucleotides to GT/AG. With lariat sequencing and splicing inhibitor assays, we show that F. borealis introns are spliced by the U2 spliceosome, which thus evolved to a different selectivity, with neither novel U1 small nuclear RNA (snRNA) types nor major remodeling of its protein and snRNA complements. This genome-wide recolonization by non-canonical introns emphasizes the importance of transposons as a resource of novel introns in a context of massive intron loss. An evolution of the spliceosome may also permit to neutralize harmful transposons through their conversion into introns.

Massive Changes of Genome Size Driven by Expansions of Non-autonomous Transposable Elements.

In eukaryotes, genome size correlates little with the number of coding genes or the level of organismal complexity (C-value paradox). The underlying causes of variations in genome size, whether adaptive or neutral, remain unclear, although several biological traits often covary with it [1-5]. Rapid increases in genome size occur mainly through whole-genome duplications or bursts in the activity of transposable elements (TEs) [6]. The very small and compact genome of Oikopleura dioica, a tunicate of the larvacean class, lacks elements of most ancient families of animal retrotransposons [7, 8]. Here, we sequenced the genomes of six other larvaceans, all of which are larger than that of Oikopleura (up to 12 times) and which increase in size with greater body length. Although no evidence was found for whole-genome duplications within the group of species, the global amount of TEs strongly correlated with genome size. Compared to other metazoans, however, the TE diversity was reduced in all species, as observed previously in O. dioica, suggesting a common ancestor with a compacted genome. Strikingly, non-autonomous elements, particularly short interspersed nuclear elements (SINEs), massively contributed to genome size variation through species-specific independent amplifications, ranging from 3% in the smallest genome up to 49% in the largest. Variations in SINE abundance explain as much as 83% of interspecific genome size variation. These data support an indirect influence of autonomous TEs on genome size via their ability to mobilize non-autonomous elements.

Prevalence of Mutation-Prone Microhomology-Mediated End Joining in a Chordate Lacking the c-NHEJ DNA Repair Pathway.

Classical non-homologous end joining (c-NHEJ), a fundamental pathway that repairs double-strand breaks in DNA, is almost universal in eukaryotes and involves multiple proteins highly conserved from yeast to human [1]. The genes encoding these proteins were not detected in the genome of Oikopleura dioica, a new model system of tunicate larvaceans known for its very compact and highly rearranged genome [2-4]. After showing their absence in the genomes of six other larvacean species, the present study examined how O. dioica oocytes and embryos repair double-strand DNA breaks (DSBs), using two approaches: the injection of linearized plasmids, which resulted in their rapid end joining, and a newly established CRISPR Cas9 technique. In both cases, end joining merged short microhomologous sequences surrounding the break (mainly 4 bp long), thus inducing deletions larger than for the tunicate ascidian Ciona intestinalis and human cells. A relatively high frequency of nucleotide insertions was also observed. Finally, a survey of genomic indels supports the involvement of microhomology-mediated repair in natural conditions. Overall, O. dioica repairs DSBs as other organisms do when their c-NHEJ pathway is experimentally rendered deficient, using another mode of end joining with the same effect as alternative NHEJ (a-NHEJ) or microhomology-mediated end joining (MMEJ) [5-7]. We discuss how the exceptional loss of c-NHEJ and its replacement by a more mutation-prone mechanism may have contributed to reshaping this genome and even been advantageous under pressure for genome compaction.

Development of the house secreting epithelium, a major innovation of tunicate larvaceans, involves multiple homeodomain transcription factors.

The mechanisms driving innovations that distinguish large taxons are poorly known and essentially accessible via a candidate gene approach. A spectacular acquisition by tunicate larvaceans is the house, a complex extracellular filtration device. Its components are secreted by the oikoplastic epithelium which covers the animal trunk. Here we describe the development of this epithelium in larvae through the formation of specific cellular territories known to produce distinct sets of house proteins (Oikosins). It involves cell divisions and morphological differentiation but very limited cell migration. A diverse set of homeobox genes, most often duplicated in the genome, are transiently and site-specifically expressed in the trunk epithelium at early larval stages. Using RNA interference, we show that two prop duplicates are involved in the differentiation of a region on and around the dorsal midline, regulating morphology and the production of a specific oikosin. Our observations favor a scenario in which multiple homeobox genes and most likely other developmental transcription factors were recruited for this innovation. Their frequent duplications probably predated, but were not required for the emergence of the house.

Evidence for a centrosome-attracting body like structure in germ-soma segregation during early development, in the urochordate Oikopleura dioica.

BACKGROUND: Germ cell formation has been investigated in sessile forms of tunicates. This process involves the release of a subset of maternal transcripts from the centrosome-attracting body (CAB) in the progenitor cells of the germ line. When germ-soma segregation is completed, CAB structures are missing from the newly formed primordial germ cells (PGCs). In free-swimming tunicates, knowledge about germ cell formation is lacking. In this investigation, comparative gene expression and electron microscopy studies were used to address germ cell formation in Oikopleura dioica (O. dioica). RESULTS: We found that the RNA localization pattern of pumilio (pum1) is similar to the pattern described for a subset of maternal transcripts marking the posterior end of ascidian embryos. Transcripts marking the posterior end are called postplasmic or posterior-end mark (PEM) transcripts. We found no localization of vasa (vas) transcripts to any sub-region within the germ-line precursor cells. Expression of vas4 was detected in the newly formed PGCs. Electron microscopy studies confirmed the presence of structures with similar morphology to CAB. In the same cytoplasmic compartment, we also identified pum1 transcripts and an epitope recognized by an antibody to histone H3 phosphorylated on serine 28. CONCLUSIONS: Our findings support that a CAB-like structure participates in the segregation of maternal pum1 transcripts during germ-soma separation in O. dioica.

Embryonic expression of endogenous retroviral RNAs in somatic tissues adjacent to the Oikopleura germline.

Selective pressure to maintain small genome size implies control of transposable elements, and most old classes of retrotransposons are indeed absent from the very compact genome of the tunicate Oikopleura dioica. Nonetheless, two families of retrotransposons are present, including the Tor elements. The gene organization within Tor elements is similar to that of LTR retrotransposons and retroviruses. In addition to gag and pol, many Tor elements carry a third gene encoding viral envelope-like proteins (Env) that may mediate infection. We show that the Tor family contains distinct classes of elements. In some classes, env mRNA is transcribed from the 5'LTR as in retroviruses. In others, env is transcribed from an additional promoter located downstream of the 5'LTR. Tor Env proteins are membrane-associated glycoproteins which exhibit some features of viral membrane fusion proteins. Whereas some elements are expressed in the adult testis, many others are specifically expressed in embryonic somatic cells adjacent to primordial germ cells. Such embryonic expression depends on determinants present in the Tor elements and not on their surrounding genomic environment. Our study shows that unusual modes of transcription and expression close to the germline may contribute to the proliferation of Tor elements.

Modification of the larval swimming behavior in Oikopleura dioica, a chordate with a miniaturized central nervous system by dsRNA injection into fertilized eggs.

Using RNA interference, we have selectively perturbed neurotransmitter-related features of the larval swimming behavior of Oikopleura dioica, a tunicate with a central nervous system comprising about 130 neurons. We injected dsRNA into fertilized eggs to knockdown the expression of the genes, respectively, encoding ChAT (choline acetyltransferase) and GAD (glutamic acid decarboxylase), enzymes critical for the biosynthesis of acetylcholine and GABA. These two neurotransmitters have conserved roles during evolution, particularly within chordate motor systems, where they mediate respectively neuromuscular and central inhibitory signals. In Oikopleura, interference with ChAT expression prevented the normal bidirectional, propagating tail movement characteristic of swimming, permitting only repeated unilateral tail bends. Proper swimming was never observed, and the resting period between episodes of activity was lengthened. This phenotype is most likely caused by the reduction of transcription observed for both the targeted ChAT gene and the VAChT gene (Vesicular Acetylcholine Transporter), both genes being transcribed from the same operon. Interference with GAD expression led to an uncoordinated version of swimming with a spiral movement trajectory, but with episodes similar in duration and cycle frequency to those of normal swimming. Our results suggest locomotor functions for ChAT and GABA that are more subtle than previously proposed for tunicates and opens the way for a genetic dissection of Oikopleura neuronal circuits, which are likely to be among the most simplified in the chordate phylum.

Oikopleura dioica culturing made easy: a low-cost facility for an emerging animal model in EvoDevo.

The genome sequencing and the development of RNAi knockdown technologies in the urochordate Oikopleura dioica are making this organism an attractive emergent model in the field of EvoDevo. To succeed as a new animal model, however, an organism needs to be easily and affordably cultured in the laboratory. Nowadays, there are only two facilities in the world capable to indefinitely maintain Oikopleura dioica, one in the SARS institute (Bergen, Norway) and the other in the Osaka University (Japan). Here, we describe the setup of a new facility in the University of Barcelona (Spain) in which we have modified previously published husbandry protocols to optimize the weekly production of thousands of embryos and hundreds of mature animals using the minimum amount of space, human resources, and technical equipment. This optimization includes novel protocols of cryopreservation and solid cultures for long-term maintenance of microalgal stocks-Chaetoceros calcitrans, Isochrysis sp., Rhinomonas reticulate, and Synechococcus sp.-needed for Oikopleura dioica feeding. Our culture system maintains partially inbred lines healthy with similar characteristics to wild animals, and it is easily expandable to satisfy on demand the needs of any laboratory that may wish to use Oikopleura dioica as a model organism.

OikoBase: a genomics and developmental transcriptomics resource for the urochordate Oikopleura dioica.

We report the development of OikoBase (http://oikoarrays.biology.uiowa.edu/Oiko/), a tiling array-based genome browser resource for Oikopleura dioica, a metazoan belonging to the urochordates, the closest extant group to vertebrates. OikoBase facilitates retrieval and mining of a variety of useful genomics information. First, it includes a genome browser which interrogates 1260 genomic sequence scaffolds and features gene, transcript and CDS annotation tracks. Second, we annotated gene models with gene ontology (GO) terms and InterPro domains which are directly accessible in the browser with links to their entries in the GO (http://www.geneontology.org/) and InterPro (http://www.ebi.ac.uk/interpro/) databases, and we provide transcript and peptide links for sequence downloads. Third, we introduce the transcriptomics of a comprehensive set of developmental stages of O. dioica at high resolution and provide downloadable gene expression data for all developmental stages. Fourth, we incorporate a BLAST tool to identify homologs of genes and proteins. Finally, we include a tutorial that describes how to use OikoBase as well as a link to detailed methods, explaining the data generation and analysis pipeline. OikoBase will provide a valuable resource for research in chordate development, genome evolution and plasticity and the molecular ecology of this important marine planktonic organism.

Conservation and divergence of chemical defense system in the tunicate Oikopleura dioica revealed by genome wide response to two xenobiotics.

BACKGROUND: Animals have developed extensive mechanisms of response to xenobiotic chemical attacks. Although recent genome surveys have suggested a broad conservation of the chemical defensome across metazoans, global gene expression responses to xenobiotics have not been well investigated in most invertebrates. Here, we performed genome survey for key defensome genes in Oikopleura dioica genome, and explored genome-wide gene expression using high density tiling arrays with over 2 million probes, in response to two model xenobiotic chemicals - the carcinogenic polycyclic aromatic hydrocarbon benzo[a]pyrene (BaP) the pharmaceutical compound Clofibrate (Clo). RESULTS: Oikopleura genome surveys for key genes of the chemical defensome suggested a reduced repertoire. Not more than 23 cytochrome P450 (CYP) genes could be identified, and neither CYP1 family genes nor their transcriptional activator AhR was detected. These two genes were present in deuterostome ancestors. As in vertebrates, the genotoxic compound BaP induced xenobiotic biotransformation and oxidative stress responsive genes. Notable exceptions were genes of the aryl hydrocarbon receptor (AhR) signaling pathway. Clo also affected the expression of many biotransformation genes and markedly repressed genes involved in energy metabolism and muscle contraction pathways. CONCLUSIONS: Oikopleura has the smallest number of CYP genes among sequenced animal genomes and lacks the AhR signaling pathway. However it appears to have basic xenobiotic inducible biotransformation genes such as a conserved genotoxic stress response gene set. Our genome survey and expression study does not support a role of AhR signaling pathway in the chemical defense of metazoans prior to the emergence of vertebrates.

Involvement of Prop1 homeobox gene in the early development of fish pituitary gland.

When mutated in mammals, paired-like homeobox Prop1 gene produces highly variable pituitary phenotypes with impaired regulation of Pit1 and eventually defective synthesis of Pit1-regulated pituitary hormones. Here we have identified fish prop1 orthologs, confirmed their pituitary-specific expression, and blocked the splicing of zebrafish prop1 transcripts using morpholino oligonucleotides. Very early steps of the gland formation seemed unaffected based on morphology and expression of early placodal marker pitx. Prop1 knock-down reduced the expression of pit1, prl (prolactin) and gh (growth hormone), as expected if the function of Prop1 is conserved throughout vertebrates. Less expectedly, lim3 was down regulated. This gene is expressed from early stages of vertebrate pituitary development but is not known to be Prop1-dependent. In situ hybridizations on prop1 morphants using probes for the pan pituitary gene pitx3 and for the hormone gene markers prl, gh and tshß, revealed abnormal shape, growth and cellular organization of the developed adenohypophysis. Strikingly, the effects of prop1 knock-down on adenohypophysis morphology and gene expression were gradually reversed during late development, despite persistent splice-blocking of transcripts. Therefore, prop1 function appears to be conserved between mammals and fish, at least for the mediation of hormonal cell type differentiation via pit1, but the existence of other fish-specific pathways downstream of prop1 are suggested by our observations.

Culture optimization for the emergent zooplanktonic model organism Oikopleura dioica.

The pan-global marine appendicularian, Oikopleura dioica, shows considerable promise as a candidate model organism for cross-disciplinary research ranging from chordate genetics and evolution to molecular ecology research. This urochordate, has a simplified anatomical organization, remains transparent throughout an exceptionally short life cycle of less than 1 week and exhibits high fecundity. At 70 Mb, the compact, sequenced genome ranks among the smallest known metazoan genomes, with both gene regulatory and intronic regions highly reduced in size. The organism occupies an important trophic role in marine ecosystems and is a significant contributor to global vertical carbon flux. Among the short list of bona fide biological model organisms, all share the property that they are amenable to long-term maintenance in laboratory cultures. Here, we tested diet regimes, spawn densities and dilutions and seawater treatment, leading to optimization of a detailed culture protocol that permits sustainable long-term maintenance of O. dioica, allowing continuous, uninterrupted production of source material for experimentation. The culture protocol can be quickly adapted in both coastal and inland laboratories and should promote rapid development of the many original research perspectives the animal offers.

Spatio-temporal expression patterns of anterior Hox genes in Atlantic salmon (Salmo salar).

Hox genes encode transcription factors that play important roles in patterning the anterior-posterior (A-P) body axis. In vertebrates, up to 14 Hox genes are physically linked in 4-13 chromosomal clusters. Their expression patterns obey spatial and temporal collinearity. Genes located at the 3' end of the clusters are expressed earlier and more anteriorly than those at the 5' end. To investigate how the expression of Hox genes has evolved after very recent ( approximately 25-100 Mya) and relatively recent ( approximately 320-350 Mya) genome duplications, we focused on three paralogous groups of salmon anterior Hox genes, in which gene duplicates have been retained. RNA-RNA whole mount in situ hybridization with gene-specific probes at early development stages showed essentially conserved expression patterns for most genes when compared to mouse and zebrafish orthologs. However, changes in spatial expression were observed for the ancient fish gene duplicate HoxB3b, while recently duplicated genes showed divergence in their expression levels.

Independent and dynamic reallocation of pitx gene expression during vertebrate evolution, with emphasis on fish pituitary development.

Several transcription regulators play key roles during pituitary morphogenesis. Well known intrinsic signals of the adenohypophysis such as the K(50)paired-like homeodomain proteins regulate commitment, proliferation, differential specification and maintenance of adenohypophyseal cells. We have cloned and successively characterized the mRNA localization of three pitx gene-pairs and three of their splice variants in salmon, pitx1alpha, pitx1beta; pitx2alpha, pitx2beta; pitx3alpha, pitx3beta; pitx1alphash, pitx1betash and pitx2alphaA. The high level of conservation between the pitx paralog-pairs indicates that they likely arose from lineage-specific genome duplication. We also report the isolation of a pitx1 gene in zebrafish. Comparative ISH studies of zebrafish, salmon and mouse pitx genes indicate both conservation and divergence of spatial expression domains in vertebrates. Significant differences were observed between the expressions of pitx orthologs during pituitary development. We suggest that the ancestral pituitary expression at early and late events of morphogenesis is preserved in different species through complementary shuffling of expression between the distinct pitx members of the family. Moreover, ISH analysis of the pitx salmon repertoire shows rapid evolution in this lineage, differences in spatio-temporal expression are observed between gene duplicates.

Embryology of a planktonic tunicate reveals traces of sessility.

A key problem in understanding deuterostome evolution has been the origin of the chordate body plan. A biphasic life cycle with a sessile adult and a free-swimming larva is traditionally considered ancestral in chordates with subsequent neotenic loss of the sessile adult stage. Molecular phylogenies challenged this view, suggesting that the primitive life cycle in chordates was entirely free-living as in modern day larvaceans. Here, we report the precise cell lineage and fate map in the normal embryo of the larvacean Oikopleura dioica, using 4D microscopy technique and transmission electron microscopy. We document the extraordinary rapidity of cleavage and morphogenetic events until hatching and demonstrate that--compared with ascidians--fate restriction occurs considerably earlier in O. dioica and that clonal organization of the cell lineage is more tightly coupled to tissue fate. We show that epidermal cells in the trunk migrate through 90 degrees, reminiscent of events in ascidian metamorphosis and that the axis of bilateral symmetry in the tail rotates in relation to the trunk. We argue that part of the tail muscle cells are ectomesodermal, because they are more closely associated with prospective epidermis than with other tissues in the cell lineage. Cladistic comparison with other deuterostomes suggests that these traits are derived within tunicates strengthening the hypothesis that the last common ancestor of tunicates had a sessile adult and thus support traditional morphology-derived scenarios. Our results allow hypothesizing that molecular developmental mechanisms known from ascidian models are restricted to fewer, yet identifiable, cells in O. dioica.

Differential evolution of the 13 Atlantic salmon Hox clusters.

Hox cluster organization represents a valuable marker to study the effects of recent genome duplication in salmonid fish (25-100 Mya). Using polymerase chain reaction amplification of cDNAs, BAC library screening, and genome walking, we reconstructed 13 Hox clusters in the Atlantic salmon containing 118 Hox genes including 8 pseudogenes. Hox paralogs resulting from the genome duplication preceding the radiation of ray-finned fish have been much better preserved in salmon than in other model teleosts. The last genome duplication in the salmon lineage has been followed by the loss of 1 of the 4 HoxA clusters. Four rounds of genome duplication after the vertebrate ancestor salmon Hox clusters display the main organizational features of vertebrate Hox clusters, with Hox genes exclusively that are densely packed in the same orientation. Recently, duplicated Hox clusters have engaged a process of divergence, with several cases of pseudogenization or asymmetrical evolution of Hox gene duplicates, and a marked erosion of identity in noncoding sequences. Strikingly, the level of divergence attained strongly depends on the Hox cluster pairs rather than on the Hox genes within each cluster. It is particularly high between both HoxBb clusters and both HoxDa clusters, whereas both HoxBa clusters remained virtually identical. Positive selection on the Hox protein-coding sequences could not be detected.

Development of the caudal nerve cord, motoneurons, and muscle innervation in the appendicularian urochordate Oikopleura dioica.

The development of the caudal nerve cord and muscle innervation in the appendicularian Oikopleura dioica was assessed using differential interference contrast and confocal microscopy, phalloidin staining of actin, and in situ hybridization for the neuronal markers tubulin and choline acetyltransferase (ChAT). The caudal nerve cord first appears as a stream of tubulin mRNA-positive neurons that extends into the tail from the caudal ganglion. By this stage a few actin-rich nerve fibers course longitudinally along the cord. As the tail lengthens, the caudal nerve cord extends and becomes more fasciculated and the neurons cluster at stereotyped longitudinal positions. The number of neurons in the nerve cord reaches a relatively stable maximum of about 29. A subset of neurons in the caudal ganglion and caudal nerve cord expresses ChAT mRNA. These putative motoneurons are distributed along nearly the full extent of the tail in numbers consistent with an independent innervation of each tail muscle cell. The longitudinal series of putative motoneurons is not aligned with the muscle cells, but peripheral nerve fibers extending to the muscle cells are, indicating that motor axons grow along the cord before exiting adjacent to their peripheral target. Muscle innervation occurs roughly coincident with the onset of ChAT mRNA expression. Our results provide the first molecular identification of motoneurons and the first developmental characterization of the motor system in an appendicularian and help set the stage for gene expression studies aimed at understanding the evolution of developmental patterning in this part of the chordate central nervous system.