Miniature genome in the marine chordate Oikopleura dioica

De novo genome assembly and comparative genomics for the colonial ascidian Botrylloides violaceus

Ascidians have the potential to reveal fundamental biological insights related to coloniality, regeneration, immune function, and the evolution of these traits. This study implements a hybrid assembly technique to produce a genome assembly and annotation for the botryllid ascidian, Botrylloides violaceus. A hybrid genome assembly was produced using Illumina, Inc. short and Oxford Nanopore Technologies long-read sequencing technologies. The resulting assembly is comprised of 831 contigs, has a total length of 121 Mbp, N50 of 1 Mbp, and a BUSCO score of 96.1%. Genome annotation identified 13 K protein-coding genes. Comparative genomic analysis with other tunicates reveals patterns of conservation and divergence within orthologous gene families even among closely related species. Characterization of the Wnt gene family, encoding signaling ligands involved in development and regeneration, reveals conserved patterns of subfamily presence and gene copy number among botryllids. This supports the use of genomic data from nonmodel organisms in the investigation of biological phenomena.

A miniaturized nigrostriatal-like circuit regulating locomotor performance in a protochordate

To gain insight into the evolution of motor control systems at the origin of vertebrates, we have investigated higher-order motor circuitry in the protochordate Oikopleura dioica. We have identified a highly miniaturized circuit in Oikopleura with a projection from a single pair of dopaminergic neurons to a small set of synaptically coupled GABAergic neurons, which in turn exert a disinhibitory descending projection onto the locomotor central pattern generator. The circuit is reminiscent of the nigrostriatopallidal system in the vertebrate basal ganglia, in which disinhibitory circuits release specific movements under the modulatory control of dopamine. We demonstrate further that dopamine is required to optimize locomotor performance in Oikopleura, mirroring its role in vertebrates. A dopamine-regulated disinhibitory locomotor control circuit reminiscent of the vertebrate nigrostriatopallidal system was thus already present at the origin of ancestral chordates and has been maintained in the face of extreme nervous system miniaturization in the urochordate lineage.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

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.

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.

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.

The filter-house of the larvacean Oikopleura dioica. A complex extracellular architecture: From fiber production to rudimentary state to inflated house

While cellulose is the most abundant macromolecule in the biosphere, most animals are unable to produce cellulose with the exception of tunicates. Some tunicates have evolved the ability to secrete a complex house containing cellulosic fibers, yet little is known about the early stages of the house building process. Here, we investigate the rudimentary house of Oikopleura dioica for the first time using complementary light and electron microscopic techniques. In addition, we digitally modeled the arrangement of chambers, nets, and filters of the functional, expanded house in three dimensions based on life-video-imaging. Combining 3D-reconstructions based on serial histological semithin-sections, confocal laser scanning microscopy, transmission electron microscopy, scanning electron microscopy (SEM), and focused ion beam (FIB)-SEM, we were able to elucidate the arrangement of structural components, including cellulosic fibers, of the rudimentary house with a focus on the food concentration filter. We developed a model for the arrangement of folded structures in the house rudiment and show it is a precisely preformed structure with identifiable components intricately correlated with specific cells. Moreover, we demonstrate that structural details of the apical surfaces of Nasse cells provide the exact locations and shapes to produce the fibers of the house and interact among each other, with Giant Fol cells, and with the fibers to arrange them in the precise positions necessary for expansion of the house rudiment into the functional state. The presented data and hypotheses advance our knowledge about the interrelation of structure and function on different biological levels and prompt investigations into this astonishing biological object.

Telomere-to-telomere assembly of the genome of an individual Oikopleura dioica from Okinawa using Nanopore-based sequencing

Background

The larvacean Oikopleura dioica is an abundant tunicate plankton with the smallest (65–70 Mbp) non-parasitic, non-extremophile animal genome identified to date. Currently, there are two genomes available for the Bergen (OdB3) and Osaka (OSKA2016) O. dioica laboratory strains. Both assemblies have full genome coverage and high sequence accuracy. However, a chromosome-scale assembly has not yet been achieved.

Results

Here, we present a chromosome-scale genome assembly (OKI2018_I69) of the Okinawan O. dioica produced using long-read Nanopore and short-read Illumina sequencing data from a single male, combined with Hi-C chromosomal conformation capture data for scaffolding. The OKI2018_I69 assembly has a total length of 64.3 Mbp distributed among 19 scaffolds. 99% of the assembly is contained within five megabase-scale scaffolds. We found telomeres on both ends of the two largest scaffolds, which represent assemblies of two fully contiguous autosomal chromosomes. Each of the other three large scaffolds have telomeres at one end only and we propose that they correspond to sex chromosomes split into a pseudo-autosomal region and X-specific or Y-specific regions. Indeed, these five scaffolds mostly correspond to equivalent linkage groups in OdB3, suggesting overall agreement in chromosomal organization between the two populations. At a more detailed level, the OKI2018_I69 assembly possesses similar genomic features in gene content and repetitive elements reported for OdB3. The Hi-C map suggests few reciprocal interactions between chromosome arms. At the sequence level, multiple genomic features such as GC content and repetitive elements are distributed differently along the short and long arms of the same chromosome.

