Germline development during embryogenesis of the larvacean, Oikopleura dioica.

Germ cells develop into eggs and sperms and represent a lineage that survives through multiple generations. Germ cell specification during embryogenesis proceeds through one of two basic modes: either the cell-autonomous mode or the inductive mode. In the cell-autonomous mode, specification of germ cell fate involves asymmetric partitioning of the specialized maternal cytoplasm, known as the germplasm. Oikopleura dioica is a larvacean (class Appendicularia) and a chordate. It is regarded as a promising animal model for studying chordate development because of its short life cycle (5 days) and small genome size (∼60 â€‹Mb). We show that their embryos possess germplasm, as observed in ascidians (class Ascidiacea). The vegetal cytoplasm shifted towards the future posterior pole before the first cleavage occurred. A bilateral pair of primordial germ cells (PGC, B11 â€‹cells) was formed at the posterior pole at the 32-cell stage through two rounds of unequal cleavage. These B11 â€‹cells did not undergo further division before hatching of the tadpole-shaped larvae. The centrosome-attracting body (CAB) is a subcellular structure that contains the germplasm and plays crucial roles in germ cell development in ascidians. The presence of CAB with germplasm was observed in the germline lineage cells of larvaceans via electron microscopy and using extracted embryos. The CAB appeared at the 8-cell stage and persisted until the middle stage of embryogenesis. The antigen for the phosphorylated histone 3 antibody was localized to the CAB and persisted in the PGC until hatching after the CAB disappeared. Maternal snail mRNA, which encodes a transcription factor, was co-localized with the antigen for the H3S28p antibody. Furthermore, we found a novel PGC-specific subcellular structure that we call the germ body (GB). This study thus highlights the conserved and non-conserved features of germline development between ascidians and larvaceans. The rapid development and short life cycle (five days) of O. dioica would open the way to genetically analyze germ cell development in the future.

Formation of the brain by stem cell divisions of large neuroblasts in Oikopleura dioica, a simple chordate.

Stem cell division contributes to the generation of various cell types during animal development, especially a diverse pool of neural cells in the nervous system. One example is reiterated unequal stem cell divisions, in which a large stem cell undergoes a series of oriented unequal divisions to produce a chain of small daughter cells that differentiate. We show that reiterated unequal stem cell divisions are involved in the formation of the brain in simple chordate appendicularians (larvaceans). Two large neuroblasts in the anterior and middle of the brain-forming region of hatched larvae were observed. They produced at least 30 neural cells out of 96 total brain cells before completion of brain formation at 10 hours after fertilization by reiterated unequal stem cell divisions. The daughter cells of the anterior neuroblast were postmitotic, and the number was at least 19. The neuroblast produced small daughter neural cells posteriorly every 20 min. The neural cells first moved toward the dorsal side, turned in the anterior direction, aligned in a single line according to their birth order, and showed collective movement to accumulate in the anterior part of the brain. The anterior neuroblast originated from the right-anterior blastomeres of the eight-cell embryos and the right a222 blastomere of the 64-cell embryo. The posterior neuroblast also showed reiterated unequal stem cell divisions, and generated at least 11 neural cells. Sequential unequal stem cell divisions without stem cell growth have been observed in protostomes, such as insects and annelids. The results provide the first examples of this kind of stem cell division during brain formation in non-vertebrate deuterostomes.

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) Ca2+ 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 Ca2+ 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 Ca2+ 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.

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.

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.

Gap junction-dependent coordination of intercellular calcium signalling in the developing appendicularian tunicate Oikopleura dioica.

