Orthologs at the Base of the Olfactores Clade.

Tunicate orthologs in the human genome comprise just 84 genes of the 19,872 protein-coding genes and 23 of the 16,528 non-coding genes, yet they stand at the base of the Olfactores clade, which radiated to generate thousands of tunicate and vertebrate species. What were the powerful drivers among these genes that enabled this process? Many of these orthologs are present in gene families. We discuss the biological role of each family and the orthologs' quantitative contribution to the family. Most important was the evolution of a second type of cadherin. This, a Type II cadherin, had the property of detaching the cell containing that cadherin from cells that expressed the Type I class. The set of such Type II cadherins could now detach and move away from their Type I neighbours, a process which would eventually evolve into the formation of the neural crest, "the fourth germ layer", providing a wide range of possibilities for further evolutionary invention. A second important contribution were key additions to the broad development of the muscle and nerve protein and visual perception toolkits. These developments in mobility and vision provided the basis for the development of the efficient predatory capabilities of the Vertebrata.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species.

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

Molecular control of cellulosic fin morphogenesis in ascidians.

BACKGROUND: The tunicates form a group of filter-feeding marine animals closely related to vertebrates. They share with them a number of features such as a notochord and a dorsal neural tube in the tadpole larvae of ascidians, one of the three groups that make tunicates. However, a number of typical chordate characters have been lost in different branches of tunicates, a diverse and fast-evolving phylum. Consequently, the tunic, a sort of exoskeleton made of extracellular material including cellulose secreted by the epidermis, is the unifying character defining the tunicate phylum. In the larva of ascidians, the tunic differentiates in the tail into a median fin (with dorsal and ventral extended blades) and a caudal fin. RESULTS: Here we have performed experiments in the ascidian Phallusia mammillata to address the molecular control of tunic 3D morphogenesis. We have demonstrated that the tail epidermis medio-lateral patterning essential for peripheral nervous system specification also controls tunic elongation into fins. More specifically, when tail epidermis midline identity was abolished by BMP signaling inhibition, or CRISPR/Cas9 inactivation of the transcription factor coding genes Msx or Klf1/2/4/17, median fin did not form. We postulated that this genetic program should regulate effectors of tunic secretion. We thus analyzed the expression and regulation in different ascidian species of two genes acquired by horizontal gene transfer (HGT) from bacteria, CesA coding for a cellulose synthase and Gh6 coding for a cellulase. We have uncovered an unexpected dynamic history of these genes in tunicates and high levels of variability in gene expression and regulation among ascidians. Although, in Phallusia, Gh6 has a regionalized expression in the epidermis compatible with an involvement in fin elongation, our functional studies indicate a minor function during caudal fin formation only. CONCLUSIONS: Our study constitutes an important step in the study of the integration of HGT-acquired genes into developmental networks and a cellulose-based morphogenesis of extracellular material in animals.

The alternative enzymes-bearing tunicates lack multiple widely distributed genes coding for peripheral OXPHOS subunits.

The respiratory chain alternative enzymes (AEs) NDX and AOX from the tunicate Ciona intestinalis (Ascidiacea) have been xenotopically expressed and characterized in human cells in culture and in the model organisms Drosophila melanogaster and mouse, with the purpose of developing bypass therapies to combat mitochondrial diseases in human patients with defective complexes I and III/IV, respectively. The fact that the genes coding for NDX and AOX have been lost from genomes of evolutionarily successful animal groups, such as vertebrates and insects, led us to investigate if the composition of the respiratory chain of Ciona and other tunicates differs significantly from that of humans and Drosophila, to accommodate the natural presence of AEs. We have failed to identify in tunicate genomes fifteen orthologous genes that code for subunits of the respiratory chain complexes; all of these putatively missing subunits are peripheral to complexes I, III and IV in mammals, and many are important for complex-complex interaction in supercomplexes (SCs), such as NDUFA11, UQCR11 and COX7A. Modeling of all respiratory chain subunit polypeptides of Ciona indicates significant structural divergence that is consistent with the lack of these fifteen clear orthologous subunits. We also provide evidence using Ciona AOX expressed in Drosophila that this AE cannot access the coenzyme Q pool reduced by complex I, but it is readily available to oxidize coenzyme Q molecules reduced by glycerophosphate oxidase, a mitochondrial inner membrane-bound dehydrogenase that is not involved in SCs. Altogether, our results suggest that Ciona AEs might have evolved in a mitochondrial inner membrane environment much different from that of mammals and insects, possibly without SCs; this correlates with the preferential functional interaction between these AEs and non-SC dehydrogenases in heterologous mammalian and insect systems. We discuss the implications of these findings for the applicability of Ciona AEs in human bypass therapies and for our understanding of the evolution of animal respiratory chain.

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.

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.