Conclusions

We show that a hybrid approach of integrating multiple sequencing technologies with chromosome conformation information results in an accurate de novo chromosome-scale assembly of O. dioica’s highly polymorphic genome. This genome assembly opens up the possibility of cross-genome comparison between O. dioica populations, as well as of studies of chromosomal evolution in this lineage.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07512-6.

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.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

Genomic basis of environmental adaptation in the leathery sea squirt (Styela clava)

Tunicates occupy the evolutionary position at the boundary of invertebrates and vertebrates. It exhibits adaptation to broad environmental conditions and is distributed globally. Despite hundreds of years of embryogenesis studies, the genetic basis of the invasive habits of ascidians remains largely unknown. The leathery sea squirt, Styela clava, is an important invasive species. We used the chromosomal-level genome and transcriptome of S. clava to explore its genomic- and molecular-network-based mechanisms of adaptation to environments. Compared with Ciona intestinalis type A (C. robusta), the size of the S. clava genome was expanded by 2-fold, although the gene number was comparable. An increase in transposon number and variation in dominant types were identified as potential expansion mechanisms. In the S. clava genome, the number of genes encoding the heat-shock protein 70 family and members of the complement system was expanded significantly, and cold-shock protein genes were transferred horizontally into the S. clava genome from bacteria. The expanded gene families potentially play roles in the adaptation of S. clava to its environments. The loss of key genes in the galactan synthesis pathway might explain the distinct tunic structure and hardness compared with the ascidian Ciona species. We demonstrated further that the integrated thyroid hormone pathway participated in the regulation of larval metamorphosis that provides S. clava with two opportunities for adapting to their environment. Thus, our report of the chromosomal-level leathery sea squirt genome provides a comprehensive genomic basis for the understanding of environmental adaptation in tunicates.

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.

H3S28P Antibody Staining of Okinawan Oikopleura dioica Suggests the Presence of Three Chromosomes

Oikopleura dioica is a ubiquitous marine zooplankton of biological interest owing to features that include dioecious reproduction, a short life cycle, conserved chordate body plan, and a compact genome. It is an important tunicate model for evolutionary and developmental research, as well as investigations into marine ecosystems. The genome of north Atlantic O. dioica comprises three chromosomes. However, comparisons with the genomes of O. dioica sampled from mainland and southern Japan revealed extensive sequence differences. Moreover, historical studies have reported widely varying chromosome counts. We recently initiated a project to study the genomes of O. dioica individuals collected from the coastline of the Ryukyu (Okinawa) Islands in southern Japan. Given the potentially large extent of genomic diversity, we employed karyological techniques to count individual animals' chromosomes in situ using centromere-specific antibodies directed against H3S28P, a prophase-metaphase cell cycle-specific marker of histone H3. Epifluorescence and confocal images were obtained of embryos and oocytes stained with two commercial anti-H3S28P antibodies (Abcam ab10543 and Thermo Fisher 07-145). The data lead us to conclude that diploid cells from Okinawan O. dioica contain three pairs of chromosomes, in line with the north Atlantic populations. The finding facilitates the telomere-to-telomere assembly of Okinawan O. dioica genome sequences and gives insight into the genomic diversity of O. dioica from different geographical locations. The data deposited in the EBI BioImage Archive provide representative images of the antibodies' staining properties for use in epifluorescent and confocal based fluorescent microscopy.

Inferring Tunicate Relationships and the Evolution of the Tunicate Hox Cluster with the Genome of Corella inflata

Tunicates, the closest living relatives of vertebrates, have served as a foundational model of early embryonic development for decades. Comparative studies of tunicate phylogeny and genome evolution provide a critical framework for analyzing chordate diversification and the emergence of vertebrates. Toward this goal, we sequenced the genome of Corella inflata (Ascidiacea, Phlebobranchia), so named for the capacity to brood self-fertilized embryos in a modified, "inflated" atrial chamber. Combining the new genome sequence for Co. inflata with publicly available tunicate data, we estimated a tunicate species phylogeny, reconstructed the ancestral Hox gene cluster at important nodes in the tunicate tree, and compared patterns of gene loss between Co. inflata and Ciona robusta, the prevailing tunicate model species. Our maximum-likelihood and Bayesian trees estimated from a concatenated 210-gene matrix were largely concordant and showed that Aplousobranchia was nested within a paraphyletic Phlebobranchia. We demonstrated that this relationship is not an artifact due to compositional heterogeneity, as had been suggested by previous studies. In addition, within Thaliacea, we recovered Doliolida as sister to the clade containing Salpida and Pyrosomatida. The Co. inflata genome provides increased resolution of the ancestral Hox clusters of key tunicate nodes, therefore expanding our understanding of the evolution of this cluster and its potential impact on tunicate morphological diversity. Our analyses of other gene families revealed that several cardiovascular associated genes (e.g., BMP10, SCL2A12, and PDE2a) absent from Ci. robusta, are present in Co. inflata. Taken together, our results help clarify tunicate relationships and the genomic content of key ancestral nodes within this phylogeny, providing critical insights into tunicate evolution.

Animal, Fungi, and Plant Genome Sequences Harbor Different Non-Canonical Splice Sites

