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

Identification of astrocyte-like cells in an adult ascidian during regeneration of the central nervous system

The mechanisms underlying regeneration of the central nervous system (CNS) following lesions have been studied extensively in both vertebrate and invertebrate models. To shed light on regeneration, ascidians, a sister group of vertebrates and with remarkable ability to regenerate their brains, constitute an appropriate model system. Glial cells have been implicated in regeneration in vertebrates; however, their role in the adult ascidian CNS regeneration is unknown. A model of degeneration and regeneration using the neurotoxin 3-acetylpyridine (3AP) in the brain of the ascidian Styela plicata was used to identify astrocyte-like cells and investigate their role. We studied the CNS of control ascidians (injected with artificial sea water) and of ascidians whose CNS was regenerating (1 and 10 days after the injection with 3AP). Our results show that the mRNA of the ortholog of glutamine synthetase (GS), a glial-cell marker in vertebrates, is increased during the early stages of regeneration. Confirming the identity of GS, the protein was identified via immunostaining in a cell population during the same regeneration stage. Last, a single ortholog of GS (GSII) is present in ascidian and amphioxus genomes, while two types exist in fungi, some invertebrates, and vertebrates, suggesting that ascidians have lost the GSI type. Taken together, our findings revealed that a cell population expressing glial-cell markers may play a role in regeneration in adult ascidians. This is the first report of astrocyte-like cells in the adult ascidian CNS, and contributes to understanding of the evolution of glial cells among metazoans.

Navigating spatio-temporal microbiome dynamics: Environmental factors and trace elements shape the symbiont community of an invasive marine species

The proliferation of marine invasive species is a mounting concern. While the role of microbial communities in invasive ascidian species is recognized, the role of seasonal shifts in microbiome composition remains largely unexplored. We sampled five individuals of the invasive ascidian Styela plicata quarterly from January 2020 to October 2021 in two harbours, examining gills, tunics, and surrounding water. By analysing Amplicon Sequence Variants (ASVs) and seawater trace elements, we found that compartment (seawater, tunic, or gills) was the primary differentiating factor, followed by harbour. Clear seasonal patterns were evident in seawater bacteria, less so in gills, and absent in tunics. We identified compartment-specific bacteria, as well as seasonal indicator ASVs and ASVs correlated with trace element concentrations. Among these bacteria, we found that Endozoicomonas, Hepatoplasma and Rhodobacteraceae species had reported functions which might be necessary for overcoming seasonality and trace element shifts. This study contributes to understanding microbiome dynamics in invasive holobiont systems, and the patterns found indicate a potential role in adaptation and invasiveness.

Spatially resolved single-cell atlas of ascidian endostyle provides insight into the origin of vertebrate pharyngeal organs

The pharyngeal endoderm, an innovation of deuterostome ancestors, contributes to pharyngeal development by influencing the patterning and differentiation of pharyngeal structures in vertebrates; however, the evolutionary origin of the pharyngeal organs in vertebrates is largely unknown. The endostyle, a distinct pharyngeal organ exclusively present in basal chordates, represents a good model for understanding pharyngeal organ origins. Using Stereo-seq and single-cell RNA sequencing, we constructed aspatially resolved single-cell atlas for the endostyle of the ascidian Styela clava. We determined the cell composition of the hemolymphoid region, which illuminates a mixed ancestral structure for the blood and lymphoid system. In addition, we discovered a cluster of hair cell-like cells in zone 3, which has transcriptomic similarity with the hair cells of the vertebrate acoustico-lateralis system. These findings reshape our understanding of the pharynx of the basal chordate and provide insights into the evolutionary origin of multiplexed pharyngeal organs.

Temporospatial hierarchy and allele-specific expression of zygotic genome activation revealed by distant interspecific urochordate hybrids

Zygotic genome activation (ZGA) is a universal process in early embryogenesis of metazoan, when the quiescent zygotic nucleus initiates global transcription. However, the mechanisms related to massive genome activation and allele-specific expression (ASE) remain not well understood. Here, we develop hybrids from two deeply diverged (120 Mya) ascidian species to symmetrically document the dynamics of ZGA. We identify two coordinated ZGA waves represent early developmental and housekeeping gene reactivation, respectively. Single-cell RNA sequencing reveals that the major expression wave exhibits spatial heterogeneity and significantly correlates with cell fate. Moreover, allele-specific expression occurs in a species- rather than parent-related manner, demonstrating the divergence of cis-regulatory elements between the two species. These findings provide insights into ZGA in chordates.