Specification of distinct cell types in a sensory-adhesive organ important for metamorphosis in tunicate larvae.

The papillae of tunicate larvae contribute sensory, adhesive, and metamorphosis-regulating functions that are crucial for the biphasic lifestyle of these marine, non-vertebrate chordates. We have identified additional molecular markers for at least 5 distinct cell types in the papillae of the model tunicate Ciona, allowing us to further study the development of these organs. Using tissue-specific CRISPR/Cas9-mediated mutagenesis and other molecular perturbations, we reveal the roles of key transcription factors and signaling pathways that are important for patterning the papilla territory into a highly organized array of different cell types and shapes. We further test the contributions of different transcription factors and cell types to the production of the adhesive glue that allows for larval attachment during settlement, and to the processes of tail retraction and body rotation during metamorphosis. With this study, we continue working towards connecting gene regulation to cellular functions that control the developmental transition between the motile larva and sessile adult of Ciona.

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi.

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Effects of color variation and physiological state on ascidian microbiomes.

Ascidians, known for their color variation, host species-specific microbial symbiont communities. Some ascidians can also transition into a nonfiltering (resting) physiological state. Recent studies suggest that the microbial symbiont communities may vary across different physiological states and color morphs of the host. The colonial ascidian, Polyclinum constellatum, which exhibits several color morphs in the Caribbean Sea, periodically ceases its filtering activity. To investigate if color variation in P. constellatum is indicative of sibling speciation, we sequenced fragments of the ribosomal 18S rRNA and the mitochondrial cytochrome oxidase subunit I genes. Additionally, we sequenced a fragment of the 16S rRNA gene to characterize the microbial communities of two common color morphs (red and green) in colonies that were either actively filtering (active) or nonfiltering (resting). Phylogenetic analyses of both ascidian genes resulted in well-supported monophyletic clades encompassing all color variants of P. constellatum. Interestingly, no significant differences were observed among the microbial communities of the green and red morphs, suggesting that color variation in this species is a result of intraspecific variation. However, the host's physiological state significantly influenced the microbial community structure. Nonfiltering (resting) colonies hosted higher relative abundances of Kiloniella (Alphaproteobacteria) and Fangia (Gammaproteobacteria), while filtering colonies hosted more Reugeria (Alphaproteobacteria) and Endozoicomonas (Gammaproteobacteria). This study demonstrates that microbial symbiont communities serve as reliable indicators of the taxonomic state of their host and are strongly influenced by the host's feeding condition.

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.

Drop it all: extraction-free detection of targeted marine species through optimized direct droplet digital PCR.

Molecular biomonitoring programs increasingly use environmental DNA (eDNA) for detecting targeted species such as marine non-indigenous species (NIS) or endangered species. However, the current molecular detection workflow is cumbersome and time-demanding, and thereby can hinder management efforts and restrict the "opportunity window" for rapid management responses. Here, we describe a direct droplet digital PCR (direct-ddPCR) approach to detect species-specific free-floating extra-cellular eDNA (free-eDNA) signals, i.e., detection of species-specific eDNA without the need for filtration or DNA extraction, with seawater samples. This first proof-of-concept aquarium study was conducted with three distinct marine species: the Mediterranean fanworm Sabella spallanzanii, the ascidian clubbed tunicate Styela clava, and the brown bryozoan Bugula neritina to evaluate the detectability of free-eDNA in seawater. The detectability of targeted free-eDNA was assessed by directly analysing aquarium marine water samples using an optimized species-specific ddPCR assay. The results demonstrated the consistent detection of S. spallanzanii and B. neritina free-eDNA when these organisms were present in high abundance. Once organisms were removed, the free-eDNA signal exponentially declined, noting that free-eDNA persisted between 24-72 h. Results indicate that organism biomass, specimen characteristics (e.g., stress and viability), and species-specific biological differences may influence free-eDNA detectability. This study represents the first step in assessing the feasibility of direct-ddPCR technology for the detection of marine species. Our results provide information that could aid in the development of new technology, such as a field development of ddPCR systems, which could allow for automated continuous monitoring of targeted marine species, enabling point-of-need detection and rapid management responses.

