The Rohde-like cells at the posterior end of the dorsal nerve cord of amphioxus (Cephalochordata).

For postmetamorphic specimens of amphioxus (Cephalochordata), serial block-face scanning electron microscopy (SBSEM) is used to describe the long-ignored Rohde-like cells (RLCs) at the extreme posterior end of the dorsal nerve cord. These cells, numbering about three dozen in all, are divisible into a group with larger diameters running near the dorsal side of the cord and a more ventral group with smaller diameters closely associated with the central canal of the neurocoel. It is possible that the smaller ventral cells might be generated at the ependymal zone of the dorsal nerve cord and later migrate to a dorsal position, although a functional reason for this remains a mystery. All the RLCs have conspicuous regions of microvilli covering as much as 40% of their surface; limited data (by others) on the more anterior bona fide Rohde cells also indicate an extensive microvillar surface. Thus, both the RLCs and the better-known Rohde cells appear to be rhabdomeric photoreceptors, although a specific function for this feature is currently unknown. Even more perplexingly, although the Rohde cells are quintessential neurons extending giant processes, each RLC comprises a perikaryon that does not bear any neurites.

Functional characterization of growth hormone releasing hormone and its receptor in amphioxus with implication for origin of hypothalamic-pituitary axis.

Growth hormone-releasing hormone (GHRH) has been widely shown to stimulate growth hormone (GH) production via binding to GHRH receptor GHRHR in various species of vertebrates, but information regarding the functional roles of GHRH and GHRHR in the protochordate amphioxus remains rather scarce. We showed here that two mature peptides, BjGHRH-1 and BjGHRH-2, encoded by BjGHRH precursor, and a single BjGHRHR protein were identified in the amphioxus Branchiostoma. japonicum. Like the distribution profiles of vertebrate GHRHs and GHRHRs, both the genes Bjghrh and Bjghrhr were widely expressed in the different tissues of amphioxus, including in the cerebral vesicle, Hatschek's pit, neural tube, gill, hepatic caecum, notochord, testis and ovary. Moreover, both BjGHRH-1 and BjGHRH-2 interacted with BjGHRHR, and triggered the cAMP/PKA signal pathway in a dose-dependent manner. Importantly, BjGHRH-1 and BjGHRH-2 were both able to activate the expression of GH-like gene in the cells of Hatschek's pit. These indicate that a functional vertebrate-like GHRH-GHRHR axis had already emerged in amphioxus, which is a seminal innovation making physiological divergence including reproduction, growth, metabolism, stress and osmoregulation possible during the early evolution of vertebrates.

An amphioxus neurula stage cell atlas supports a complex scenario for the emergence of vertebrate head mesoderm.

The emergence of new structures can often be linked to the evolution of novel cell types that follows the rewiring of developmental gene regulatory subnetworks. Vertebrates are characterized by a complex body plan compared to the other chordate clades and the question remains of whether and how the emergence of vertebrate morphological innovations can be related to the appearance of new embryonic cell populations. We previously proposed, by studying mesoderm development in the cephalochordate amphioxus, a scenario for the evolution of the vertebrate head mesoderm. To further test this scenario at the cell population level, we used scRNA-seq to construct a cell atlas of the amphioxus neurula, stage at which the main mesodermal compartments are specified. Our data allowed us to validate the presence of a prechordal-plate like territory in amphioxus. Additionally, the transcriptomic profile of somite cell populations supports the homology between specific territories of amphioxus somites and vertebrate cranial/pharyngeal and lateral plate mesoderm. Finally, our work provides evidence that the appearance of the specific mesodermal structures of the vertebrate head was associated to both segregation of pre-existing cell populations, and co-option of new genes for the control of myogenesis.

A novel C-type lectin protein (BjCTL5) interacts with apoptosis stimulating proteins of p53 (ASPP) to activate NF-κB signaling pathway in primitive chordate.

