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.

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.