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

Differentiation of endostyle cells by Nkx2-1 and FoxE in the ascidian Ciona intestinalis type A: insights into shared gene regulation in glandular- and thyroid-equivalent elements of the chordate endostyle

Due to similarities in iodine concentrations and peroxidase activities, the thyroid in vertebrates is considered to originate from the endostyle of invertebrate chordates even though it is a glandular (mucus-producing) organ for aquatic suspension feeding. Among chordates with an endostyle, urochordates are useful evolutionary research models for the study of vertebrate traits. The ascidian Ciona intestinalis forms an endostyle with specific components of glandular- and thyroid-related elements, and molecular markers have been identified for these components. Since we previously examined a simple endostyle in the larvacean Oikopleura dioica, the expression of the thyroid-related transcription factor genes, Ciona Nkx2-1 and FoxE, was perturbed by TALEN-mediated gene knockout in the present study to elucidate the shared and/or divergent features of a complex ascidian endostyle. The knockout of Ciona Nkx2-1 and FoxE exerted different effects on the morphology of the developing endostyle. The knockout of Nkx2-1 eliminated the expression of both glandular and thyroidal differentiation marker genes, e.g., vWFL1, vWFL2, CiEnds1, TPO, and Duox, while that of FoxE eliminated the expression of the differentiation marker genes, TPO and CiEnds1. The supporting element-related expression of Pax2/5/8a, Pax2/5/8b, FoxQ1, and β-tubulin persisted in the hypoplastic endostyles of Nkx2-1- and FoxE-knockout juveniles. Although the gene regulation of ascidian-specific CiEnds1 remains unclear, these results provide insights into the evolution of the vertebrate thyroid as well as the urochordate endostyle.

Encyclopaedia of eukaryotic DNA methylation: from patterns to mechanisms and functions

DNA methylation is an epigenetic modification with a very long evolutionary history. However, DNA methylation evolves surprisingly rapidly across eukaryotes. The genome-wide distribution of methylation diversifies rapidly in different lineages, and DNA methylation is lost altogether surprisingly frequently. The growing availability of genomic and epigenomic sequencing across organisms highlights this diversity but also illuminates potential factors that could explain why both the DNA methylation machinery and its genome-wide distribution evolve so rapidly. Key to this are new discoveries about the fitness costs associated with DNA methylation, and new theories about how the fundamental biochemical mechanisms of DNA methylation introduction and maintenance could explain how new genome-wide patterns of methylation evolve.

Metazoans and Intrinsic Apoptosis: An Evolutionary Analysis of the Bcl-2 Family

The B-cell lymphoma-2 (Bcl-2) family is a group of genes regulating intrinsic apoptosis, a process controlling events such as development, homeostasis and the innate and adaptive immune responses in metazoans. In higher organisms, Bcl-2 proteins coordinate intrinsic apoptosis through their regulation of the integrity of the mitochondrial outer membrane; this function appears to have originated in the basal metazoans. Bcl-2 genes predate the cnidarian-bilaterian split and have been identified in porifera, placozoans and cnidarians but not ctenophores and some nematodes. The Bcl-2 family is composed of two groups of proteins, one with an α-helical Bcl-2 fold that has been identified in porifera, placozoans, cnidarians, and almost all higher bilaterians. The second group of proteins, the BH3-only group, has little sequence conservation and less well-defined structures and is found in cnidarians and most bilaterians, but not porifera or placozoans. Here we examine the evolutionary relationships between Bcl-2 proteins. We show that the structures of the Bcl-2-fold proteins are highly conserved over evolutionary time. Some metazoans such as the urochordate Oikopleura dioica have lost all Bcl-2 family members. This gene loss indicates that Bcl-2 regulated apoptosis is not an absolute requirement in metazoans, a finding mirrored in recent gene deletion studies in mice. Sequence analysis suggests that at least some Bcl-2 proteins lack the ability to bind BH3-only antagonists and therefore potentially have other non-apoptotic functions. By examining the foundations of the Bcl-2 regulated apoptosis, functional relationships may be clarified that allow us to understand the role of specific Bcl-2 proteins in evolution and disease.

Evolution of CDK1 Paralog Specializations in a Lineage With Fast Developing Planktonic Embryos