C-type lectin proteins (CTLs), a group of pattern recognition receptors (PRRs), play pivotal roles in immune responses. However, the signal transduction and regulation of CTLs in cephalochordates have yet to be explored. In this study, we examined the composition of CTLs in Branchiostoma japonicum, identifying a total of 272 CTLs. These CTLs underwent further analysis concerning domain arrangement, tandem and segmental duplication events. A multidomain C-type lectin gene, designated as BjCTL5, encompassing CLECT, KR, CUB, MAM, and SR domains, was the focal point of our investigation. BjCTL5 exhibits ubiquitous expression across all detected tissues and is responsive to stimulation by LPS, mannose, and poly (I:C). The recombinant protein of BjCTL5 can bind to Escherichia coli and Staphylococcus aureus, inducing their agglutination and inhibiting the proliferation of S. aureus. Yeast two-hybrid, CoIP, and confocal immunofluorescence experiments revealed the interaction between BjCTL5 and apoptosis-stimulating proteins of p53, BjASPP. Intriguingly, BjCTL5 was observed to induce the luciferase activity of the NF-κB promoter in HEK293T cells. These results suggested a potential interaction between BjCTL5 and BjASPP, implicating that they involve in the activation of the NF-κB signaling pathway, which provides an evolutionary viewpoint on NF-κB signaling pathway in primitive chordate.

Molecular and functional characterization of a type-1 cystatin in amphioxus (Branchiostoma japonicum).

Cystatins comprise a vast superfamily of evolutionary conserved proteins, predominantly recognized for their roles as endogenous inhibitors by regulating the activity of cysteine proteases. Emerging lines of research evidence also provides insight into their alternative roles in a spectrum of biological and pathological processes, including neurodegenerative disorders, tumor progression, inflammatory diseases, and immune response. Nowadays, various type-1 cystatins (stefins) have been demonstrated among a variety of discovered vertebrate groups, while little is known about the related homologue in cephalochordate amphioxus, which are repositioned at the base of the chordate phylum. In the present study, a single type-1 cystatin homologue in Branchiostoma japonicum was first successfully cloned and designated as Bjcystatin-1. The deduced Bjcystatin-1 protein is structurally characterized by the presence of typical wedge-shaped cystatin features, including the 'QxVxG' and 'Px' motif, as well as the conserved N-terminal glycine residue. Phylogenomic analyses utilizing different cystatin counterparts affirmed the close evolutionary relationship of Bjcystatin-1 and type-1 cystatin homologue. Bjcystatin-1 was predominantly expressed in the gills and hind-gut in a tissue-specific pattern, and its expression was remarkably up-regulated in response to challenge with bacteria or their signature molecules LPS and LTA, suggesting the involvement in immune response. Additionally, the recombinant Bjcystatin-1 (rBjcystatin-1) protein showed significant inhibitory activity towards papain and binding ability to LPS and LTA, indicating its hypothesized role as a pattern recognition receptor in immune response. Subcellular localization results also showed that Bjcystatin-1 was located in the cytoplasm and nucleus, and its overexpression could attenuate the activation of LPS-induced nuclear transcription factors NF-κB. Taken together, our study suggests that amphioxus Bjcystatin-1 acts as a dual role in protease inhibitor and an immunocompetent factor, providing new insights into the immune defense effect of type-1 cystatin in amphioxus.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus.

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

scRNA-seq analysis of cells comprising the amphioxus notochord.