Genomic Resources and Annotations for a Colonial Ascidian, the Light-Bulb Sea Squirt Clavelina lepadiformis

Ascidian embryos have been studied since the birth of experimental embryology at the end of the 19th century. They represent textbook examples of mosaic development characterized by a fast development with very few cells and invariant cleavage patterns and lineages. Ascidians belong to tunicates, the vertebrate sister group, and their study is essential to shed light on the emergence of vertebrates. Importantly, deciphering developmental gene regulatory networks has been carried out mostly in two of the three ascidian orders, Phlebobranchia and Stolidobranchia. To infer ancestral developmental programs in ascidians, it is thus essential to carry out molecular embryology in the third ascidian order, the Aplousobranchia. Here, we present genomic resources for the colonial aplousobranch Clavelina lepadiformis: a transcriptome produced from various embryonic stages, and an annotated genome. The assembly consists of 184 contigs making a total of 233.6 Mb with a N50 of 8.5 Mb and a L50 of 11. The 32,318 predicted genes capture 96.3% of BUSCO orthologs. We further show that these resources are suitable to study developmental gene expression and regulation in a comparative framework within ascidians. Additionally, they will prove valuable for evolutionary and ecological studies.

Establishment of a developmental atlas and transgenetic tools in the ascidian Styela clava

The ascidian Styela clava is an ecologically important species that is distributed along coastal regions worldwide. It has a long history as a model animal for evolutionary and developmental biology research owing to its phylogenetic position between vertebrates and invertebrates, and its classical mosaic expression patterns. However, the standard developmental atlas and protocols and tools for molecular manipulation of this organism are inadequate. In this study, we established a standard developmental table and provided a web-based digital image resource for S. clava embryogenesis at each developmental stage from fertilized eggs to hatching larvae by utilizing confocal laser microscopy and 3D reconstruction images. It takes around 10 h for fertilized eggs to develop into swimming larvae and 20-30 min to complete the tail regression processes at the metamorphic stage. We observed that the notochord cells in S. clava embryos did not produce an extracellular lumen like Ciona robusta, but showed polarized elongation behaviors, providing us an ideal comparative model to study tissue morphogenesis. In addition, we established a chemical-washing procedure to remove the chorion easily from the fertilized eggs. Based on the dechorionation technique, we further realized transgenic manipulation by electroporation and successfully applied tissue-specific fluorescent labeling in S. clava embryos. Our work provides a standard imaging atlas and powerful genetic tools for investigating embryogenesis and evolution using S. clava as a model organism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00200-2.

Ciona spp. and ascidians as bioindicator organisms for evaluating effects of endocrine disrupting chemicals: A discussion paper

In context of testing, screening and monitoring of endocrine-disrupting (ED) type of environmental pollutants, tunicates could possibly represent a particularly interesting group of bioindicator organisms. These primitive chordates are already important model organisms within developmental and genomics research due to their central position in evolution and close relationship to vertebrates. The solitary ascidians, such as the genus Ciona spp. (vase tunicates), could possibly be extra feasible as ED bioindicators. They have a free-swimming, tadpole-like larval stage that develops extremely quickly (<20 h under favorable conditions), has a short life cycle (typically 2-3 months), are relatively easy to maintain in laboratory culture, have fully sequenced genomes, and transgenic embryos with 3D course data of the embryo ontogeny are available. In this article, we discuss possible roles of Ciona spp. (and other solitary ascidians) as ecotoxicological bioindicator organisms in general but perhaps especially for effect studies of contaminants with presumed endocrine disrupting modes of action.

Ancestral glycoprotein hormone and its cognate receptor present in primitive chordate ascidian: Molecular identification and functional characterization