The active site of the essential CDK1 kinase is generated by core structural elements, among which the PSTAIRE motif in the critical αC-helix, is universally conserved in the single CDK1 ortholog of all metazoans. We report serial CDK1 duplications in the chordate, Oikopleura. Paralog diversifications in the PSTAIRE, activation loop substrate binding platform, ATP entrance site, hinge region, and main Cyclin binding interface, have undergone positive selection to subdivide ancestral CDK1 functions along the S-M phase cell cycle axis. Apparent coevolution of an exclusive CDK1d:Cyclin Ba/b pairing is required for oogenic meiosis and early embryogenesis, a period during which, unusually, CDK1d, rather than Cyclin Ba/b levels, oscillate, to drive very rapid cell cycles. Strikingly, the modified PSTAIRE of odCDK1d shows convergence over great evolutionary distance with plant CDKB, and in both cases, these variants exhibit increased specialization to M-phase.

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.

The Genome of the "Sea Vomit" Didemnum vexillum

Tunicates are the sister group of vertebrates and thus occupy a key position for investigations into vertebrate innovations as well as into the consequences of the vertebrate-specific genome duplications. Nevertheless, tunicate genomes have not been studied extensively in the past, and comparative studies of tunicate genomes have remained scarce. The carpet sea squirt Didemnum vexillum, commonly known as “sea vomit”, is a colonial tunicate considered an invasive species with substantial ecological and economical risk. We report the assembly of the D. vexillum genome using a hybrid approach that combines 28.5 Gb Illumina and 12.35 Gb of PacBio data. The new hybrid scaffolded assembly has a total size of 517.55 Mb that increases contig length about eightfold compared to previous, Illumina-only assembly. As a consequence of an unusually high genetic diversity of the colonies and the moderate length of the PacBio reads, presumably caused by the unusually acidic milieu of the tunic, the assembly is highly fragmented (L50 = 25,284, N50 = 6539). It is sufficient, however, for comprehensive annotations of both protein-coding genes and non-coding RNAs. Despite its shortcomings, the draft assembly of the “sea vomit” genome provides a valuable resource for comparative tunicate genomics and for the study of the specific properties of colonial ascidians.

The pericardium of Oikopleura dioica (Tunicata, Appendicularia) contains two distinct cell types and is rotated by 90 degrees to the left

The planktonic Oikopleura dioica belongs to Tunicata, the probable sister taxon to Craniota, and might show plesiomorphic characters, conserved from the common lineage of Tunicata and Craniota. In O. dioica a pericardium in a position similar to other chordates but also to the heart and pericardium of craniates is found. Surprisingly, little is known about the ultrastructure of the pericardium in O. dioica. Here, we show based on electron microscopy that the pericardium is completely lined by a single layer of 16 epithelial cells: 6 epithelial myocardial cells on the left side of the pericardium and 10 peritoneal cells constituting the right side. One of the peritoneal cells, situated at the ventral border between peritoneal cells and myocardial cells has an extension that anchors the pericardium to the basal lamina beneath the latero-ventral epidermis. The primary body cavity of O. dioica appears quite uniformly clear in electron microscopic aspect but several sheets, resembling the basal lamina of the pericardium cross the larger spaces of the body cavity and connect to the pericardial basal lamina. This is the first detailed description of two distinct cell types in the epithelial lining of the pericardium of O. dioica. In comparison with other chordates, we conclude that two cell types can be reconstructed for the last common ancestor of Chordata at least. The position of the pericardium at the intersection of trunk and tail in combination with the basal-lamina like sheets spanning the hemocoel is probably of importance for the function of the circulation of the hemocoelic fluid. Similar to the tail, the axis of the pericardium is shifted through 90 degrees to the left as compared to the main body axis of the trunk and we infer that this shift is an apomorphic character of Appendicularia.

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.

The Evolutionary History of Ephs and Ephrins: Toward Multicellular Organisms