Cephalochordates occupy a key phylogenetic position for deciphering the origin and evolution of chordates, since they diverged earlier than urochordates and vertebrates. The notochord is the most prominent feature of chordates. The amphioxus notochord features coin-shaped cells bearing myofibrils. Notochord-derived hedgehog signaling contributes to patterning of the dorsal nerve cord, as in vertebrates. However, properties of constituent notochord cells remain unknown at the single-cell level. We examined these properties using Iso-seq analysis, single-cell RNA-seq analysis, and in situ hybridization (ISH). Gene expression profiles broadly categorize notochordal cells into myofibrillar cells and non-myofibrillar cells. Myofibrillar cells occupy most of the central portion of the notochord, and some cells extend the notochordal horn to both sides of the ventral nerve cord. Some notochord myofibrillar genes are not expressed in myotomes, suggesting an occurrence of myofibrillar genes that are preferentially expressed in notochord. On the other hand, non-myofibrillar cells contain dorsal, lateral, and ventral Müller cells, and all three express both hedgehog and Brachyury. This was confirmed by ISH, although expression of hedgehog in ventral Müller cells was minimal. In addition, dorsal Müller cells express neural transmission-related genes, suggesting an interaction with nerve cord. Lateral Müller cells express hedgehog and other signaling-related genes, suggesting an interaction with myotomes positioned lateral to the notochord. Ventral Müller cells also expressed genes for FGF- and EGF-related signaling, which may be associated with development of endoderm, ventral to the notochord. Lateral Müller cells were intermediate between dorsal/ventral Müller cells. Since vertebrate notochord contributes to patterning and differentiation of ectoderm (nerve cord), mesoderm (somite), and endoderm, this investigation provides evidence that an ancestral or original form of vertebrate notochord is present in extant cephalochordates.

Serial block-face scanning electron microscopy of the tail tip of post-metamorphic amphioxus finds novel myomeres with odd shapes and unusually prominent sclerocoels.

Serial block-face scanning electron microscopy of the tail tip of post-metamorphic amphioxus (Branchiostoma floridae) revealed some terminal myomeres never been seen before with other techniques. The morphology of these myomeres differed markedly from the chevron shapes of their more anterior counterparts. Histologically, these odd-shaped myomeres ranged from empty vesicles bordered by undifferentiated cells to ventral sacs composed of well-developed myotome, dermatome, and sclerotome. Strikingly, several of these ventral sacs gave rise to a nipple-like dorsal projection composed either entirely of sclerotome or a mixture of sclerotome and myotome. Considered as a whole, from posterior to anterior, these odd-shaped posterior myomeres suggested that their more substantial ventral part may represent the ventral limb of a chevron, while the delicate projection represents a nascent dorsal limb. This scenario contrasts with formation of chevron-shaped myomeres along most of the antero-posterior axis. Although typical chevron formation in amphioxus is surprisingly poorly studied, it seems to be attained by a dorso-ventral extension of the myomere accompanied by the assumption of a V-shape; this is similar to what happens (at least superficially) in developing fishes. Another unusual feature of the odd-shaped posterior myomeres of amphioxus is their especially distended sclerocoels. One possible function for these might be to protect the posterior end of the central nervous system from trauma when the animals burrow into the substratum.

Comparison of Evolutionary Relationships between Branchiostoma floridae, Ciona intestinalis, and Homo sapiens Globins Provide Evidence of Gene Co-Option and Convergent Evolution.

Globins have been studied as model proteins to elucidate the principles of protein evolution. This was achieved by understanding the relationship between amino acid sequence, three-dimensional structure, physicochemical properties, and physiological function. Previous molecular phylogenies of chordate globin genes revealed the monophyletic evolution of urochordate globins and suggested convergent evolution. However, to provide evidence of convergent evolution, it is necessary to determine the physicochemical and functional similarities between vertebrates and urochordate globins. In this study, we determined the expression patterns of Ciona globin genes using real-time RT-PCR. Two genes (Gb-1 and Gb-2) were predominantly expressed in the branchial sac, heart, and hemocytes and were induced under hypoxia. Combined with the sequence analysis, our findings suggest that Gb-1/-2 correspond to vertebrate hemoglobin-α/-ß. However, we did not find a robust similarity between Gb-3, Gb-4, and vertebrate globins. These results suggested that, even though Ciona globins obtained their unique functions differently from vertebrate globins, the two of them shared some physicochemical features and physiological functions. Our findings offer a good example for understanding the molecular mechanisms underlying gene co-option and convergence, which could lead to evolutionary innovations.

Cephalochordate Hemocytes: First Demonstration for Asymmetron lucayanum (Bahamas Lancelet) Plus Augmented Description for Branchiostoma floridae (Florida Amphioxus).