Most protein-encoding genes in eukaryotes contain introns, which are interwoven with exons. Introns need to be removed from initial transcripts in order to generate the final messenger RNA (mRNA), which can be translated into an amino acid sequence. Precise excision of introns by the spliceosome requires conserved dinucleotides, which mark the splice sites. However, there are variations of the highly conserved combination of GT at the 5' end and AG at the 3' end of an intron in the genome. GC-AG and AT-AC are two major non-canonical splice site combinations, which have been known for years. Recently, various minor non-canonical splice site combinations were detected with numerous dinucleotide permutations. Here, we expand systematic investigations of non-canonical splice site combinations in plants across eukaryotes by analyzing fungal and animal genome sequences. Comparisons of splice site combinations between these three kingdoms revealed several differences, such as an apparently increased CT-AC frequency in fungal genome sequences. Canonical GT-AG splice site combinations in antisense transcripts are a likely explanation for this observation, thus indicating annotation errors. In addition, high numbers of GA-AG splice site combinations were observed in Eurytemoraaffinis and Oikopleuradioica. A variant in one U1 small nuclear RNA (snRNA) isoform might allow the recognition of GA as a 5' splice site. In depth investigation of splice site usage based on RNA-Seq read mappings indicates a generally higher flexibility of the 3' splice site compared to the 5' splice site across animals, fungi, and plants.

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.

Exon 3 of the NUMB Gene Emerged in the Chordate Lineage Coopting the NUMB Protein to the Regulation of MDM2

MDM2 regulates a variety of cellular processes through its dual protein:protein interaction and ubiquitin ligase activities. One major function of MDM2 is to bind and ubiquitinate P53, thereby regulating its proteasomal degradation. This function is in turn controlled by the cell fate determinant NUMB, which binds to and inhibits MDM2 via a short stretch of 11 amino acids, contained in its phosphotyrosine-binding (PTB) domain, encoded by exon 3 of the NUMB gene. The NUMB-MDM2-P53 circuitry is relevant to the specification of the stem cell fate and its subversion has been shown to be causal in breast cancer leading to the emergence of cancer stem cells. While extensive work on the evolutionary aspects of the MDM2/P53 circuitry has provided hints as to how these two proteins have evolved together to maintain conserved and linked functions, little is known about the evolution of the NUMB gene and, in particular, how it developed the ability to regulate MDM2 function. Here, we show that NUMB is a metazoan gene, which acquired exon 3 in the common ancestor of the Chordate lineage, first being present in the Cephalochordate and Tunicate subphyla, but absent in invertebrates. We provide experimental evidence showing that since its emergence, exon 3 conferred to the PTB domain of NUMB the ability to bind and to regulate MDM2 functions.

Switching of INCENP paralogs controls transitions in mitotic chromosomal passenger complex functions

A single inner centromere protein (INCENP) found throughout eukaryotes modulates Aurora B kinase activity and chromosomal passenger complex (CPC) localization, which is essential for timely mitotic progression. It has been proposed that INCENP might act as a rheostat to regulate Aurora B activity through mitosis, with successively higher activity threshold levels for chromosome alignment, the spindle checkpoint, anaphase spindle transfer and finally spindle elongation and cytokinesis. It remains mechanistically unclear how this would be achieved. Here, we reveal that the urochordate, Oikopleura dioica, possesses two INCENP paralogs, which display distinct localizations and subfunctionalization in order to complete M-phase. INCENPa was localized on chromosome arms and centromeres by prometaphase, and modulated Aurora B activity to mediate H3S10/S28 phosphorylation, chromosome condensation, spindle assembly and transfer of the CPC to the central spindle. Polo-like kinase (Plk1) recruitment to CDK1 phosphorylated INCENPa was crucial for INCENPa-Aurora B enrichment on centromeres. The second paralog, INCENPb was enriched on centromeres from prometaphase, and relocated to the central spindle at anaphase onset. In the absence of INCENPa, meiotic spindles failed to form, and homologous chromosomes did not segregate. INCENPb was not required for early to mid M-phase events but became essential for the activity and localization of Aurora B on the central spindle and midbody during cytokinesis in order to allow abscission to occur. Together, our results demonstrate that INCENP paralog switching on centromeres modulates Aurora B kinase localization, thus chronologically regulating CPC functions during fast embryonic divisions in the urochordate O. dioica. Abbreviations: CCAN: constitutive centromere-associated network; CENPs: centromere proteins; cmRNA: capped messenger RNA; CPC: chromosomal passenger complex; INCENP: inner centromere protein; Plk1: polo-like kinase 1; PP1: protein phosphatase 1; PP2A: protein phosphatase 2A; SAC: spindle assembly checkpoint; SAH: single α-helix domain.

Regulation and evolution of muscle development in tunicates

For more than a century, studies on tunicate muscle formation have revealed many principles of cell fate specification, gene regulation, morphogenesis, and evolution. Here, we review the key studies that have probed the development of all the various muscle cell types in a wide variety of tunicate species. We seize this occasion to explore the implications and questions raised by these findings in the broader context of muscle evolution in chordates.

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.