Eph receptor (Eph) and ephrin signaling regulate fundamental developmental processes through both forward and reverse signaling triggered upon cell–cell contact. In vertebrates, they are both classified into classes A and B, and some representatives have been identified in many metazoan groups, where their expression and functions have been well studied. We have extended previous phylogenetic analyses and examined the presence of Eph and ephrins in the tree of life to determine their origin and evolution. We have found that 1) premetazoan choanoflagellates may already have rudimental Eph/ephrin signaling as they have an Eph-/ephrin-like pair and homologs of downstream-signaling genes; 2) both forward- and reverse-downstream signaling might already occur in Porifera since sponges have most genes involved in these types of signaling; 3) the nonvertebrate metazoan Eph is a type-B receptor that can bind ephrins regardless of their membrane-anchoring structure, glycosylphosphatidylinositol, or transmembrane; 4) Eph/ephrin cross-class binding is specific to Gnathostomata; and 5) kinase-dead Eph receptors can be traced back to Gnathostomata. We conclude that Eph/ephrin signaling is of older origin than previously believed. We also examined the presence of protein domains associated with functional characteristics and the appearance and conservation of downstream-signaling pathways to understand the original and derived functions of Ephs and ephrins. We find that the evolutionary history of these gene families points to an ancestral function in cell–cell interactions that could contribute to the emergence of multicellularity and, in particular, to the required segregation of cell populations.

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.

Phylogenetic Analyses of Glycosyl Hydrolase Family 6 Genes in Tunicates: Possible Horizontal Transfer

Horizontal gene transfer (HGT) is the movement of genetic material between different species. Although HGT is less frequent in eukaryotes than in bacteria, several instances of HGT have apparently shaped animal evolution. One well-known example is the tunicate cellulose synthase gene, CesA, in which a gene, probably transferred from bacteria, greatly impacted tunicate evolution. A Glycosyl Hydrolase Family 6 (GH6) hydrolase-like domain exists at the C-terminus of tunicate CesA, but not in cellulose synthases of other organisms. The recent discovery of another GH6 hydrolase-like gene (GH6-1) in tunicate genomes further raises the question of how tunicates acquired GH6. To examine the probable origin of these genes, we analyzed the phylogenetic relationship of GH6 proteins in tunicates and other organisms. Our analyses show that tunicate GH6s, the GH6-1 gene, and the GH6 part of the CesA gene, form two independent, monophyletic gene groups. We also compared their sequence signatures and exon splice sites. All tunicate species examined have shared splice sites in GH6-containing genes, implying ancient intron acquisitions. It is likely that the tunicate CesA and GH6-1 genes existed in the common ancestor of all extant tunicates.

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.

Protein phosphatase 2A is essential to maintain meiotic arrest, and to prevent Ca2+ burst at spawning and eventual parthenogenesis in the larvacean Oikopleura dioica

Unfertilized eggs of most animals are arrested at a certain point in the meiotic cell cycles. Reinitiation of meiosis and the start of embryogenesis are triggered by fertilization. This arrest is essential for preventing parthenogenetic activation and for promoting proper initiation of development by fertilization. In the larvacean Oikopleura dioica, which is a simple model organism for studies of chordate development, the unfertilized egg is arrested at metaphase of meiosis I. We show here that protein phosphatase 2A (PP2A) is essential for maintenance of meiotic arrest after spawning of oocytes. Knockdown (KD) of the maternal PP2A catalytic subunit, which was found in functional screening of maternal factors, caused unfertilized eggs to spontaneously release polar bodies after spawning, and then start pseudo-cleavages without fertilization, namely, parthenogenesis. Parthenogenetic embryos failed to undergo proper mitosis and cytokinesis because of lack of a centrosome, which is to be brought into the egg by a sperm. Activation of the KD oocytes was triggered by possible rise of ambient and intracellular pH upon their release from the gonad into seawater at spawning. Live recording of intracellular calcium level of the KD oocytes indicated that the pH rise caused an aberrant Ca²⁺ burst, which mimicked the Ca²⁺ burst that occurs at fertilization. Then, the aberrant Ca²⁺ burst triggered meiosis resumption through Calcium/calmodulin-dependent protein kinase (CaMK II). Therefore, PP2A is essential for maintenance of meiotic arrest and prevention of parthenogenesis by suppressing the aberrant Ca²⁺ burst at spawning.

Widespread use of the "ascidian" mitochondrial genetic code in tunicates

