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

Transcriptome sequencing of seven deep marine invertebrates

We present 4k video and whole transcriptome data for seven deep-sea invertebrate animals collected in the Eastern Pacific Ocean during a research expedition onboard the Schmidt Ocean Institute's R/V Falkor in August of 2021. The animals include one jellyfish (Atolla sp.), three siphonophores (Apolemia sp., Praya sp., and Halistemma sp.), one larvacean (Bathochordaeus mcnutti), one tunicate (Pyrosomatidae sp.), and one ctenophore (Lampocteis sp.). Four of the animals were sequenced with long-read RNA sequencing technology, such that the reads themselves define a reference assembly for those animals. The larvacean tissues were successfully preserved in situ and has paired long-read reference data and short read quantitative transcriptomic data for within-specimen analyses of gene expression. Additionally, for three animals we provide quantitative image data, and a 3D model for one siphonophore. The paired image and transcriptomic data can be used for species identification, species description, and reference genetic data for these deep-sea animals.

An in situ digital synthesis strategy for the discovery and description of ocean life

Revolutionary advancements in underwater imaging, robotics, and genomic sequencing have reshaped marine exploration. We present and demonstrate an interdisciplinary approach that uses emerging quantitative imaging technologies, an innovative robotic encapsulation system with in situ RNA preservation and next-generation genomic sequencing to gain comprehensive biological, biophysical, and genomic data from deep-sea organisms. The synthesis of these data provides rich morphological and genetic information for species description, surpassing traditional passive observation methods and preserved specimens, particularly for gelatinous zooplankton. Our approach enhances our ability to study delicate mid-water animals, improving research in the world's oceans.

Gene expression patterns of Salpa thompsoni reveal remarkable differences in metabolism and reproduction near the Antarctic Polar Front

Salpa thompsoni is an important grazer in the Southern Ocean and most abundant in the Antarctic Polar Front (APF) region. During recent decades, their distribution expanded southwards. However, it is unclear whether salps can maintain their populations in the high Antarctic regions throughout the year owing to a poor understanding of their physiological responses to changing environmental conditions. We examined gene expression signatures of salps collected in two geographically close regions south of the APF that differed in water mass composition and productivity. The observed differences in the expression of genes related to reproductive, cellular and metabolic processes reflect variations in water temperature and food supply between the two regions studied here. Our study contributes to a better understanding of the physiological responses of S. thompsoni to changing environmental conditions, and how the species may adapt to a changing environment through potential geographical population shifts under future climate change scenarios.

Salpa genome and developmental transcriptome analyses reveal molecular flexibility enabling reproductive success in a rapidly changing environment

Ocean warming favors pelagic tunicates, such as salps, that exhibit increasingly frequent and rapid population blooms, impacting trophic dynamics and composition and human marine-dependent activities. Salp blooms are a result of their successful reproductive life history, alternating seasonally between asexual and sexual protogynous (i.e. sequential) hermaphroditic stages. While predicting future salp bloom frequency and intensity relies on an understanding of the transitions during the sexual stage from female through parturition and subsequent sex change to male, these transitions have not been explored at the molecular level. Here we report the development of the first complete genome of S. thompsoni and the North Atlantic sister species S. aspera. Genome and comparative analyses reveal an abundance of repeats and G-quadruplex (G4) motifs, a highly stable secondary structure, distributed throughout both salp genomes, a feature shared with other tunicates that perform alternating sexual-asexual reproductive strategies. Transcriptional analyses across sexual reproductive stages for S. thompsoni revealed genes associated with male sex differentiation and spermatogenesis are expressed as early as birth and before parturition, inconsistent with previous descriptions of sequential sexual differentiation in salps. Our findings suggest salp are poised for reproductive success at birth, increasing the potential for bloom formation as ocean temperatures rise.

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.

WOS, Steinberg DK, Swadling KM, Tarling GA, Thorpe SE, Veytia D, Ward P, Weldrick CK, Yang G. Status, Change, and Futures of Zooplankton in the Southern Ocean

In the Southern Ocean, several zooplankton taxonomic groups, euphausiids, copepods, salps and pteropods, are notable because of their biomass and abundance and their roles in maintaining food webs and ecosystem structure and function, including the provision of globally important ecosystem services. These groups are consumers of microbes, primary and secondary producers, and are prey for fishes, cephalopods, seabirds, and marine mammals. In providing the link between microbes, primary production, and higher trophic levels these taxa influence energy flows, biological production and biomass, biogeochemical cycles, carbon flux and food web interactions thereby modulating the structure and functioning of ecosystems. Additionally, Antarctic krill (Euphausia superba) and various fish species are harvested by international fisheries. Global and local drivers of change are expected to affect the dynamics of key zooplankton species, which may have potentially profound and wide-ranging implications for Southern Ocean ecosystems and the services they provide. Here we assess the current understanding of the dominant metazoan zooplankton within the Southern Ocean, including Antarctic krill and other key euphausiid, copepod, salp and pteropod species. We provide a systematic overview of observed and potential future responses of these taxa to a changing Southern Ocean and the functional relationships by which drivers may impact them. To support future ecosystem assessments and conservation and management strategies, we also identify priorities for Southern Ocean zooplankton research.

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.

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.

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.

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.

nrDNA:mtDNA copy number ratios as a comparative metric for evolutionary and conservation genetics