We characterized spontaneous Ca2+ signals in Oikopleura dioica embryos from pre-fertilization to gastrula stages following injection of GCaMP6 mRNA into unfertilized eggs. The unfertilized egg exhibited regular, transient elevations in intracellular Ca2+ concentration with an average duration of 4-6 s and an average frequency of about 1 every 2.5 min. Fertilization was accompanied by a longer Ca2+ transient that lasted several minutes. Thereafter, regular Ca2+ transients were reinstated that spread within seconds among blastomeres and gradually increased in duration (by about 50%) and decreased in frequency (by about 20%) by gastrulation. Peak amplitudes also exhibited a dynamic, with a transitory drop occurring at about the 4-cell stage and a subsequent rise. Each peak was preceded by about 15 s by a smaller and shorter Ca2+ increase (about 5% of the main peak amplitude, average duration 3 s), which we term the "minipeak". By gastrulation, Ca2+ transients exhibited a stereotyped initiation site on either side of the 32-64-cell embryo, likely in the nascent muscle precursor cells, and spread thereafter symmetrically in a stereotyped spatial pattern that engaged blastomeres giving rise to all the major tissue lineages. The rapid spread of the transients relative to the intertransient interval created a coordinated wave that, on a coarse time scale, could be considered an approximate synchronization. Treatment with the divalent cations Ni2+ or Cd2+ gradually diminished peak amplitudes, had only moderate effects on wave frequency, but markedly disrupted wave synchronization and normal development. The T-type Ca2+ channel blocker mibefradil similarly disrupted normal development, and eliminated the minipeaks, but did not affect wave synchronization. To assess the role of gap junctions in calcium wave spread and coordination, we first characterized the expression of two Oikopleura connexins, Od-CxA and Od-CxB, both of which are expressed during pre-gastrulation and gastrula stages, and then co-injected double-stranded inhibitory RNAs together with CGaMP6 to suppress connexin expression. Connexin mRNA knockdown led to a gradual increase in Ca2+ transient peak width, a decrease of interpeak interval and a marked disruption of wave synchronization. As seen with divalent cations and mibefradil, this desynchronization was accompanied by a disruption of normal development.

Mouth opening is mediated by separation of dorsal and ventral daughter cells of the lip precursor cells in the larvacean, Oikopleura dioica.

Mouth formation involves the processes of mouth opening, formation of the oral cavity, and the development of associated sensory organs. In deuterostomes, the surface ectoderm and the anterior part of the archenteron are reconfigured and reconnected to make a mouth opening. This study of the larval development of the larvacean, Oikopleura dioica, investigates the cellular organization of the oral region, the developmental processes of the mouth, and the formation of associated sensory cells. O. dioica is a simple chordate whose larvae are transparent and have a small number of constituent cells. It completes organ morphogenesis in 7 h, between hatching 3 h after fertilization and the juvenile stage at 10 h, when it attains adult form and starts to feed. It has two types of mechanosensory cell embedded in the oral epithelium, which is a single layer of cells. There are twenty coronal sensory cells in the circumoral nerve ring and two dorsal sensory organ cells. Two bilateral lip precursor cells (LPCs), facing the anterior surface, divide dorsoventrally and make a wedge-shaped cleft between the two daughter cells named the dorsal lip cell (DLC) and the ventral lip cell (VLC). Eventually, the DLC and VLC become detached and separated into dorsal and ventral lips, triggering mouth opening. This is an intriguing example of cell division itself contributing to morphogenesis. The boundary between the ectoderm and endoderm is present between the lip cells and coronal sensory cells. All oral sensory cells, including dorsal sensory organ cells, were of endodermal origin and were not derived from the ectodermal placode. These observations on mouth formation provide a cellular basis for further studies at a molecular level, in this simple chordate.

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.

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.

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.

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.

Tracing Homopolymers in Oikopleura dioica's Mitogenome.

Oikopleura dioica is a planktonic tunicate (Appendicularia class) found extensively across the marine waters of the globe. The genome of a single male individual collected from Okinawa, Japan was sequenced using the single-molecule PacBio Hi-Fi method and assembled with NOVOLoci. The mitogenome is 39,268 bp long, featuring a large control region of around 22,000 bp. We annotated the proteins atp6, cob, cox1, cox2, cox3, nad1, nad4, and nad5, and found one more open reading frame that did not match any known gene. This study marks the first complete mitogenome assembly for an appendicularian, and reveals that A and T homopolymers cumulatively account for nearly half of its length. This reference sequence will be an asset for environmental DNA and phylogenetic studies.

Diffusion tubes: a method for the mass culture of ctenophores and other pelagic marine invertebrates.

The culture of pelagic marine invertebrates, especially the ctenophore Mnemiopsis leidyi, has been demonstrated in past studies dating back to the 1960s; however, the mass culture of delicate pelagic invertebrates has remained elusive. By using a pair of acrylic tubes and enabling water diffusion between them, we have been able to reliably and cost effectively mass culture several genera of ctenophores (Pleurobrachia, Hormiphora, Bolinopsis, Mnemiopsis and Leucothea), one species of siphonophore (Nanomia) and one species of larvacean (Oikopleura). The simple, compact method is effective enough to support two permanent exhibits of ctenophores at the Monterey Bay Aquarium while minimizing live food culture requirements with the potential to support further investigation of pelagic marine invertebrate ontogeny, ecology and genomics.

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.