The glycoprotein hormone (GPH) system is fundamentally significant in regulating the physiology of chordates, such as thyroid activity and gonadal function. However, the knowledge of the GPH system in the primitive chordate ascidian species is largely lacking. Here, we reported an ancestral GPH system in the ascidian (Styela clava), which consists of GPH α subunit (Sc-GPA2), GPH β subunit (Sc-GPB5), and the cognate leucine-rich repeat-containing G protein-coupled receptor (Sc-GPHR). Comparative structure analysis revealed that distinct from vertebrate GPH β subunits, Sc-GPB5 was less conserved, showing an atypical N-terminal sequence with a type II transmembrane domain instead of a typical signal peptide. By investigating the presence of recombinant Sc-GPA2 and Sc-GPB5 in cell lysates and culture media of HEK293T cells, we confirmed that these two subunits could be secreted out of the cells via distinct secretory pathways. The deglycosylation experiments demonstrated that N-linked glycosylation only occurred on the conserved cysteine residue (N78) of Sc-GPA2, whereas Sc-GPB5 was non-glycosylated. Although Sc-GPB5 exhibited distinct topology and biochemical properties in contrast to its chordate counterparts, it could still interact with Sc-GPA2 to form a heterodimer. The Sc-GPHR was then confirmed to be activated by tethered Sc-GPA2/GPB5 heterodimer on the Gs-cAMP pathway, suggesting that Sc-GPA2/GPB5 heterodimer-initiated Gs-cAMP signaling pathway is evolutionarily conserved in chordates. Furthermore, in situ hybridization and RT-PCR results revealed the co-expression patterns of Sc-GPA2 and Sc-GPB5 with Sc-GPHR transcripts, respectively in ascidian larvae and adults, highlighting the potential functions of Sc-GPA2/GPB5 heterodimer as an autocrine/paracrine neurohormone in regulating metamorphosis of larvae and physiological functions of adults. Our study systematically investigated the GPA2/GPB5-GPHR system in ascidian for the first time, which offers insights into understanding the function and evolution of the GPH system within the chordate lineage.

Transcriptional Analysis of the Endostyle Reveals Pharyngeal Organ Functions in Ascidian

The endostyle is a pharyngeal organ with an opening groove and cilia in invertebrate chordates (amphioxus and ascidian) and cyclostomate (lamprey), serving as a filter-feeding tract and thyroid-secreting location. Emerging evidence implies its complex cellular composition and potentially versatile functions. Multiple cell types in the endostyle have been thought to be progenitors of complex organs in advanced vertebrates. To describe the expression profile and the potential functions, bulk RNA sequencing on the endostyle in ascidian Styela clava was conducted and distinct markers were selected by multileveled comparative analysis. Transcriptional data assay and qRT-PCR-verified results showed the regional expression patterns of Hox genes in the longitudinal axis. Organ-specific markers of the endostyle was proposed by comparing expression with the main organs of the ascidian. A cross-species transcriptional profile projection between the endostyle and organs from Danio rerio and Homo sapiens indicates a robust homogenous relationship to the thyroid and digestive system of the endostyle. The high similarity between the endostyle and the head kidney in zebrafish/the bone marrow in human implies uniquely profound functions of the pharyngeal organ in proto-vertebrates. Our result revealed that the transcriptional profile of the human parathyroid gland was similar to the ascidian endostyle, indicating the evolutionary origin of vertebrate hormone secretion organs.

Thyroid-like hormone signaling in invertebrates and its potential role in initial screening of thyroid hormone system disrupting chemicals

This review examines the presence and evolution of thyroid-like systems in selected aquatic invertebrates to determine the potential use of these organisms in screens for vertebrate thyroid hormone axis disrupting chemicals (THADCs). Such a screen might support the phasing out of some vertebrate testing. Although arthropods including crustaceans do not contain a functional thyroid signaling system, elements of such a system exist in the aquatic phyla mollusks, echinoderms, tunicates, and cephalochordates. These phyla can synthesize thyroid hormone, which has been demonstrated in some groups to induce the nuclear thyroid hormone receptor (THR). Thyroid hormone may act in these phyla through interaction with a membrane integrin receptor. Thyroid hormone regulates inter alia metamorphosis but, unlike in vertebrates, this does not occur via receptor activation by the ligands triiodothyronine (T3) and thyroxine (T4). Instead, the unliganded nuclear receptor itself controls metamorphosis in mollusks, echinoderms, and tunicates, whereas the T3 derivative tri-iodothyroacetic acid (TRIAC) acts as a THR ligand in cephalochordates. In view of this, it may be possible to develop an invertebrate-based screen that is sensitive to vertebrate THADCs that interfere with thyroid hormone synthesis or metabolism along with interaction with membrane receptors. The review makes some recommendations for the need to develop an appropriate test method. Integr Environ Assess Manag 2023;19:63-82. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Molecular Characterization and Expression Analysis of Putative Class C (Glutamate Family) G Protein-Coupled Receptors in Ascidian Styela clava

