Decelerated genome evolution in modern vertebrates revealed by analysis of multiple lancelet genomes

Characterization of Myxovirus resistance (Mx) gene from Chinese seabass Lateolabrax maculatus: Insights into the evolution and function of Mx genes

Chinese seabass (Lateolabrax maculatus) stands out as one of the most sought-after and economically significant species in aquaculture within China. Diseases of L. maculatus occur frequently due to the degradation of the germplasm, the aggravation of environmental pollution of water, and the reproduction of pathogenic microorganisms, inflicting considerable economic losses on the Chinese seabass industry. The Myxovirus resistance (Mx) gene plays pivotal roles in the antiviral immune response ranging from mammals to fish. However, the function of the Mx gene in L. maculatus is still unknown. Firstly, the origin and evolutionary history of Mx proteins was elucidated in this study. Subsequently, an Mx gene from L. maculatus (designed as LmMxA gene) was identified, and its functions in combating antiviral and antibacterial threats were investigated. Remarkably, our findings suggested that while Mx group genes were present in chordates, DYN group genes were present in everything from single-celled animals to humans. Furthermore, our investigation revealed that the LmMxA mRNA level increased in the kidney, spleen and liver subsequent to Vibrio anguillarum and poly(I:C) challenged. Immunofluorescence analysis indicated that LmMxA is predominantly localization in the nucleus and the cytoplasm. Notably, the expression of MAVS, IFN1 and Mx1 increased when LmMxA was overexpression within the EPC cells. Moreover, through assessment via cytopathic effect (CPE), virus titer, and antibacterial activity, it becomes evident that LmMxA exerts a dual role in bolstering both antiviral and antibacterial immune responses. These compelling findings laid the foundation for further exploring the mechanism of LmMxA in response to innate immunity of L. maculatus.

FishTEDB 2.0: an update fish transposable element (TE) database with new functions to facilitate TE research

We launched the initial version of FishTEDB in 2018, which aimed to establish an open-source, user-friendly, data-rich transposable element (TE) database. Over the past 5 years, FishTEDB 1.0 has gained approximately 10 000 users, accumulating more than 450 000 interactions. With the unveiling of extensive fish genome data and the increasing emphasis on TE research, FishTEDB needs to extend the richness of data and functions. To achieve the above goals, we introduced 33 new fish species to FishTEDB 2.0, encompassing a wide array of fish belonging to 48 orders. To make the updated database more functional, we added a genome browser to visualize the positional relationship between TEs and genes and the estimated TE insertion time in different species. In conclusion, we released a new version of the fish TE database, FishTEDB 2.0, designed to assist researchers in the future study of TE functions and promote the progress of biological theories related to TEs. Database URL: https://www.fishtedb.com/.

The IL-17 pathway intertwines with neurotrophin and TLR/IL-1R pathways since its domain shuffling origin

The IL-17 pathway displays remarkably diverse functional modes between different subphyla, classes, and even orders, yet its driving factors remains elusive. Here, we demonstrate that the IL-17 pathway originated through domain shuffling between a Toll-like receptor (TLR)/IL-1R pathway and a neurotrophin-RTK (receptor-tyrosine-kinase) pathway (a Trunk-Torso pathway). Unlike other new pathways that evolve independently, the IL-17 pathway remains intertwined with its donor pathways throughout later evolution. This intertwining not only influenced the gains and losses of domains and components in the pathway but also drove the diversification of the pathway's functional modes among animal lineages. For instance, we reveal that the crustacean female sex hormone, a neurotrophin inducing sex differentiation, could interact with IL-17Rs and thus be classified as true IL-17s. Additionally, the insect prothoracicotropic hormone, a neurotrophin initiating ecdysis in Drosophila by binding to Torso, could bind to IL-17Rs in other insects. Furthermore, IL-17R and TLR/IL-1R pathways maintain crosstalk in amphioxus and zebrafish. Moreover, the loss of the Death domain in the pathway adaptor connection to IκB kinase and stress-activated protein kinase (CIKSs) dramatically reduced their abilities to activate nuclear factor-kappaB (NF-κB) and activator protein 1 (AP-1) in amphioxus and zebrafish. Reinstating this Death domain not only enhanced NF-κB/AP-1 activation but also strengthened anti-bacterial immunity in zebrafish larvae. This could explain why the mammalian IL-17 pathway, whose CIKS also lacks Death, is considered a weak signaling activator, relying on synergies with other pathways. Our findings provide insights into the functional diversity of the IL-17 pathway and unveil evolutionary principles that could govern the pathway and be used to redesign and manipulate it.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Genome-Wide Comparative Analysis of SRCR Gene Superfamily in Invertebrates Reveals Massive and Independent Gene Expansions in the Sponge and Sea Urchin

Without general adaptative immunity, invertebrates evolved a vast number of heterogeneous non-self recognition strategies. One of those well-known adaptations is the expansion of the immune receptor gene superfamily coding for scavenger receptor cysteine-rich domain containing proteins (SRCR) in a few invertebrates. Here, we investigated the evolutionary history of the SRCR gene superfamily (SRCR-SF) across 29 metazoan species with an emphasis on invertebrates. We analyzed their domain architectures, genome locations and phylogenetic distribution. Our analysis shows extensive genome-wide duplications of the SRCR-SFs in Amphimedon queenslandica and Strongylocentrotus purpuratus. Further molecular evolution study reveals various patterns of conserved cysteines in the sponge and sea urchin SRCR-SFs, indicating independent and convergent evolution of SRCR-SF expansion during invertebrate evolution. In the case of the sponge SRCR-SFs, a novel motif with seven conserved cysteines was identified. Exon-intron structure analysis suggests the rapid evolution of SRCR-SFs during gene duplications in both the sponge and the sea urchin. Our findings across nine representative metazoans also underscore a heightened expression of SRCR-SFs in immune-related tissues, notably the digestive glands. This observation indicates the potential role of SRCR-SFs in reinforcing distinct immune functions in these invertebrates. Collectively, our results reveal that gene duplication, motif structure variation, and exon-intron divergence might lead to the convergent evolution of SRCR-SF expansions in the genomes of the sponge and sea urchin. Our study also suggests that the utilization of SRCR-SF receptor duplication may be a general and basal strategy to increase immune diversity and tissue specificity for the invertebrates.

Identification and expression analysis of transcripts involved in taurine biosynthesis during early ontogeny of tropical gar Atractosteus tropicus