Genome-wide analysis of early vascular tunic repair and regeneration for Botrylloides digenesis reveals striking similarities to human wound healing.

The early stages of regeneration after injury are similar to those of wound healing. The ascidian Botrylloides diegensis can regenerate an entire adult from a small fragment of vascular tunic following the removal of all zooids in an injury-induced regeneration model. We investigated the molecular and cellular changes following injury to determine the differences between the healing process and the initiation of whole-body regeneration (WBR). We conducted transcriptome analysis at specific time points during regeneration and wound healing to identify differentially expressed genes (DEGs) and the unique biological processes associated with each state. Our findings revealed 296 DEGs at 10 h post-injury (hpi), with 71 highly expressed in healed tissue and 225 expressed during the WBR process. These DEGs were predicted to play roles in tissue reorganization, integrin signaling, extracellular matrix organization, and the innate immune system. Pathway analysis of the upregulated genes in the healed tunic indicated functional enrichment related to tissue repair, as has been observed in other species. Additionally, we examined the cell types in the tunic and ampullae in both tissue states using histology and in situ hybridization for six genes identified by transcriptome analysis. We observed strong mRNA expression in cells within the WBR tunic, and in small RNA-positive granules near the tunic edge. We hypothesized that many of these genes function in the compaction of the ampullae tunic, which is a pivotal process for WBR and dormancy in B. diegensis, and in an immune response. These findings establish surprising similarities between ascidian regeneration and human wound healing, emphasizing the potential for future investigations into human regenerative and repair mechanisms. This study provides valuable insights into the gene sets specifically activated during regeneration compared to wound healing, shedding light on the divergent activities of these processes.

The transcription factor AP2 and downstream genes shared by asexual reproduction and zooidal regeneration in the tunicate, Polyandrocarpa misakiensis.

Epithelial outpocketing, tunic softening, mesenchymal cell death, dedifferentiation/transdifferentiation, and resistance to environmental stress are major events that occur during asexual reproduction by budding in the tunicate, Polyandrocarpa misakiensis. To identify the molecules underlying these events and compare them with those operating in regeneration, differential gene expression profiles were developed in buds and zooids. Among approximately 40,000 contigs, 21 genes were identified as potentially being involved in asexual reproduction. Genes related to tunic softening, phagocytosis-stimulating opsonin, and stress resistance were activated in the very early stage of budding. At the later stage of budding when buds separated from the parent and entered the developmental stage, genes for cell adhesion, cell death, and differentiation were activated. The transcription factor AP2 was spatio-temporally expressed in a similar pattern to the tunic-softening gene endoglucanase (EndoG). AP2 mRNA activated EndoG when introduced into zooids by electroporation. Eight out of 21 budding-related genes were significantly activated by AP2 mRNA. Polyandrocarpa zooids possess regenerative potential other than budding. Zooidal regeneration accompanied cell death/phagocytosis, cell-cell adhesion/communication, and dedifferentiation/redifferentiation. Consistent with morphological features, eight related genes including SP8 transcription factor were activated during zooidal regeneration. Most of these genes were identical to those induced by AP2 mRNA, indicating that asexual reproduction in P. misakiensis shares AP2-regulated downstream genes with zooidal regeneration. The present results suggest that SP8 may be indispensable for both budding and regeneration and that the potential dedifferentiation-related gene SOXB1 plays a minor role in zooidal regeneration.

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.

Stem cell-mediated development, regeneration, chimerism, and aging in the colonial chordate Botryllus schlosseri.

Stem cells are units of biological organization, responsible for tissue and organ development and regeneration. I study stem cell biology, aging, and the evolution of immunity using the colonial chordate Botryllus schlosseri as a model system. This organism is uniquely suited for this study because it is closely related to vertebrates, undergoes weekly cycles of stem cell mediated regeneration, is long lived and has a recognition system and robust immune system. I have led the Botryllus genome project and developed a novel method to obtain a synthetic long read sequence, identified Botryllus stem cells and stem cell niches, isolated the gene that controls self/non self-recognition and characterized its immune system on the cellular and molecular levels. Recently, I led the Botryllus atlas project to characterize the two developmental pathways, embryogenesis (sexual) and blastogenesis (asexual), revealing the unique molecular landscapes for each developmental mode and investigated the molecular clock and neurodegeneration pathways in young and old colonies and investigated the molecular clock and neurodegeneration pathways in young and old colonies. These results and the resources we developed are used by my lab and others to further study stem cell and immune cell properties during development, regeneration, transplantation, and aging.