Background: Ascidians, a tunicate class, use a mitochondrial genetic code that is distinct from vertebrates and other invertebrates. Though it has been used to translate the coding sequences from other tunicate species on a case-by-case basis, it is has not been investigated whether this can be done systematically. This is an important because a) some tunicate mitochondrial sequences are currently translated with the invertebrate code by repositories such as NCBI GenBank, and b) uncertainties about the genetic code to use can complicate or introduce errors in phylogenetic studies based on translated mitochondrial protein sequences.

Methods: We collected publicly available nucleotide sequences for non-ascidian tunicates including appendicularians such as Oikopleura dioica, translated them using the ascidian mitochondrial code, and built multiple sequence alignments covering all tunicate classes.

Results: All tunicates studied here appear to translate AGR codons to glycine instead of serine (invertebrates) or as a stop codon (vertebrates), as initially described in ascidians. Among Oikopleuridae, we suggest further possible changes in the use of the ATA (Ile → Met) and TGA (Trp → Arg) codons.

Conclusions: We recommend using the ascidian mitochondrial code in automatic translation pipelines of mitochondrial sequences for all tunicates. Further investigation is required for additional species-specific differences.

Trans-splicing of mRNAs links gene transcription to translational control regulated by mTOR

BACKGROUND

In phylogenetically diverse organisms, the 5' ends of a subset of mRNAs are trans-spliced with a spliced leader (SL) RNA. The functions of SL trans-splicing, however, remain largely enigmatic.

RESULTS

We quantified translation genome-wide in the marine chordate, Oikopleura dioica, under inhibition of mTOR, a central growth regulator. Translation of trans-spliced TOP mRNAs was suppressed, consistent with a role of the SL sequence in nutrient-dependent translational control of growth-related mRNAs. Under crowded, nutrient-limiting conditions, O. dioica continued to filter-feed, but arrested growth until favorable conditions returned. Upon release from unfavorable conditions, initial recovery was independent of nutrient-responsive, trans-spliced genes, suggesting animal density sensing as a first trigger for resumption of development.

CONCLUSION

Our results are consistent with a proposed role of trans-splicing in the coordinated translational down-regulation of nutrient-responsive genes under growth-limiting conditions.

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.

Spermatogenesis as a tool for staging gonad development in the gonochoric appendicularian Oikopleura dioica Fol 1872

Oikopleura dioica, the only gonochoric species among appendicularians, has a spematozoon with a mid-piece and a conspicuous acrosome that, during fertilisation, undergoes a reaction forming an acrosomal process. To provide more insight into the spermatogenesis of a holoplanktonic tunicate species that completes its life cycle in three to five days, changes in the testis during individual growth have been examined. Spermatogenesis has been subdivided into seven stages based on ultrastructural features during the formation and organisation of the male gonad and the relationships between its macroscopic anatomy and the events of sperm differentiation. Gametes undergo highly synchronised differentiation due to the presence of widespread syncytial structures. Both meiosis and spermiogenesis are brief, and the passage from spermatocytes to spermatids involves a progressive segregation of the germ cells from the syncytial mass with the formation of large cytoplasmic bridges and volume reduction for nucleus compacting and cytoplasmic material changing. The nucleus is small and penetrated anteriorly by a complex acrosome and posteriorly by the distal centriole and part of the flagellum. In spermatids, the single, large mitochondrion appears laterally to the nucleus, and finally, in spermatozoa, it migrates into the mid-piece, wrapping the proximal portion of the axoneme. Because this mitochondrial position is reached only in the late phases of spermatogenesis, it suggests that appendicularians have derived oligopyrenic sperms in which the small nucleus results from adaptation to the assembly of numerous spermatozoa inside the narrow space of the testis compacted in the genital cavity. The formulation of a staging system of gonad development in a model tunicate species known for having the most compacted genome in chordates led to a comparison of histological observations with recent molecular data, improving the characterisation of its biology and life cycle in light of evolutionary implications.

ORTHOSCOPE Analysis Reveals the Presence of the Cellulose Synthase Gene in All Tunicate Genomes but Not in Other Animal Genomes