Developmental system drift in motor ganglion patterning between distantly related tunicates

BACKGROUND

The larval nervous system of the solitary tunicate Ciona is a simple model for the study of chordate neurodevelopment. The development and connectivity of the Ciona motor ganglion have been studied in fine detail, but how this important structure develops in other tunicates is not well known.

METHODS AND RESULTS

By comparing gene expression patterns in the developing MG of the distantly related tunicate Molgula occidentalis, we found that its patterning is highly conserved compared to the Ciona MG. MG neuronal subtypes in Molgula were specified in the exact same positions as in Ciona, though the timing of subtype-specific gene expression onset was slightly shifted to begin earlier, relative to mitotic exit and differentiation. In transgenic Molgula embryos electroporated with Dmbx reporter plasmids, we were also able to characterize the morphology of the lone pair of descending decussating neurons (ddNs) in Molgula, revealing the same unique contralateral projection seen in Ciona ddNs and their putative vertebrate homologs the Mauthner cells. Although Dmbx expression labels the ddNs in both species, cross-species transgenic assays revealed significant changes to the regulatory logic underlying Dmbx transcription. We found that Dmbx cis-regulatory DNAs from Ciona can drive highly specific reporter gene expression in Molgula ddNs, but Molgula sequences are not active in Ciona ddNs.

CONCLUSIONS

This acute divergence in the molecular mechanisms that underlie otherwise functionally conserved cis-regulatory DNAs supports the recently proposed idea that the extreme genetic plasticity observed in tunicates may be attributed to the extreme rigidity of the spatial organization of their embryonic cell lineages.

Polymorphic adaptations in metazoans to establish and maintain photosymbioses

Mutualistic symbioses are common throughout the animal kingdom. Rather unusual is a form of symbiosis, photosymbiosis, where animals are symbiotic with photoautotrophic organisms. Photosymbiosis is found among sponges, cnidarians, flatworms, molluscs, ascidians and even some amphibians. Generally the animal host harbours a phototrophic partner, usually a cyanobacteria or a unicellular alga. An exception to this rule is found in some sea slugs, which only retain the chloroplasts of the algal food source and maintain them photosynthetically active in their own cytosol - a phenomenon called 'functional kleptoplasty'. Research has focused largely on the biodiversity of photosymbiotic species across a range of taxa. However, many questions with regard to the evolution of the ability to establish and maintain a photosymbiosis are still unanswered. To date, attempts to understand genome adaptations which could potentially lead to the evolution of photosymbioses have only been performed in cnidarians. This knowledge gap for other systems is mainly due to a lack of genetic information, both for non-symbiotic and symbiotic species. Considering non-photosymbiotic species is, however, important to understand the factors that make symbiotic species so unique. Herein we provide an overview of the diversity of photosymbioses across the animal kingdom and discuss potential scenarios for the evolution of this association in different lineages. We stress that the evolution of photosymbiosis is probably based on genome adaptations, which (i) lead to recognition of the symbiont to establish the symbiosis, and (ii) are needed to maintain the symbiosis. We hope to stimulate research involving sequencing the genomes of various key taxa to increase the genomic resources needed to understand the most fundamental question: how have animals evolved the ability to establish and maintain a photosymbiosis?

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

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

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.

Dynamics of genomic innovation in the unicellular ancestry of animals

Which genomic innovations underpinned the origin of multicellular animals is still an open debate. Here, we investigate this question by reconstructing the genome architecture and gene family diversity of ancestral premetazoans, aiming to date the emergence of animal-like traits. Our comparative analysis involves genomes from animals and their closest unicellular relatives (the Holozoa), including four new genomes: three Ichthyosporea and Corallochytrium limacisporum. Here, we show that the earliest animals were shaped by dynamic changes in genome architecture before the emergence of multicellularity: an early burst of gene diversity in the ancestor of Holozoa, enriched in transcription factors and cell adhesion machinery, was followed by multiple and differently-timed episodes of synteny disruption, intron gain and genome expansions. Thus, the foundations of animal genome architecture were laid before the origin of complex multicellularity – highlighting the necessity of a unicellular perspective to understand early animal evolution.

DOI:http://dx.doi.org/10.7554/eLife.26036.001

Hundreds of millions of years ago, some single-celled organisms gained the ability to work together and form multicellular organisms. This transition was a major step in evolution and took place at separate times in several parts of the tree of life, including in animals, plants, fungi and algae.

Animals are some of the most complex organisms on Earth. Their single-celled ancestors were also quite genetically complex themselves and their genomes (the complete set of the organism’s DNA) already contained many genes that now coordinate the activity of the cells in a multicellular organism.

The genome of an animal typically has certain features: it is large, diverse and contains many segments (called introns) that are not genes. By seeing if the single-celled relatives of animals share these traits, it is possible to learn more about when specific genetic features first evolved, and whether they are linked to the origin of animals.

Now, Grau-Bové et al. have studied the genomes of several of the animal kingdom’s closest single-celled relatives using a technique called whole genome sequencing. This revealed that there was a period of rapid genetic change in the single-celled ancestors of animals during which their genes became much more diverse. Another ‘explosion’ of diversity happened after animals had evolved. Furthermore, the overall amount of the genomic content inside cells and the number of introns found in the genome rapidly increased in separate, independent events in both animals and their single-celled ancestors.