In fishes, the availability of taurine is regulated during ontogenetic development, where its endogenous synthesis capacity is species dependent. Thus, different pathways and involved enzymes have been described: pathway I (cysteine sulfinate-dependent pathway), cysteine dioxygenase type 1 (cdo1) and cysteine sulfinic acid decarboxylase (csad); pathway II (cysteic acid pathway), cdo1 and glutamic acid decarboxylase (gad); and pathway III (cysteamine pathway), 2-aminoethanethiol dioxygenase (ado); whereas taurine transporter (taut) is responsible for taurine entry into cells on the cell membrane and the mitochondria. This study determined if the tropical gar (Atractosteus tropicus), an ancient holostean fish model, has the molecular mechanism to synthesize taurine through the identification and analysis expression of transcripts coding for proteins involved in its biosynthesis and transportation, at different embryo-larvae stages and in different organs of juveniles (31 dah). We observed a fluctuating expression of all transcripts involved in the three pathways at all analyzed stages. All transcripts are expressed during the beginning of larval development; however, ado and taut show a peak expression at 9 dah, and all transcripts but csad decreased at 23 dah, when the organism ended the larval period. Furthermore, at 31 dah, we observed taut expression in all examined organs. The transcripts involved in pathways I and III are expressed differently across all organs, whereas pathway II was only observed in the brain, eye, and skin. The results suggested that taurine biosynthesis in tropical gar is regulated during its early development before first feeding, and the pathway might also be organ-type dependent.

Ancient fish lineages illuminate toll-like receptor diversification in early vertebrate evolution

Since its initial discovery over 50 years ago, understanding the evolution of the vertebrate RAG- mediated adaptive immune response has been a major area of research focus for comparative geneticists. However, how the evolutionary novelty of an adaptive immune response impacted the diversity of receptors associated with the innate immune response has received considerably less attention until recently. Here, we investigate the diversification of vertebrate toll-like receptors (TLRs), one of the most ancient and well conserved innate immune receptor families found across the Tree of Life, integrating genomic data that represent all major vertebrate lineages with new transcriptomic data from Polypteriformes, the earliest diverging ray-finned fish lineage. Our analyses reveal TLR sequences that reflect the 6 major TLR subfamilies, TLR1, TLR3, TLR4, TLR5, TLR7, and TLR11, and also currently unnamed, yet phylogenetically distinct TLR clades. We additionally recover evidence for a pulse of gene gain coincident with the rise of the RAG-mediated adaptive immune response in jawed vertebrates, followed by a period of rapid gene loss during the Cretaceous. These gene losses are primarily concentrated in marine teleost fish and synchronous with the mid Cretaceous anoxic event, a period of rapid extinction for marine species. Finally, we reveal a mismatch between phylogenetic placement and gene nomenclature for up to 50% of TLRs found in clades such as ray-finned fishes, cyclostomes, amphibians, and elasmobranchs. Collectively, these results provide an unparalleled perspective of TLR diversity and offer a ready framework for testing gene annotations in non-model species.

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.

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.

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 B. 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 B. 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.

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.

'Cannibalism' of exogenous DNA sequences: The ancestral form of adaptive immunity which entails recognition of danger

Adaptive immunity is a sophisticated form of immune response capable of retaining the molecular memory of a very great diversity of target antigens (epitopes) as non-self. It is capable of reactivating itself upon a second encounter with an immunoglobulin or T-cell receptor antigen-binding site with a known epitope that had previously primed the host immune system. It has long been considered that adaptive immunity is a highly evolved form of non-self recognition that appeared quite late in speciation and complemented a more generalist response called innate immunity. Innate immunity offers a relatively non-specific defense (although mediated by sensors that could specifically recognize virus or bacteria compounds) and which does not retain a memory of the danger. But this notion of recent acquisition of adaptive immunity is challenged by the fact that another form of specific recognition mechanisms already existed in prokaryotes that may be able to specifically auto-protect against external danger. This recognition mechanism can be considered a primitive form of specific (adaptive) non-self recognition. It is based on the fact that many archaea and bacteria use a genome editing system that confers the ability to appropriate viral DNA sequences allowing prokaryotes to prevent host damage through a mechanism very similar to adaptive immunity. This is indistinctly called, 'endogenization of foreign DNA' or 'viral DNA predation' or, more pictorially 'DNA cannibalism'. For several years evidence has been accumulating, highlighting the crucial role of endogenization of foreign DNA in the fundamental processes related to adaptive immunity and leading to a change in the dogma that adaptive immunity appeared late in speciation.

Superconserved receptors expressed in the brain: Expression, function, motifs and evolution of an orphan receptor family

GPR27, GPR85 and GPR173 constitute a small family of G protein-coupled receptors (GPCR) that share the distinctive characteristics of being highly conserved throughout vertebrate evolution and predominantly expressed in the brain. Accordingly, they have been coined as "Superconserved Receptors Expressed in the Brain" (SREB), although their expression profile is more complex than what was originally thought. SREBs have no known validated endogenous ligands and are thus labeled as "orphan" receptors. The investigation of this particular category of uncharacterized receptors holds great promise both in terms of physiology and drug development. In the largest GPCR family, the Rhodopsin-like or Class A, around 100 receptors are considered orphans. Because GPCRs are the most successful source of drug targets, the discovery of a novel function or ligand most likely will lead to significant breakthroughs for the discovery of innovative therapies. The high level of conservation is one of the characteristic features of the SREBs. We propose herein a detailed analysis of the putative evolutionary origin of this family. We highlight the properties that distinguish SREBs from other rhodopsin-like GPCRs. We present the current evidence for these receptors downstream signaling pathways and functions. We discuss the pharmacological challenge for the identification of natural or synthetic ligands of orphan receptors like SREBs. The different SREB-related scientific questions are presented with a highlight on what should be addressed in the near future, including the confirmation of published evidence and their validation as drug targets. In particular, we discuss in which pathological conditions these receptors may be of great relevance to solve unmet medical needs.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications

BACKGROUND

Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates.

RESULTS

We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements.

CONCLUSIONS

In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

Two novel mollusk short-form ApeC-containing proteins act as pattern recognition proteins for peptidoglycan

The Apextrin C-terminal (ApeC) domain is a new protein domain largely specific to aquatic invertebrates. In amphioxus, a short-form ApeC-containing protein (ACP) family is capable of binding peptidoglycan (PGN) and agglutinating bacteria via its ApeC domain. However, the functions of ApeC in other phyla remain unknown. Here we examined 130 ACPs from gastropods and bivalves, the first and second biggest mollusk classes. They were classified into nine groups based on their phylogenetics and architectures, including three groups of short-form ACPs, one group of apextrins and two groups of ACPs of complex architectures. No groups have orthologs in other phyla and only four groups have members in both gastropods and bivalves, suggesting that mollusk ACPs are highly diversified. We selected one bivalve ACP (CgACP1; from the oyster Crossostrea gigas) and one gastropod ACP (BgACP1; from the snail Biomphalaria glabrata) for functional experiments. Both are highly-expressed, secreted short-form ACPs and hence comparable to the amphioxus ACPs previously reported. We found that recombinant CgACP1 and BgACP1 bound with yeasts and several bacteria with different affinities. They also agglutinated these microbes, but showed no inhibiting or killing effects. Further analyses show that both ACPs had high affinities to the Lys-type PGN from S. aureus but weak or no affinities to the DAP-type PGN from Bacillus subtilis. Both recombinant ACPs displayed weak or no affinities to other microbial cell wall components, including lipopolysaccharide (LPS), lipoteichoic acid (LTA), zymosan A, chitin, chitosan and cellulose, as well as to several PGN moieties, including muramyl dipeptide (MDP), N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc). Besides, CgACP1 had the highest expression in the gill and could be greatly up-regulated quickly after bacterial challenge. This is reminiscent of the amphioxus ACP1/2 which serve as essential mucus lectins in the gill. Taken together, the current findings from mollusk and amphioxus ACPs suggest several basic common traits for the ApeC domains, including the high affinity to Lys-type PGN, the bacterial binding and agglutinating capacity, and the role as mucus proteins to protect the mucosal surface.