AbstractWithin phylum Chordata, the subphylum Cephalochordata (amphioxus and lancelets) has figured large in considerations of the evolutionary origin of the vertebrates. To date, these discussions have been predominantly based on knowledge of a single cephalochordate genus (Branchiostoma), almost to the exclusion of the other two genera (Asymmetron and Epigonichthys). This uneven pattern is illustrated by cephalochordate hematology, until now known entirely from work done on Branchiostoma. The main part of the present study is to describe hemocytes in the dorsal aorta of a species of Asymmetron by serial block-face scanning electron microscopy. This technique, which demonstrates three-dimensional fine structure, showed that the hemocytes have a relatively uniform morphology characterized by an oval shape and scanty cytoplasm. Ancillary information is also included for Branchiostoma hemocytes, known from previous studies to have relatively abundant cytoplasm; our serial block-face scanning electron microscopy provides more comprehensive views of the highly variable shapes of these cells, which typically extend one or several pseudopodium-like protrusions. The marked difference in hemocyte morphology found between Asymmetron and Branchiostoma was unexpected and directs attention to investigating comparable cells in the genus Epigonichthys. A broader knowledge of the hemocytes in all three cephalochordate genera would provide more balanced insights into the evolution of vertebrate hematopoiesis.

Amphioxus as a model to study the evolution of development in chordates.

Cephalochordates and tunicates represent the only two groups of invertebrate chordates, and extant cephalochordates - commonly known as amphioxus or lancelets - are considered the best proxy for the chordate ancestor, from which they split around 520 million years ago. Amphioxus has been an important organism in the fields of zoology and embryology since the 18th century, and the morphological and genomic simplicity of cephalochordates (compared to vertebrates) makes amphioxus an attractive model for studying chordate biology at the cellular and molecular levels. Here we describe the life cycle of amphioxus, and discuss the natural histories and habitats of the different species of amphioxus. We also describe their use as laboratory animal models, and discuss the techniques that have been developed to study different aspects of amphioxus.

Amphioxus adenosine-to-inosine tRNA-editing enzyme that can perform C-to-U and A-to-I deamination of DNA.

Adenosine-to-inosine tRNA-editing enzyme has been identified for more than two decades, but the study on its DNA editing activity is rather scarce. We show that amphioxus (Branchiostoma japonicum) ADAT2 (BjADAT2) contains the active site 'HxE-PCxxC' and the key residues for target-base-binding, and amphioxus ADAT3 (BjADAT3) harbors both the N-terminal positively charged region and the C-terminal pseudo-catalytic domain important for recognition of substrates. The sequencing of BjADAT2-transformed Escherichia coli genome suggests that BjADAT2 has the potential to target E. coli DNA and can deaminate at TCG and GAA sites in the E. coli genome. Biochemical analyses further demonstrate that BjADAT2, in complex with BjADAT3, can perform A-to-I editing of tRNA and convert C-to-U and A-to-I deamination of DNA. We also show that BjADAT2 preferentially deaminates adenosines and cytidines in the loop of DNA hairpin structures of substrates, and BjADAT3 also affects the type of DNA substrate targeted by BjADAT2. Finally, we find that C89, N113, C148 and Y156 play critical roles in the DNA editing activity of BjADAT2. Collectively, our study indicates that BjADAT2/3 is the sole naturally occurring deaminase with both tRNA and DNA editing capacity identified so far in Metazoa.

Genome-wide transcriptomics and microRNAomics analyses uncover multi-faceted mechanisms to cope with copper stress in ancient macrobenthos amphioxus.

The mechanisms underlying the toxicity of environmental stress are unclear for marine macrobenthos. Copper/Cu has posed the most serious threats to amphioxus, an ancient and model benthic cephalochordate. Herein, a dynamic change in the physiological parameters (GR, SOD, ATP, and MDA) was detected with ROS accumulation in Branchiostoma belcheri exposed to 0.3 mg·L-1 Cu. Transcriptomes and microRNAomes of B. belcheri were generated to investigate the molecular mechanisms by which this amphioxus copes with Cu exposure. Time-specific genes identified at different time points after exposure were involved in the stimulus and immune response, detoxification and ionic homeostasis, aging and the nervous system, sequentially, with prolongation of exposure time, forming a dynamic process of molecular response to Cu stress. In total, 57 differentially expressed miRNAs were identified under Cu stress. Transcriptomics-miRNAomics analyses indicate that these miRNAs targeted genes associated with many key biological processes such as xenobiotics degradation, oxidative stress, and energy metabolism. The constructed miRNA-mRNA-pathway network uncovered a broad post-transcriptional regulatory mechanism in B. belcheri to cope with Cu stress. Overall, this integrated analyses show that enhanced defense response, accelerated ROS elimination, and repressed ATP production constitute a comprehensive strategy to cope with Cu toxicity in the ancient macrobenthos.