Identifying genetic cues of functional relevance is key to understanding the drivers of evolution and increasingly important for the conservation of biodiversity. This study introduces nuclear ribosomal DNA (nrDNA) to mitochondrial DNA (mtDNA) copy number ratios as a metric with which to screen for this functional genetic variation prior to more extensive omics analyses. To illustrate the metric, quantitative PCR was used to estimate nrDNA (18S) to mtDNA (16S) copy number ratios in muscle tissue from samples of two zooplankton species: Salpa thompsoni caught near Elephant Island (Southern Ocean) and S. fusiformis sampled off Gough Island (South Atlantic). Average 18S:16S ratios in these samples were 9:1 and 3:1, respectively. nrDNA 45S arrays and mitochondrial genomes were then deep sequenced to uncover the sources of intra-individual genetic variation underlying these 18S:16S copy number differences. The deep sequencing profiles obtained were consistent with genetic changes resulting from adaptive processes, including an expansion of nrDNA and damage to mtDNA in S. thompsoni, potentially in response to the polar environment. Beyond this example from zooplankton, nrDNA:mtDNA copy number ratios offer a promising metric to help identify genetic variation of functional relevance in animals more broadly.

A phylogenomic framework and timescale for comparative studies of tunicates

Background

Tunicates are the closest relatives of vertebrates and are widely used as models to study the evolutionary developmental biology of chordates. Their phylogeny, however, remains poorly understood, and to date, only the 18S rRNA nuclear gene and mitogenomes have been used to delineate the major groups of tunicates. To resolve their evolutionary relationships and provide a first estimate of their divergence times, we used a transcriptomic approach to build a phylogenomic dataset including all major tunicate lineages, consisting of 258 evolutionarily conserved orthologous genes from representative species.

Results

Phylogenetic analyses using site-heterogeneous CAT mixture models of amino acid sequence evolution resulted in a strongly supported tree topology resolving the relationships among four major tunicate clades: (1) Appendicularia, (2) Thaliacea + Phlebobranchia + Aplousobranchia, (3) Molgulidae, and (4) Styelidae + Pyuridae. Notably, the morphologically derived Thaliacea are confirmed as the sister group of the clade uniting Phlebobranchia + Aplousobranchia within which the precise position of the model ascidian genus Ciona remains uncertain. Relaxed molecular clock analyses accommodating the accelerated evolutionary rate of tunicates reveal ancient diversification (~ 450–350 million years ago) among the major groups and allow one to compare their evolutionary age with respect to the major vertebrate model lineages.

Conclusions

Our study represents the most comprehensive phylogenomic dataset for the main tunicate lineages. It offers a reference phylogenetic framework and first tentative timescale for tunicates, allowing a direct comparison with vertebrate model species in comparative genomics and evolutionary developmental biology studies.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0499-2) contains supplementary material, which is available to authorized users.

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.

ANISEED 2017: extending the integrated ascidian database to the exploration and evolutionary comparison of genome-scale datasets

ANISEED (www.aniseed.cnrs.fr) is the main model organism database for tunicates, the sister-group of vertebrates. This release gives access to annotated genomes, gene expression patterns, and anatomical descriptions for nine ascidian species. It provides increased integration with external molecular and taxonomy databases, better support for epigenomics datasets, in particular RNA-seq, ChIP-seq and SELEX-seq, and features novel interactive interfaces for existing and novel datatypes. In particular, the cross-species navigation and comparison is enhanced through a novel taxonomy section describing each represented species and through the implementation of interactive phylogenetic gene trees for 60% of tunicate genes. The gene expression section displays the results of RNA-seq experiments for the three major model species of solitary ascidians. Gene expression is controlled by the binding of transcription factors to cis-regulatory sequences. A high-resolution description of the DNA-binding specificity for 131 Ciona robusta (formerly C. intestinalis type A) transcription factors by SELEX-seq is provided and used to map candidate binding sites across the Ciona robusta and Phallusia mammillata genomes. Finally, use of a WashU Epigenome browser enhances genome navigation, while a Genomicus server was set up to explore microsynteny relationships within tunicates and with vertebrates, Amphioxus, echinoderms and hemichordates.

Nonprotein-Coding RNAs as Regulators of Development in Tunicates

Tunicates, or urochordates, are a group of small marine organisms that are found widely throughout the seas of the world. As most plausible sister group of the vertebrates, they are of utmost importance for a comprehensive understanding of chordate evolution; hence, they have served as model organisms for many aspects of the developmental biology. Current genomic analysis of tunicates indicates that their genomes evolved with a fast rate not only at the level of nucleotide substitutions but also in terms of genomic organization. The latter involves genome reduction, rearrangements, as well as the loss of some important coding and noncoding RNA (ncRNAs) elements and even entire genomic regions that are otherwise well conserved. These observations are largely based on evidence from comparative genomics resulting from the analysis of well-studied gene families such as the Hox genes and their noncoding elements. In this chapter, the focus lies on the ncRNA complement of tunicates, with a particular emphasis on microRNAs, which have already been studied extensively for other animal clades. MicroRNAs are known as important regulators of key genes in animal development, and they are intimately related to the increase morphological complexity in higher metazoans. Here we review the discovery, evolution, and genome organization of the miRNA repertoire, which has been drastically reduced and restructured in tunicates compared to the chordate ancestor. Known functions of microRNAs as regulators of development in tunicates are a central topic. For instance, we consider the role of miRNAs as regulators of the muscle development and their importance in the regulation of the differential expression during the oral siphon regeneration. Beyond microRNAs, we touch upon the functions of some other ncRNAs such as yellow crescent RNA, moRNAs, RMST lncRNAs, or spliced-leader (SL) RNAs, which have diverse functions associated with the embryonic development, neurogenesis, and mediation of mRNA stability in general.

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.