Encyclopaedia of eukaryotic DNA methylation: from patterns to mechanisms and functions

DNA methylation is an epigenetic modification with a very long evolutionary history. However, DNA methylation evolves surprisingly rapidly across eukaryotes. The genome-wide distribution of methylation diversifies rapidly in different lineages, and DNA methylation is lost altogether surprisingly frequently. The growing availability of genomic and epigenomic sequencing across organisms highlights this diversity but also illuminates potential factors that could explain why both the DNA methylation machinery and its genome-wide distribution evolve so rapidly. Key to this are new discoveries about the fitness costs associated with DNA methylation, and new theories about how the fundamental biochemical mechanisms of DNA methylation introduction and maintenance could explain how new genome-wide patterns of methylation evolve.

Deep whole-genome resequencing sheds light on the distribution and effect of amphioxus SNPs

Background

Amphioxus is a model organism for vertebrate evolutionary research. The significant contrast between morphological phenotypic similarity and high-level genetic polymorphism among amphioxus populations has aroused scientists' attention. Here we resequenced 21 amphioxus genomes to over 100X depth and mapped them to a haploid reference.

Results

More than 11.5 million common SNPs were detected in the amphioxus population, which mainly affect genes enriched in ion transport, signal transduction and cell adhesion, while protein structure analysis via AlphaFold2 revealed that these SNPs fail to bring effective structural variants.

Conclusions

Our work provides explanation for “amphioxus polymorphism paradox” in a micro view, and generates an enhanced genomic dataset for amphioxus research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-022-01038-w.

Ciliary photoreceptors in sea urchin larvae indicate pan-deuterostome cell type conservation

BACKGROUND

The evolutionary history of cell types provides insights into how morphological and functional complexity arose during animal evolution. Photoreceptor cell types are particularly broadly distributed throughout Bilateria; however, their evolutionary relationship is so far unresolved. Previous studies indicate that ciliary photoreceptors are homologous at least within chordates, and here, we present evidence that a related form of this cell type is also present in echinoderm larvae.

RESULTS

Larvae of the purple sea urchin Strongylocentrotus purpuratus have photoreceptors that are positioned bilaterally in the oral/anterior apical neurogenic ectoderm. Here, we show that these photoreceptors express the transcription factor Rx, which is commonly expressed in ciliary photoreceptors, together with an atypical opsin of the G_O family, opsin3.2, which localizes in particular to the cilia on the cell surface of photoreceptors. We show that these ciliary photoreceptors express the neuronal marker synaptotagmin and are located in proximity to pigment cells. Furthermore, we systematically identified additional transcription factors expressed in these larval photoreceptors and found that a majority are orthologous to transcription factors expressed in vertebrate ciliary photoreceptors, including Otx, Six3, Tbx2/3, and Rx. Based on the developmental expression of rx, these photoreceptors derive from the anterior apical neurogenic ectoderm. However, genes typically involved in eye development in bilateria, including pax6, six1/2, eya, and dac, are not expressed in sea urchin larval photoreceptors but are instead co-expressed in the hydropore canal.

CONCLUSIONS

Based on transcription factor expression, location, and developmental origin, we conclude that the sea urchin larval photoreceptors constitute a cell type that is likely homologous to the ciliary photoreceptors present in chordates.

A Potential Method for Rapid Screening of Amphioxus Founder Harboring Germline Mutation and Transgene

Amphioxus is a promising model organism for understanding the origin and evolution of vertebrates due to its basal phylogenetic position among chordates. We here compared the mutation efficacy and mutation type of tail tips and gametes of amphioxus founders injected with Cas9 protein and six different sgRNAs targeting five distinct genes, and revealed a strong correlation for mutation efficacy and a mild correlation for mutation type among the two tissues. In addition, we also observed a positive relationship between gene insertions observed in tail tips and gametes of amphioxus founders injected with Tol2 transposase and two different transgenic constructs. Finally, we showed that amphioxus larvae which had their tail tips cut at the 3-4 gill-slit stage were able to recover within 6 days and developed a normal number of gonads at the adult stage, and that F0 larvae carry similar mutation efficacy and type in the posterior end to that in the tail tips after their metamorphosis. Together, these findings suggest a great potential for obtaining valid amphioxus founders with desired mutations and transgenes at as early as the early larval stage, which will certainly speed up the generation of amphioxus mutants and transgenes and make it more cost- and labor-effective.

Two Amphioxus ApeC-Containing Proteins Bind to Microbes and Inhibit the TRAF6 Pathway

The apextrin C-terminal (ApeC) domain is a class of newly discovered protein domains with an origin dating back to prokaryotes. ApeC-containing proteins (ACPs) have been found in various marine and aquatic invertebrates, but their functions and the underlying mechanisms are largely unknown. Early studies suggested that amphioxus ACP1 and ACP2 bind to bacterial cell walls and have a role in immunity. Here we identified another two amphioxus ACPs (ACP3 and ACP5), which belong to the same phylogenetic clade with ACP1/2, but show distinct expression patterns and sequence divergence (40-50% sequence identities). Both ACP3 and ACP5 were mainly expressed in the intestine and hepatic cecum, and could be up-regulated after bacterial challenge. Both prokaryotic-expressed recombinant ACP3 and ACP5 could bind with several species of bacteria and yeasts, showing agglutinating activity but no microbicidal activity. ELISA assays suggested that their ApeC domains could interact with peptidoglycan (PGN), but not with lipoteichoic acid (LTA), lipopolysaccharides (LPS) and zymosan A. Furthermore, they can only bind to Lys-type PGN from Staphylococcus aureus, but not to DAP-type PGN from Bacillus subtilis and not to moieties of PGN such as MDPs, NAMs and NAGs. This recognition spectrum is different from that of ACP1/2. We also found that when expressed in mammalian cells, ACP3 could interact with TRAF6 via a conserved non-ApeC region, which inhibited the ubiquitination of TRAF6 and hence suppressed downstream NF-κB activation. This work helped define a novel subfamily of ACPs, which have conserved structures, and have related yet diversified molecular functions. Its members have dual roles, with ApeC as a lectin and a conserved unknown region as a signal transduction regulator. These findings expand our understanding of the ACP functions and may guide future research on the role of ACPs in different animal clades.

A Preliminary Single-Cell RNA-Seq Analysis of Embryonic Cells That Express Brachyury in the Amphioxus, Branchiostoma japonicum

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.