Future research is needed to investigate whether other multicellular life forms – such as plants, fungi and algae – originated in the same way as animal life. Understanding how the genetic material of animals evolved also helps us to understand the genetic structures that affect our health. For example, genes that coordinate the behavior of cells (and so are important for multicellular organisms) also play a role in cancer, where cells break free of this regulation to divide uncontrollably.

DOI:http://dx.doi.org/10.7554/eLife.26036.002

Structure and ultrastructure of eyes and brains of Thalia democratica (Thaliacea, Tunicata, Chordata)

Salps are marine planktonic chordates that possess an obligatory alternation of reproductive modes in subsequent generations. Within tunicates, salps represent a derived life cycle and are of interest in considerations of the evolutionary origin of complex anatomical structures and life history strategies. In the present study, the eyes and brains of both the sexual, aggregate blastozooid and the asexual, solitary oozooid stage of Thalia democratica (Forskål, ) were digitally reconstructed in detail based on serial sectioning for light and transmission electron microscopy. The blastozooid stage of T. democratica possesses three pigment cup eyes, situated in the anterior ventral part of the brain. The eyes are arranged in a way that the optical axes of each eye point toward different directions. Each eye is an inverse eye that consists of two different cell types: pigment cells (pigc) and rhabdomeric photoreceptor cells (prcs). The oozooid stage of T. democratica is equipped with a single horseshoe-shaped eye, positioned in the anterior dorsal part of the brain. The opening of the horseshoe-shaped eye points anteriorly. Similar to the eyes of the blastozooid, the eye of the oozooid consists of pigment cells and rhabdomeric photoreceptor cells. The rhabdomeric photoreceptor cells possess apical microvilli that form a densely packed presumably photosensitive receptor part adjacent to the concave side of the pigc. We suggest correspondences of the individual eyes in the blastozooid stage to respective parts of the single horseshoe-shaped eye in the oozooid stage and hypothesize that the differences in visual structures and brain anatomies evolved as a result of the aggregate life style of the blastozooid as opposed to the solitary life style of the oozooid.

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.

Studying the genetic basis of speciation in high gene flow marine invertebrates

A growing number of genes responsible for reproductive incompatibilities between species (barrier loci) exhibit the signals of positive selection. However, the possibility that genes experiencing positive selection diverge early in speciation and commonly cause reproductive incompatibilities has not been systematically investigated on a genome-wide scale. Here, I outline a research program for studying the genetic basis of speciation in broadcast spawning marine invertebrates that uses a priori genome-wide information on a large, unbiased sample of genes tested for positive selection. A targeted sequence capture approach is proposed that scores single-nucleotide polymorphisms (SNPs) in widely separated species populations at an early stage of allopatric divergence. The targeted capture of both coding and non-coding sequences enables SNPs to be characterized at known locations across the genome and at genes with known selective or neutral histories. The neutral coding and non-coding SNPs provide robust background distributions for identifying FST-outliers within genes that can, in principle, identify specific mutations experiencing diversifying selection. If natural hybridization occurs between species, the neutral coding and non-coding SNPs can provide a neutral admixture model for genomic clines analyses aimed at finding genes exhibiting strong blocks to introgression. Strongylocentrotid sea urchins are used as a model system to outline the approach but it can be used for any group that has a complete reference genome available.

Independent evolution of genomic characters during major metazoan transitions

Metazoan evolution encompasses a vast evolutionary time scale spanning over 600 million years. Our ability to infer ancestral metazoan characters, both morphological and functional, is limited by our understanding of the nature and evolutionary dynamics of the underlying regulatory networks. Increasing coverage of metazoan genomes enables us to identify the evolutionary changes of the relevant genomic characters such as the loss or gain of coding sequences, gene duplications, micro- and macro-synteny, and non-coding element evolution in different lineages. In this review we describe recent advances in our understanding of ancestral metazoan coding and non-coding features, as deduced from genomic comparisons. Some genomic changes such as innovations in gene and linkage content occur at different rates across metazoan clades, suggesting some level of independence among genomic characters. While their contribution to biological innovation remains largely unclear, we review recent literature about certain genomic changes that do correlate with changes to specific developmental pathways and metazoan innovations. In particular, we discuss the origins of the recently described pharyngeal cluster which is conserved across deuterostome genomes, and highlight different genomic features that have contributed to the evolution of this group. We also assess our current capacity to infer ancestral metazoan states from gene models and comparative genomics tools and elaborate on the future directions of metazoan comparative genomics relevant to evo-devo studies.

Rapid Evolutionary Rates and Unique Genomic Signatures Discovered in the First Reference Genome for the Southern Ocean Salp, Salpa thompsoni (Urochordata, Thaliacea)