Tunicates or urochordates—comprising ascidians, larvaceans, and salps—are the only metazoans that can synthesize cellulose, a biological function usually associated with bacteria and plants but not animals. Tunicate cellulose or tunicine is a major component of the outer acellular coverage (tunic) of the entire body of these organisms. Previous studies have suggested that the prokaryotic cellulose synthase gene (CesA) was horizontally transferred into the genome of a tunicate ancestor. However, no convenient tools have been devised to determine whether only tunicates harbor CesA. ORTHOSCOPE is a recently developed tool used to identify orthologous genes and to examine the phylogenic relationship of molecules within major metazoan taxa. The present analysis with this tool revealed the presence of CesA orthologs in all sequenced tunicate genomes but an absence in other metazoan genomes. This supports an evolutionary origin of animal cellulose and provides insights into the evolution of this animal taxon.

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.

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.

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.

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.

Electronic supplementary material

The online version of this article (10.1186/s12861-018-0165-5) contains supplementary material, which is available to authorized users.

Distinct core promoter codes drive transcription initiation at key developmental transitions in a marine chordate

Background

Development is largely driven by transitions between transcriptional programs. The initiation of transcription at appropriate sites in the genome is a key component of this and yet few rules governing selection are known. Here, we used cap analysis of gene expression (CAGE) to generate bp-resolution maps of transcription start sites (TSSs) across the genome of Oikopleura dioica, a member of the closest living relatives to vertebrates.

Results

Our TSS maps revealed promoter features in common with vertebrates, as well as striking differences, and uncovered key roles for core promoter elements in the regulation of development. During spermatogenesis there is a genome-wide shift in mode of transcription initiation characterized by a novel core promoter element. This element was associated with > 70% of male-specific transcription, including the use of cryptic internal promoters within operons. In many cases this led to the exclusion of trans-splice sites, revealing a novel mechanism for regulating which mRNAs receive the spliced leader. Binding of the cell cycle regulator, E2F1, is enriched at the TSS of maternal genes in endocycling nurse nuclei. In addition, maternal promoters lack the TATA-like element found in zebrafish and have broad, rather than sharp, architectures with ordered nucleosomes. Promoters of ribosomal protein genes lack the highly conserved TCT initiator. We also report an association between DNA methylation on transcribed gene bodies and the TATA-box.

Conclusions

Our results reveal that distinct functional promoter classes and overlapping promoter codes are present in protochordates like in vertebrates, but show extraordinary lineage-specific innovations. Furthermore, we uncover a genome-wide, developmental stage-specific shift in the mode of TSS selection. Our results provide a rich resource for the study of promoter structure and evolution in Metazoa.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4504-5) contains supplementary material, which is available to authorized users.

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.

Modified whole-mount in situ hybridisation and immunohistochemistry protocols without removal of the vitelline membrane in the appendicularian Oikopleura dioica

The appendicularian Oikopleura dioica is a planktonic chordate that retains a tadpole shape throughout its life. Its simple and transparent body, invariant cell lineages, fast development and available genome and transcriptome resources make it a promising model organism for research in developmental biology. However, large-scale analysis of gene expression in O. dioica is limited owing to the laborious and time-consuming process of manual removal of the vitelline membrane, because devitellinisation of pre-hatching embryos causes failure of normal development. Therefore, in this study, modified procedures were developed for whole-mount in situ hybridisation (WISH) and immunohistochemistry (WIHC). This protocol enables rapid mRNA or protein detection without a manual devitellination step for each specimen. The critical procedure is brief treatment of the vitelline membrane of living embryos with 0.05% actinase E before fixation. Two minutes of treatment was optimal for the penetration of antisense RNA probes and antibodies through the vitelline membrane. This WISH protocol was applicable for chromogenic and fluorescent tyramide signal amplification reactions. Using the new protocol, we found eight genes with tissue-specific expression in the tail muscle, trunk epidermis, heart, pharynx, oesophagus, stomach or gill openings of developing larvae. This procedure also allowed for the detection of exogenous FLAG-tagged histone-enhanced green fluorescent protein by WIHC using anti-FLAG antibody. This study provides a useful and convenient tool for studying spatial and temporal gene expression patterns in this simple chordate model and should facilitate handling large amounts of genetic data from transcriptome-based approaches and other techniques such as treatments with chemical inhibitors.