Transposable Elements and Stress in Vertebrates: An Overview

Since their identification as genomic regulatory elements, Transposable Elements (TEs) were considered, at first, molecular parasites and later as an important source of genetic diversity and regulatory innovations. In vertebrates in particular, TEs have been recognized as playing an important role in major evolutionary transitions and biodiversity. Moreover, in the last decade, a significant number of papers has been published highlighting a correlation between TE activity and exposition to environmental stresses and dietary factors. In this review we present an overview of the impact of TEs in vertebrate genomes, report the silencing mechanisms adopted by host genomes to regulate TE activity, and finally we explore the effects of environmental and dietary factor exposures on TE activity in mammals, which is the most studied group among vertebrates. The studies here reported evidence that several factors can induce changes in the epigenetic status of TEs and silencing mechanisms leading to their activation with consequent effects on the host genome. The study of TE can represent a future challenge for research for developing effective markers able to detect precocious epigenetic changes and prevent human diseases.

Conservative route to genome compaction in a miniature annelid

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

Genome size evolution: towards new model systems for old questions

Genome size (GS) variation is a fundamental biological characteristic; however, its evolutionary causes and consequences are the topic of ongoing debate. Whether GS is a neutral trait or one subject to selective pressures, and how strong these selective pressures are, may remain open questions. Fundamentally, the genomic sequences responsible for this variation directly impact the potential evolutionary outcomes and, equally, are the targets of different evolutionary pressures. For example, duplications and deletions of genic regions (large or small) can have immediate and drastic phenotypic effects, while an expansion or contraction of non-coding DNA is less likely to cause catastrophic phenotypic effects. However, in the long term, the accumulation or deletion of ncDNA is likely to have larger effects. Modern sequencing technologies are allowing for the dissection of these proximate causes, but a combination of these new technologies with more traditional evolutionary experiments and approaches could revolutionize this debate and potentially resolve many of these arguments. Here, I discuss an ambitious way forward for GS research, putting it in context of historical debates, theories and sometimes contradictory evidence, and highlighting the promise of combining new sequencing technologies and analytical developments with more traditional experimental evolution approaches.

Molecular Fingerprint of Amphioxus Frontal Eye Illuminates the Evolution of Homologous Cell Types in the Chordate Retina

The evolution of the vertebrate eye remains so far unresolved. Amphioxus frontal eye pigment cells and photoreceptors were proposed to be homologous to vertebrate photoreceptors and retinal pigmented epithelium, based on ultrastructural morphology and gene expression analysis in B. floridae. Here, we present comparative molecular data using two additional amphioxus species, a closely related B. lanceolatum, and the most divergent A. lucayanum. Taking advantage of a unique set of specific antibodies we characterized photoreceptors and putative interneurons of the frontal eye and investigated its neuronal circuitry. Our results corroborate generally conserved molecular fingerprint among cephalochordate species. Furthermore, we performed pharmacological perturbations and found that the Notch signaling pathway, a key regulator of retina development in vertebrates, is required for correct ratios among frontal eye cell types. In summary, our study provides a valuable insight into cell-type relationships in chordate visual organs and strengthens the previously proposed homology between amphioxus frontal eye and vertebrate eyes.