A preliminary genome sequence has been assembled for the Southern Ocean salp, Salpa thompsoni (Urochordata, Thaliacea). Despite the ecological importance of this species in Antarctic pelagic food webs and its potential role as an indicator of changing Southern Ocean ecosystems in response to climate change, no genomic resources are available for S. thompsoni or any closely related urochordate species. Using a multiple-platform, multiple-individual approach, we have produced a 318,767,936-bp genome sequence, covering >50% of the estimated 602 Mb (±173 Mb) genome size for S. thompsoni. Using a nonredundant set of predicted proteins, >50% (16,823) of sequences showed significant homology to known proteins and ∼38% (12,151) of the total protein predictions were associated with Gene Ontology functional information. We have generated 109,958 SNP variant and 9,782 indel predictions for this species, serving as a resource for future phylogenomic and population genetic studies. Comparing the salp genome to available assemblies for four other urochordates, Botryllus schlosseri, Ciona intestinalis, Ciona savignyi and Oikopleura dioica, we found that S. thompsoni shares the previously estimated rapid rates of evolution for these species. High mutation rates are thus independent of genome size, suggesting that rates of evolution >1.5 times that observed for vertebrates are a broad taxonomic characteristic of urochordates. Tests for positive selection implemented in PAML revealed a small number of genes with sites undergoing rapid evolution, including genes involved in ribosome biogenesis and metabolic and immune process that may be reflective of both adaptation to polar, planktonic environments as well as the complex life history of the salps. Finally, we performed an initial survey of small RNAs, revealing the presence of known, conserved miRNAs, as well as novel miRNA genes; unique piRNAs; and mature miRNA signatures for varying developmental stages. Collectively, these resources provide a genomic foundation supporting S. thompsoni as a model species for further examination of the exceptional rates and patterns of genomic evolution shown by urochordates. Additionally, genomic data will allow for the development of molecular indicators of key life history events and processes and afford new understandings and predictions of impacts of climate change on this key species of Antarctic pelagic ecosystems.

Automated detection of ncRNAs in the draft genome sequence of a colonial tunicate: the carpet sea squirt Didemnum vexillum

BACKGROUND

The colonial ascidian Didemnum vexillum, sea carpet squirt, is not only a key marine organism to study morphological ancestral patterns of chordates evolution but it is also of great ecological importance due to its status as a major invasive species. Non-coding RNAs, in particular microRNAs (miRNAs), are important regulatory genes that impact development and environmental adaptation. Beyond miRNAs, not much in known about tunicate ncRNAs.

RESULTS

We provide here a comprehensive homology-based annotation of non-coding RNAs in the recently sequenced genome of D. vexillum. To this end we employed a combination of several computational approaches, including blast searches with a wide range of parameters, and secondary structured centered survey with infernal. The resulting candidate set was curated extensively to produce a high-quality ncRNA annotation of the first draft of the D. vexillum genome. It comprises 57 miRNA families, 4 families of ribosomal RNAs, 22 isoacceptor classes of tRNAs (of which more than 72 % of loci are pseudogenes), 13 snRNAs, 12 snoRNAs, and 1 other RNA family. Additionally, 21 families of mitochondrial tRNAs and 2 of mitochondrial ribosomal RNAs and 1 long non-coding RNA.

CONCLUSIONS

The comprehensive annotation of the D. vexillum non-coding RNAs provides a starting point towards a better understanding of the restructuring of the small RNA system in ascidians. Furthermore it provides a valuable research for efforts to establish detailed non-coding RNA annotations for other recently published and recently sequences in tunicate genomes.

Bilaterian-like promoters in the highly compact Amphimedon queenslandica genome

The regulatory systems underlying animal development must have evolved prior to the emergence of eumetazoans (cnidarians and bilaterians). Although representatives of earlier-branching animals - sponges ctenophores and placozoans - possess most of the developmental transcription factor families present in eumetazoans, the DNA regulatory elements that these transcription factors target remain uncharted. Here we characterise the core promoter sequences, U1 snRNP-binding sites (5' splice sites; 5'SSs) and polyadenylation sites (PASs) in the sponge Amphimedon queenslandica. Similar to unicellular opisthokonts, Amphimedon's genes are tightly packed in the genome and have small introns. In contrast, its genes possess metazoan-like core promoters populated with binding motifs previously deemed to be specific to vertebrates, including Nrf-1 and Krüppel-like elements. Also as in vertebrates, Amphimedon's PASs and 5'SSs are depleted downstream and upstream of transcription start sites, respectively, consistent with non-elongating transcripts being short-lived; PASs and 5'SSs are more evenly distributed in bidirectional promoters in Amphimedon. The presence of bilaterian-like regulatory DNAs in sponges is consistent with these being early and essential innovations of the metazoan gene regulatory repertoire.

Evolutionary loss of peroxisomes--not limited to parasites

Background

Peroxisomes are ubiquitous eukaryotic organelles that compartmentalize a variety of metabolic pathways that are primarily related to the oxidative metabolism of lipids and the detoxification of reactive oxygen species. The importance of peroxisomes is underscored by serious human diseases, which are caused by disorders in peroxisomal functions. Some eukaryotic lineages, however, lost peroxisomes. These organisms are mainly anaerobic protists and some parasitic lineages including Plasmodium and parasitic platyhelminths. Here we performed a systematic in-silico analysis of peroxisomal markers among metazoans to assess presence of peroxisomes and peroxisomal enzymes.

Results