Diversity and evolution of TIR-domain-containing proteins in bivalves and Metazoa: New insights from comparative genomics

The Toll/interleukin-1 receptor (TIR) domain has a fundamental role in the innate defence response of plants, vertebrate and invertebrate animals. Mostly found in the cytosolic side of membrane-bound receptor proteins, it mediates the intracellular signalling upon pathogen recognition via heterotypic interactions. Although a number of TIR-domain-containing (TIR-DC) proteins have been characterized in vertebrates, their evolutionary relationships and functional role in protostomes are still largely unknown. Due to the high abundance and diversity of TIR-DC proteins in bivalve molluscs, we investigated this class of marine invertebrates as a case study. The analysis of the available genomic and transcriptomic data allowed the identification of over 400 full-length sequences and their classification in protein families based on sequence homology and domain organization. In addition to TLRs and MyD88 adaptors, bivalves possess a surprisingly large repertoire of intracellular TIR-DC proteins, which are conserved across a broad range of metazoan taxa. Overall, we report the expansion and diversification of TIR-DC proteins in several invertebrate lineages and the identification of many novel protein families possibly involved in both immune-related signalling and embryonic development.

Genome-wide survey of miRNAs and their evolutionary history in the ascidian, Halocynthia roretzi

BACKGROUND

miRNAs play essential roles in the modulation of cellular functions via degradation and/or translation attenuation of target mRNAs. They have been surveyed in a single ascidian genus, Ciona. Recently, an annotated draft genome sequence for a distantly related ascidian, Halocynthia roretzi, has become available, but miRNAs in H. roretzi have not been previously studied.

RESULTS

We report the prediction of 319 candidate H. roretzi miRNAs, obtained through three complementary methods. Experimental validation suggests that more than half of these candidate miRNAs are expressed during embryogenesis. The majority of predicted H. roretzi miRNAs appear specific to ascidians or tunicates, and only 32 candidates, belonging to 25 families, are widely conserved across metazoans.

CONCLUSION

Our study presents a comprehensive identification of candidate H. roretzi miRNAs. This resource will facilitate the study of the mechanisms for miRNA-controlled gene regulatory networks during ascidian development. Further, our analysis suggests that the majority of Halocynthia miRNAs are specific to ascidian or tunicates, with only a small number of widely conserved miRNAs. This result is consistent with the general notion that animal miRNAs are less conserved between taxa than plant ones.

DNA interference-mediated screening of maternal factors in the chordate Oikopleura dioica

The maternal contribution to the oocyte cytoplasm plays an important role during embryogenesis because it is involved in early cell fate specification and embryonic axis establishment. However, screening projects targeting maternal factors have only been conducted in a limited number of animal models, such as nematodes, fruit flies, and zebrafish, while few maternal genes have been analysed because of difficulties encountered in inhibiting gene products already expressed in the ovaries. Therefore, simple and efficient methods for large-scale maternal screening are necessary. The appendicularian Oikopleura dioica is a planktonic tunicate member of the chordates. Gonadal microinjection and a novel gene knockdown method, DNA interference (DNAi), have been developed for use in this animal with the aim of inhibiting gene functions during oogenesis within the gonad. In this study, we adapted these methods for large-scale maternal factor screening, and observed malformation phenotypes related to some maternal factors. Approximately 2000 (56.9%) ovary-enriched gene products were screened, of which the knockdown of seven encoding genes resulted in various abnormalities during embryonic development. Most of these were related to microtubules and cell adhesion-related proteins. We conclude that DNAi is a potentially powerful screening tool for the identification of novel maternal factors in chordates.

Sex-specific chromatin landscapes in an ultra-compact chordate genome

Background

In multicellular organisms, epigenome dynamics are associated with transitions in the cell cycle, development, germline specification, gametogenesis and inheritance. Evolutionarily, regulatory space has increased in complex metazoans to accommodate these functions. In tunicates, the sister lineage to vertebrates, we examine epigenome adaptations to strong secondary genome compaction, sex chromosome evolution and cell cycle modes.

Results