In this study, we performed the genome-wide domain analysis and sequence alignment on the genome of Styela clava, and obtained a repertoire of 204 putative GPCRs, which exhibited a highly reduced gene number compared to vertebrates and cephalochordates. In this repertoire, six Class C GPCRs, including four metabotropic glutamate receptors (Sc-GRMs), one calcium-sensing receptor (Sc-CaSR), and one gamma-aminobutyric acid (GABA) type B receptor 2-like (Sc-GABABR2-like) were identified, with the absence of type 1 taste and vomeronasal receptors. All the Sc-GRMs and Sc-CaSR contained the typical "Venus flytrap" and cysteine-rich domains required for ligand binding and subsequent propagation of conformational changes. In swimming larvae, Sc-grm3 and Sc-casr were mainly expressed at the junction of the sensory vesicle and tail nerve cord while the transcripts of Sc-grm4, Sc-grm7a, and Sc-grm7b appeared at the anterior trunk, which suggested their important functions in neurotransmission. The high expression of these Class C receptors at tail-regression and metamorphic juvenile stages hinted at their potential involvement in regulating metamorphosis. In adults, the transcripts were highly expressed in several peripheral tissues, raising the possibility that S. clava Class C GPCRs might function as neurotransmission modulators peripherally after metamorphosis. Our study systematically characterized the ancestral chordate Class C GPCRs to provide insights into the origin and evolution of these receptors in chordates and their roles in regulating physiological and morphogenetic changes relevant to the development and environmental adaption.

Cultivation of gut microorganisms of the marine ascidian Halocynthia roretzi reveals their potential roles in the environmental adaptation of their host

It has long been known that abundant symbiotic bacteria exist in the tunic and gut of marine ascidians, and that these play crucial roles in host development, physiological metabolism, and environmental adaptation. However, the identity, roles and functions of these symbiotic bacteria are known for only a few strains. In this study, we isolated and cultivated 263 strains of microorganisms from the intestine of the marine ascidian Halocynthia roretzi through a combination of aerobic and anaerobic culture approaches. Most cultivated species, both aerobic and anaerobic, from ascidian stool samples belonged to the genus Bacillus based on 16S rDNA sequencing identification and phylogenetic assays. The distribution of cultured bacteria varied with seasonal changes in environmental conditions. To explore the functions of cultured bacteria, we screened out a strain of Serratia sp. whose extracts showed high antibacterial activity against aquatic pathogens. These findings revealed the potential roles of gut microorganisms in ascidian defense and environmental adaptation, thus providing insights into the interaction and co-evolution between gut bacteria and their hosts.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00131-4.

The Genome of the "Sea Vomit" Didemnum vexillum

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

The HSP/co-chaperone network in environmental cold adaptation of Chilo suppressalis

Winter cold is one of the major environmental stresses for ectotherm species. Chilo suppressalis, a notorious lepidopteran pest of rice, has a wide geographic region that includes temperate zones with severe environmental conditions. Although C. suppressalis exhibits remarkable cold tolerance, its cold-adaptation mechanisms remain unclear. Here, we used bioinformatics approaches to evaluate transcript levels of genes comprising the C. suppressalis heat shock protein (HSP)/co-chaperone network in response to cold-induced stress. Using all such genes identified in the C. suppressalis genome, we experimentally examined the corresponding transcript levels under cold-acclimation or intermittent cold-shock stresses in diapause and non-diapausing larvae. In total, we identified 19 HSPs and 8 HSP co-chaperones in the C. suppressalis genome. Nine (hsp90, hsp75, hsp70, hsp40, small hsp, activator of 90 kDa heat shock protein ATPase-like, heat shock factor, heat shock factor binding protein 1-like and HSPB1-associated protein 1) were highly cold-inducible and likely comprise the principal cold-response HSP/co-chaperone network in C. suppressalis. We also found that transcriptional regulation of the HSP/co-chaperone networks response differs between cold-acclimation and short-term cold-shock. Moreover, activation of the HSP/co-chaperone network depends on the diapause state of overwintering larvae and cold acclimation may further increase larval cold tolerance. These results provide key new insights in the cold-adaptation mechanisms in C. suppressalis.