Our analyses reveal an obvious loss of peroxisomes in all tested flukes, tapeworms, and parasitic roundworms of the order Trichocephalida. Intriguingly, peroxisomal markers are absent from the genome of the free-living tunicate Oikopleura dioica, which inhabits oxygen-containing niches of sea waters. We further map the presence and predicted subcellular localization of putative peroxisomal enzymes, showing that in organisms without the peroxisomal markers the set of these enzymes is highly reduced and none of them contains a predicted peroxisomal targeting signal.

Conclusions

We have shown that several lineages of metazoans independently lost peroxisomes and that the loss of peroxisomes was not exclusively associated with adaptation to anaerobic habitats and a parasitic lifestyle. Although the reason for the loss of peroxisomes from O. dioica is unclear, organisms lacking peroxisomes, including the free-living O. dioica, share certain typical r-selected traits: high fecundity, limited ontogenesis and relatively low complexity of the gene content. We hypothesize that peroxisomes are generally the first compartment to be lost during evolutionary reductions of the eukaryotic cell.

Reviewers

This article was reviewed by Michael Gray and Nick Lane.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0101-6) contains supplementary material, which is available to authorized users.

Transcriptional control analyses of the Xiphophorus melanoma oncogene

Melanoma development in interspecific hybrids of Xiphophorus is induced by the overexpression of the mutationally activated receptor tyrosine kinase Xmrk in pigment cells. Based on the melanocyte specificity of the transcriptional upregulation, a pigment cell-specific promoter region was postulated for xmrk, the activity of which is controlled in healthy purebred fish by the molecularly still unidentified regulator locus R. However, as yet the xmrk promoter region is still poorly characterized. In order to contribute to a better understanding of xmrk expression regulation, we performed a functional analysis of the entire putative gene regulatory region of the oncogene using conventional plasmid-based reporter systems as well as a newly established method employing BAC-derived luciferase reporter constructs in melanoma and non-melanoma cell lines. Using the melanocyte-specific mitfa promoter as control, we could demonstrate that our in vitro system is able to reliably monitor regulation of transcription through cell type-specific regulatory sequences. We found that sequences within 200kb flanking the xmrk oncogene do not lead to any specific transcriptional activation in melanoma compared to control cells. Hence, xmrk reporter constructs fail to faithfully reproduce the endogenous transcriptional regulation of the oncogene. Our data therefore strongly indicate that the melanocyte-specific transcription of xmrk is not the consequence of pigment cell-specific cis-regulatory elements in the promoter region. This hints at additional regulatory mechanisms involved in transcriptional control of the oncogene, thereby suggesting a key role for epigenetic mechanisms in oncogenic xmrk overexpression and thereby in tumor development in Xiphophorus.

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.

The transcriptome of an amphioxus, Asymmetron lucayanum, from the Bahamas: a window into chordate evolution

Cephalochordates, the sister group of tunicates plus vertebrates, have been called “living fossils” due to their resemblance to fossil chordates from Cambrian strata. The genome of the cephalochordate Branchiostoma floridae shares remarkable synteny with vertebrates and is free from whole-genome duplication. We performed RNA sequencing from larvae and adults of Asymmetron lucayanum, a cephalochordate distantly related to B. floridae. Comparisons of about 430 orthologous gene groups among both cephalochordates and 10 vertebrates using an echinoderm, a hemichordate, and a mollusk as outgroups showed that cephalochordates are evolving more slowly than the slowest evolving vertebrate known (the elephant shark), with A. lucayanum evolving even more slowly than B. floridae. Against this background of slow evolution, some genes, notably several involved in innate immunity, stand out as evolving relatively quickly. This may be due to the lack of an adaptive immune system and the relatively high levels of bacteria in the inshore waters cephalochordates inhabit. Molecular dating analysis including several time constraints revealed a divergence time of ∼120 Ma for A. lucayanum and B. floridae. The divisions between cephalochordates and vertebrates, and that between chordates and the hemichordate plus echinoderm clade likely occurred before the Cambrian.

The Global Invertebrate Genomics Alliance (GIGA): developing community resources to study diverse invertebrate genomes