Characterization of GSDME in amphioxus provides insights into the functional evolution of GSDM-mediated pyroptosis.

Members of the gasdermin (GSDM) family are pore-forming effectors that cause membrane permeabilization and pyroptosis, a lytic proinflammatory type of cell death. To reveal the functional evolution of GSDM-mediated pyroptosis at the transition from invertebrates to vertebrates, we conducted functional characterization of amphioxus GSDME (BbGSDME) and found that it can be cleaved by distinct caspase homologs, yielding the N253 and N304 termini with distinct functions. The N253 fragment binds to cell membrane, triggers pyroptosis, and inhibits bacterial growth, while the N304 performs negative regulation of N253-mediated cell death. Moreover, BbGSDME is associated with bacteria-induced tissue necrosis and transcriptionally regulated by BbIRF1/8 in amphioxus. Interestingly, several amino acids that are evolutionarily conserved were found to be important for the function of both BbGSDME and HsGSDME, shedding new lights on the functional regulation of GSDM-mediated inflammation.

In Amphioxus Embryos, Some Neural Tube Cells Resemble Differentiating Coronet Cells of Fishes and Tunicates.

AbstractFor neurula embryos of amphioxus (chordate subphylum Cephalochordata), the anterior region of the neural tube was studied with transmission electron microscopy. This survey demonstrated previously unreported cells, each characterized by a cilium bearing on its shaft a protruding lateral bubble packed with vesicles. Such cilia resemble those known from immature coronet cells in other chordates-namely, fishes in the Vertebrata and ascidians and appendicularians in the Tunicata. This wide occurrence of coronet-like cells raises questions about their possible homologies within the phylum Chordata. When considered at the level of the whole cell, such homology is not well supported. For example, the fish cells are generally thought to be glia, while the tunicate cells are considered to be neurons; moreover, cytoplasmic smooth endoplasmic reticulum, which is predominant in the former, is undetectable in the latter. In contrast, a more convincing case for homology can be made by limiting comparisons to the cell apices with their modified cilia. In addition to the fine-structural similarities between fishes and tunicates already mentioned, nonvisual opsins have been found associated with the vesicles in the modified cilia of both groups. Such opsins are thought to link photoreception to endocrine output controlling behavior. Further work would be needed to test the idea that the amphioxus diencephalic cells with lateral bubble cilia might similarly be opsin rich and could provide insights into the evolutionary history of the coronet cells within the phylum Chordata.

An ancient defense mechanism: Conservation of gasdermin-mediated pyroptosis.

The gasdermins are a family of pore-forming proteins involved in various cellular processes such as cell death and inflammation. A new study in PLOS Biology explores the evolutionary history of gasdermins across metazoans, highlighting the conservation and divergence of gasdermin E.

Genomic Analysis of Amphioxus Reveals a Wide Range of Fragments Homologous to Viral Sequences.

Amphioxus species are considered living fossils and are important in the evolutionary study of chordates and vertebrates. To explore viral homologous sequences, a high-quality annotated genome of the Beihai amphioxus (Branchiostoma belcheri beihai) was examined using virus sequence queries. In this study, 347 homologous fragments (HFs) of viruses were identified in the genome of B. belcheri beihai, of which most were observed on 21 genome assembly scaffolds. HFs were preferentially located within protein-coding genes, particularly in their CDS regions and promoters. A range of amphioxus genes with a high frequency of HFs is proposed, including histone-related genes that are homologous to the Histone or Histone H2B domains of viruses. Together, this comprehensive analysis of viral HFs provides insights into the neglected role of viral integration in the evolution of amphioxus.