Across the 70 MB Oikopleura dioica genome, we profiled 19 histone modifications, and RNA polymerase II, CTCF and p300 occupancies, to define chromatin states within two homogeneous tissues with distinct cell cycle modes: ovarian endocycling nurse nuclei and mitotically proliferating germ nuclei in testes. Nurse nuclei had active chromatin states similar to other metazoan epigenomes, with large domains of operon-associated transcription, a general lack of heterochromatin, and a possible role of Polycomb PRC2 in dosage compensation. Testis chromatin states reflected transcriptional activity linked to spermatogenesis and epigenetic marks that have been associated with establishment of transgenerational inheritance in other organisms. We also uncovered an unusual chromatin state specific to the Y-chromosome, which combined active and heterochromatic histone modifications on specific transposable elements classes, perhaps involved in regulating their activity.

Conclusions

Compacted regulatory space in this tunicate genome is accompanied by reduced heterochromatin and chromatin state domain widths. Enhancers, promoters and protein-coding genes have conserved epigenomic features, with adaptations to the organization of a proportion of genes in operon units. We further identified features specific to sex chromosomes, cell cycle modes, germline identity and dosage compensation, and unusual combinations of histone PTMs with opposing consensus functions.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-016-0110-4) contains supplementary material, which is available to authorized users.

DNA interference: DNA-induced gene silencing in the appendicularian Oikopleura dioica

RNA interference is widely employed as a gene-silencing system in eukaryotes for host defence against invading nucleic acids. In response to invading double-stranded RNA (dsRNA), mRNA is degraded in sequence-specific manner. So far, however, DNA interference (DNAi) has been reported only in plants, ciliates and archaea, and has not been explored in Metazoa. Here, we demonstrate that linear double-stranded DNA promotes both sequence-specific transcription blocking and mRNA degradation in developing embryos of the appendicularian Oikopleura dioica. Introduced polymerase chain reaction (PCR) products or linearized plasmids encoding Brachyury induced tail malformation and mRNA degradation. This malformation was also promoted by DNA fragments of the putative 5'-flanking region and intron without the coding region. PCR products encoding Zic-like1 and acetylcholine esterase also induced loss of sensory organ and muscle acetylcholinesterase activity, respectively. Co-injection of mRNA encoding EGFP and mCherry, and PCR products encoding these fluorescent proteins, induced sequence-specific decrease in the green or red fluorescence, respectively. These results suggest that O. dioica possesses a defence system against exogenous DNA and RNA, and that DNA fragment-induced gene silencing would be mediated through transcription blocking as well as mRNA degradation. This is the first report of DNAi in Metazoa.

Maternal and zygotic transcriptomes in the appendicularian, Oikopleura dioica: novel protein-encoding genes, intra-species sequence variations, and trans-spliced RNA leader

RNA sequencing analysis was carried out to characterize egg and larval transcriptomes in the appendicularian, Oikopleura dioica, a planktonic chordate, which is characterized by rapid development and short life cycle of 5 days, using a Japanese population of the organism. De novo transcriptome assembly matched with 16,423 proteins corresponding to 95.4% of the protein-encoding genes deposited in the OikoBase, the genome database of the Norwegian population. Nucleotide and amino acid sequence identities between the Japanese and Norwegian O. dioica were estimated to be around 91.0 and 94.8%, respectively. We discovered 175 novel protein-encoding genes: 144 unigenes were common to both the Japanese and Norwegian populations, whereas 31 unigenes were not found in the OikoBase genome reference. Among the total 12,311 unigenes, approximately 63% were detected in egg-stage RNAs, whereas 99% were detected in larval stage RNAs; 3772 genes were up-regulated, and 1336 genes were down-regulated more than four-fold in the larvae. Gene ontology analyses characterized gene activities in these two developmental stages. We found a messenger RNA (mRNA) 5' trans-spliced leader, which was observed in 40.8% of the total unique transcripts. It showed preferential linkage to adenine at the 5' ends of the downstream exons. Trans-splicing was observed more frequently in egg mRNAs compared with larva-specific mRNAs.

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.

Trans-splicing and operons in metazoans: translational control in maternally regulated development and recovery from growth arrest