Over 95% of all metazoan (animal) species comprise the "invertebrates," but very few genomes from these organisms have been sequenced. We have, therefore, formed a "Global Invertebrate Genomics Alliance" (GIGA). Our intent is to build a collaborative network of diverse scientists to tackle major challenges (e.g., species selection, sample collection and storage, sequence assembly, annotation, analytical tools) associated with genome/transcriptome sequencing across a large taxonomic spectrum. We aim to promote standards that will facilitate comparative approaches to invertebrate genomics and collaborations across the international scientific community. Candidate study taxa include species from Porifera, Ctenophora, Cnidaria, Placozoa, Mollusca, Arthropoda, Echinodermata, Annelida, Bryozoa, and Platyhelminthes, among others. GIGA will target 7000 noninsect/nonnematode species, with an emphasis on marine taxa because of the unrivaled phyletic diversity in the oceans. Priorities for selecting invertebrates for sequencing will include, but are not restricted to, their phylogenetic placement; relevance to organismal, ecological, and conservation research; and their importance to fisheries and human health. We highlight benefits of sequencing both whole genomes (DNA) and transcriptomes and also suggest policies for genomic-level data access and sharing based on transparency and inclusiveness. The GIGA Web site (http://giga.nova.edu) has been launched to facilitate this collaborative venture.

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.

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.

Genetic and genomic toolbox of the chordate Ciona intestinalis

The experimental malleability and unique phylogenetic position of the sea squirt Ciona intestinalis as part of the sister group to the vertebrates have helped establish these marine chordates as model organisms for the study of developmental genetics and evolution. Here we summarize the tools, techniques, and resources available to the Ciona geneticist, citing examples of studies that employed such strategies in the elucidation of gene function in Ciona. Genetic screens, germline transgenesis, electroporation of plasmid DNA, and microinjection of morpholinos are all routinely employed, and in the near future we expect these to be complemented by targeted mutagenesis, homologous recombination, and RNAi. The genomic resources available will continue to support the design and interpretation of genetic experiments and allow for increasingly sophisticated approaches on a high-throughput, whole-genome scale.

The Ruegeria pomeroyi acuI gene has a role in DMSP catabolism and resembles yhdH of E. coli and other bacteria in conferring resistance to acrylate

The Escherichia coli YhdH polypeptide is in the MDR012 sub-group of medium chain reductase/dehydrogenases, but its biological function was unknown and no phenotypes of YhdH⁻ mutants had been described. We found that an E. coli strain with an insertional mutation in yhdH was hyper-sensitive to inhibitory effects of acrylate, and, to a lesser extent, to those of 3-hydroxypropionate. Close homologues of YhdH occur in many Bacterial taxa and at least two animals. The acrylate sensitivity of YhdH⁻ mutants was corrected by the corresponding, cloned homologues from several bacteria. One such homologue is acuI, which has a role in acrylate degradation in marine bacteria that catabolise dimethylsulfoniopropionate (DMSP) an abundant anti-stress compound made by marine phytoplankton. The acuI genes of such bacteria are often linked to ddd genes that encode enzymes that cleave DMSP into acrylate plus dimethyl sulfide (DMS), even though these are in different polypeptide families, in unrelated bacteria. Furthermore, most strains of Roseobacters, a clade of abundant marine bacteria, cleave DMSP into acrylate plus DMS, and can also demethylate it, using DMSP demethylase. In most Roseobacters, the corresponding gene, dmdA, lies immediately upstream of acuI and in the model Roseobacter strain Ruegeria pomeroyi DSS-3, dmdA-acuI were co-regulated in response to the co-inducer, acrylate. These observations, together with findings by others that AcuI has acryloyl-CoA reductase activity, lead us to suggest that YdhH/AcuI enzymes protect cells against damaging effects of intracellular acryloyl-CoA, formed endogenously, and/or via catabolising exogenous acrylate. To provide “added protection” for bacteria that form acrylate from DMSP, acuI was recruited into clusters of genes involved in this conversion and, in the case of acuI and dmdA in the Roseobacters, their co-expression may underpin an interaction between the two routes of DMSP catabolism, whereby the acrylate product of DMSP lyases is a co-inducer for the demethylation pathway.

Origin and evolution of spliceosomal introns

Evolution of exon-intron structure of eukaryotic genes has been a matter of long-standing, intensive debate. The introns-early concept, later rebranded ‘introns first’ held that protein-coding genes were interrupted by numerous introns even at the earliest stages of life's evolution and that introns played a major role in the origin of proteins by facilitating recombination of sequences coding for small protein/peptide modules. The introns-late concept held that introns emerged only in eukaryotes and new introns have been accumulating continuously throughout eukaryotic evolution. Analysis of orthologous genes from completely sequenced eukaryotic genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists, suggesting that many ancestral introns have persisted since the last eukaryotic common ancestor (LECA). Reconstructions of intron gain and loss using the growing collection of genomes of diverse eukaryotes and increasingly advanced probabilistic models convincingly show that the LECA and the ancestors of each eukaryotic supergroup had intron-rich genes, with intron densities comparable to those in the most intron-rich modern genomes such as those of vertebrates. The subsequent evolution in most lineages of eukaryotes involved primarily loss of introns, with only a few episodes of substantial intron gain that might have accompanied major evolutionary innovations such as the origin of metazoa. The original invasion of self-splicing Group II introns, presumably originating from the mitochondrial endosymbiont, into the genome of the emerging eukaryote might have been a key factor of eukaryogenesis that in particular triggered the origin of endomembranes and the nucleus. Conversely, splicing errors gave rise to alternative splicing, a major contribution to the biological complexity of multicellular eukaryotes. There is no indication that any prokaryote has ever possessed a spliceosome or introns in protein-coding genes, other than relatively rare mobile self-splicing introns. Thus, the introns-first scenario is not supported by any evidence but exon-intron structure of protein-coding genes appears to have evolved concomitantly with the eukaryotic cell, and introns were a major factor of evolution throughout the history of eukaryotes. This article was reviewed by I. King Jordan, Manuel Irimia (nominated by Anthony Poole), Tobias Mourier (nominated by Anthony Poole), and Fyodor Kondrashov. For the complete reports, see the Reviewers’ Reports section.