Whole-genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole-genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long-term population dynamics of Branchiostoma. The inferred ancestral N e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10-9) or third glaciation (mutation rate = 9 × 10-9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long-term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.

Deeply conserved synteny resolves early events in vertebrate evolution

Although it is widely believed that early vertebrate evolution was shaped by ancient whole-genome duplications, the number, timing and mechanism of these events remain elusive. Here, we infer the history of vertebrates through genomic comparisons with a new chromosome-scale sequence of the invertebrate chordate amphioxus. We show how the karyotypes of amphioxus and diverse vertebrates are derived from 17 ancestral chordate linkage groups (and 19 ancestral bilaterian groups) by fusion, rearrangement and duplication. We resolve two distinct ancient duplications based on patterns of chromosomal conserved synteny. All extant vertebrates share the first duplication, which occurred in the mid/late Cambrian by autotetraploidization (that is, direct genome doubling). In contrast, the second duplication is found only in jawed vertebrates and occurred in the mid-late Ordovician by allotetraploidization (that is, genome duplication following interspecific hybridization) from two now-extinct progenitors. This complex genomic history parallels the diversification of vertebrate lineages in the fossil record.

Identification of RAG-like transposons in protostomes suggests their ancient bilaterian origin

Background

V(D) J recombination is essential for adaptive immunity in jawed vertebrates and is initiated by the RAG1-RAG2 endonuclease. The RAG1 and RAG2 genes are thought to have evolved from a RAGL (RAG-like) transposon containing convergently-oriented RAG1-like (RAG1L) and RAG2-like (RAG2L) genes. Elements resembling this presumptive evolutionary precursor have thus far only been detected convincingly in deuterostomes, leading to the model that the RAGL transposon first appeared in an early deuterostome.

Results

We have identified numerous RAGL transposons in the genomes of protostomes, including oysters and mussels (phylum Mollusca) and a ribbon worm (phylum Nemertea), and in the genomes of several cnidarians. Phylogenetic analyses are consistent with vertical evolution of RAGL transposons within the Bilateria clade and with its presence in the bilaterian ancestor. Many of the RAGL transposons identified in protostomes are intact elements containing convergently oriented RAG1L and RAG2L genes flanked by terminal inverted repeats (TIRs) and target site duplications with striking similarities with the corresponding elements in deuterostomes. In addition, protostome genomes contain numerous intact RAG1L-RAG2L adjacent gene pairs that lack detectable flanking TIRs. Domains and critical active site and structural amino acids needed for endonuclease and transposase activity are present and conserved in many of the predicted RAG1L and RAG2L proteins encoded in protostome genomes.

Conclusions

Active RAGL transposons were present in multiple protostome lineages and many were likely transmitted vertically during protostome evolution. It appears that RAGL transposons were broadly active during bilaterian evolution, undergoing multiple duplication and loss/fossilization events, with the RAGL genes that persist in present day protostomes perhaps constituting both active RAGL transposons and domesticated RAGL genes. Our findings raise the possibility that the RAGL transposon arose earlier in evolution than previously thought, either in an early bilaterian or prior to the divergence of bilaterians and non-bilaterians, and alter our understanding of the evolutionary history of this important group of transposons.

Co-evolution of mutagenic genome editors and vertebrate adaptive immunity

The adaptive immune systems of all vertebrates rely on self-DNA mutating enzymes to assemble their antigen receptors in lymphocytes of their two principal lineages. In jawed vertebrates, the RAG1/2 recombinase directs V(D)J recombination of B cell and T cell receptor genes, whereas the activation-induced cytidine deaminase AID engages in their secondary modification. The recombination activating genes (RAG) 1 and 2 evolved from an ancient transposon-encoded genome modifier into a self-DNA mutator serving adaptive immunity; this was possible as a result of domestication, involving several changes in RAG1 and RAG2 proteins suppressing transposition and instead facilitating-coupled cleavage and recombination. By contrast, recent evidence supports the notion that the antigen receptors of T-like and B-like cells of jawless vertebrates, designated variable lymphocyte receptors (VLRs), are somatically assembled through a process akin to gene conversion that is believed to be dependent on the activities of distant relatives of AID, the cytidine deaminases CDA1 and CDA2, respectively. It appears, therefore, that the precursors of AID and CDAs underwent a domestication process that changed their target range from foreign nucleic acids to self-DNA; this multi-step evolutionary process ensured that the threat to host genome integrity was minimized. Here, we review recent findings illuminating the evolutionary steps associated with the domestication of the two groups of genome editors, RAG1/2 and cytidine deaminases, indicating how they became the driving forces underlying the emergence of vertebrate adaptive immune systems.

Chordate PIAS proteins act as conserved repressors of the TRAF6 self-polyubiquitination

In mammals, PIAS proteins are important SUMO E3 ligases and act as versatile regulators of over sixty different proteins, including components from the NF-κB pathways. But the PIAS functions are not well-understood due to complicated molecular mechanisms and multiple gene paralogs with overlapping roles, which is especially true in lower vertebrates where dedicated studies are scarce. As a basal chordate with a single PIAS gene, amphioxus is a convenient model to study PIAS from the evolutionary perspective. TRAF6 is a critical adaptor of the NF-κB pathways but it is not known whether TRAF6 is regulated by PIAS. Here we discover that in mammalian cells, amphioxus PIAS inhibited NF-κB activation by co-localizing and binding with TRAF6. The interaction relied on the N-terminal SAP and PINIT domains of PIAS. TRAF6 is an E3 ubiquitin ligase, which initiates downstream NF-κB signaling by promoting its self-ubiquitination. Both amphioxus SUMO1 and Ubc9 (SUMO E2 ligase) could suppress TRAF6 self-ubiquitination and NF-κB activation, suggesting that the SUMOylation activity competed away the ubiquitination activity of TRAF6. However, we show that the wild-type PIAS and the mutant PIAS without SUMO E3 ligase activity both could inhibit TRAF6-mediated NF-κB activation by reducing TRAF6 self-ubiquitination. This implies that SUMO ligase activity is not the only mechanism for PIAS to negatively regulate TRAF6. Finally, we tested the interactions between human PIAS1-4 and TRAF6. It reveals that human PIAS1, 3 and 4, but not 2, were able to repress NF-κB activation by reducing TRAF6 self-ubiquitination. Taken together, our study discovers a conserved regulatory interaction between chordate PIAS and TRAF6. It therefore sheds light on the complicated role of PIAS in immune regulation, and may help to understand the PIAS functions in other lower chordate taxa, such as jawless and jawed fishes.

Embryo bioassays with aquatic animals for toxicity testing and hazard assessment of emerging pollutants: A review

This review article gathers the available information on the use of embryo-tests as high-throughput tools for toxicity screening, hazard assessment and prioritization of new and existing chemical compounds. The approach is contextualized considering the new legal trends for animal experimentation, fostering the 3R policy, with reduction of experimental animals, addressing the potential of embryo-tests as high-throughput toxicity screening and prioritizing tools. Further, the current test guidelines, such as the ones provided by OECD and EPA, focus mainly in a limited number of animal lineages, particularly vertebrates and arthropods. To extrapolate hazard assessment to the ecosystem scale, a larger diversity of taxa should be tested. The use of new experimental animal models in toxicity testing, from a representative set of taxa, was thoroughly revised and discussed in this review. Here, we critically review current tools and the main advantages and drawbacks of different animal models and set researcher priorities.

Whole-Genome Resequencing of Twenty Branchiostoma belcheri Individuals Provides a Brand-New Variant Dataset for Branchiostoma

As the extant representatives of the basal chordate lineage, amphioxi (including the genera Branchiostoma, Asymmetron and Epigonichthys) play important roles in tracing the state of chordate ancestry. Previous studies have reported that members of the Branchiostoma species have similar morphological phenotypic characteristics, but in contrast, there are high levels of genetic polymorphisms in the populations. Here, we resequenced 20 Branchiostomabelcheri genomes to an average depth of approximately 12.5X using the Illumina HiSeq 2000 platform. In this study, over 52 million variations (~12% of the total genome) were detected in the B. belcheri population, and an average of 12.8 million variations (~3% of the total genome) were detected in each individual, confirming that Branchiostoma is one of the most genetically diverse species sequenced to date. Demographic inference analysis highlighted the role of historical global temperature in the long-term population dynamics of Branchiostoma, and revealed a population expansion at the Greenlandian stage of the current geological epoch. We detected 594 Single nucleotide polymorphism and 148 Indels in the Branchiostoma mitochondrial genome, and further analyzed their genetic mutations. A recent study found that the epithelial cells of the digestive tract in Branchiostoma can directly phagocytize food particles and convert them into absorbable nontoxic nutrients using powerful digestive and immune gene groups. In this study, we predicted all potential mutations in intracellular digestion-associated genes. The results showed that most “probably damaging” mutations were related to rare variants (MAF<0.05) involved in strengthening or weakening the intracellular digestive capacity of Branchiostoma. Due to the extremely high number of polymorphisms in the Branchiostoma genome, our analysis with a depth of approximately 12.5X can only be considered a preliminary analysis. However, the novel variant dataset provided here is a valuable resource for further investigation of phagocytic intracellular digestion in Branchiostoma and determination of the phenotypic and genotypic features of Branchiostoma.

LanceletDB: an integrated genome database for lancelet, comparing domain types and combination in orthologues among lancelet and other species

Lancelet (amphioxus) represents the most basally divergent extant chordate (cephalochordates) that diverged from the other two chordate lineages (urochordates and vertebrates) more than half a billion years ago. As it occupies a key position in evolution, it is considered as one of the best proxies for understanding the chordate ancestral state. Thus, the construction of a database with multiple lancelet genomes and gene annotation data, including protein domains, is urgently needed to investigate the loss and gain of domains in orthologues among species, especially ancient domain types (non-vertebrate-specific domains) and novel domain combination, which is helpful for providing new insight into the chordate ancestral state and vertebrate evolution. Here, we present an integrated genome database for lancelet, LanceletDB, which provides reference haploid genome sequence and annotation data for lancelet (Branchiostoma belcheri), including gene models and annotation, protein domain types, gene expression pattern in embryogenesis, different expression sequence tag sets and alternative polyadenylation (APA) sites profiled by the sequencing APA sites method. Especially, LanceletDB allows comparison of domain types and combination in orthologues among type species so as to decode the ancient domain types and novel domain combination during evolution. We also integrated the released diploid lancelet genome annotation data (Branchiostoma floridae) to expand LanceletDB and extend its usefulness. These data are available through the search and analysis page, basic local alignment search tool page and genome browser to provide an integrated display.

Functional Variation of IL-1R-Associated Kinases in the Conserved MyD88-TRAF6 Pathway during Evolution

IL-1R-associated kinases (IRAK) are important regulators in the TLR/IL-1R pathways, but their function appears inconsistent between Drosophila, bony fishes, and vertebrates. This causes a difficulty to understand the IRAK functions. As a step to reveal the evolution of IRAKs, in this study, we performed comparative and functional analysis of IRAKs by exploiting the amphioxus, a pivotal taxon connecting invertebrates and vertebrates. Sequence and phylogenetic analysis indicated three major IRAK lineages: IRAK1/2/3 is a vertebrate-specific lineage, IRAK4 is an ancient lineage conserved between invertebrate and vertebrates, and Pelle is another ancient lineage that is preserved in protostomes and invertebrate deuterostomes but lost in vertebrate deuterostomes. Pelle is closer neither to IRAK4 nor to IRAK1/2/3, hence suggesting no clear functional analogs to IRAK1/2/3 in nonvertebrates. Functional analysis showed that both amphioxus IRAK4 and Pelle could suppress NF-κB activation induced by MyD88 and TRAF6, which are unlike mammalian and Drosophila IRAKs, but, surprisingly, similar to bony fish IRAK4. Also unlike Drosophila IRAKs, no interaction was detected between amphioxus IRAK4 and Pelle, although both of them were shown capable of binding MyD88. These findings, together with previous reports, show that unlike other signal transducers in the TLR/IL-1R pathways, such as MyD88 and TRAF6, the functions of IRAKs are highly variable during evolution and very specialized in different major animal taxa. Indeed, we suggest that the functional variability of IRAKs might confer plasticity to the signal transduction of the TLR/IL-1R pathways, which in return helps the species to evolve against the pathogens.

The family of amphioxus chitin synthases offers insight into the evolution of chitin formation in chordates

Chitin is a very important and widely-used biopolymer in fungi and lower metazoans, but mysteriously disappears in mammals. Recent studies reveal that at least lower vertebrates have chitin synthases (CS) and use them to synthesize endogenous chitin. Amphioxus, a basal chordate, therefore becomes critical to understand the evolution of CS, as it occupies the transitional position from invertebrates to vertebrates, and is considered as a good proxy to the chordate ancestor. Here, by exploiting multiple genome assemblies, high-depth RNA-seq data and synteny relations, we identify 11-12 CS genes for each amphioxus species. It represents the largest CS gene pool ever found in eukaryotes so far. As comparison, most metazoans have one or two CSs. Amphioxus is the only chordate that has both the very ancient type-I CS family and the more broadly distributed type-II CS family. Specifically, amphioxus has only one type-II CS but 10-11 type-I CSs, which means that amphioxus is the only metazoan with a greatly expanded type-I CS family. Further analysis suggests that the chordate ancestor have at least one type-II CS and an expanded of type-I CS family. We hypothesize that: these ancient CSs are mostly retained in amphioxus; but the whole type-I CS family was lost in urochordates and vertebrates; the type-II CS was later duplicated into two lineages in vertebrates and followed by stochastic losses, till all type-II CSs were eventually lost in birds and mammals. Finally, our expression profiling and preliminary gene knockout analysis suggest that amphioxus CSs could have highly diverse but mildly overlapping functions in various tissues and organs. Taken together, these findings not only provide insights into the evolution of chordate CSs, lay a foundation for further functional study of the chordate CSs. After all, it is mysterious that our chordate ancestor needed so many isoenzymes for chitin formation.

Genome-Wide Identification and Transcriptomic Analysis of MicroRNAs Across Various Amphioxus Organs Using Deep Sequencing

Amphioxus is the closest living invertebrate proxy of the vertebrate ancestor. Systematic gene identification and expression profile analysis of amphioxus organs are thus important for clarifying the molecular mechanisms of organ function formation and further understanding the evolutionary origin of organs and genes in vertebrates. The precise regulation of microRNAs (miRNAs) is crucial for the functional specification and differentiation of organs. In particular, those miRNAs that are expressed specifically in organs (OSMs) play key roles in organ identity, differentiation, and function. In this study, the genome-wide miRNA transcriptome was analyzed in eight organs of adult amphioxus Branchiostoma belcheri using deep sequencing. A total of 167 known miRNAs and 23 novel miRNAs (named novel_mir), including 139 conserved miRNAs, were discovered, and 79 of these were identified as OSMs. Additionally, analyses of the expression patterns of eight randomly selected known miRNAs demonstrated the accuracy of the miRNA deep sequencing that was used in this study. Furthermore, potentially OSM-regulated genes were predicted for each organ type. Functional enrichment of these predicted targets, as well as further functional analyses of known OSMs, was conducted. We found that the OSMs were potentially to be involved in organ-specific functions, such as epidermis development, gonad development, muscle cell development, proteolysis, lipid metabolism, and generation of neurons. Moreover, OSMs with non-organ-specific functions were detected and primarily include those related to innate immunity and response to stimuli. These findings provide insights into the regulatory roles of OSMs in various amphioxus organs.

DDX23, an Evolutionary Conserved dsRNA Sensor, Participates in Innate Antiviral Responses by Pairing With TRIF or MAVS

DExD/H-box helicases play essential roles in RNA metabolism, and emerging data suggests that they have additional functions in antiviral immunity across species. However, little is known about this evolutionarily conserved family in antiviral responses in lower species. Here, through isolation of poly(I:C)-binding proteins in amphioxus, an extant basal chordate, we found that DExD/H-box helicases DHX9, DHX15, and DDX23 are responsible for cytoplasmic dsRNA detection in amphioxus. Since the antiviral roles of DDX23 have not been characterized in mammals, we performed further poly(I:C) pull-down assays and found that human DDX23 binds to LMW poly(I:C) through its N-terminal region, suggesting that DDX23 is an evolutionarily conserved dsRNA sensor. Knockdown of human DDX23 enhanced the replication of VSV and reduced the activation of the NF-κB and IRF3. Moreover, when stimulated with poly(I:C) or VSV, human DDX23 translocated from the nucleus to the cytoplasm and formed complexes with TRIF or MAVS to initiate downstream signaling. Collectively, this comparative immunological study not only defined DDX23 as an emerging nuclear pattern recognition receptor (PRR) for the innate sensing of an RNA virus, but also extended the essential role of the DExD/H helicase family in viral RNA sensing from mammals to basal chordates.

Using Amphioxus as a Basal Chordate Model to Study BMP Signaling Pathway

The BMP signaling pathway has been shown to be involved in different aspects of embryonic development across diverse metazoan phyla. Comparative studies on the roles of the BMP signaling pathway provide crucial insights into the evolution of the animal body plans. In this chapter, we present the general workflow on how to investigate the roles of BMP signaling pathway during amphioxus embryonic development. As amphioxus are basal invertebrate chordates, studies on the BMP signaling pathway in amphioxus could elucidate the functional evolution of BMP pathway in the chordate group. Here, we describe methods for animal husbandry, spawning induction, and manipulation of the BMP signaling pathway during embryonic development through drug inhibitors and recombinant proteins. We also introduce an efficient method of using mesh baskets to handle amphioxus embryos for fluorescence immunostaining and multicolor fluorescence in situ hybridization and to assay the effects of manipulating BMP signaling pathway during amphioxus embryogenesis.

Molecular and evolutionary aspects of the protochordate digestive system

The digestive system is a functional unit consisting of an endodermal tubular structure (alimentary canal) and accessory organs that function in nutrition processing in most triploblastic animals. Various morphologies and apparatuses are formed depending on the phylogenetical relationship and food habits of the specific species. Nutrition processing and morphogenesis of the alimentary canal and accessory organs have both been investigated in vertebrates, mainly humans and mammals. When attempting to understand the evolutionary processes that led to the vertebrate digestive system, however, it is useful to examine other chordates, specifically protochordates, which share fundamental functional and morphogenetic molecules with vertebrates, which also possess non-duplicated genomes. In protochordates, basic anatomical and physiological studies have mainly described the characteristic traits of suspension feeders. Recent progress in genome sequencing has allowed researchers to comprehensively detail protochordate genes and has compared the genetic backgrounds among chordate nutrition processing and alimentary canal/accessory organ systems based on genomic information. Gene expression analyses have revealed spatiotemporal gene expression profiles in protochordate alimentary canals. Additionally, to investigate the basis of morphological diversity in the chordate alimentary canal and accessory organs, evolutionary developmental research has examined developmental transcription factors related to morphogenesis and anterior-posterior pattering of the alimentary canal and accessory organs. In this review, we summarize the current knowledge of molecules involved in nutrition processing and the development of the alimentary canal and accessory organs with innate immune and endocrine roles in protochordates and we explore the molecular basis for understanding the evolution of the chordate digestive system.

Phagocytic intracellular digestion in amphioxus (Branchiostoma)

The digestive methods employed by amphioxus (Branchiostoma)—both intracellular phagocytic digestion and extracellular digestion—have been discussed since 1937. Recent studies also show that epithelial cells lining the Branchiostoma digestive tract can express many immune genes. Here, in Branchiostoma belcheri, using a special tissue fixation method, we show that some epithelial cells, especially those lining the large diverticulum protruding from the gut tube, phagocytize food particles directly, and Branchiostoma can rely on this kind of phagocytic intracellular digestion to obtain energy throughout all stages of its life. Gene expression profiles suggest that diverticulum epithelial cells have functional features of both digestive cells and phagocytes. In starved Branchiostoma, these cells accumulate endogenous digestive and hydrolytic enzymes, whereas, when sated, they express many kinds of immune genes in response to stimulation by phagocytized food particles. We also found that the distal hindgut epithelium can phagocytize food particles, but not as many. These results illustrate phagocytic intercellular digestion in Branchiostoma, explain why Branchiostoma digestive tract epithelial cells express typical immune genes and suggest that the main physiological function of the Branchiostoma diverticulum is different from that of the vertebrate liver.

Platanus-allee is a de novo haplotype assembler enabling a comprehensive access to divergent heterozygous regions

The ultimate goal for diploid genome determination is to completely decode homologous chromosomes independently, and several phasing programs from consensus sequences have been developed. These methods work well for lowly heterozygous genomes, but the manifold species have high heterozygosity. Additionally, there are highly divergent regions (HDRs), where the haplotype sequences differ considerably. Because HDRs are likely to direct various interesting biological phenomena, many genomic analysis targets fall within these regions. However, they cannot be accessed by existing phasing methods, and we have to adopt costly traditional methods. Here, we develop a de novo haplotype assembler, Platanus-allee ( http://platanus.bio.titech.ac.jp/platanus2 ), which initially constructs each haplotype sequence and then untangles the assembly graphs utilizing sequence links and synteny information. A comprehensive benchmark analysis reveals that Platanus-allee exhibits high recall and precision, particularly for HDRs. Using this approach, previously unknown HDRs are detected in the human genome, which may uncover novel aspects of genome variability.

Genome-wide gene expression analysis reveals novel insights into the response to nitrite stress in gills of Branchiostoma belcheri

Amphioxus has been widely used as a model for the comparative immunology of vertebrates. Studies have reported that gene expression changes in the amphioxus gill in response to biotic stress, such as microbial and their mimic challenge, but little is known about how gene expression is affected by abiotic stress in the marine environment, such as nitrite. A lack of information regarding gene expression response to abiotic stress hinders a comprehensive understanding of gill defense response in amphioxus. Here, RNA sequencing was used to carry out gene expression profiling analyses of Branchiostoma belcheri gills under nitrite stress. Six libraries were created for the control and treatment groups, including three biological replicates. In total, 2416 differently expressed genes (DEGs) were detected in response to nitrite stress, of which 1522 DEGs were up-regulated in the treatment group in comparison to the control, while the remaining 894 DEGs were down-regulated genes. Functional enrichment revealed that these DEGs are primarily involved in disease, innate immunity, xenobiotic biodegradation and metabolism, and biomolecular processes and apoptosis. We screened 11 key nitrite-responsive DEGs to detect their expression responses to nitrite stress at different time points, and validate the sequencing data using real time quantitative PCR. The results indicated that the expression of gene encoding CYP3A, POD, CASPR1, GST, MAO, DDH, and XDH/XO were induced, while those encoding MRC, GT, DNASE1L, and RIPK5 were reduced, to participate in the anti-nitrite response. This study provides a useful resource for research of molecular toxicology in amphioxus under environmental stress.

Identification of circular RNAs and their altered expression under poly(I:C) challenge in key antiviral immune pathways in amphioxus

Amphioxus is a key model for studying comparative immunity of vertebrates. Circular RNA (circRNA), as RNAs with a circular structure, has received little attention until recently, where several studies have reported that circRNA expression changes are involved in the immune response in animals. However, circRNA and its immune role in amphioxus have not been previously studied. Here, circRNAs in Chinese amphioxus (Branchiostoma belcheri) were sequenced, and 1859 circRNAs were identified using two algorithms (find_circ and CIRI). The analysis of miRNA target sites on circRNAs showed that 332 circRNAs may function as miRNA sponges. Furthermore, we identified circRNAs that were conserved between B. belcheri and vertebrates, tracing the origin of these circRNAs within chordates. Additionally, in combination with several key antiviral immune (poly(I:C), pIC) pathways identified in our previous B. belcheri studies, nine circRNAs potentially involved in these pathways were identified using bioinformatic predictions. Among these nine circRNAs, eight were selected to examine their expression response in B. belcheri challenged by pIC in comparison to control using real-time quantitative PCR. The results showed that four circRNAs were induced as part of the antiviral response against pIC, while expression of two circRNAs was decreased, and the expression levels of the remaining two were not significantly altered after pIC challenge. This work is the first to identify circRNAs and reveal their antiviral role in amphioxus. Therefore, it opens a new window to explore the comparative immunology of circRNAs in chordates and the regulatory roles of circRNAs in antiviral immunity in amphioxus.

Origin and evolution of plexins, semaphorins, and Met receptor tyrosine kinases

The transition from unicellular to multicellular organisms poses the question as to when genes that regulate cell-cell interactions emerged during evolution. The receptor and ligand pairing of plexins and semaphorins regulates cellular interactions in a wide range of developmental and physiological contexts. We surveyed here genomes of unicellular eukaryotes and of non-bilaterian and bilaterian Metazoa and performed phylogenetic analyses to gain insight into the evolution of plexin and semaphorin families. Remarkably, we detected plexins and semaphorins in unicellular choanoflagellates, indicating their evolutionary origin in a common ancestor of Choanoflagellida and Metazoa. The plexin domain structure is conserved throughout all clades; in contrast, semaphorins are structurally diverse. Choanoflagellate semaphorins are transmembrane proteins with multiple fibronectin type III domains following the N-terminal Sema domain (termed Sema-FN). Other previously not yet described semaphorin classes include semaphorins of Ctenophora with tandem immunoglobulin domains (Sema-IG) and secreted semaphorins of Echinoderamata (Sema-SP, Sema-SI). Our study also identified Met receptor tyrosine kinases (RTKs), which carry a truncated plexin extracellular domain, in several bilaterian clades, indicating evolutionary origin in a common ancestor of Bilateria. In addition, a novel type of Met-like RTK with a complete plexin extracellular domain was detected in Lophotrochozoa and Echinodermata (termed Met-LP RTK). Our findings are consistent with an ancient function of plexins and semaphorins in regulating cytoskeletal dynamics and cell adhesion that predates their role as axon guidance molecules.

The evolution of transcriptional repressors in the Notch signaling pathway: a computational analysis

Background

The Notch signaling pathway governs the specification of different cell types in flies, nematodes and vertebrates alike. Principal components of the pathway that activate Notch target genes are highly conserved throughout the animal kingdom. Despite the impact on development and disease, repression mechanisms are less well studied. Repressors are known from arthropods and vertebrates that differ strikingly by mode of action: whereas Drosophila Hairless assembles repressor complexes with CSL transcription factors, competition between activator and repressors occurs in vertebrates (for example SHARP/MINT and KyoT2). This divergence raises questions on the evolution: Are there common ancestors throughout the animal kingdom?

Results

Available genome databases representing all animal clades were searched for homologues of Hairless, SHARP and KyoT2. The most distant species with convincing Hairless orthologs belong to Myriapoda, indicating its emergence after the Mandibulata-Chelicarata radiation about 500 million years ago. SHARP shares motifs with SPEN and SPENITO proteins, present throughout the animal kingdom. The CSL interacting domain of SHARP, however, is specific to vertebrates separated by roughly 600 million years of evolution. KyoT2 bears a C-terminal CSL interaction domain (CID), present only in placental mammals but highly diverged already in marsupials, suggesting introduction roughly 100 million years ago. Based on the LIM-domains that characterize KyoT2, homologues can be found in Drosophila melanogaster (Limpet) and Hydra vulgaris (Prickle 3 like). These lack the CID of KyoT2, however, contain a PET and additional LIM domains. Conservation of intron/exon boundaries underscores the phylogenetic relationship between KyoT2, Limpet and Prickle. Most strikingly, Limpet and Prickle proteins carry a tetra-peptide motif resembling that of several CSL interactors. Overall, KyoT2 may have evolved from prickle and Limpet to a Notch repressor in mammals.

Conclusions

Notch repressors appear to be specific to either chordates or arthropods. Orthologues of experimentally validated repressors were not found outside the phylogenetic group they have been originally identified. However, the data provide a hypothesis on the evolution of mammalian KyoT2 from Prickle like ancestors. The finding of a potential CSL interacting domain in Prickle homologues points to a novel, very ancestral CSL interactor present in the entire animal kingdom.

Electronic supplementary material

The online version of this article (10.1186/s41065-019-0081-0) contains supplementary material, which is available to authorized users.

Genomic resources for the study of echinoderm development and evolution

Echinoderms are important research models for a wide range of biological questions. In particular, echinoderm embryos are exemplary models for dissecting the molecular and cellular processes that drive development and testing how these processes can be modified through evolution to produce the extensive morphological diversity observed in the phylum. Modern attempts to characterize these processes depend on some level of genomic analysis; from querying annotated gene sets to functional genomics experiments to identify candidate cis-regulatory sequences. Given how essential these data have become, it is important that researchers using available datasets or performing their own genome-scale experiments understand the nature and limitations of echinoderm genomic analyses. In this chapter we highlight the current state of echinoderm genomic data and provide methodological considerations for common approaches, including analysis of transcriptome and functional genomics datasets.

Evolution of cis-regulatory modules for the head organizer gene goosecoid in chordates: comparisons between Branchiostoma and Xenopus

BACKGROUND

In cephalochordates (amphioxus), the notochord runs along the dorsal to the anterior tip of the body. In contrast, the vertebrate head is formed anterior to the notochord, as a result of head organizer formation in anterior mesoderm during early development. A key gene for the vertebrate head organizer, goosecoid (gsc), is broadly expressed in the dorsal mesoderm of amphioxus gastrula. Amphioxus gsc expression subsequently becomes restricted to the posterior notochord from the early neurula. This has prompted the hypothesis that a change in expression patterns of gsc led to development of the vertebrate head during chordate evolution. However, molecular mechanisms of head organizer evolution involving gsc have never been elucidated.

RESULTS

To address this question, we compared cis-regulatory modules of vertebrate organizer genes between amphioxus, Branchiostoma japonicum, and frogs, Xenopus laevis and Xenopus tropicalis. Here we show conservation and diversification of gene regulatory mechanisms through cis-regulatory modules for gsc, lim1/lhx1, and chordin in Branchiostoma and Xenopus. Reporter analysis using Xenopus embryos demonstrates that activation of gsc by Nodal/FoxH1 signal through the 5' upstream region, that of lim1 by Nodal/FoxH1 signal through the first intron, and that of chordin by Lim1 through the second intron, are conserved between amphioxus and Xenopus. However, activation of gsc by Lim1 and Otx through the 5' upstream region in Xenopus are not conserved in amphioxus. Furthermore, the 5' region of amphioxus gsc recapitulated the amphioxus-like posterior mesoderm expression of the reporter gene in transgenic Xenopus embryos.

CONCLUSIONS

On the basis of this study, we propose a model, in which the gsc gene acquired the cis-regulatory module bound with Lim1 and Otx at its 5' upstream region to be activated persistently in anterior mesoderm, in the vertebrate lineage. Because Gsc globally represses trunk (notochord) genes in the vertebrate head organizer, this cooption of gsc in vertebrates appears to have resulted in inhibition of trunk genes and acquisition of the head organizer and its derivative prechordal plate.