Identification of a novel amphioxus leucine-rich repeat receptor involved in phagocytosis reveals a role for Slit2-N-type LRR in bacterial elimination.

The basal chordate amphioxus is a model for tracing the origin and evolution of vertebrate immunity. To explore the evolution of immunoreceptor signaling pathways, we searched the associated receptors of the amphioxus Branchiostoma belcheri (Bb) homolog of immunoreceptor signaling adaptor protein Grb2. Mass-spectrum analysis of BbGrb2 immunoprecipitates from B. belcheri intestine lysates revealed a folate receptor (FR) domain- and leucine-rich repeat (LRR)-containing protein (FrLRR). Sequence and structural analysis showed that FrLRR is a membrane protein with a predicted curved solenoid structure. The N-terminal Fr domain contains very few folate-binding sites; the following LRR region is a Slit2-type LRR, and a GPI-anchored site was predicted at the C-terminus. RT-PCR analysis showed FrLRR is a transcription-mediated fusion gene of BbFR-like and BbSlit2-N-like genes. Genomic DNA structure analysis implied the B. belcheri FrLRR gene locus and the corresponding locus in Branchiostoma floridae might be generated by exon shuffling of a Slit2-N-like gene into an FR gene. RT-qPCR, immunostaining, and immunoblot results showed that FrLRR was primarily distributed in B. belcheri intestinal tissue. We further demonstrated that FrLRR localized to the cell membrane and lysosomes. Functionally, FrLRR mediated and promoted bacteria-binding and phagocytosis, and FrLRR antibody blocking or Grb2 knockdown inhibited FrLRR-mediated phagocytosis. Interestingly, we found that human Slit2-N (hSlit2-N) also mediated direct bacteria-binding and phagocytosis which was inhibited by Slit2-N antibody blocking or Grb2 knockdown. Together, these results indicate FrLRR and hSlit2-N may function as phagocytotic-receptors to promote phagocytosis through Grb2, implying the Slit2-N-type-LRR-containing proteins play a role in bacterial binding and elimination.

SUMO and PIAS repress NF-κB activation in a basal chordate.

Small ubiquitin-like modifier (SUMO) regulates various biological processes, including the MyD88/TICAMs-IRAKs-TRAF6-NF-κB pathway, one of the core immune pathways. However, its functions are inconsistent between invertebrates and vertebrates and have rarely been investigated in lower chordates, including amphioxus and fishes. Here, we investigated the SUMOylation gene system in the amphioxus, a living basal chordate. We found that amphioxus has a SUMOylation system that has a complete set of genes and preserves several ancestral traits. We proceeded to study their molecular functions using the mammal cell lines. Both amphioxus SUMO1 and SUMO2 were shown to be able to attach to NF-κB Rel and to inhibit NF-κB activation by 50-75% in a dose-dependent fashion. The inhibition by SUMO2 could be further enhanced by the addition of the SUMO E2 ligase UBC9. In comparison, while human SUMO2 inhibited RelA, human SUMO1 slightly activated RelA. We also showed that, similar to human PIAS1-4, amphioxus PIAS could serve as a SUMO E3 ligase and promote its self-SUMOylation. This suggests that amphioxus PIAS is functionally compatible in human cells. Moreover, we showed that amphioxus PIAS is not only able to inhibit NF-κB activation induced by MyD88, TICAM-like, TRAF6 and IRAK4 but also able to suppress NF-κB Rel completely in the presence of SUMO1/2 in a dose-insensitive manner. This suggests that PIAS could effectively block Rel by promoting Rel SUMOylation. In comparison, in humans, only PIAS3, but not PIAS1/2/4, has been reported to promote NF-κB SUMOylation. Taken together, the findings from amphioxus, together with those from mammals and other species, not only offer insights into the functional volatility of the animal SUMO system, but also shed light on its evolutionary transitions from amphioxus to fish, and ultimately to humans.

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates.

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.