Genome-Wide Identification, Comparison, and Expression Analysis of Transcription Factors in Ascidian Styela clava

Tunicates include diverse species, as they are model animals for evolutionary developmental biology study. The embryonic development of tunicates is known to be extensively regulated by transcription factors (TFs). Styela clava, the globally distributed invasive tunicate, exhibits a strong capacity for environmental adaptation. However, the TFs were not systematically identified and analyzed. In this study, we reported 553 TFs categorized into 60 families from S. clava, based on the whole genome data. Comparison of TFs analysis among the tunicate species revealed that the gene number in the zinc finger superfamily displayed the most significant discrepancy, indicating this family was under the highly evolutionary selection and might be related to species differentiation and environmental adaptation. The greatest number of TFs was discovered in the Cys2His2-type zinc finger protein (zf-C2H2) family in S. clava. From the point of temporal view, more than half the TFs were expressed at the early embryonic stage. The expression correlation analysis revealed the existence of a transition for TFs expression from early embryogenesis to the later larval development in S. clava. Eight Hox genes were identified to be located on one chromosome, exhibiting different arrangement and expression patterns, compared to Ciona robusta (C. intestinalis type A). In addition, a total of 23 forkhead box (fox) genes were identified in S. clava, and their expression profiles referred to their potential roles in neurodevelopment and sensory organ development. Our data, thus, provides crucial clues to the potential functions of TFs in development and environmental adaptation in the leathery sea squirt.

Chromosome-level genome assembly of the Arctic fox (Vulpes lagopus) using PacBio sequencing and Hi-C technology

The Arctic fox (Vulpes lagopus) is the only fox species occurring in the Arctic and has adapted to its extreme climatic conditions. Currently, the molecular basis of its adaptation to the extreme climate has not been characterized. Here, we applied PacBio sequencing and chromosome structure capture technique to assemble the first V. lagopus genome assembly, which is assembled into chromosome fragments. The genome assembly has a total length of 2.345 Gb with a contig N50 of 31.848 Mb and a scaffold N50 of 131.537 Mb, consisting of 25 pseudochromosomal scaffolds. The V. lagopus genome had approximately 32.33% repeat sequences. In total, 21,278 protein-coding genes were predicted, of which 99.14% were functionally annotated. Compared with 12 other mammals, V. lagopus was most closely related to V. Vulpes with an estimated divergence time of ~7.1 Ma. The expanded gene families and positively selected genes potentially play roles in the adaptation of V. lagopus to Arctic extreme environment. This high-quality assembled genome will not only promote future studies of genetic diversity and evolution in foxes and other canids but also provide important resources for conservation of Arctic species.

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.

Seasonal dynamics and starvation impact on the gut microbiome of urochordate ascidian Halocynthia roretzi

Background

Gut microbiota plays important roles in host animal metabolism, homeostasis and environmental adaptation. However, the interplay between the gut microbiome and urochordate ascidian, the most closet relative of vertebrate, remains less explored. In this study, we characterized the gut microbial communities of urochordate ascidian (Halocynthia roretzi) across the changes of season and starvation stress using a comprehensive set of omic approaches including 16S rRNA gene amplicon sequencing, shotgun metagenomics, metabolomic profiling, and transcriptome sequencing.

Results

The 16S rRNA gene amplicon profiling revealed that ascidians harbor indigenous gut microbiota distinctly different to the marine microbial community and significant variations in composition and abundance of gut bacteria, with predominant bacterial orders representing each season. Depressed alpha-diversities of gut microbiota were observed across starvation stress when compared to the communities in aquafarm condition. Synechococcales involving photosynthesis and its related biosynthesis was reduced in abundance while the enrichments of Xanthomonadales and Legionellales may facilitate bile acid biosynthesis during starvation. Metabolomics analysis found that long chain fatty acids, linolenic acid, cyanoamino acid, and pigments derived from gut bacteria were upregulated, suggesting a beneficial contribution of the gut microbiome to the ascidian under starvation stress.

Conclusions

Our findings revealed seasonal variation of ascidian gut microbiota. Defense and energy-associated metabolites derived from gut microbiome may provide an adaptive interplay between gut microbiome and ascidian host that maintains a beneficial metabolic system across season and starvation stress. The diversity-generating metabolisms from both microbiota and host might lead to the co-evolution and environmental adaptation.

Graphical Abstract