Polycistronic mRNAs transcribed from operons are resolved via the trans-splicing of a spliced-leader (SL) RNA. Trans-splicing also occurs at monocistronic transcripts. The phlyogenetically sporadic appearance of trans-splicing and operons has made the driving force(s) for their evolution in metazoans unclear. Previous work has proposed that germline expression drives operon organization in Caenorhabditis elegans, and a recent hypothesis proposes that operons provide an evolutionary advantage via the conservation of transcriptional machinery during recovery from growth arrested states. Using a modified cap analysis of gene expression protocol we mapped sites of SL trans-splicing genome-wide in the marine chordate Oikopleura dioica. Tiled microarrays revealed the expression dynamics of trans-spliced genes across development and during recovery from growth arrest. Operons did not facilitate recovery from growth arrest in O. dioica. Instead, we found that trans-spliced transcripts were predominantly maternal. We then analyzed data from C. elegans and Ciona intestinalis and found that an enrichment of trans-splicing and operon gene expression in maternal mRNA is shared between all three species, suggesting that this may be a driving force for operon evolution in metazoans. Furthermore, we found that the majority of known terminal oligopyrimidine (TOP) mRNAs are trans-spliced in O. dioica and that the SL contains a TOP-like motif. This suggests that the SL in O. dioica confers nutrient-dependent translational control to trans-spliced mRNAs via the TOR-signaling pathway. We hypothesize that SL-trans-splicing provides an evolutionary advantage in species that depend on translational control for regulating early embryogenesis, growth and oocyte production in response to nutrient levels.

Retroduplication and loss of parental genes is a mechanism for the generation of intronless genes in Ciona intestinalis and Ciona savignyi

Tunicates, the sister clade of vertebrates, have miniature genomes and numerous intronless genes compared to other animals. It is still unclear how the tunicates acquired such a large number of intronless genes. Here, we analyzed sequences and intron-exon organizations of homologous genes from two closely related tunicates, Ciona intestinalis and Ciona savignyi. We found seven cases in which ancestral introns of a gene were completely lost in a species after their divergence. In four cases, both the intronless copy and the intron-containing copy were present in the genome, indicating that the intronless copy was generated by retroduplication. In the other three cases, the intron-containing copy was absent, implying it was lost after retroduplication. This result suggests that retroduplication and loss of parental genes is a major mechanism for the accumulation of intronless genes in tunicates.

Advances in biotechnology and informatics to link variation in the genome to phenotypes in plants and animals

Advances in our understanding of genome structure provide consistent evidence for the existence of a core genome representing species classically defined by phenotype, as well as conditionally dispensable components of the genome that shows extensive variation between individuals of a given species. Generally, conservation of phenotypic features between species reflects conserved features of the genome; however, this is evidently not necessarily always the case as demonstrated by the analysis of the tunicate chordate Oikopleura dioica. In both plants and animals, the methylation activity of DNA and histones continues to present new variables for modifying (eventually) the phenotype of an organism and provides for structural variation that builds on the point mutations, rearrangements, indels, and amplification of retrotransposable elements traditionally considered. The translation of the advances in the structure/function analysis of the genome to industry is facilitated through the capture of research outputs in "toolboxes" that remain accessible in the public domain.

The 2013 Nucleic Acids Research Database Issue and the online molecular biology database collection

The 20th annual Database Issue of Nucleic Acids Research includes 176 articles, half of which describe new online molecular biology databases and the other half provide updates on the databases previously featured in NAR and other journals. This year’s highlights include two databases of DNA repeat elements; several databases of transcriptional factors and transcriptional factor-binding sites; databases on various aspects of protein structure and protein–protein interactions; databases for metagenomic and rRNA sequence analysis; and four databases specifically dedicated to Escherichia coli. The increased emphasis on using the genome data to improve human health is reflected in the development of the databases of genomic structural variation (NCBI’s dbVar and EBI’s DGVa), the NIH Genetic Testing Registry and several other databases centered on the genetic basis of human disease, potential drugs, their targets and the mechanisms of protein–ligand binding. Two new databases present genomic and RNAseq data for monkeys, providing wealth of data on our closest relatives for comparative genomics purposes. The NAR online Molecular Biology Database Collection, available at http://www.oxfordjournals.org/nar/database/a/, has been updated and currently lists 1512 online databases. The full content of the Database Issue is freely available online on the Nucleic Acids Research website (http://nar.oxfordjournals.org/).